序列卷积和运算及其可视化呈现

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