pi=3.1416;m=1;N=2^9;t=linspace(-10,10,N);
u=exp(-1/2*t.^(2*m));
u1=abs(u);r=52.3;
s=u.*exp(j*r*u1.^2);
s1=fft(s,N+1);
ts=t(2)-t(1);
ws=pi/ts;
w=ws*(-N/2:N/2)/N;
plot(w/(2*pi),fftshift(abs(s1).^2)/5);
程序运行完没有图像,怎么回事?
| 欢迎光临 光电工程师社区 (http://bbs.oecr.com/) | Powered by Discuz! X3.2 |