全部評論(1條)
-
- lsbb198801 2016-11-12 04:41:32
- 我也在學(xué)習(xí)濾波,有空交流。 一個離散系統(tǒng)可以看作是一個濾波器,系統(tǒng)的輸出就是輸入經(jīng)過濾波器濾波的結(jié)果。那么,MATLAB實(shí)現(xiàn)這個離散系統(tǒng)的過程,也就是濾波器的實(shí)現(xiàn)過程了。 而差分方程可以表示一個系統(tǒng),其系統(tǒng)函數(shù)可以輕松求得。如,差分方程y(n)-0.8y(n-1)-0.5y(n-2)=0.7x(n)+0.3x(n-1),其脈沖響應(yīng)可用filter函數(shù)或者impz函數(shù)來求 MATLAB代碼: pulse=[1,zeros(1,63)];%定義一個單位脈沖信號 b=[0.7 0.3];%系統(tǒng)函數(shù)的分子項(xiàng)系數(shù) a=[1 -0.8 -0.5];%系統(tǒng)函數(shù)的分母項(xiàng) h1=filter(b,a,pulse);%濾波,輸入為單位脈沖,輸出為脈沖響應(yīng)。 stem(h1)
-
贊(3)
回復(fù)(0)
熱門問答
- matlab做一個濾波器
- 如題 隨便怎樣的 Z好有一定的文字說明 不要ctrl+c ctrl+v的 謝謝了
2016-11-11 10:14:37
414
1
- matlab fir濾波器
2011-09-08 05:04:49
363
3
- 如何用MATLAB做一個低通濾波器?
- 幅度|H|=1,0=2009-12-02 11:02:16 441 3
- 可以幫我做一個濾波器嗎?
2018-12-06 23:25:46
395
0
- matlab中設(shè)計(jì)fir濾波器
2011-11-30 15:44:48
376
3
- 用matlab設(shè)計(jì)一個濾波器
- 數(shù)字濾波器設(shè)計(jì)及在語音信號分析中的應(yīng)用。 步驟: 1、語音信號采集 錄制一段課程設(shè)計(jì)學(xué)生的語音信號并保存為文件,要求長度不小于10秒,并對錄制的信號進(jìn)行采樣;錄制時可以使用Windows自帶的錄音機(jī),或者使用其它專業(yè)的錄音軟件,錄制時需要配備錄音硬件... 數(shù)字濾波器設(shè)計(jì)及在語音信號分析中的應(yīng)用。 步驟: 1、語音信號采集 錄制一段課程設(shè)計(jì)學(xué)生的語音信號并保存為文件,要求長度不小于10秒,并對錄制的信號進(jìn)行采樣;錄制時可以使用Windows自帶的錄音機(jī),或者使用其它專業(yè)的錄音軟件,錄制時需要配備錄音硬件(如麥克風(fēng)),為便于比較,需要在安靜、干擾小的環(huán)境下錄音。 2、語音信號分析 使用MATLAB繪出采樣后的語音信號的時域波形和頻譜圖。根據(jù)頻譜圖求出其帶寬,并說明語音信號的采樣頻率不能低于多少赫茲。 3、含噪語音信號合成 在MATLAB軟件平臺下,給原始的語音信號疊加上噪聲,噪聲類型分為如下幾種:(1)白噪聲;(2)單頻噪色(正弦干擾);(3)多頻噪聲(多正弦干擾);(4)其它干擾,可設(shè)置為低頻、高頻、帶限噪聲,或Chirp干擾、沖激干擾。繪出疊加噪聲后的語音信號時域和頻譜圖,在視覺上與原始語音信號圖形對比,也可通過Windows播放軟件從聽覺上進(jìn)行對比,分析并體會含噪語音信號頻譜和時域波形的改變。 4、數(shù)字濾波器設(shè)計(jì)及濾波,完成以下題目中的一個 給定濾波器的規(guī)一化性能指標(biāo)(參考指標(biāo),實(shí)際中依據(jù)每個同學(xué)所疊加噪聲情況而定)例如:通帶截止頻率wp=0.25*pi, 阻通帶截止頻率ws=0.3*pi; 通帶Z大衰減Rp=1 dB; 阻帶Z小衰減Rs=15 dB,每個題目至少設(shè)計(jì)出5個用不同方法的不同類型濾波器。 題目(1):采用窗函數(shù)法與等波紋法分別設(shè)計(jì)各型FIR濾波器(低通、高通、帶通、帶阻中的至少3種類型)來對疊加噪聲前后的語音信號進(jìn)行濾波處理,繪出濾波器的頻域響應(yīng),繪出濾波后信號的時域波形和頻譜,并對濾波前后的信號進(jìn)行對比,分析信號的變化;在相同的性能指標(biāo)下比較各方法的濾波效果,并從理論上進(jìn)行分析(或解釋)。 題目(2):采用雙線性變換法與脈沖響應(yīng)不變法,分別利用不同的原型低通濾波器(Butterworth型與切比雪夫I型)來設(shè)計(jì)各型IIR濾波器(低通、高通、帶通、帶阻中的至少3種類型),繪出濾波器的頻域響應(yīng);并用這些數(shù)字濾波器對含噪語音信號分別進(jìn)行濾波處理,比較不同方法下設(shè)計(jì)出來的數(shù)字濾波器的濾波效果,并從理論上進(jìn)行分析(或解釋)。 5、回放語音信號 對濾波后的語音信號進(jìn)行回放,感覺濾波前后語音信號的變化。 6、設(shè)計(jì)一個語音信號分析與處理系統(tǒng)界面(選作) 利用MATLAB的界面設(shè)計(jì)功能,設(shè)計(jì)一個MATLAB環(huán)境下的人機(jī)交互界面,完成上述濾波器的設(shè)計(jì)及濾波過程,要求能夠接收輸入的設(shè)計(jì)參數(shù),查看設(shè)計(jì)結(jié)果及濾波結(jié)果。我的郵箱是471147021@QQ.com 哪位高手知道這方面的問題,給我點(diǎn)提示嘛,或者將答案發(fā)到我的郵箱里,謝謝了!?。? 展開
2009-06-24 07:58:04
718
2
- MATLAB中FIR濾波器如何設(shè)計(jì)
2017-09-03 07:45:45
434
1
- matlab怎么調(diào)用fdatool的濾波器
- 就是我用fdatool設(shè)計(jì)了一個濾波器 但是在matlab的主界面上怎么使用它 就是比如說有一個函數(shù)y=sin(x) 加入噪聲之后的函數(shù)叫A,怎么讓A通過這個濾波器得到我想要的波形 答得好加分
2017-09-26 04:37:28
434
2
- 如何用MATLAB設(shè)計(jì)希爾伯特濾波器
- 要求把頻率150K的正弦信號移相90度,幅頻特性為1,采樣率每周期12點(diǎn)。
2012-04-30 03:07:59
384
2
- matlab怎么調(diào)用fir濾波器的函授
2016-07-04 13:40:23
373
1
- 求matlab設(shè)計(jì)濾波器的大神解答?
- Rp=0.7; Rs=20; Wp=950*2/Fs; Ws=860*2/Fs; [n,Wn]=ellipord(Wp,Ws,Rp,Rs,'s'); [b,a]=ellip(n,Rp,Rs,Wp,'high') [H,w]=freqz(b,a,512); 本人設(shè)計(jì)的是橢圓高通數(shù)字濾波器。請問一下:Wp=950*2/Fs; Ws=860*2/Fs; 為什么不用乘以π。wp和ws代表什么?
2017-12-16 02:20:07
328
1
- 急求,matlab自適應(yīng)濾波器的程序
2007-06-11 18:54:50
583
1
- matlab中怎么讓信號通過濾波器
- %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 巴特沃夫?yàn)V波器 wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30; [N,wc]=buttord(wp,ws,Rp,As,'s')%計(jì)算率波器的階數(shù)和3dB截止頻率 [B,A]=butter(N,wc,'s');%計(jì)算濾波器系統(tǒng)函數(shù)分子分母多項(xiàng)式 fk=0:800/512:8000; wk=2*pi*f... %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 巴特沃夫?yàn)V波器 wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30; [N,wc]=buttord(wp,ws,Rp,As,'s')%計(jì)算率波器的階數(shù)和3dB截止頻率 [B,A]=butter(N,wc,'s');%計(jì)算濾波器系統(tǒng)函數(shù)分子分母多項(xiàng)式 fk=0:800/512:8000; wk=2*pi*fk; Hk=freqs(B,A,wk); figure plot(fk/1000,20*log10(abs(Hk))); grid on,xlabel('頻率(kHz)'),ylabel('幅度(dB)') title('低通濾波器') axis([0,4,-35,5]) %%%%%%%%%%%%%%%%%%%%% 信號為m=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)].*cos(2*pi*fc.*t).*cos(2*pi*fc.*t); cos(2*pi*fc.*t)相干載波 展開
2016-05-20 01:23:03
473
2
- 用matlab設(shè)計(jì)一個簡單的濾波器
- 我有一串離散的數(shù)據(jù),圖像比較混亂,其表達(dá)式可能為y=a1*sin(x+n1)+...an*sin(100x+n100) 但是x前面的系數(shù)超過20了我就不考慮了,所以我想把這個離散的數(shù)據(jù)點(diǎn)變的平滑一些,把x前面的系數(shù)超過20的全部濾掉。隔行如隔山,這方面知識確實(shí)不懂,望大神給寫個mat... 我有一串離散的數(shù)據(jù),圖像比較混亂,其表達(dá)式可能為y=a1*sin(x+n1)+...an*sin(100x+n100) 但是x前面的系數(shù)超過20了我就不考慮了,所以我想把這個離散的數(shù)據(jù)點(diǎn)變的平滑一些,把x前面的系數(shù)超過20的全部濾掉。隔行如隔山,這方面知識確實(shí)不懂,望大神給寫個matlab程序! 展開
2014-03-02 09:09:44
370
3
- matlab中fdesign濾波器filter design有什么區(qū)別
2017-03-19 11:04:38
467
1
- matlab怎么用以設(shè)計(jì)好濾波器進(jìn)行濾波
2014-12-28 15:48:19
434
2
- 如何用matlab做非因果fir濾波器
2016-07-08 05:01:46
443
1
- matlab中想用filter函數(shù)實(shí)現(xiàn)IIR濾波器
2010-12-09 08:31:47
423
2
- 在Matlab中已知頻率響應(yīng)曲線如何設(shè)計(jì)濾波器
- 如題,我已經(jīng)知道了一個頻率響應(yīng),存在數(shù)組h中,怎樣把這個設(shè)計(jì)成濾波器。
2009-04-09 19:53:24
403
1
- 濾波器matlab中wp=2*pi*fpass/fs意思
- % IIR Lowpass Use Butterworth clear; fs=20; fpass=4; fstop=5; Ap=0.5; As=10; wp=2*pi*fpass/fs; ws=2*pi*fstop/fs; omegap=tan(wp/2); omegas=tan(ws/2); ep=sqrt(10^(Ap/10)-1); es=sqrt(10^(As/10)-1); N=ceil(log(es/e... % IIR Lowpass Use Butterworth clear; fs=20; fpass=4; fstop=5; Ap=0.5; As=10; wp=2*pi*fpass/fs; ws=2*pi*fstop/fs; omegap=tan(wp/2); omegas=tan(ws/2); ep=sqrt(10^(Ap/10)-1); es=sqrt(10^(As/10)-1); N=ceil(log(es/ep)/log(omegas/omegap)); omega0=omegap/ep^(1/N); K=floor(N/2); for i=1:K theta(i)=pi*(N-1+2*i)/(2*N); end for i=1:K G(i)=omega0^2/(1-2*omega0*cos(theta(i))+omega0^2); end for i=1:K a1(i)=2*(omega0^2-1)/(1-2*omega0*cos(theta(i))+omega0^2); end for i=1:K a2(i)=(1+2*omega0*cos(theta(i))+omega0^2)/(1-2*omega0*cos(theta(i))+omega0^2); end if K<(N/2) G0=omega0/(omega0+1); a0=(omega0-1)/(omega0+1); end w=0:pi/300:pi; Hw2=1./(1+(tan(w/2)/omega0).^(2*N)); plot(w/pi,Hw2); grid; 展開
2014-03-24 05:54:13
679
2
4月突出貢獻(xiàn)榜
推薦主頁
最新話題





參與評論
登錄后參與評論