《matlab用于超短脉冲中啁啾与色散概念的理解》由会员分享，可在线阅读，更多相关《matlab用于超短脉冲中啁啾与色散概念的理解（13页珍藏版）》请在金锄头文库上搜索。
1、Matlab用于超短脉冲中啁啾与色散的理解超短脉冲是谐振腔中增益与损耗达到平衡、色散与非线性达到平衡时所产生的。超快光学中两个重要的概念就是啁啾与色散,而利用matlab的绘图功能可以更加形象的理解这两个概念的物理意义。一、啁啾的概念1、超短脉冲的光场在时域可以表示为：E(r,t)，A(r,t)-expitQ(t)kr(1)0其中A(r,t)为光场振幅，为中心频率，(t)为时域相位与啁啾相关，kr为空间位置0带来的相位变化。由此光场的瞬时频率可表示为：d(t)，+和丿(2)0dt若为常数，则不随时间而变化，这种脉冲称为傅里叶变换极限脉冲(transformlimitedpulse)，又称为无啁
2、啾脉冲。若与时间有关，则随时间而变化，称为啁啾脉冲(chirppulse)。当按时间的一次函数增加，则称为线性上啁啾(linearup-chirp)。当按时间的一次函数减小,则称为线性下啁啾(lineardown-chirp)。啁啾只影响光场的相位，表现为光场振荡频率的变化，而不影响光场的强度分布。2、以高斯脉冲为例，用matlab程序验证以上说明的正确性。2.1 Z=0位置处无啁啾高斯脉冲可以表示为：12E(t)，A-exp(-)-exp(it)(3)020p其中是与脉宽有关的量。其光场实部、强度及频率随时间变化图如下：(程序见附录)pHO44曲图1.1无啁啾高斯脉冲光场实部、强度及频率随时
3、间的变化2.2 Z=0位置处线性上啁啾高斯脉冲可以表示为:E(t)=Aexp(-2)expi(t+Q12)0,202pe为正的二阶相位量，其光场实部、强度及频率随时间变化图如下:2图1.2线性上啁啾高斯脉冲光场实部、强度及频率随时间的变化2.3 Z=0位置处线性上啁啾高斯脉冲可以表示为：E(t)=Aexp(-2)expi(t+Q12)0,202pe为负的二阶相位量，其光场实部、强度及频率随时间变化图如下:2图1.3线性下啁啾高斯脉冲光场实部、强度及频率随时间的变化3、讨论从图1.1、1.2、1.3中可以看出响啁啾只影响光场的相位，表现为光场振荡频率的变化，而不影响光场的强度分布。线性上啁啾和线
4、性下啁啾分别代表者光场振荡频率随时间逐渐增加和减小。二、色散的概念当光波通过介质时不同频率的光具有不同的折射率，从而具有不同的速度，传播相同的距离时间不同，由此出射的光波中不同频率成分就会分开，使脉宽发生变化，这种现象称为色散。由于超短脉冲光谱较宽，色散作用较明显，所以即使通过一个很薄的光学器件都必须要考虑色散的影响。1、脉冲在色散介质中传播的一般规律将(1)式进行傅里叶变换得：E(z,)，A(z,)expiQ()(6)其中申()，-k()z，将其展开成泰勒级数得弟()，弟()弟|(-)+士弟|(-)2+士弟|(-)30002!003!000,0,0分别称为群延时时间(groupdelay)、
5、群延迟色散(groupdelaydispersion,GDD)、三阶色散(thirdorderdispersion,TOD)。Z=0处脉冲(1)通过色散介质后光场变为E(z,t-(p)，1ei0(t-0)+0)+u/42兀0卜A(t)e(t)e-i(t-t)2/(20)dtg(8)下面对无啁啾和线性啁啾高斯脉冲进行模拟。2、无啁啾高斯脉冲通过色散介质122.1无啁啾高斯脉冲E(t)，Aexp(-2ln2)exp(it)通过色散介质(考虑到二阶0T20P色散)后场强变为：(9)E(z,t-0)，A(20)-1/2(a2+p2)-i/4expi(t-0)+0-0/2000xexp-ap2(a2+p
6、2)-112exp-iap2(a2+p2)-112其中a，(2Jln2/t)2/2,P，1/(20),0，arctan(P/a).传播后的脉宽pTp,out，1+a02/(T)41/2tpp(10)其中a，16(ln2)2，从上式可以看出不管介质是正色散还是负色散，传播后的脉冲都会加宽。出射后的相位Q(t)=(0/2)(02+T4/a)-1t2+(t-0)+0-0/2(11)p00由式子可以看出相位与时间成二次函数，则频率是时间的一次函数，产生了线性啁啾。2.2模拟验证的变比光场强屋丽时冋的菟建叮沪光场竟本陌盯间的SHEVb光埼兵卽補盯间的雯比0.B/-EQB/-S04Q2J-0yJ光IS豊度
7、越时问的旻住Lira1J7BII17BI-202先埼顛车俺时问wait*计*202Vb遥陽荽訓帕时间的喪it10-0朗|=|:3览叮搭0出-2a21血打严光1=115）的更化BJ-2QQQKID*202110-*J0Us6d2ao.Q0tTcJ.T-图2.1无啁啾高斯脉冲及通过正负色散介质后脉冲强度及频率随时间的变化2.3讨论从图2.1可以看出无啁啾高斯脉冲通过正色散介质后产生上啁啾，通过负色散介质后产生下啁啾。并且不管是正色散还是负色散介质，脉宽都会展开，峰值功率降低。3、线性啁啾高斯脉冲通过色散介质12o3.1线性啁啾高斯脉冲E(t)=Aexp(21n2)expi(/T)t2exp(i,t
8、)通过色散介0T2p0P质(考虑到二阶色散)后场强变为：E(z,t0)=A(2)-i/2g+(bp)2i/4expi,(t0)p+0/20/2000xexpB2a2(bB)2-112expiBa2b(b+B)a2(bB)2-112a=(2ln2/t)2/2,B=1/(20),b=/t,0=arctan(b/a),其中pp0=arctan(a2b2)bB/(aB)传播后的脉宽T=(120/T)2a02/(T)41/2t(10)p,outppp其中a=16(ln2)2，出射后的相位(t)=(1+20/t)(5/t)+a02/2T4/(1+20/t)2pppp+a0S/t4t2+,(t0)0+0/2
9、0/2p003.2模拟验证：线性上啁啾脉冲通过色散介质后电场与强度的变化：(11)光垢哀注耐问由莊代-1012WhIO1光场实宜逐时间的莊优D口QnoL+LIl&r光情营虞財吋问的恋it-2-I0123比hID*光场营膜闍吋鬥的娈代F-41-QQQQ2”1113跑HID14光埼冬営屈町间的兗优口-2-1013讯KID1*E42Dn-Q让址1中益4k.jTK=A-liWIrH杓F-2123H1O-*OHl陡邑呂山览J-4iIQ)tE.-2-10123旳HIO-*.5L量求KiaM图3.1线性上啁啾脉冲通过色散介质后电场与强度的变化:线性下啁啾脉冲通过色散介质后电场与强度的变化：图3.1线性下啁啾
10、脉冲通过色散介质后电场与强度的变化：3.3讨论从图3.1、3.2可以看出，啁啾脉冲通过色散介质后其脉冲可能增加也可能减小。当上啁啾脉冲通过正色散介质或下啁啾脉冲通过负色散介质时，其脉宽增加，色散更加严重；当上啁啾脉冲通过负色散介质或下啁啾脉冲通过正色散介质时，其脉宽减小，色散得到缓解。所以在超快激光光路中提供适当的色散可以压缩脉宽，提高峰值功率。这也是获得超快激光所常用的手段，比如利用棱镜对、光栅对或啁啾镜等压缩脉宽。附录：文章中所用到的matlab程序%啁啾%无啁啾高斯脉冲clearclcA0=1;lambda=800e-9;c=3e8;w0=2*pi*c/lambda;taup=5*lam
11、bda/c;t=(-3:0.01:3)*taup;%清除内存%清屏%振幅归一化%光波长选择800nm%光速%中心角频率%脉宽相关量选择为脉冲周期的5倍%选择时间范围Et=AO*exp(-(t/taup).2).*exp(i*w0*t);%无啁啾高斯脉冲subplot(1,3,1),plot(t,Et);xlabel(t/s);ylabel(E的实部);title(光场实部随时间的变化)；subplot(l,3,2),plot(t,abs(Et).2);xlabel(t/s);ylabel(光场强度)；title(光场强度随时间的变化)；w=w0;subplot(1,3,3),plot(t,w)
12、;xlabel(t/s);ylabel(频率);title(光场频率随时间的变化);%线性上啁啾高斯脉冲clearclcA0=1;lambda=800e-9;c=3e8;w0=2*pi*c/lambda;taup=5*lambda/c;t=(-3:0.01:3)*taup;phi2=(w0/10).2;%清除内存%清屏%振幅归一化%光波长选择800nm%光速%中心角频率%脉宽相关量选择为脉冲周期的5倍%选择时间范围%选择二阶相位量的大小Et=A0*exp(-(t/taup).2).*exp(i*(w0*t+phi2*t.2);%线性上啁啾高斯脉冲subplot(1,3,1),plot(t,Et
13、);xlabel(t/s);ylabel(E的实部);title(光场实部随时间的变化)；subplot(l,3,2),plot(t,abs(Et).2);xlabel(t/s);ylabel(光场强度);title(光场强度随时间的变化)；w=w0+2*phi2*t;subplot(1,3,3),plot(t,w);xlabel(t/s);ylabel(频率);title(光场频率随时间的变化);%线性下啁啾高斯脉冲clearclcA0=1;lambda=800e-9;c=3e8;w0=2*pi*c/lambda;taup=5*lambda/c;t=(-3:0.01:3)*taup;phi2=(w0/10).2;%清除内存%清屏%振幅归一化%光波长选择800nm%光速%中心角频率%脉宽相关量选择为脉冲周期的5倍%选择时间范围%选择二阶相位量的大小Et=AO*exp(-(t/taup).2).*exp(i*(w0*t+phi2*t.2);%线性下啁啾高斯脉冲subplot(1,3,1),plot(t,Et);xlabel(t/s);ylabel(E的实部);title(光场实部随时间的变化)；subplot(l,3,2),plot(t,abs(Et).2);xlabel(t/s);yl
