第4章matlab在信號(hào)處理中的應(yīng)用-wenkub.com
2024-09-24 16:08 本頁(yè)面
【正文】 h(n)39。xlabel(39。stem(B,39。 A=[ones(1,16),zeros(1,N15)]。 第 4 章 MATLAB在信號(hào)處理中的應(yīng)用 [例 447] 試用頻率抽樣法設(shè)計(jì)一個(gè) FIR低通濾波器 ， 該濾波器的截止頻率為 ， 頻率抽樣點(diǎn)數(shù)為 33。 第 4 章 MATLAB在信號(hào)處理中的應(yīng)用 1. 頻率抽樣法的基本原理 對(duì)所期望的濾波器的頻率響應(yīng) Hd(ejw)進(jìn)行等間隔采樣獲得 H(k)， 利用h(n)=IDFT[H(k)]求得 FIR的單位沖激響應(yīng) 。 title(39。w/pi39。)。),ylabel(39。)。 % 矩形窗 hn1=hdn.*wn139。 nr=nr+eps*((nr)==0)。 MATLAB源程序?yàn)椋? M=21。 在 MATLAB信號(hào)處理工具箱中為用戶提供了 Boxcar (矩形 )、 Bartlet(巴特利特 )、 Hanning(漢寧 )等窗函數(shù) ， 這些窗函數(shù)的調(diào)用格式相同 。 第 4 章 MATLAB在信號(hào)處理中的應(yīng)用 FIR數(shù)字濾波器設(shè)計(jì) 格式： w = boxcar(M) 功能：返回 M點(diǎn)矩形窗序列。)。 %繪出頻率響應(yīng)曲線 plot(W*Fs/(2*pi),abs(H))。 %求濾波器階數(shù) [num,den]=ellip(N,rp,rs,wc,39。ws=[ws1,ws2]/(Fs/2)。rp=。MATLAB源程序設(shè)計(jì)如下： wp1=650。 ylabel(39。grid。 %實(shí)現(xiàn)低通向高通的轉(zhuǎn)變 [num1,den1]=ss2tf(AT,BT,CT,DT)。)。rs=40。 [例 439] 設(shè)計(jì)一個(gè)數(shù)字信號(hào)處理系統(tǒng) ， 它的采樣率為 Fs＝ 100Hz， 希望在該系統(tǒng)中設(shè)計(jì)一個(gè) Butterworth型高通數(shù)字濾波器 ， 使其通帶中允許的最小衰減為 ， 阻帶內(nèi)的最小衰減為 40dB， 通帶上限臨界頻率為 30Hz， 阻帶下限臨界頻率為 40Hz。 雙線性變換利用頻率變換關(guān)系： 11112)()(??????zzTsasHzH[例 438] MATLAB源程序如下： num=1。 采樣頻率 Fs的默認(rèn)值為 Fs=1。 Hk=fft(hn,L)。 %對(duì)序列 x(n)賦值 , M=15 hn= .^(1:20)。 2． 一維快速逆傅里葉變換函數(shù) ifft 格式： x=ifft(X, N) 功能：采用 FFT算法計(jì)算序列向量 X的 N點(diǎn) IDFT變換 。 ，求該系統(tǒng)對(duì)信號(hào) 的響應(yīng)。 y=filter(B,A,x)。 [例 433] 設(shè)系統(tǒng)差分方程為 )()1()( nxnyny ???)()( 32 nRnx n?MATLAB源程序?yàn)椋? B=1。 freqz(B,A) 該程序運(yùn)行所繪出的幅頻與 相頻性曲線如圖 。 W=0::2*pi*5。 [例 431] 已知某模擬濾波器的系統(tǒng)函數(shù) 1)( 234 ????? sssssH a求該模擬濾波器的頻率響應(yīng) 。 3. 求連續(xù)系統(tǒng)的單位階躍響應(yīng)函數(shù) step( ) 格式： [Y,T] = step(sys) 功能：返回系統(tǒng)的單位階躍響應(yīng) Y和仿真所用的時(shí)間向量 T，自動(dòng)選擇仿真的時(shí)間范圍。 ( ) ( ) / ( )g z nu m z de n z?第 4 章 MATLAB在信號(hào)處理中的應(yīng)用 時(shí)域響應(yīng)的其它函數(shù) 1. 求連續(xù) LTI系統(tǒng)的零輸入響應(yīng)函數(shù) initial( ) 格式： [y,t,x]=initial(a,b,c,d,x0) 功能：計(jì)算出連續(xù)時(shí)間 LTI系統(tǒng)由于初始狀態(tài) x0所引起的零輸入響應(yīng) y。 [ 1 ] [ ] [ ][ ] [ ] [ ]x n a x n b u ny n c x n d u n? ? ??????對(duì)輸入序列 u的響應(yīng) y和狀態(tài)記錄 x。 %系統(tǒng) 2 [nums,dens]=series(num1,den1,num2,den2) %實(shí)現(xiàn)兩個(gè)系統(tǒng)級(jí)聯(lián) [nump,denp]=parallel(num1,den1,num2,den2) %實(shí)現(xiàn)兩個(gè)系統(tǒng)并聯(lián) [numf,denf]=feedback(num1,den1,num2,den2) %實(shí)現(xiàn)兩個(gè)系統(tǒng)反饋 程序運(yùn)行結(jié)果為 ： nums = 0 0 2 ； dens = 1 3 2 nump = 0 3 4 ； denp = 1 3 2 numf = 0 1 2 ； denf = 1 3 4 因此 ， 各系統(tǒng)的傳遞函數(shù)分別為： 2322 ?? ss23432 ???sss 23 22 ???ss s第 4 章 MATLAB在信號(hào)處理中的應(yīng)用 線性時(shí)不變系統(tǒng)的響應(yīng) 線性時(shí)不變系統(tǒng)的時(shí)域響應(yīng) 1．連續(xù) LTI系統(tǒng)的響應(yīng) 2．離散 LTI系統(tǒng)的響應(yīng) ? ? ? ???? ??????? dthxthtxtxTty )()()(*)()()(用 MATLAB中的卷積函數(shù) conv( )來(lái)實(shí)現(xiàn)。 MATLAB源程序?yàn)椋? num1=1。 den=[1,1]。 其中 sos表示二次分式 ， g為比例系數(shù) ， sos為 L 6的矩陣 ， 即 01 11 21 11 210 1 2 1 211L L L L Lb b b a asosb b b a a??????? （ 4－ 15） 1． ss2tf函數(shù) 格式： [num, den]=ss2tf(A,B,C,D,iu) 功能：將指定輸入量 iu的線性系統(tǒng) (A,B,C,D)轉(zhuǎn)換為傳遞函數(shù)模型 [num,den]。 數(shù)字定義： 2. 信號(hào)功率 MATLAB實(shí)現(xiàn) : P=sum(x.*conj(x))/N。 兩序列的相關(guān)運(yùn)算 兩序列相關(guān)運(yùn)算： ? ??nmnxnxmy )()()( 21 。 MATLAB實(shí)現(xiàn)： y=x。 5．正 (余 )弦序列 )c o s ()( ???? nnx直接實(shí)現(xiàn)： n=[ns:nf]。 第 4 章 MATLAB在信號(hào)處理中的應(yīng)用 3．實(shí)指數(shù)序列 nanx ?)(n? Ra?直接實(shí)現(xiàn)： n=[ns:nf]。 圖 離散時(shí)間信號(hào)圖形 第 4 章 MATLAB在信號(hào)處理中的應(yīng)用 1．單位脈沖序列 直接實(shí)現(xiàn)： x=zeros(1,N)。)
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1