i tried to fix the problem i didnt know how :( please help

1 次查看（过去 30 天）

abed 2013-12-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/108553-i-tried-to-fix-the-problem-i-didnt-know-how-please-help

编辑： Friedrich 2013-12-4

clear all;
distance=15;
beta2=1;
N=1;
mshape=input('m=0 for seech, m>0 for super-gaussian= ');
chirp0=0;
nt=1024;tmax=32;
stepnum=round(20*distance*N.^2);
deltaz=distance/stepnum;
dtau=(2*tmax)/nt;
tau=(-nt/2:nt/2-1)*dtau;
omega=(pi/tmax)*(0:nt/2-1) (-nt/2:-1);
if mshape==0
    uu=sech(tau).*exp(-0.5i*chirp0*tau.^2);
else
    uu=exp(-0.5*(1+1i*chirp0).*tau.^(2*mshape));
end
%%%%%
temp=fftshift(ifft(uu)).*(nt*dtau)/sqrt(2*pi);
figure;
subplot(2,1,1);
plot(tau,abs(uu).^2,'--k');hold on;
axis([-5 5 0 inf]);
xlabel('time');
ylabel('power');
title('input and output pulse shape and spectrum');
subplot(2,1,2);
plot(fftshift(omega)/(2*pi),abs(temp).^2,'--k');hold on;
axis([-.5 .5 0 inf]);
xlabel('frequency');
ylabel('power');
%%%%%
dispersion=exp(0.5*1i*beta2*omega.^2*delta2);
hhz=1i*N^2*deltaz;
%%%%%
temp=uu.*exp(abs(uu).^2.*hhz/2);
for n=1:stepnum
    ftemp=ifft(temp).*dispersion;
    uu=fft(f_temp);
    temp=uu.*exp(abs(uu).^2.*hhz);
end
uu=temp.*exp(-abs(uu).^2.*hhz);
temp=fftshift(ifft(uu)).*(nt*dtau)/sqrt(2*pi);
subplot(2,1,1)
plot(tau,abs(uu).^2,'-k')
subplot(2,1,2)
plot(fftshift(omega)/(2*pi),abs(temp).^2,'-k')

2 个评论
显示无隐藏无

Andreas Goser 2013-12-4

What's the problem?

dpb 2013-12-4

And format the code--insert a couple of blanks before the first line of code and observe the preview window...if it isn't legible yet, read the Help...

请先登录，再进行评论。

请先登录，再回答此问题。