| 查看: 578 | 回复: 1 | |||
jiasen1980金虫 (正式写手)
|
[交流]
【求助】那位大侠帮看下哪里出问题了啊
|
|
Po；%输入光强，单位W alpha；%光纤损耗值，单位dB/km gamma；%光纤非线性参数 to；%初始脉冲宽度，单位秒 C；%第一次计算输入的啁啾参数 b2；%波数的倒数 cputime=0; tic; ln=1; i=sqrt(1); pi=3.1415926535; alph=alpha/(4.343); Ld=(to^2)/(abs(b2)); %扩散长度，单位是m Ao=sqrt(Po); %光振幅 tau =4096e12: 1e12: 4095e12; dt=1e12; h=1000;%步长 for ii=0.1:0.1:1.5 %不同的光纤长度不同，这个量可变 z=ii*Ld; u=Ao*exp((( 1+i*(C))/ 2)*(tau/to).^2); figure(1) plot(abs(u),'r'); title('Input Pulse'); xlabel('Time'); ylabel('Amplitude'); grid on; hold on; l=max(size(u)); fwhm1=find(abs(u)>abs(max(u)/2)); fwhm1=length(fwhm1); dw=1/l/dt*2*pi; w=(1* l/2:1:l/21)* dw; u=fftshift(u); %零延迟对中的谱 w=fftshift(w); %零延迟对中的谱 spectrum=fft(fftshift(u)); %快速离散傅立叶变换 for jj=h:h:z spectrum=spectrum.*exp(g1) ; %g1为线性算符e的指数表达式 f=ifft(spectrum); %快速离散反傅立叶变换 f=f.*exp(g2);%g2为非线性算符e的指数表达式 spectrum=fft(f); %快速离散傅立叶变换 spectrum=spectrum.*exp(g1) ; end f=ifft(spectrum); %快速离散反傅立叶变换 op_pulse(ln,  =abs(f);%保存在所有间隔点上的输出脉冲fwhm=find(abs(f)>abs(max(f)/2)); fwhm=length(fwhm); ratio=fwhm/fwhm1; pbratio(ln)=ratio; dd=atand((abs(imag(f)))/(abs(real(f)))); phadisp(ln)=dd;%保存脉冲相位 ln=ln+1; end toc; cputime=toc; figure(2); mesh(op_pulse(1:1:ln1, );title('Pulse Evolution'); xlabel('Time'); ylabel('distance'); zlabel('amplitude'); figure(3) plot(pbratio(1:1:ln1),' k'); xlabel('Number of steps'); ylabel('Pulse broadening ratio'); grid on; hold on; figure(5) plot(phadisp(1:1:ln1),' k'); xlabel('distance travelled'); ylabel('phase change'); grid on; hold on; disp('CPU time:'), disp(cputime); ??? Undefined function or variable 'Po' |
回复此楼» 猜你喜欢
280求调剂
已经有7人回复
-
环境领域全国重点实验室招收博士1-2名
已经有5人回复
-
303求调剂
已经有6人回复
-
290求调剂
已经有10人回复
-
0703化学求调剂 总分331
已经有3人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有8人回复
-
0703化学调剂
已经有3人回复
-
材料，纺织，生物（0856、0710），化学招生啦
已经有9人回复
-
一志愿中国海洋大学，生物学，301分，求调剂
已经有3人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有9人回复
jiasen1980
金虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 1391.2
- 散金: 1600
- 红花: 7
- 帖子: 459
- 在线: 494.7小时
- 虫号: 580859
- 注册: 2008-07-18
- 专业: 光学
2楼2010-09-09 14:16:39
