用MATLAB语言编写一个描述时域卷积运算过程的演示软件

实验一：用MATLAB语言编写一个描述时域卷积运算过程的演示软

件

1.1实验目的:

深入理解卷积的运算过程。

1.2实验内容:

以矩形方波的自卷积计算为例，用MATLAB语言在计算机屏幕上演示出卷积运算过程中变量t的替换、反折、平移和积分运算，并画出卷积运算结果曲线。

1.3实验过程：

在信号与系统时域分析中，对于线性时不变系统，系统的零状态响应y(t)或y[n]就是系统的激励x(t)或x[n]与系统的单位冲激响应h(t)或h[n]的的卷积。

在分析线性系统过程中知道卷积运算将输入信号、输出信号以及冲激响应三者之间联系起来,主要有时域和频域两种计算方法。在时域法中涉及到信号的平移、反转运算以及冲激响应的概念;而在频域法中,涉及到信号的傅里叶变换和频率响应的概念。可以看出,卷积运算将各种物理概念、分析方法以及知识点有机地衔接起来,是培养学生综合能力的一个重要方面

T = 1;

A = 1;

left = -3/2;

right = 3/2;

delta = 1e-3;

t = left:delta:right;

x = zeros(size(t));

head = round((-T/2-left)/delta);

tail = round((T/2-left)/delta);

for i=head:tail

x(i) = A;

end

y = x;

figure(1);

title('hanshu1');

plot(t, x);

xlabel('tau');

ylabel('x(t)');

figure(2);

plot(t, y);

xlabel('tau');

ylabel('y(t)');

z = zeros(size(x));

figure(3);

axis([-3*T/2 3*T/2 0 A*A*2]);

for s = -T:delta:T

j = round((s-left)/delta);

if(s<0)

z(j) = A*A*(s+T);

else

z(j) = A*A*(T-s);

end

figure(3);

axis([-3*T/2 3*T/2 0 A*A*2]);

plot(t(j), z(j));

hold on;

pause(1e-3);

End

1.4实习总结：

通过本次实习，我基本熟悉了MATLAB语言，对matlab语言的使用也有了一定的认识。通过老师和同学们的帮助基本可以用MATLAB语言在计算机屏幕上演示出卷积运算过程中变量t的替换、反折、平移和积分运算，并能画出卷积运算结果曲线。