myconv.m

来自「卷积实验代码 实现多种功能 包括线性卷积和圆周卷积的动态演示」· M 代码 · 共 45 行

%x1=[1,1,1,1,0,0,1,1,1,1,0,0];x2=[0,1,2,1,0,0,0,1,2,1,0,0];n0=8;
function yn0=myconv( x1,x2,n0 )
%y(n0)
%y(n)=x1(n)*x2(n)

n1=0:length(x1)-1;
n2=0:length(x2)-1;

subplot(2,3,1);
stem(n1,x1);
xlabel('n');
ylabel('x1(n)');
title('x1(n)');


subplot(2,3,2);
stem(n2,x2);
xlabel('n');
ylabel('x2(n)');
title('x2(n)');

[x2,n2] = sigfold(x2,n2);
subplot(2,3,3);
stem(n2,x2);
xlabel('n');
ylabel('x2(-n)');
title('x2(-n)');

[x2,n2]=sigshift(x2,n2,n0);
subplot(2,3,4);
stem(n2,x2);
xlabel('n');
string1=['x2(-n+',num2str(n0),')'];
ylabel(string1);
title(string1);

[y,n]=sigmult(x1,n1,x2,n2);
subplot(2,3,5);
stem(n,y);
xlabel('n');
string2=['x1(n).x2(-n+',num2str(n0),')'];
ylabel(string2);
title(string2);

yn0=y*ones(length(n),1);