关于数字滤波器的幅频相频曲线图

wxm27811

题目如图所示

然后只能求出来差分方程，还不是编程求的。不知道能否编程求呢？另外，幅频相频图是不是需要采样频率呢？这个好像没有给出来啊。多谢解答！

ChaChing

若没记错, 看看Control System Toolbox !! :)

wxm27811

若没记错, 看看Control System Toolbox !!
ChaChing 发表于 2011-8-2 11:23

谢谢！怎么看呢？可以具体一下么~

wxm27811

我的代码是这样的。

1. clear;clc
   
2. B=[1 0 0 0 1];
   
3. A=[1 0 0 0 -0.8145];
   
4. H=freqz(B,A);
   
5. figure;
   
6. plot(abs(H))
   
7. figure;
   
8. plot(angle(H))
   
9. n=1:200;
   
10. x=sin(pi*n/2)+5*cos(pi*n);
    
11. y=filter(B,A,x);
    
12. figure;
    
13. stem(x)
    
14. hold on
    
15. stem(y,'r')

复制代码

但是，一看图就好像不对啊。这哪里是滤波器嘛。

wxm27811

好吧。还是解决了。
作为这个问题的结果，看看这个吧，或许将来有参考的地方。

1. clear;
   
2. B=[1,0,0,0,1];
   
3. A=[1,0,0,0,-0.8145];
   
4. figure;
   
5. freqz(B,A); %绘制幅频相频图
   
6. n=1:200;
   
7. x=sin(pi*n/2)+5*cos(pi*n);
   
8. y=filter(B,A,x);%求取滤波后的输出信号
   
9. figure;
   
10. subplot(2,1,1);
    
11. plot(x);%绘出输入信号
    
12. subplot(2,1,2);
    
13. plot(n,y,'g');%绘出输出信号
    
14. figure;
    
15. plot(n,x,'bo-',n,y,'r.-')%绘制输入输出对比图
    
16. xlim([180 200]);%仅选取稳态部分

复制代码