function?[ynny]?=?my_convu（xnnxhnnh）
%UNtitleD2?此處顯示有關此函數的摘要
%?整個過程分為四個部分：分別為xn的顯示、hn的顯示、hn的翻轉移位、輸出yn的顯示

%???此處顯示詳細說明
N=length（hn）;%計算序列長度
M=length（xn）;
a1=min（nx-N-（nh+N-1））;%計算擴列范圍
m1=min（a1nh）;
a2=max（（nx+M-1）+N-1-nh）;
m2=max（a2nh+N-1）;
n1=nx+nh;%計算輸出序列的位置
n2=（nx+M-1）+（nh+N-1）;
m=m1:m2;%擴列的最終范圍

%圖一：擴列顯示第一個序列
xn=[zeros（1nx-m1）?xn?zeros（1m2-（nx+M-1））];%輸入xn的擴列
subplot（221）;stem（mxn‘r‘‘f‘）;grid?on;%畫圖
ylabel（‘x（n）‘）;xlabel（‘n‘）;title（‘輸入序列x（n）‘）;
%?pause（1）;

%圖二：擴列顯示第二個序列
hn1=[zeros（1nh-m1）?hn?zeros（1m2-（nh+N-1））];%hn1為擴列后的輸入序列hn
subplot（222）;stem（mhn1‘r‘‘f‘）;grid?on;%畫圖
ylabel（‘h（n）‘）xlabel（‘n‘）title（‘輸入序列h（n）‘）
%?pause（1）;

%圖三：先翻轉第二個序列并顯示，然后移位
hn2=fliplr（hn1）;%hn2為擴列后的hn1的翻轉序列
s