畢業(yè)設(shè)計-數(shù)字信號課程設(shè)計報告(已改無錯字)
2023-04-23 09:52:52 本頁面
【正文】 subplot(2,1,1)。plot(x1) title(39。加單頻正弦函數(shù)干擾后的信號時域圖39。)。ylabel(39。幅值39。)。xlabel(39。序列 n39。)。subplot(2,1,2)。plot(f,abs(y1(1:512)))。title(39。加單頻正弦函數(shù)干擾后的信號頻譜圖39。)。ylabel(39。幅值39。)。xlabel(39。Hz39。)。t=0:1/fs:1/fs*(length(x)1)。z=*chirp(t,6000,1,9000)。 x1=x+z39。sound(x1,fs)。y1=fft(x1,1024)。 %對加噪音后的信號取1024個點的fft變換figure(1)subplot(2,1,1)。plot(x1) title(39。加chirp干擾信號后的信號時域圖39。)。ylabel(39。幅值39。)。xlabel(39。序列 n39。)。subplot(2,1,2)。plot(f,abs(y1(1:512)))。title(39。加chirp干擾信號后的信號頻譜圖39。)。ylabel(39。幅值39。)。xlabel(39。Hz39。)。t=0:1/fs:(size(x)1)/fs。 %將所加噪聲信號的點數(shù)調(diào)整到與原始信號相同a=[*sin(2*pi*2000*t)]39。 b=[*sin(2*pi*3000*t)]39。 c=[*sin(2*pi*3500*t)]39。 x2=x+a+b+c。sound(x2,fs)。y2=fft(x2,1024)。 %對加噪音后的信號取1024個點的fft變換figure(1)subplot(2,1,1)。plot(x2) title(39。加多個正弦干擾信號后的信號時域圖39。)。ylabel(39。幅值39。)。xlabel(39。序列 n39。)。subplot(2,1,2)。plot(f,abs(y2(1:512)))。title(39。加多個正弦干擾信號后的信號頻譜圖39。)。ylabel(39。幅值39。)。xlabel(39。Hz39。)。%雙線性變換法設(shè)計Butterworth低通濾波器%Fp=1000hz。Fs=1600hz。wp=*pi。ws=*pi。Rp=1。Rs=30。Ts=1/fs。wp1=2/Ts*tan(wp/2)。 %將模擬指標(biāo)轉(zhuǎn)換成數(shù)字指標(biāo)ws1=2/Ts*tan(ws/2)。 [N,Wn]=buttord(wp1,ws1,Rp,Rs,39。s39。)。 %選擇濾波器的最小階數(shù)[b,a]=butter(N,Wn,39。s39。)。 [bz,az]=bilinear(b,a,fs)。 %用雙線性變換法實現(xiàn)模擬濾波 %創(chuàng)建butterworth模擬濾波器[H,W]=freqz(bz,az)。 %繪制幅頻特性響應(yīng)曲線figure(3)。plot(W*fs/(2*pi),abs(H))xlabel(39。頻率／Hz39。)ylabel(39。頻率響應(yīng)幅度39。)title(39。Butterworth39。)f1=filter(bz,az,x1)。figure(4)subplot(2,1,1)plot(f1)。 %畫出濾波后的時域圖title(39。濾波后的時域波形39。)。sound(f1,fs)。 %播放濾波后的信號F0=fft(f1,1024)。f=fs*(0:511)/1024。subplot(2,1,2)F1=plot(f,abs(F0(1:512)))。 %畫出濾波后的頻譜圖title(39。濾波后的頻譜39。)xlabel(39。Hz39。)。ylabel(39。fuzhi39。)。%用脈沖響應(yīng)不變法，設(shè)計切比雪夫低通濾波器；wp=*pi。ws=*pi。Rp=1。Rs=30。Ts=1/fs。wp1=wp/Ts。ws1=ws/Ts。 %將模擬指標(biāo)轉(zhuǎn)換成數(shù)字指標(biāo)[N,Wn]=cheb1ord(wp1,ws1,Rp,Rs,39。s39。)。 %選擇濾波器的最小階數(shù)[b,a]=cheby1(N,Rp,Wn,39。s39。)。 [bz,az]=bilinear(b,a,fs)。 %用雙線性變換法實現(xiàn)模擬濾波 %創(chuàng)建butterworth模擬濾波器[H,W]=freqz(bz,az)。 figure(2)。plo
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1