现代信号处理课程设计报告(6)

xlabel('N=8');ylabel('X1(k)'); m=0:15;

y2=fft(x1,16);16点fft

subplot(3,2,3);stem(m,abs(y2),'b.'); xlabel('N=16');ylabel('X2(k)'); n=0:7;

x2=[1,1,1,1,1,1,1,1];%矩形窗序列表示 figure;

subplot(3,1,1);stem(n,x2); m=0:7;

y1=fft(x2,8);%矩形窗序列的8点fft subplot(3,1,2);stem(m,abs(y1),'b.'); xlabel('N=8');ylabel('X2(k)'); m=0:15;

y2=fft(x2,16);点fft

subplot(3,1,3);stem(m,abs(y2),'b.'); xlabel('N=16');ylabel('X2(k)')

26

这是冲激序列的8点和16点fft.因为冲激序列只是在x=0处为1，所以N的点数对FFT的影响并不大。

27

这是矩形序列的8点和16点fft

4.4滤波器的设计 4.4.1、IIRDF设计

f=1000;

wp=2*20/f;ws=2*30/f;;ap=1;as=20 [n,w1]=buttord(wp,ws,ap,as); [b,a]=butter(n,w1); [H,W]=freqz(b,a,f);

subplot(2,1,1);plot(W*f/(2*pi),abs(H)); xlabel('频率');ylabel('幅值'); angH=angle(H);

subplot(2,1,2);plot(W*f/(2*pi),angH/pi); xlabel('频率');ylabel('相位');

28

这是该低通滤波器的幅频和相频特性

4.4.2、FIR滤波器的设计

As=70; ws=0.2* pi; wp=0.3* pi; tr_width=wp-ws;

M=ceil((As-7.95)*2*pi/14.36./tr_width+1)+1; disp(['滤波器的长度为',num2str(M)]); beta=0.1102*(As-8.7); n=[0:1:M-1];

disp(['线性相位斜率为',num2str(beta)])

29

delta_w=2*pi/1000;

w_kai=(kaiser(M,beta)); %调用窗函数 wc=(ws+wp)/2; r=(M-1)/2; n=[0:1:(M-1)]; m=n-r+eps;

hd=sin(wc*m)./(pi*m);%求理想脉冲响应 h=hd.*w_kai';%加窗 stem(n,w_kai); title('凯泽窗 '); axis([0 M-1 0 1.1]); ylabel('w(n)'); figure; freqz(h,1); fs=1000; t=0:1/20000:2;

x=sin(2*pi*2000*t)+sin(2*pi*5000*t)+sin(2*pi*8000*t); x1=filter(h,2,x);%对信号进行滤波 figure; n=5000:5100; subplot(211);

30