數(shù)字信號(hào)處理課設(shè)報(bào)告基于matlab測(cè)速儀設(shè)計(jì)-文庫(kù)吧資料
2025-02-28 08:48本頁(yè)面
【正文】 字濾波器組，并對(duì)濾波器組加窗降低副瓣，通過(guò)頻域方法測(cè)量速度。) %% 驗(yàn)證指標(biāo) Rs1=min(20*log10(abs(H(1:ceil(wpl/(pi/M))))))。 title(39。w(n)39。)。xlabel(39。grid subplot(325)。相頻響應(yīng) 39。)。ylabel(39。w pi*rad/s39。plot(w/pi,angle(H))。grid。YTick39。manual39。YTickMode39。,[ ])。,39。,39。 set(gca,39。幅頻響應(yīng) 39。)。ylabel(39。w pi*rad/s39。plot(w/pi,20*log10(abs(H)))。)。title(39。幅值 39。)。xlabel(39。 subplot(322)。hd(n)39。)。ylabel(39。n39。stem(n,hd)。 %% 繪圖 figure。 %% 求 FIR濾波器的頻響 [H,w]=freqz(h,1)。 %% 用兩個(gè)低通濾波器相減得到帶通濾波器 hd=(sin(m*wch)sin(m*wcl))./(pi*m)。 m=na+eps。 a=(N1)/2。 wcl=*(wpl+wstl)。 Bw=wsthwph。 wsth=*pi。 wph=*pi。 wpl=*pi。 wstl=*pi。close all。 開始 根據(jù)題目要求確定濾波器參數(shù) 采用 blackman 窗，計(jì)算濾波器階 數(shù)數(shù) 數(shù) 低通濾波器相減得到帶通濾波器 計(jì)算系統(tǒng)單位脈沖相應(yīng) hd，加窗點(diǎn)乘 分別畫出 h、 wn，不加窗濾波器和加窗濾波器 驗(yàn)證指標(biāo) 結(jié)束 ? 繪出實(shí)驗(yàn)得到的圖形 ? 附上相應(yīng)的源程序 %{ 基本題三 用窗函數(shù)設(shè)計(jì)法設(shè)計(jì)第一類線性相位 FIR 數(shù)字帶通濾波器，指標(biāo)為 低阻帶：ω 1s=， As=60dB 低通帶：ω 1p=， Rp=1dB 高通帶：ω 2p=， Rp=1dB 高阻帶：ω 2s=， As=60dB 畫出濾波器的頻率響應(yīng)曲線和單位脈沖響應(yīng)曲線，給出窗長(zhǎng)度，實(shí)際通帶衰減和阻帶衰減。) grid on 題目三： 利用窗函數(shù)法設(shè)計(jì)數(shù)字帶通濾波器，指標(biāo)為 低阻帶：ω 1s=， As=60dB 低通帶：ω 1p=， Rp=1dB 高通帶：ω 2p=， Rp=1dB 高阻帶：ω 2s=， As=60dB 分析每一種窗的時(shí)域情況和頻域抑制旁瓣情況，選擇一種合適的窗函數(shù)，畫出濾波器的頻率響應(yīng)曲線和單位脈沖響應(yīng)曲線，給出實(shí)際窗長(zhǎng)度、通 帶衰減和阻帶衰減。 subplot(2,2,4) stem(t,hn) title(39。 模擬濾波器 h(n)39。) grid on axis([0 30 0]) %******************************************** % 單位沖激響應(yīng) h(n) [hn,t] = impz(bz,az,50)。 ) ylabel(39。) xlabel(39。 subplot(2,2,2) plot(w/pi,magHzdB) title(39。 %%用沖激響應(yīng)不變法將模擬濾波器變換成數(shù)字濾波器 [Hz,w]=freqz(bz,az)。dB39。頻率 w/pi39。模擬濾波器 20log|H(jw)|39。 magHadB=20*log10(magHa)。 %將系統(tǒng)函數(shù)的零極點(diǎn)轉(zhuǎn)化為系統(tǒng)函數(shù)一般形式的系數(shù)，比如得到了系統(tǒng)函數(shù) H(s)的零極點(diǎn)的值， z 為零點(diǎn)的值， p為極點(diǎn)的值， k 為系數(shù) [Ha,w]=freqs(b,a)。 %%模擬濾波器的極點(diǎn) %% %%模擬濾波器 H(s) K=(wp^N)/e/2^(N1)。 %%求ξ i=(N+1):2*N p=sin((2*i1)*pi/2/N)*sinh(phi)+j*cos((2*i1)*pi/2/N)*cosh(phi) p=wp*p。 N=ceil(N)。 lamuda=ws/wp。 e=sqrt(10^(Rp/10)1)。 wp=*pi。 %} %% clear all clc clf Rp=1。討論 T在沖激響應(yīng)不變法中的影響。重疊保留法 N點(diǎn)圓周卷積 39。 yi(i,:)=ifft(Yi(i,:)) % 反變換 end y1=yi(:,M:L)39。 for i=1:K+1 Xi(i,:)=fft(xi(i,:))。 IFFT hn=[hn zeros(1,LM)]。 xi(i,:)=xseg。 % 后面需要補(bǔ)零的個(gè)數(shù) xn=[zeros(1,M1),xn,zeros(1,bzero)]。 % 求分組數(shù)，向下取整 leftlenx=NK*(L(M1))。 % hn 的長(zhǎng)度 L=8。 %% 重疊保留法 N=length(xn)。線性卷積 39。 stem(yn)。 figure。 %% 線性卷積 hn=ones(1,5)。100點(diǎn) DFT的模 39。|X(k)|39。)。 xlabel(39。 subplot(326)。 k=0:length(xn)1。100點(diǎn)離散時(shí)間序列 39。)。ylabel(39。n39。 stem(0:length(xn)1,xn)。 xn=cos(*pi*n)+cos(*pi*n)。) grid on。) title(39。ylabel(39。k39。 stem(k,X)。 X=abs(fft(xn))。)。 title(39。x(n)39。)。 xlabel(39。 subplot(323)。) grid on。) title(39。ylabel(39。k39。 stem(k,X)。 X=abs(fft(xn))。)。 title(39。x(n)39。)。 xlabel(39。 subplot(321)。 n=0:9。 Ts=。clc。補(bǔ)零沒(méi)有增加
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1