數(shù)字信號處理課設(shè)報告基于matlab測速儀設(shè)計-資料下載頁
2025-02-20 08:48本頁面
【導(dǎo)讀】1)用桿形圖畫出離散時間序列；2)取x，求x的10點FFT，并繪出圖形；6)為了能分辨出信號的頻譜成分，所需的頻率分辨率是多少？至少應(yīng)記錄多長時。7)用FFT作出離散譜，離散譜線與實際頻率如何對應(yīng)？點，并從每一段產(chǎn)生的輸出中取出m個樣點，再連接在一起以得到所要求的輸出序列。計濾波算法程序。試用h={1,1,1,1,1}驗證你的程序，并。繪出線性卷積輸出結(jié)果和通過FFT重疊保留法的輸出結(jié)果?；卮饐栴}并寫明計算依據(jù)對結(jié)果進(jìn)行分析比較?？偨Y(jié)實驗所得的主要結(jié)論。本實驗表明，時域原序列補(bǔ)零不能提高譜分析的分辨率F。是在記錄時間內(nèi)的時域抽樣點數(shù)。補(bǔ)零沒有增加時域“有效”數(shù)據(jù)，所以不會提高譜分析精
【正文】 寫出設(shè)計過程中所遇難點及解決的方法 五、 結(jié)論和分析 給出設(shè)計的主要結(jié)論和仿真結(jié)果，對結(jié)論進(jìn)行分析和對比 設(shè)計窗函數(shù)并求頻響 設(shè)計帶通濾波器并加窗 判斷速度，顯示結(jié)果 信號譜分析 設(shè)計濾波器組并進(jìn)行旁瓣抑制 結(jié) 束 六、 代碼和結(jié)果 %{ 提高題 基于 Matlab 測速儀設(shè)計 ? 目的 利用離散付里葉變換進(jìn)行頻譜分析，設(shè)計數(shù)字濾波器組，并對濾波器組加窗降低副瓣，通過 頻 域方法測量速度。 ? 要求 測速精度： 5km/h ? 步驟 1. 按技術(shù)指標(biāo)要求確定抽樣頻率，獲得抽樣輸入信號。 2. 按技術(shù)指標(biāo)要求確定傅立葉變換的點數(shù)，對輸入信號進(jìn)行頻域分析。 3. 設(shè)計窄帶濾波器組。 4. 設(shè)計加權(quán)窗以降低副瓣。 5. 采用門限等方法確定目標(biāo)速度。 測速范圍： 精度： 70km/h 150km/h 5km/h 換算： fd: 606Hz 1300Hz %} %% %%設(shè)計回波信號 clc。clear all。 v1=input(39。輸入速度 v1=?39。)。 fs=6000。 %%fs為 4~5倍 fd Ts=1/fs。 N=1200。%% N 如何取值？？？ n=0:N1。 t=n*Ts。 %%v=75/。 v=(v1*1000)./3600。 fd=(2*v)./(*)。%%計算多普勒頻移 y=sin(2*pi*fd*t)。%% subplot(4,2,1)。 y1=y(1:100)。 plot(y1)。 grid。title(39。(a)回波信號時域波形 39。) result=(fd606)/。 result=ceil(result)。 %% %%blackman窗 Bw=2*pi*150/fs。 %%確定數(shù)字濾波器過渡帶寬 M=ceil(11*pi/Bw)。 %%窗函數(shù)的圖像 wn=(blackman(M))39。 subplot(4,2,2)。 plot(wn)。 grid。title(39。(b) 窗函數(shù)圖形 39。) %% %%求窗函數(shù)的頻響 Hw=abs(fft(wn,1024))。 f=[0:1023]*fs/1024。 subplot(4,2,3)。 plot(f,20*log10(Hw/max(Hw)))。 axis([0,fs/10,100,20])。 grid。title(39。(c) 窗函數(shù)幅頻特性 39。) %% %%設(shè)計帶通濾波器 wcl=2*pi*。 wch=2*pi*。 a=(M1)/2。 n=0:M1。 m=na+eps。 %通過加任意小值 eps避免 %了 n=a 時零比零的出現(xiàn)。 hd=(sin(m*wch)sin(m*wcl))./(pi*m)。%% 用兩個低通濾波器相減得到帶通濾波器 %% %%加窗 h=hd.*(blackman(M))39。 %% %%求 FIR帶通濾波器的頻響 Ha=abs(fft(hd,1024))。 f=[0:1023]*fs/1024。 subplot(4,2,4)。 plot(f,20*log10(Ha/max(Ha)))。 axis([0,fs/3,60,20])。 grid。title(39。(d) FlR 帶通濾波器的幅頻響應(yīng) 39。) %% %%設(shè)計濾波器組 L=16。 for i=0:L wcl=2*pi*(+*i)/fs。 wch=2*pi*(+*i)/fs。 h=(sin(m*wch)sin(m*wcl))./(pi*m)。 Ha=abs(fft(h,1024))。 subplot(4,2,5)。 plot(f,20*log10(Ha/max(Ha)))。 axis([0,fs/3,60,20])。 grid。title(39。(e) FlR 帶通濾波器組的幅頻響應(yīng) 39。) hold on。 h1=h.*(blackman(M))39。 Hb=abs(fft(h1,1024))。 subplot(4,2,6)。 plot(f,20*log10(Hb/max(Hb)))。 axis([0,fs/3,60,20])。 hold on。 grid。title(39。(f) FlR 帶通濾波器組加窗后的幅頻響應(yīng) 39。)。 Hb=abs(fft(h1,1024))。 subplot(4,2,7)。 plot(f,20*log10(Hb/max(Hb)))。 axis([500,fs/3,60,20])。 hold on。 end %% %%對回波信號做 FFT Y=abs(fft(y))。 subplot(4,2,7)。 f=0:fs/N:fsfs/N。 plot(f,20*log10(Y/max(Y)),39。red39。)。 axis([0,fs/3,60,20])。 fprintf(39。速度為 :%d\n39。,v1)。 fprintf(39。對應(yīng)濾波第 %d\n39。,result)。 運行結(jié)果： 六、 拓展問題 長度太短會對什么造成影響？ 頻域加權(quán)各自的用途。各種窗形式的衰減 、主瓣寬度等參數(shù)及用途。 FFT對模擬信號作譜分析是一種近似的譜分析。取 FFT的點數(shù)多（信號補(bǔ)零）一些會得到怎樣的結(jié)果？ 第三部分： 課設(shè)總結(jié)
法律信息相關(guān)推薦
鄂ICP備17016276號-1