关于FFT相位角的问题

francexuqu

我在做fft变换的时候，发现相位谱得到的比较怪异，跳跃很大，估计出了问题，请各位大虾给与指点，小妹先谢了！

clear all;
fid=fopen('ft-ey-trace');                              %打开数据文件
m=fscanf(fid,'%e',[6,inf]);
fclose(fid);

subplot(3,1,1);
plot(m(1,:),m(5,:));
xlabel('t');
ylabel('E(t)');

%%%%%fft%%%%%%%%%%%%%%
N=length(m(1,:));
tp=20;                         %unit is laser period;
Ts=tp./N;              % time interval;
Ws=1./Ts;

F=fft(m(5,:));  
Fp=F(1:N./2+1)*Ts;
Ip=abs(Fp).*abs(Fp);

W=Ws.*(0:N./2)./N;
theta=angle(F(1:N./2+1));      %估计这里出现了问题 (可以unwrap)

subplot(3,1,2);
semilogy(W,Ip,'r');
xlabel('\omega/\omega_0');
ylabel('Intensity(任意单位)');

subplot(3,1,3);
plot(W,theta,'g');
xlabel('\omega/\omega_0');
ylabel('angle');
请问利用unwrap后，对应第二张图，计算的对不？

[ 本帖最后由 francexuqu 于 2007-12-7 00:07 编辑 ]

VibrationMaster

就是这样子,通常相位的导数有意义,变动小

francexuqu

这位老师你现在在线吧 ，我有问题想请教请教你，我的QQ：30146028
msn：liulifluency@hotmail.com
可以联系你吗？

VibrationMaster

发邮件到 ChenKuiFu@hotmail.com