數(shù)字信號處理大作業(yè)

1、數(shù)字信號大作業(yè) 數(shù)字信號大作業(yè)題目一：一、設(shè)采樣頻率為 fs=2Hz，信號 ，f1=0.2Hz, 噪聲 ) 3 2 cos( 7 ) 2 sin( 10 ) ( 1 11 n ff n ff n xs s? ? ? ?為幅度為 10 的白噪聲(randn 指令產(chǎn)生)，輸入信號為 ,試用 FFT 對信號 ) (n s ? ? ? ? ) ( 1 n s n x n x ? ?作頻譜分析（幅頻特性、相頻特性） 。 ) ( ), ( ),

2、( 1 n s n x n xFFT 階次：N=256,如果對 作 FIR 低通線性相位濾波，濾波器通帶范圍（ ） ，階 ? ? n x ? 2 . 0 , 0次 33，阻帶衰減 50db,求濾波后輸出波形及幅頻特性（ ） 。 ) ( ) ( ) ( n x n h n y ? ?程序：N=256;n=0:N-1;fs=2;f1=0.2;w=2*pi/fs;x1=10*sin(w*f1*(0:N-1))+7*cos(w*3*f1*(0:

3、N-1));%x1信號sn=10*randn(1,N);%白噪聲信號xn=x1+sn;%xn信號fx=-0.5:1/N:0.5-1/N;%設(shè)定橫坐標范圍figure(1);X=fft(xn,N);%快速傅立葉變換subplot(321);plot(fx,abs(X));title('X的幅頻');subplot(322);plot(fx,angle(X(1:N))*180/pi);title('X的相頻'

4、);X1=fft(x1,N);subplot(323);plot(fx,abs(X1));title('X1的幅頻');subplot(324);plot(fx,angle(X1(1:N))*180/pi);title('X1的相頻');S=fft(sn,N);subplot(325);plot(fx,abs(S));title('S的幅頻');subplot(326);plot(fx,a

5、ngle(S(1:N))*180/pi);title('S的相頻');%FIR低通線性相位濾波，通帶范圍（0,0.2pi） ，階次33，阻帶衰減50dB%因為阻帶衰減為50dB，查表得選海明窗figure(2);%以下是低通濾波N1=33;w0=0.2*pi;%通帶范圍hn=fir1(N1-1,w0/pi);%默認是hamming窗freqz(hn,1);%濾波器的頻率響應(yīng)figure(3);Hw=filter(hn,1

6、,xn)%通過濾波后得到的信號subplot(311);plot(Hw);title('時域信號');Xa=fft(Hw);subplot(312);-0.5 0 0.505001000X信 信 信-0.5 0 0.5-2000200X信 信 信-0.5 0 0.505001000X1信 信 信-0.5 0 0.5-2000200X1信 信 信-0.5 0 0.50200400S信 信 信-0.5 0 0.5-20002

7、00S信 信 信低通濾波器海明窗（hamming）：0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1-1500-1000-5000Normalized Frequency (?? rad/sample)Phase (degrees)0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1-100-50050Normalized Frequency (?? rad/sample)Ma