matlab做相干卷积信号处理

t=1:10;          %采集10个数据
ts=t*pi/4/10;    %在0~pi/4区间等间距取10个值
y1=exp(-ts/2).*sin(2*ts);%y1函数
y2=sin(2*ts);    %y2函数
subplot(411);
stem(ts,y1);     %做出y1图形
title('y1序列');
subplot(412);
stem(ts,y2);     %做出y2图像
title('y2序列');
z1=conv(y1,y2);  %y1与y2线性卷积
subplot(413);
stem(z1);        %做出线性卷积图形
title('y1与y2线性卷积');
axis([0 20 0.95*min(z1) 1.1*max(z1)]);
z2=cconv(y1,y2,10);%y1与y2循环卷积
subplot(414);
stem(z2);       %做出循环卷积图形
title('y1与y2循环卷积');
axis([0 10 0.95*min(z2) 1.1*max(z2)]);

t=1:10;          %采集10个数据
ts=t*pi/4/10;    %在0~pi/4区间等间距取10个值
y1=exp(-ts/2).*sin(2*ts);%y1函数
y2=sin(2*ts);    %y2函数
x1=xcorr(y1,y2); %y1与y2线性相关
subplot(311);
stem(x1);        %画出线相关图
title('y1与y2线性相关');
V=a3(y2);        %调用a3循环左移函数