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