基于matlab的離散系統(tǒng)時(shí)域分析課程設(shè)計(jì)

1、<p>　　課程名稱: 數(shù)字信號(hào)處理</p><p>　　系 別：信息工程系 </p><p>　　專業(yè)班級(jí)：通信工程專業(yè)</p><p>　　課程設(shè)計(jì)題目：基于MATLAB的離散系統(tǒng)的時(shí)域分析及實(shí)現(xiàn)</p><p><b>　　一.引言</b></p><p>　　MATLAB是

2、由美國mathworks公司發(fā)布的主要面對(duì)科學(xué)計(jì)算、可視化以及交互式程序設(shè)計(jì)的高科技計(jì)算環(huán)境。它將數(shù)值分析、矩陣計(jì)算、科學(xué)數(shù)據(jù)可視化以及非線性動(dòng)態(tài)系統(tǒng)的建模和仿真等諸多強(qiáng)大功能集成在一個(gè)易于使用的視窗環(huán)境中，為科學(xué)研究、工程設(shè)計(jì)以及必須進(jìn)行有效數(shù)值計(jì)算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計(jì)語言（如C、Fortran）的編輯模式，代表了當(dāng)今國際科學(xué)計(jì)算軟件的先進(jìn)水平?！?lt;/p>&l

3、t;p><b>　　二.設(shè)計(jì)目的</b></p><p>　　掌握離散時(shí)間序列卷積和MATLAB實(shí)現(xiàn)的方法。</p><p>　　掌握離散系統(tǒng)的單位響應(yīng)及其MATLAB的實(shí)現(xiàn)方法。</p><p>　　掌握用MATLAB表示離散系統(tǒng)的時(shí)域響應(yīng)。</p><p><b>　　二.設(shè)計(jì)題目及要求</b&

4、gt;</p><p>　　題目：利用MATLAB實(shí)現(xiàn)對(duì)離散系統(tǒng)時(shí)域的分析</p><p><b>　　要求：</b></p><p>　　1.熟悉離散時(shí)間序列的卷積和原理</p><p>　　2.熟悉編程實(shí)現(xiàn)離散時(shí)間序列的卷積和</p><p>　　3.熟悉離散系統(tǒng)的單位響應(yīng)及其常用函數(shù)的調(diào)用格

5、式</p><p><b>　　三.設(shè)計(jì)原理</b></p><p>　　1. （1）單位序列（單位脈沖序列）MATLAB描述。</p><p>　　MATLAB函數(shù)可寫為</p><p>　　k=[k1:k2]; %k1,k2為時(shí)間序列的起始及終止時(shí)間序號(hào)</p><p&g

6、t;　　fk=[(k-k0)==0]; %k0為§(k)在時(shí)間軸上的位移量</p><p>　　stem(k,fk)</p><p>　?。?）單位階躍序列MATLAB描述。</p><p>　　MARLAB函數(shù)可寫為</p><p>　　k=[k1:k2] %k1,k2為時(shí)間序列的起始及終止時(shí)

7、間序號(hào)</p><p>　　fk=[(k-k0)>=0] %k0為（k）在時(shí)間軸上的位移量</p><p>　　stem(k,fk)</p><p>　　2. 離散系統(tǒng)的時(shí)域響應(yīng)</p><p>　　MATLAB提供了求離散系統(tǒng)響應(yīng)的專用函數(shù)filter()</p><p>　　filter函數(shù)能求出差

8、分方程描述的離散系統(tǒng)在指定時(shí)間范圍內(nèi)的輸入序列所產(chǎn)生的響應(yīng)序列的數(shù)值解。其調(diào)用格式如下：</p><p>　　filter(b,a,x)</p><p>　　其中，b和a是由描述系統(tǒng)的差分方程的系數(shù)決定的、表示離散系統(tǒng)的兩個(gè)行向量；x是包含輸入序列非零樣點(diǎn)的行向量。上述命令將求出系統(tǒng)在與x的取樣時(shí)間點(diǎn)相同的輸出序列樣值。</p><p>　　3.離散時(shí)間序列的卷積&

9、lt;/p><p>　　設(shè)序列f1(k)在～范圍內(nèi)非零，f2(k)在～范圍內(nèi)非零，則f1(k)時(shí)域?qū)挾葹長1=-+1，f2(k)的時(shí)域?qū)挾葹長2=-+1。由卷積和的定義可得，序列y(k)=f1(k)*f2(k)的時(shí)域?qū)挾葹長=L1+L2-1，且只在(+)～(++L1+L2-2)范圍內(nèi)非零。因此，對(duì)于f1(k)和f2(k)均為有限期間非零的情況，只要計(jì)算序列y(k)在 (+)～(++L1+L2-2)范圍內(nèi)的序列值，便可

10、以表征序列y(k)。</p><p>　　MATLAB的conv()函數(shù)可以用來快速求出兩個(gè)離散序列的卷積和。conv()函數(shù)的調(diào)用格式為</p><p>　　y=conv(f1,f2)</p><p>　　式中f1和f2為待卷積兩序列的向量表示；y為卷積結(jié)果。</p><p><b>　　四.程序?qū)崿F(xiàn)</b></

11、p><p>　　1.畫出單位序列(k)在-3k3區(qū)間的波形。</p><p>　　程序：function impseq(k1,k2,k0)</p><p>　　k1=-3;k2=3;k0=0;</p><p>　　k=[k1:k2];</p><p>　　fk=[(k==k0)];</p><p>

12、　　stem(k,fk,'fill')</p><p>　　axis([k1,k2,0,1.5])</p><p><b>　　圖示：</b></p><p>　　2.畫出單位階躍序列(k)在-2k5區(qū)間的波形。</p><p>　　程序：k1=-2;k2=5;k0=0;</p><p

13、>　　k=[k1:k2];</p><p>　　fk=[(k>=k0)];</p><p>　　stem(k,fk,'fill')</p><p>　　axis([k1,k2,0,1.5])</p><p>　　title('單位階躍')</p><p><b>　

14、　圖示：</b></p><p>　　3.利用函數(shù)filter()，把描述離散系統(tǒng)的差分方程為：</p><p>　　y(k)-0.25y(k-1)+0.5y(k-2)=f(k)+f(k-1)</p><p>　　且該系統(tǒng)輸入序列為f(k)=(k),</p><p>　　用MATLAB 求系統(tǒng)的單位序列響應(yīng) 和系統(tǒng)的響應(yīng)y(k)。

15、</p><p>　　程序如下：function y=filter(b,a,x)</p><p>　　k=-100:100;</p><p>　　uk=[zeros(1,100),ones(1,101)];</p><p><b>　　b=[1,1];</b></p><p>　　a=[1,-0.

16、25,0.5];</p><p>　　x=((1/2).^k).*uk;</p><p>　　y=filter(b,a,x);</p><p>　　h=impz(b,a,k);</p><p>　　stem(k,h);</p><p>　　xlabel('k');</p><p>

17、;　　title('單位序列響應(yīng)')</p><p>　　figure(2);stem(k,y);</p><p>　　title('系統(tǒng)響應(yīng)')</p><p>　　xlabel('k')</p><p><b>　　圖示：</b></p><p>

18、;　　4.已知離散信號(hào)f1(k)和f2(k)，求y(k)=f1(k)*f2(k)</p><p>　　f1(k)= f2(k)=</p><p>　　用前面介紹的conv()函數(shù)來完成。并說明f1(k)和f2(k)的時(shí)域?qū)挾扰c序列y(k)的時(shí)域?qū)挾鹊年P(guān)系。</p><p>　　程序：function y=conv(f1,f2)</p><

19、;p>　　f1=[-1,0,1,2,3];</p><p>　　f2=[-1,0,1,2,3,4];</p><p>　　y=conv(f1,f2);</p><p>　　N=length(y)</p><p>　　stem(0:N-1,y)</p><p>　　結(jié)果： f1和f2的時(shí)域?qū)挾确謩e為5，6，y的時(shí)域

20、寬度為5+6-1=10.</p><p><b>　　圖示：</b></p><p><b>　　五、設(shè)計(jì)心得</b></p><p>　　通過本次設(shè)計(jì)，我感覺受益匪淺。首先，通過自己親手查閱資料、編程、調(diào)試，培養(yǎng)了我的實(shí)踐能力和獨(dú)立分析問題、解決問題的能力，進(jìn)一步鞏固和加深了對(duì)課程內(nèi)容的理解和綜合運(yùn)用。</p>

21、;<p>　　同時(shí)，我對(duì)于MATLAB的應(yīng)用更加熟練了，通過運(yùn)用本軟件，我對(duì)于離散系統(tǒng)的時(shí)域的特點(diǎn)掌握的更加牢固了。</p><p><b>　　六、參考文獻(xiàn)</b></p><p>　　【1】梁虹. 信號(hào)與系統(tǒng)分析及MATLAB實(shí)現(xiàn) 電子工業(yè)出版社 2002.02第12章</p><p>　　【2】王立寧等.MATLAB與通信