【正文】 加上一些附加的零值。)。axis([4 16 0 3])。其調(diào)用格式如下：y=impulse(sys,t)y=step(sys,t)其中：t表示計算系統(tǒng)響應(yīng)的時間抽樣點向量，sys表示線性時域系統(tǒng)模型。title(39。兩個序列的加減乘除是對應(yīng)離散樣點值的加減乘除，因此，可通過MATLAB點乘和點除、序列移位和反折來實現(xiàn)，與連續(xù)時間信號處理方法基本一樣。)。stem(n3,x,39。)。n39。)。stem(nx,x,39。xlabel(39。y(n)=x(n)*h(n)39。grid on。subplot(2,2,2)。axis([4 4 ])。axis([10 10 *pi])。a=[1 1 ]。filter函數(shù)求出的單位沖激響應(yīng)39。axis([5 95 ])。 n2=0:length(hn2)1。.39。subplot(2,1,1)stem(n,hn),xlabel(39。第一天，在網(wǎng)上尋找程序的時候，發(fā)現(xiàn)我的這個題目在網(wǎng)上找程序十分的困難，別人的程序都可以直接找到做好的，我的只能一點一點找，再自己一點一點改，這讓我感到十分的不公平，不過還是熬過來了?？芍罆r不變系統(tǒng)的變換關(guān)系不隨時間變化而變化，即系統(tǒng)的輸出隨輸入的移位而相應(yīng)移位但形狀不變。a=[1]。xn=[1,zeros(1,30)]。hn1=filter(B1,A1,xn)。grid。h(n)39。MATLAB程序如下：clear all。*W)。)。k=N:1:N。t1=4:dt:4。xlabel(39。stem(nh,h,39。y=conv(x,h)。impz(b,a,30)。fill39。grid on。xlabel(39。n4=n1。title(39。)。axis([4 16 0 3])。)。)。MATLAB求離散時間信號的卷積和命令為conv，其調(diào)用格式如下：y=conv(x,h)其中：x與h表示離散時間信號值的向量；y為卷積結(jié)果。線性時不變系統(tǒng)既滿足疊加原理又具有時不變特性，它可以用單位脈沖響應(yīng)來表示。f(t)系統(tǒng)y(t)，Af(t)系統(tǒng)Ay(t)（2）、疊加性若激勵f1(t)與f2(t)產(chǎn)生的響應(yīng)分別為y1(t)， y2(t)，則激勵f1(t)+f2(t)產(chǎn)生的響應(yīng)即為y1(t)+y2(t)，此性質(zhì)稱為疊加性。其實現(xiàn)的MATLAB程序代碼如下：clear all。title(39。subplot(3,1,3)。已知一連續(xù)時域系統(tǒng)的微分方程如下：下面用MATLAB命令繪出0 ≤t ≤4 范圍內(nèi)系統(tǒng)的沖激響應(yīng) h(t)和階躍響應(yīng) g(t)。沖激響應(yīng)39。下面列舉例子來說明。grid on。fill39。程序在運行過程中，其源代碼如下：function y=uDT(n)y=n=0。)。其仿真波形如下圖：圖4 單位取樣響應(yīng)(3)卷積和運算由于系統(tǒng)的零狀態(tài)響應(yīng)是激勵與系統(tǒng)的單位取樣響應(yīng)的卷積，因此，卷積運算在離散時間信號處理領(lǐng)域被廣泛應(yīng)用。fill39。n39。)。axis([4 4 ])。plot(W,abs(Fw))。xlabel(39。xlabel(39。n1=0:89。)。set(gcf,39。subplot(2,1,1)stem(n2,hn2,39。), title(39。n39。)由線性系統(tǒng)的仿真波形可以看出線性系統(tǒng)的輸入與輸出滿足線性疊加定理，由非線性系統(tǒng)的仿真波形可以看出它不滿足線性疊加定理。系統(tǒng)函數(shù)：H(z)=z/()b=[1,]。A3=[,]。xn=[1,zeros(1,30)]。)。ylabel(39。系統(tǒng)由差分方程描述： 分別用filter函數(shù)和impz函數(shù)繪出其單位沖激響應(yīng)h(n)的線圖。Fsw=Ts*fst*exp(j*t239。:39。N=500。dt=。grid on。subplot(3,1,2)。h=.^nh.*(uDT(nh)uDT(nh8))。n=0:30。stem(n,y,39。)。axis([15 15 0 1])。