求救

chy19840208

谁有基于MATLAB的FIR数字滤波器的实例程序啊

happy

西安电子第二版的第四个实验

1. %实验4:用窗函数法设计FIR数字滤波器
   
2. clc;
   
3. N=input('输入窗函数长度(默认值15,选择范围[10-64]):\nN=','s');
   
4. 
   
5. N=str2num(N);
   
6. if N<1
   
7.     N=15;
   
8. else
   
9.     if N<10 | N>64
   
10.         N=15;
    
11.     end
    
12. end
    
13. w=input('输入截止频率(默认值pi/4,选择范围[0.5-pi]):\nwc=','s');
    
14. 
    
15. if w<1
    
16.     w=pi/4;
    
17. else
    
18.     if w<0.5 | w>pi
    
19.         w=pi/4;
    
20.     end
    
21. end
    
22. b=1;
    
23. close all;
    
24. 
    
25. i=0;
    
26. while(b)
    
27.     n=[0:(N-1)];
    
28. 
    
29. %产生理想低通滤波器单位脉冲响应hd(n)
    
30. alpha=(N-1)/2;
    
31. n=[0:(N-1)];
    
32. Q=n-alpha+eps;
    
33. hd=sin(w*Q)./(pi*Q);
    
34. 
    
35.     k=input('选择窗口类型(默认1):\n[1]矩形窗\n[2]Hamming窗\n[3]Hanning窗\n[4]Blackman\n选择:','s');
    
36.     switch(k)
    
37.         case{'1','2','3','4'}
    
38.             k=str2num(k);
    
39.         otherwise
    
40.             k=1;
    
41.     end
    
42.     switch(k)
    
43.         case 1
    
44.             B=boxcar(N);
    
45.             string='Boxcar';
    
46.         case 2
    
47.             B=hamming(N);
    
48.             string='Hamming';
    
49.         case 3
    
50.             B=hanning(N);
    
51.             string='Hanning';
    
52.         case 4
    
53.             B=blackman(N);
    
54.             string='Blackman';
    
55.     end
    
56.     string=[string,' N=',num2str(N)];
    
57.    
    
58.     h=hd.*(B)';
    
59.     [H,m]=freqz(h,[1],1024,'whole');
    
60.     mag=abs(H);
    
61.     db=20*log10((mag+eps)/max(mag));
    
62.     pha=angle(H);
    
63.     i=i+1;
    
64.    
    
65.     figure(i);
    
66.     subplot(2,2,1);
    
67.     n=0:N-1;
    
68.     stem(n,h,'.');
    
69.     axis([0 N-1 -0.1 0.3]);
    
70.     hold on;
    
71.     n=0:N-1;
    
72.     x=zeros(N);
    
73.     plot(n,x,'-');
    
74.     xlabel('n');
    
75.     ylabel('h(n)');
    
76.     title('实际低通滤波器的h(n)');
    
77.     text((0.3*N),0.27,string);
    
78.     hold off;
    
79.    
    
80.     subplot(2,2,2);
    
81.     plot(m/pi,db);
    
82.     axis([0 1 -100 0]);
    
83.     xlabel('w/pi');
    
84.     ylabel('20lg|H(e^j^w)|');
    
85.     title('衰减特性/dB');
    
86.     grid;
    
87.    
    
88.     subplot(2,2,3);
    
89.     plot(m,pha);
    
90.     hold on;
    
91.     n=0:7;
    
92.     x=zeros(8);
    
93.     plot(n,x,'-');
    
94.     axis([0 3.15 -4 4]);
    
95.     xlabel('w');
    
96.     ylabel('相位/rad');
    
97.     title('相频特性');
    
98.    
    
99.     subplot(2,2,4);
    
100.     plot(m,mag);
     
101.     axis([0 3.15 0 1.5]);
     
102.     xlabel('w');
     
103.     ylabel('|H(e^j^w)|');
     
104.     title('幅频特性');
     
105.     text(0.9,1.2,string);
     
106.    
     
107.     b=input('继续实验吗?(默认1)\n[1]继续\n[0]退出\n选择:','s');
     
108.     switch(b)
     
109.         case{'0','1'}
     
110.             b=str2num(b);
     
111.         otherwise
     
112.             b=1;
     
113.     end
     
114.     if b==1
     
115.         clc;
     
116.         N=input('输入窗函数长度(默认值15,选择范围[10-64]):\nN=','s');
     
117.         N=str2num(N);
     
118.         if length(N)<1
     
119.             N=15;
     
120.         else
     
121.             if N<10 | N>64
     
122.                 N=15;
     
123.             end
     
124.         end
     
125.         w=input('输入截止频率(默认值pi/4,选择范围[0.5-pi]):\nwc=','s');
     
126.         w=str2num(w);
     
127.         if length(w)<1
     
128.             w=pi/4;
     
129.         else
     
130.             if w<0.5 | w>pi
     
131.                 w=pi/4;
     
132.             end
     
133.         end
     
134.     end
     
135. end

复制代码