mycirconv.m

卷积实验代码 实现多种功能 包括线性卷积和圆周卷积的动态演示

%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;N=23;
function yn0=mycirconv( x1,x2,n0,N )
%y(n0)
%y(n)=x1(n)(N)x2(n)
%N=length(x2);

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)');

n2=[0:1:N-1];
x2=[x2 zeros(1,N-length(x2))];
x2=x2(mod(-n2,N)+1);
subplot(2,3,3);
stem(n2,x2);
xlabel('n');
string=['x2((-n))',num2str(N)];
ylabel(string);
title(string);

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

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

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