數(shù)字信號處理課設報告基于matlab測速儀設計-資料下載頁
2025-02-20 08:48本頁面
【導讀】1)用桿形圖畫出離散時間序列；2)取x，求x的10點FFT，并繪出圖形；6)為了能分辨出信號的頻譜成分，所需的頻率分辨率是多少？至少應記錄多長時。7)用FFT作出離散譜，離散譜線與實際頻率如何對應？點，并從每一段產(chǎn)生的輸出中取出m個樣點，再連接在一起以得到所要求的輸出序列。計濾波算法程序。試用h={1,1,1,1,1}驗證你的程序，并。繪出線性卷積輸出結果和通過FFT重疊保留法的輸出結果?；卮饐栴}并寫明計算依據(jù)對結果進行分析比較?？偨Y實驗所得的主要結論。本實驗表明，時域原序列補零不能提高譜分析的分辨率F。是在記錄時間內(nèi)的時域抽樣點數(shù)。補零沒有增加時域“有效”數(shù)據(jù)，所以不會提高譜分析精
【正文】 寫出設計過程中所遇難點及解決的方法 五、 結論和分析 給出設計的主要結論和仿真結果，對結論進行分析和對比 設計窗函數(shù)并求頻響 設計帶通濾波器并加窗 判斷速度，顯示結果 信號譜分析 設計濾波器組并進行旁瓣抑制 結 束 六、 代碼和結果 %{ 提高題 基于 Matlab 測速儀設計 ? 目的 利用離散付里葉變換進行頻譜分析，設計數(shù)字濾波器組，并對濾波器組加窗降低副瓣，通過 頻 域方法測量速度。 ? 要求 測速精度： 5km/h ? 步驟 1. 按技術指標要求確定抽樣頻率，獲得抽樣輸入信號。 2. 按技術指標要求確定傅立葉變換的點數(shù)，對輸入信號進行頻域分析。 3. 設計窄帶濾波器組。 4. 設計加權窗以降低副瓣。 5. 采用門限等方法確定目標速度。 測速范圍： 精度： 70km/h 150km/h 5km/h 換算： fd: 606Hz 1300Hz %} %% %%設計回波信號 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。) %% %%設計帶通濾波器 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 帶通濾波器的幅頻響應 39。) %% %%設計濾波器組 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 帶通濾波器組的幅頻響應 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 帶通濾波器組加窗后的幅頻響應 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。對應濾波第 %d\n39。,result)。 運行結果： 六、 拓展問題 長度太短會對什么造成影響？ 頻域加權各自的用途。各種窗形式的衰減 、主瓣寬度等參數(shù)及用途。 FFT對模擬信號作譜分析是一種近似的譜分析。取 FFT的點數(shù)多（信號補零）一些會得到怎樣的結果？ 第三部分： 課設總結
法律信息相關推薦
鄂ICP備17016276號-1