matlab用于超短脉冲中啁啾与色散概念的理解

(1+2*phi2*delta/taup)*(delta/taup)+a*phi2^2/(2*taup^4))/((1+2*phi2*delta/taup)^2+...a*phi2^2/taup^4))*t+w0;subplot(3,3,6),plot(t,wout);xlabel('t/s');ylabel('频率');title('光场频率随时间的变化');%线性上啁啾高斯脉冲通过负色散介质phi2>0phi0=0;phi1=0;phi2=(4/w0).^2;%选择二阶相位量的大小alpha=(2*sqrt(log(2))/taup).^2/2;beta=1/(2*phi2);theta=atan(beta/alpha);Eout=A0*(2*phi2)^(-0.5)*(alpha^2+(b+beta)^2)^(-0.25)*exp(i*(w0*(t-phi1)+phi0+theta2/2-theta/2))....*exp(-alpha*beta^2*(alpha^2+(b+beta)^2)^-1*t.^2).*exp(-i*beta*(alpha^2+b*(b+beta))*(alpha^2+(b+beta)^2)^-1*t.^2);a=16*(log(2))^2;taupuot=(1+a*phi2^2/(taup)^4)^0.5*taup;phiout=-((phi2/2)*(phi2^2+taup^4/a)^-1)*t.^2+w0*phi1+phi0-theta/2;subplot(3,3,7),plot(t,Eout);xlabel('t/s');ylabel('负色散介质Eout的实部');title('光场实部随时间的变化');subplot(3,3,8),plot(t,abs(Eout).^2);xlabel('t/s');ylabel('光场强度');