下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、FFT算法和低通濾波器去除噪聲作用的比較我們可以把FFT簡(jiǎn)單地看作一個(gè)變換器,輸入N+1個(gè)數(shù),輸出N+1個(gè)數(shù),但他們對(duì)應(yīng)的意義不同,如果把輸入當(dāng)作時(shí)域,則輸出為頻域,表征了其對(duì)應(yīng)域的變化快慢。 假設(shè)輸入信號(hào)本身的頻率為fc(或者說頻帶寬為fc),被頻率為fs的沖擊 串采樣(由采樣定理,fs >= 2*fc),則變換前的N1個(gè)數(shù)字對(duì)應(yīng)的x軸為t0,t1,tN=0,Ts,2*Ts,.,N*Ts (其中Ts為1/fs,為采樣周期)則變換后的N+1個(gè)數(shù)對(duì)應(yīng)的x軸變?yōu)轭l率,范圍為0fs,以fs/N為間隔,既為頻率點(diǎn)0,fs/N,2*fs/N,fs,在matlab中如果用fftshift(fft(d
2、ata),則變換后對(duì)應(yīng)x軸為-fs/2fs/2,如果滿足采樣定理的話,信號(hào)頻帶-fcfc就包含在轉(zhuǎn)換后的頻譜里面了,就不會(huì)有失真??紤]信號(hào):y(t)=exp-(sin2t+2cos4t+0.4sintsin50t)在0,2pi的情況,首先利用FFT算法對(duì)其進(jìn)行去噪聲處理,即使高頻部分信號(hào)為零,利用MATLAB進(jìn)行仿真,其程序及仿真結(jié)果如下:t=linspace(0,2*pi,28); %discretizes 0,2pi into 256 nodesy=exp(-(cos(t).2).*(sin(2*t)+2*cos(4*t)+0.4*sin(t).*sin(50*t);fy=fft(y);
3、%compute fft of yfilterfy=fy(1:6) zeros(1,28-12) fy(28-5:28); %sets fftfiltery=ifft(filterfy); %computes inverse fft of the filtered ffthold onplot(t,y,'r*') %generates the graph of the original signalplot(t,filtery) %generates the plot of the compressed signaltitle('filtered signal with
4、 FFT * unfilered signal')hold off原信號(hào)波形:FFT仿真結(jié)果與原信號(hào)的比較:下面設(shè)計(jì)一個(gè)底通濾波器實(shí)現(xiàn)消噪的功能:t=linspace(0,2*pi,28); %時(shí)間向量ts=2*pi/(28); %抽樣間隔fs=1/ts;df=1/2*pi; %頻率分辨率y=exp(-(cos(t).2).*(sin(2*t)+2*cos(4*t)+0.4*sin(t).*sin(50*t);Y,y,df1=fftseq(y,ts,df); %對(duì)信號(hào)進(jìn)行傅立葉變換 Y=Y/fs;%設(shè)計(jì)一個(gè)低通濾波器f_cutoff=4/2*pi; %低通濾波器截止頻率n_cutoff
5、=floor(f_cutoff/df1);f=0:df1:df1*(length(y)-1)-fs/2;H=zeros(size(f);H(1:n_cutoff)=ones(1,n_cutoff);H(length(f)-n_cutoff+1:length(f)=ones(1,n_cutoff);DEM=H.*Y; %經(jīng)過低通濾波器后的信號(hào)頻譜signal=real(ifft(DEM)*fs; %經(jīng)過傅立葉反變換后的信號(hào)hold onplot(t,y,'g*')plot(t,signal(1:length(t)title('filtered signal with LP
6、F unfiltered signal')hold off調(diào)用的fftseq子函數(shù):function M,m,df=fftseq(m,ts,df) %M,m,df=fftseq(m,ts,df)%M,m,df=fftseq(m,ts)%FFTSEQGenerates M, the FFT of the sequence m.%The sequence is zero padded to meet the required frequency resolution df.%ts is the sampling interval. The output df is the final fr
7、equency resolution.%Output m is the zero padded version of input m. M is the FFT.fs=1/ts;if nargin = 2 n1=0;else n1=fs/df;endn2=length(m);n=2(max(nextpow2(n1),nextpow2(n2);M=fft(m,n);m=m,zeros(1,n-n2);df=fs/n;經(jīng)過低通濾波器后的信號(hào)波形:LPF仿真結(jié)果與原信號(hào)的比較:由于有用信號(hào)的頻率往往比較低,而噪聲信號(hào)往往是高頻信號(hào),所以我們可以使高頻部分為零,來(lái)達(dá)到去除噪聲的目的。通過以上的比較,我們可以發(fā)現(xiàn),利用F
溫馨提示
- 1. 本站所有資源如無(wú)特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 人教版第13章第三節(jié)比熱容教案
- Unit 3 My father is a writer. Lesson 14(教學(xué)設(shè)計(jì))-2024-2025學(xué)年人教精通版英語(yǔ)五年級(jí)上冊(cè)2
- 建立正向激勵(lì)機(jī)制的步驟計(jì)劃
- 主管工作總結(jié)的目標(biāo)任務(wù)規(guī)劃計(jì)劃
- 2025年高考數(shù)學(xué)總復(fù)習(xí) 38 第五章 第一節(jié) 平面向量的概念與線性運(yùn)算
- 常見班級(jí)管理誤區(qū)及糾正計(jì)劃
- 品牌形象重塑的成功案例分析計(jì)劃
- 廣東省興寧一中2021-2022學(xué)年物理高一下期末達(dá)標(biāo)測(cè)試試題含解析
- 廣東省廣州越秀區(qū)培正中學(xué)2021-2022學(xué)年高一物理第二學(xué)期期末學(xué)業(yè)質(zhì)量監(jiān)測(cè)試題含解析
- 甘肅天水一中2022年物理高一第二學(xué)期期末聯(lián)考試題含解析
- 危險(xiǎn)化學(xué)品事故救援和應(yīng)急處理課件
- 物品寄存管理制度(3篇)
- 商業(yè)經(jīng)濟(jì)學(xué)課件祝合良
- 保險(xiǎn)合同糾紛審理若干問題課件
- 九年級(jí)第二次教師會(huì)議課件
- KET詞匯表(英文中文完整版)
- 《基于項(xiàng)目式學(xué)習(xí)的初中語(yǔ)文散文教學(xué)探究(論文)》
- 混凝土樓板拆除施工方案(可編輯)
- 新老物業(yè)移交表格(全套)
- 五年級(jí)上冊(cè)數(shù)學(xué)課件-8.2 用字母表示數(shù)丨蘇教版 (共27張PPT)
- 微波消解儀校準(zhǔn)規(guī)范試驗(yàn)報(bào)告
評(píng)論
0/150
提交評(píng)論