Matlab中幅频函数和相频函数的误差比较

随心就动

用tfe函数求系统频率响应函数时，相频函数的误差要比幅频函数的误差大很多，大家讨论一下原因。注：tfe已被TFESTIMATE函数取代。
给定一个系统（传递函数），首先绘制幅频和相频函数；然后给系统一个脉冲输入，求得系统脉冲响应，再用tfe求系统幅频和相频函数；可以对比幅频函数和相频函数的误差。
程序

1. N=4*400;
   
2. fs=400;
   
3. n=0:N-1;
   
4. t=n/fs;
   
5. x=zeros(1,N);
   
6. x(101)=fs;
   
7. x(102)=fs;
   
8. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 生成脉冲输入信号
   
9. b=[1];
   
10. a=[1,100];
    
11. w=linspace(0,200*2*pi,401);
    
12. [mag,ph,w]=bode(b,a,w);
    
13. w_hz=(w/(2*pi))';
    
14. H_s = tf(b,a);
    
15. y=lsim(H_s,x,t);
    
16. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 生成脉冲输出信号
    
17. subplot(221);
    
18. plot(w_hz,mag);
    
19. title('幅频函数图');
    
20. grid;
    
21. subplot(222);
    
22. plot(w_hz,ph);
    
23. title('相频函数图');
    
24. grid;
    
25. 
    
26. N_fft=256;
    
27. window=hanning(N_fft);
    
28. noverlap=N_fft/2;
    
29. dflag='none';
    
30. [Txy,f]=tfe(x,y,N_fft,fs,window,noverlap,dflag);
    
31. subplot(223);
    
32. plot(f,abs(Txy));
    
33. xlabel('Frequency / Hz')
    
34. ylabel('Magnitude')
    
35. title('幅频函数');
    
36. grid;
    
37. subplot(224);
    
38. plot(f,angle(Txy)/pi*180);
    
39. xlabel('Frequency / Hz')
    
40. ylabel('Phase / Degree')
    
41. title('相频函数');
    
42. grid;

复制代码

ChaChing

建议直接贴代码, 不然我这边下载打开后会是乱码

随心就动

ChaChing 发表于 2013-6-14 16:31

                               
登录/注册后可看大图

建议直接贴代码, 不然我这边下载打开后会是乱码

没有啊，txt的文档在百度网盘上可以直接查看，我试了一下也可以下载的

ChaChing

随心就动 发表于 2013-6-14 18:40

                               
登录/注册后可看大图

没有啊，txt的文档在百度网盘上可以直接查看，我试了一下也可以下载的

LS可能不清楚我是使用繁体的!
原先的txt档下载后, 就是看到许多乱码! 原因我也不清楚
反正好像也没几行, 直接贴上好像比较好! 不是吗?

相频函数工作上较少接触, 个人也没研究学习过, 真还不清楚LS差异缘由
直觉上, 感觉应该不全然tfestimate造成!? 或许是lsim造成的!?

同待高人路过! 或许信号版块有人知道!
若有心得亦请分享下, 让我学习下

随心就动

ChaChing 发表于 2013-6-15 17:55

                               
登录/注册后可看大图

LS可能不清楚我是使用繁体的!
原先的txt档下载后, 就是看到许多乱码! 原因我也不清楚
反正好像也没几行 ...

发帖子用附件习惯了，请见谅。

用tfestimate造和lsim得到的幅频函数有很好精度，应该不是lsim造成的，并且我对比过lsim和用simulink得到的数据，两者吻合的很好。

我发现tfestimate处理脉冲响应得到的相频函数在最后一个频率相位总是pi（即-pi），表现为相频函数发生跳变，我用测试的脉冲激励信号同样出现了这一跳变现象，原因暂时还没有找到。