MATLAB实验--连续时间信号分析.doc
文本预览下载声明
实验 1 连续时间信号分析
一、实验目的
（一）掌握使用Matlab表示连续时间信号
1、学会运用Matlab表示常用连续时间信号的方法
2、观察并熟悉常用信号的波形和特性
（二）掌握使用Matlab进行连续时间信号的相关运算
1、学会运用Matlab进行连续时间信号的时移、反褶和尺度变换
2、学会运用Matlab进行连续时间信号微分、积分运算
3、学会运用Matlab进行连续时间信号相加、相乘运算
4、学会运用Matlab进行连续时间信号卷积运算
二、实验条件
Win7系统，MATLAB R2015a
三、实验内容
1、利用Matlab命令画出下列连续信号的波形图。
（1）
（2）
（3）
2、利用Matlab命令画出复信号的实部、虚部、模和辐角。
绝对值为2 ang纵轴上下限为-pi ~ pi
3、已知信号的波形如下图所示：
试用Matlab命令画出的波形图。
4、使用微分命令求关于变量x的一阶导数；使用积分命令计算不定积分 ，定积分。
5、已知，使用命令画出两信号和及两信号乘积的波形图。其中，
6、用Matlab命令绘出下列信号的卷积积分的时域波形图。
四、实验结论和讨论
记录几个重点：
1、自变量范围问题，matlab plot绘图为描点绘图，t = x1:0.01:x2
就是将t在x1到x2之间按0.01为间隔分成若干个数值。即t还是离散的，只是数值很多。
2、图像显示时的axis以及间隔要合适，才能作出能反应函数实际情况的图形。
3、卷积函数conv在实验室环境下不好用，不知为何
五、实验思考
在进行门函数u（t）与其他函数相乘时，会遇到错误Inner matrix dimensions must agree.
经查，原因是参与运算的矩阵维数不匹配，matlab中sin（x），x，e^x的数据形式是不同的。在workspace中可以观察到。解决方案是将 * 换成 .* 即对两个矩阵对应维数进行相乘。
7
代码：
A = 2; w = 3; phi = pi/4;
t = 0:0.01:3;
f = A*cos(w*t + phi);
plot(t, f), grid on;
axis([0, 3, -2, 2]);
title(2cos(3t + PI/4))
代码：
t = -5:0.01:5;
f = (2 - exp(-t)).*uCT(t);
plot(t, f), grid on;
axis([-5, 5, 0, 2.5]);
title((2-e^(-t))u(t))
代码：
A = 1; w = pi; phi = 0;
t = -3:0.01:3;
f = (1 + A*cos(w*t + phi)).*(uCT(t) - uCT(t - 2));
plot(t, f), grid on;
axis([-3, 3, 0, 2]);
title([1+cos(t*pi)][u(t)-u(t-2)])
代码：
% Ke^( a + bi)t
K = 2;a = 0;
t = -2:0.01:3;b = (t + pi/4);
f = K * exp(a + 1i*b);
subplot(2, 2, 1);plot(t, real(f)), grid on;
axis([-2, 3, -2, 2]);title(Real);
subplot(2, 2, 2);plot(t, imag(f)), grid on;
axis([-2, 3, -2, 2]);title(imag);
subplot(2, 2, 3);plot(t, abs(f)), grid on;
axis([-2, 3, 0, 3]);title(Abs);
subplot(2, 2, 4);plot(t, angle(f)), grid on;
axis([-2, 3, -pi, pi]);title(Ang);
function u = funct(t)
u = uCT(t + 2) - uCT(t) + (-t+1).*(uCT(t) - uCT(t - 1));
end
代码：
t = -2:0.01:4;
ft1 = funct(t - 2);
ft2 = funct(3 * t);
ft3 = funct(-t);
ft4 = funct(-3 * t - 2);
subplot(2, 2, 1);plot(t, ft1), grid on;
axis([-2, 4, -1, 2]);title(f(t - 2));
subplot(2, 2, 2);plot(t, ft2), grid on;
axis([-2, 4, -1, 2]);title(f(3t))
显示全部