離散時間系統(tǒng)時域分析與仿真-免費(fèi)閱讀
【正文】 在接下來的幾天里，我對尋找到的程序進(jìn)行一點(diǎn)一點(diǎn)的糾錯和修改，在同學(xué)的幫助下成功完成了程序。系統(tǒng)函數(shù)H(z)=z/()的零極點(diǎn)分布圖39。),ylabel(39。用MATLAB語言編程計(jì)算一個IIR線性時不變系統(tǒng)沖激響應(yīng)的絕對值的和，驗(yàn)證穩(wěn)定特性。原始序列響應(yīng)39。)% the original serialsB3=[1,]。39。A1=[,]。color39。h(n)39。axis([5 95 ])。)。x=[1zeros(1,89)]。其波形如下圖：信號抽樣波形圖離散系統(tǒng)的仿真已知某濾波器的差分方程：，求沖激響應(yīng)。\omega39。hold off。Time(sec)39。plot(t1,ft,39。grid on。)。xlabel(39。Ts=1。axis([4 16 0 3])。)。)。axis([4 16 0 3])。)。 %h(n)向量顯示范圍（添加了附加的零值）x=uDT(nx)uDT(nx4)。MATLAB求離散時間信號的卷積和命令為conv，其調(diào)用格式如下：y=conv(x,h)其中：x與h表示離散時間信號值的向量；y為卷積結(jié)果。b=[1 2]。title(39。y=filter(b,a,x)。其仿真波形如下圖：離散序列波形圖線性離散系統(tǒng)分析(1)零狀態(tài)響應(yīng)離散線性時域系統(tǒng)可用下式表示，即MATLAB中的函數(shù)filter可對上式差分方程在指定時間范圍內(nèi)的輸入序列所產(chǎn)生的響應(yīng)進(jìn)行求解。fill39。)。grid on。axis([15 15 0 1])。n2=n1+3。用MATLAB命令繪出下列離散時間序列的波形圖。g(t)39。)。xlabel(39。其實(shí)現(xiàn)的MATLAB 代碼如下：clear all。y(n)=x(n)*h(n)39。stem(ny,y,39。xlabel(39。x(n)39。stem(nx,x,39。nx=1:5。下面將介紹運(yùn)用 MATLAB 計(jì)算卷積和沖激響應(yīng)。（3）、線性若激勵f1(t)與f2(t)產(chǎn)生的響應(yīng)分別為y1(t)， y2(t)，則激勵A(yù)1f1(t)+A2f2(t)產(chǎn)生的響應(yīng)即為A1y1(t)+A2y2(t)，此性質(zhì)稱為線性。時不變系統(tǒng)就是系統(tǒng)的參數(shù)不隨時間而變化，即不管輸入信號作用的時間先后，輸出信號響應(yīng)的形狀均相同，僅是從出現(xiàn)的時間不同。用數(shù)學(xué)表示為T[x(n)]=y[n]則 T[x(nn0)]=y[nn0]，這說明序列x(n)先移位后進(jìn)行變換與它先進(jìn)行變換后再移位是等效的。（4）、時不變性若激勵f(t)產(chǎn)生的響應(yīng)為y(t)，則激勵f(tt0)產(chǎn)生的響應(yīng)即為y(tt0)，此性質(zhì)稱為不變性，也稱定常性或延遲性。(1)卷積求解：由于系統(tǒng)的零狀態(tài)響應(yīng)是激勵與系統(tǒng)的單位取樣響應(yīng)的卷積，因此，卷積運(yùn)算在離散時間信號處理領(lǐng)域被廣泛應(yīng)用。 %x(n)向量顯示范圍（添加了附加的零值）nh=2:10。fill39。)。n39。fill39。)。t=0::4。時間/s39。subplot(2,1,2)。)。設(shè)a=,N=8。n3=n12。xlabel(39。axis([15 15 0 1])。grid on。)。函數(shù)filter的調(diào)用格式如：y=filter(b,a,x)其中：x為輸入的離散序列；y為輸出的離散序列；y的長度和x的長度一樣；b和a分別為差分方程右端與左端的系數(shù)向量。stem(n,y,39。系統(tǒng)零狀態(tài)響應(yīng)y(n)39。n=0:30。用MATLAB進(jìn)行卷積運(yùn)算時，無法實(shí)現(xiàn)無限的累加，只能計(jì)算時限信號的卷積。h=.^nh.*(uDT(nh)uDT(nh8))。grid on。subplot(3,1,2)。title(39。grid on。其波形如下圖：卷積和求解波形圖 (4) 信號抽樣信號抽樣是連續(xù)時間信號分析向離散時間信號分析、連續(xù)信號處理向數(shù)字信號處理的第一步，廣泛應(yīng)用于實(shí)際的各類系統(tǒng)中。dt=。Time(sec)39。N=500。axis([10 10 *pi])。:39。)。Fs
試題試卷相關(guān)推薦
鄂ICP備17016276號-1