关于FFT相位角的问题
我在做fft变换的时候，发现相位谱得到的比较怪异，跳跃很大，估计出了问题，请各位大虾给与指点，小妹先谢了！clear all;
fid=fopen('ft-ey-trace'); %打开数据文件
m=fscanf(fid,'%e',);
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 编辑 ] 就是这样子,通常相位的导数有意义,变动小
回复 #2 VibrationMaster 的帖子
这位老师你现在在线吧 ，我有问题想请教请教你，我的QQ：30146028msn：liulifluency@hotmail.com
可以联系你吗？ 发邮件到 ChenKuiFu@hotmail.com
页:
[1]