畢業(yè)設(shè)計(jì)-數(shù)字信號(hào)課程設(shè)計(jì)報(bào)告-閱讀頁(yè)
2025-04-07 09:52本頁(yè)面
【正文】 9。 %將采樣信號(hào)與噪音進(jìn)行疊加sound(x1,fs)。 %對(duì)加噪音后的信號(hào)取1024個(gè)點(diǎn)的fft變換figure(1)subplot(2,1,1)。加單頻正弦函數(shù)干擾后的信號(hào)時(shí)域圖39。ylabel(39。)。序列 n39。subplot(2,1,2)。title(39。)。幅值39。xlabel(39。)。z=*chirp(t,6000,1,9000)。sound(x1,fs)。 %對(duì)加噪音后的信號(hào)取1024個(gè)點(diǎn)的fft變換figure(1)subplot(2,1,1)。加chirp干擾信號(hào)后的信號(hào)時(shí)域圖39。ylabel(39。)。序列 n39。subplot(2,1,2)。title(39。)。幅值39。xlabel(39。)。 %將所加噪聲信號(hào)的點(diǎn)數(shù)調(diào)整到與原始信號(hào)相同a=[*sin(2*pi*2000*t)]39。 c=[*sin(2*pi*3500*t)]39。sound(x2,fs)。 %對(duì)加噪音后的信號(hào)取1024個(gè)點(diǎn)的fft變換figure(1)subplot(2,1,1)。加多個(gè)正弦干擾信號(hào)后的信號(hào)時(shí)域圖39。ylabel(39。)。序列 n39。subplot(2,1,2)。title(39。)。幅值39。xlabel(39。)。Fs=1600hz。ws=*pi。Rs=30。wp1=2/Ts*tan(wp/2)。 [N,Wn]=buttord(wp1,ws1,Rp,Rs,39。)。s39。 [bz,az]=bilinear(b,a,fs)。 %繪制幅頻特性響應(yīng)曲線figure(3)。頻率／Hz39。頻率響應(yīng)幅度39。Butterworth39。figure(4)subplot(2,1,1)plot(f1)。濾波后的時(shí)域波形39。sound(f1,fs)。f=fs*(0:511)/1024。 %畫(huà)出濾波后的頻譜圖title(39。)xlabel(39。)。fuzhi39。%用脈沖響應(yīng)不變法，設(shè)計(jì)切比雪夫低通濾波器；wp=*pi。Rp=1。Ts=1/fs。ws1=ws/Ts。s39。 %選擇濾波器的最小階數(shù)[b,a]=cheby1(N,Rp,Wn,39。)。 %用雙線性變換法實(shí)現(xiàn)模擬濾波 %創(chuàng)建butterworth模擬濾波器[H,W]=freqz(bz,az)。plot(W*fs/(2*pi),abs(H))xlabel(39。)ylabel(39。)title(39。)f1=filter(bz,az,x1)。cheb1ord39。 %畫(huà)出濾波后的時(shí)域圖title(39。)。 %播放濾波后的信號(hào)F0=fft(f1,1024)。subplot(2,1,2)。 %畫(huà)出濾波后的頻譜圖title(39。)xlabel(39。)。fuzhi39。%用窗函數(shù)法設(shè)計(jì)帶阻濾波器wlp=*pi。wus=*pi。B=wlswlp。wp=[(wls+wlp)/2/pi,(wus+wup)/2/pi]。stop39。[H,w]=freqz(hn,1024)。title(39。)。頻率 Hz39。幅度39。%sound(x2,fs)。加入chirp噪聲后的時(shí)域波形圖39。subplot(2,1,2)plot(f,abs(y3(1:512)))。加入chirp噪聲后的頻譜39。頻率 Hz39。ylabel(39。)。plot(x2)。濾波后的時(shí)域波形39。y4=fft(x2,1024)。title(39。)xlabel(39。)。幅值39。%用等波紋法設(shè)計(jì)帶阻濾波器f1=[1200,1500,3500,3800]。rp=1。dat1=(10^(rp/20)1)/(10^(rp/20)+1)。rip=[dat2,dat1,dat2]。M=M+1。x3=filter(hn,1,x2)。 y2=fft(x2,2048)。figure(1)subplot(2,1,1)plot(t,x2)。加多個(gè)正弦函數(shù)干擾信號(hào)的時(shí)域波形圖39。subplot(2,1,2)plot(f,abs(y2(1:1024)))。加多個(gè)正弦函數(shù)干擾信號(hào)的頻譜圖39。xlabel(39。)。幅值39。figure(2)subplot(2,1,1)。title(39。)。序列 n39。ylabel(39。)。plot(f,abs(y3(1:1024)))。濾波后的頻譜圖39。xlabel(39。)。幅值39。figure(3)[H,W]=freqz(hn)。xlabel(39。)。幅值39。title(
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1