【文章內(nèi)容簡介】 零的方式成為同維數(shù)的二序列向量，因而在調(diào)用該函數(shù)時，要進行相加運算的二序列向量維數(shù)可以不同：function[f,k]=lsxj(f1,f2,k1,k2) %實現(xiàn)f(k)=f1(k)+f2(k)， ffkk2是參加運算的二離散序列及其對應的時間序列向量，f和k為返回的和序列及其的時間序列向量k=min(min(k1),min(k2)):max(max(k1),max(k2))。 %構(gòu)造和序列的長度s1=zeros(1,length(k))。s2=s1。 %初始化新向量s1(find((k=min(k1))amp。(k=max(k1))==1))=f1。 %將f1中在和序列范圍內(nèi)但又無定義的點賦值為零s2(find((k=min(k2))amp。(k=max(k2))==1))=f2。 %將f2中在和序列范圍內(nèi)但又無定義的點賦值為零f=s1+s2。 %兩長度相等序列求和stem(k,f,’filled’) axis([(min(min(k1),min(k2))1),(max(max(k1),max(k2))+1),(min(f)),(max(f)+)]) %坐標軸顯示范圍三、上機練習將實驗原理中提到的例子在計算機上全部練習一遍；已知信號，畫出的波形；解。3. 已知兩個連續(xù)信號，用MATLAB繪出下列信號的波形：1） 2）3） 4）解。離散序列，用MATLAB繪出下列序列的波形。1） 2）3） 4）解：（1）（2）（3）(4)四、思考題 在對信號進行平移、反折和尺度變換時，運算順序?qū)Y(jié)果是否有影響？在運算中應該注意什么？實驗四 離散序列卷積和的MATLAB實現(xiàn)一、實驗目的1. 掌握卷積和的計算機編程方法，利用MATLAB實現(xiàn)兩個離散序列的卷積和；2. 利用卷積和求離散系統(tǒng)的響應，觀察、分析系統(tǒng)的時域特性。二、實驗原理 兩個離散序列卷積和的定義為 （）定義式可以看作是：將序列的時間軸反并將其移位個樣本，然后將移位后的乘以并在上將所得到的乘積序列相加。這種說法直接來自離散時間系統(tǒng)的線性和時不變性質(zhì)。信號可以看成是由延時和加權(quán)脈沖的線性疊加所構(gòu)成，因為一個LTI系統(tǒng)能夠用它對單個脈沖的響應來表示，那么一個LTI系統(tǒng)的輸出就應該相對于系統(tǒng)對構(gòu)成的每一個延時和加權(quán)脈沖的疊加。在數(shù)學上，這個結(jié)果就是卷積和。在離散信號與系統(tǒng)的分析過程中，我們有兩個與卷積和相關的重要結(jié)論，這就是：1），即離散序列可分解為一系列幅度由決定的單位序列及其移位序列之和。2）對線性時不變系統(tǒng)，設其輸入序列為，單位響應為，其零狀態(tài)響應為，則有：可見，離散序列卷積和的計算對進行離散信號與系統(tǒng)的分析具有非常重要的意義。設序列在區(qū)間~非零，在區(qū)間~非零，則的時域?qū)挾葹?，的時域?qū)挾葹?。由卷積和的定義可得，序列的時域?qū)挾葹?，且只在區(qū)間（+） ~++（） 2非零。因此，對于和均為有限期間非零的情況，我們只需要計算序列在區(qū)間（+）~++（） 2的序列值，便可以表征整個序列。 MATLAB的conv（ ）函數(shù)可以幫助我們快速求出兩個離散序列的卷積和。conv（ ）函數(shù)的調(diào)用格式為：f=conv (f1,f2 ) 其中為包含序列的非零樣值點的行向量，為包含序列的非零樣值點的行向量，向量f則返回序列的所有非零樣值點行向量。 下面是利用MATLAB計算兩離散序列卷積和的實用函數(shù)dconv()程序，該程序在計算出卷積和的同時，還繪出序列、及的時域波形圖，并返回的非零樣值點的對應向量。function [f,k]=dconv(f1,f2,k1,k2)%The function of putef=f1*f2%f 卷積和序列f(k)對應的非零樣值向量%k 序列f(k)的對應序號向量%f1 序列f1(k)非零樣值向量%f2 序列f2(k)非零樣值向量%k1 序列f1(k)的對應序號向量%k2 序列f2(k)的對應序號向量f=conv(f1,f2) %計算序列f1與f2的卷積和fk0=k1(1)+k2(1)。 %計算序列f非零樣值的起點位置k3=length(f1)+length(f2)2。 %計算卷積和f非零樣值的寬度k=k0:k0+k3。 %確定卷積和f非零樣值的序號向量subplot(2,2,1)。stem(k1,f1) %在子圖1繪序列f1(k)的波形title(39。f1(k)39。)xlabel(39。k39。)ylabel(39。f1(k)39。)subplot(2,2,2)。stem(k2,f2) %在子圖2繪序列f2(k)的波形title(39。f1(k)39。)xlabel(39。k39。)ylabel(39。f2(k)39。)subplot(2,2,3)。stem(k,f)。 %在子圖3繪序列f(k)的波形title(39。f1(k)與f2(k)的卷積和f(k)39。)xlabel(39。k39。)ylabel(39。f(k)39。)h=get(gca,39。position39。)。h(3)=*h(3)。set(gca,39。position39。