用matlab怎样表示褶积,用MATLAB实现信号的时域卷积.doc

用MATLAB实现信号的时域卷积

用MATLAB实现信号的时域卷积

摘要：

利用抽样方法绘制连续信号和离散信号。编制函数计算连续时间信号，离散时间信号的卷积。

关键字：

MATLAB；离散卷积；连续卷积

实验原理：

卷积运算在信号与系统分析中具有非常重要的意义，是信号与系统分析的基本方法之一。线性是不变系统的零状态响应就可以通过求解驶入激励信号与系统单位冲击响应的卷积运算获得。

在连续时间系统中，卷积运算成为卷积积分，表示为

在离散时间系统中，卷积运算称为卷积和，表示为：

实验任务及实现方法：

⑴ 编制计算离散时间序列卷积核的M函数,该函数课计算两离散序列和的卷积和。程序计算出卷积和的同时，还绘制出序列、和的时域波形图，并返回的非零样值点的对应向量。

任务(1)程序：

function dconv(x1,x2,y1,y2)

%任意两序列卷积

x1=input('请输入x1=');

x2=input('请输入x2=');

y1=input('请输入y1=');

y2=input('请输入y2=');

subplot(131),stem(x1,y1),grid on,title('y1'),axis([(min(x1)-1),(max(x1)+1),(min(y1)-1),(max(y1)+1)])

subplot(132),stem(x2,y2),grid on,title('y2'),axis([(min(x2)-1),(max(x2)+1),(min(y2)-1),(max(y2)+1)])

%_____________________________________________

xmin1=min(x1);

xmax1=max(x1);

xmin2=min(x2);

xmax2=max(x2);

t=(xmax2+xmax1)-(xmin2+xmin1)+1;

xx=0:(t-1);

ft=conv(y1,y2)

subplot(133),stem(xx,ft),grid on,title('*')

axis([(min(xx)+0.5),(max(xx)+0.5),(min(ft)-0.5),(max(ft)+0.5)])

⑵ 调用来完成下面两离散序列的卷积和运算，并绘制图形：

任务(2)程序：

function f=ddconv(m1,m2)

%-----------------------------

%离散序列1

t1=m1-5;

t2=m1+5;

t=t1:t2; %横坐标

n=length(t);

tt=t1:m1; %起始到跳变

n1=length(tt);

f1=[zeros(1,n1-1),ones(1,n-n1+1)];

%subplot(141),stem(t,f1),grid on,title('阶跃函数序列 u(n-m1) '),axis([t1 t2 -0.2 1.5])

%-------------------------------

%离散序列2

p1=m1-5;

p2=m1+5;

p=p1:p2; %横坐标

z=length(p);

pp=p1:m2; %起始到跳变

z1=length(pp);

f2=[zeros(1,z1-1),ones(1,z-z1+1)];

%subplot(142),stem(p,f2),grid on,title('离散阶跃序列 u(n-m2) '),axis([p1 p2 -0.2 1.5])

%------------------------------------

%两离散序列相减

k=min(min(t),min(p)):max(max(t),max(p)); %共同横坐标

s1=zeros(1,length(k));s2=s1;

s1(find((k>=min(t))&(k<=max(t))==1))=f1;

s2(find((k>=min(p))&(k<=max(p))==1))=f2;

f=s1-s2;

subplot(131),stem(k,f),grid on,title('x(n)=u(n)-u(n-4)')

axis([(min(min(t),min(p))-1),(max(max(t),max(p))+1),(min(f)-0.5),(max(f)+0.5)])

subplot(132),stem(k,f),grid on,title('h(n)=u(n)-u(n-4)')

axis([(min(min(t),min(p))-1),(max(max(t),max(p))+1),(min(f)-0.5),(max(f)+0.5)])

%----------------------------

xmin1=min(k);

xm