| 查看: 636 | 回复: 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' |
回复此楼» 猜你喜欢
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
今天务委会开完了，明天出结果吗
已经有13人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有8人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
jiasen1980
金虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 1391.2
- 散金: 1600
- 红花: 7
- 帖子: 459
- 在线: 494.7小时
- 虫号: 580859
- 注册: 2008-07-18
- 专业: 光学
2楼2010-09-09 14:16:39
28