《数字信号处理》用matlab实现序列的卷积与相关

实验一：用convwthn求解x(n)=[1,2,3,-1,-2],nx=[-1,3]与h(n)=[2,2,1,-1,4,-2],nh=[-3,2]的卷积。

1.实验代码

%用conv函数编出有位置矢量ny的输出y(n)的convwthnHS.m程序

function [y,ny]=convwthnHS(x,nx,h,nh)
y=conv(x,h);
ny1=nx(1)+nh(1);
ny2=nx(end)+nh(end);    %end表示最后一个下标
ny=[ny1:ny2];
代码2
%用convwthn求解x(n)=[1,2,3,-1,-2],nx=[-1,3],
%与h(n)=[2,2,1,-1,4,-2],nh=[-3,2]的卷积。

x=[1,2,3,-1,-2],nx=-1:3;
h=[2,2,1,-1,4,-2];nh=-3:2;
[y,ny]=convwthnHS(x,nx,h,nh)
stem(ny,y,'b.');
title('y(n)的图形');
xlabel('n');
ylabel('y(n)');grid;

2.实验结果

实验二：计算下列两个序列的互相关：rxy(n),x(n)=[2,1,3,2,1,5,1],y(n)=[2,1,3,4],采用rxy(n)=x(m)*y(-m)公式，即用卷积代替相关。

1.实验代码

%用convwthn函数求互相关

x=[2,1,3,2,1,5,1],y=[2,1,3,4]