**例2.4 **分別調用conv()函數和自己編寫程序計算f1(t)=cos(t)[u(t)-u(t-10)]和f2(t)=[exp(t)+exp(2t)][u(t)-u(t-10)]的卷積，比較兩種計算的結果。

T=0.1; %%%時間步長，這里可以理解為抽樣函數，

對函數進行抽樣，T越小函數越精細

t1=0:T:10; %%%時間序列，建立一個數組

f1=cos(t1); %%%信號f1

t2=t1;

f2=exp(t2)+exp(-2*t2);%%%信號f2

f=T*conv(f1,f2); %%%計算卷積，注意這里乘步長

k0=t1(1)+t2(1); %%%卷積輸出序列的起始

k3=length(f1)+length(f2)-2;

t=k0:T:(k0+T*k3); %%%卷積結果對應的時間向subplot(3,1,1); %%%繪制信號f1，subplot函數的用法是一次可以繪制多個函數圖片。subplot(m,n,p),m是輸出圖片的行數，n是輸出函數的列數，p是該函數在第幾行出現

plot(t1,f1,'linewidth',2);

title('f1(t)');

subplot(3,1,2); %%%繪制信號f2

plot(t2,f2,'linewidth',2);

title('f2(t)');

subplot(3,1,3); %%%繪制卷積結果

plot(t,f,'linewidth',2);

title('convolution of f1(t)and f2(t)');

MATLAB輸出的卷積結果

用自己編寫函數的方法

for k=1:lf1+lf2-1

y(k)=0;                                        %%%y賦初始值

for ii=max(1,k-(lf2-1)):min(k,lf1)

    y(k)=y(k)+f1(ii)*f2(k-ii+1);      %%%信號相乘和求和

end

yzsappr(k)=T*y(k);                    %%%用乘和加運算來近似積分運算

end

自己編寫函數其實就是實現matlab中conv()函數。這個函數的功能類似于圖形法求函數卷積。拿多項式求卷積為例，假設p=[1,2,3],q=[1,1]。那么其實conv（p，q）的結果就是對應的多項式乘積。

換個函數求卷積試試

按照書上的例題折騰完，換一個函數求一下卷積試一下。

隨便弄了個二次函數上去，小編學藝不精，還要繼續學習。