MATLAB实验报告（信号与线性系统分析） - 图文(2)

4. 绘制ε(k-3)，-1≤k≤5及ε（k+1），-5≤k≤2。 解：M文件：

结果：（1）： jyxulie(-1,5,-3)

（2）jyxulie(-5,2,1)

5. 考虑下面3个信号：

f1(n)=cos(2πn/N)+2cos(3πn/N) f2(n)=2cos(2n/N)+2cos(3n/N) f3(n)=cos(2πn/N)+3cos(5πn/N)

假设对每个信号N=6。试确定上述信号是否是周期的。如果是，则确定信号的周期，并画图表示出该信号的两个周期；如果不是周期的，在[0，4N]的范围内画出该信号，并说明原因。

解：代码：k=0:40;

subplot(3,1,1)

stem(k,cos(2*pi*k/6)+2*cos(3*pi*k/6),'filled') title('cos(2*pi*k/6)+2*cos(3*pi*k/6)')

subplot(3,1,2)

stem(k,2*cos(2*k/6)+2*cos(3*k/6),'filled') title('2*cos(2*k/6)+2*cos(3*k/6))') subplot(3,1,3)

stem(k,cos(2*pi*k/6)+3*cos(5*pi*k/6),'filled') title('cos(2*pi*k/6)+3*cos(5*pi*k/6)') 结果：

可得1，3是周期信号 2是非周期信号 【四】 思考题

1. 观察例2-1的执行结果，当把时间间隔p取得更小（例如为0.02）时，观察执行结果，比较两结果用何不同，为什么？

2. 例2-3中的hold on和hole off命令的作用是什么？

实验三 信号的时域运算、时域变换及MATLAB实现

【一】实验目的

1. 掌握信号时域运算的MATLAB实现方法 2. 掌握信号时域变换的MATLAB实现方法

【二】实验原理

信号的时域运算包括信号的相加、相乘，信号的时域变换包括信号的平移、反折、倒相及信号的尺度变换。

一. 连续信号的时域运算与时域变换

MATLAB可以有两种方法来表示连续信号。用这两种方法均可实现连续信号的时域运算和变换，但用符号运算的方法则较为简便。 二. 离散序列的时域运算及时域变化

对于离散序列来说，序列相加、相乘是将两序列对应时间序号的值逐项相加或相乘，平移、反折及倒相变化与连续信号的定义完全相同，这里就不再赘述。但需要注意，与连续信号不同的是，在MATLAB中，离散序列的时域运算和变换不能用符号运算来实现，而必须用向量表示的方法，即在MATLAB中离散序列的相加、相乘需表示成两个向量的相加、相乘，因而参加运算的两序列向量必须具有相同的维数。

【三】上机练习

将实验原理中提到的例子在计算机上全部练习一遍；

?1?4(t?1)?4?t?0?1. 已知信号

f(t)??10?t?2，画出f(t)及f(?2t?4)的波形；

?0其它??解：代码：syms t f

f=sym('1/4*(t+1)*(heaviside(t+4)-heaviside(t))+(heaviside(t-2)-heavis

ide(t))')

subplot(2,1,1),ezplot(f,[-4,5]) y=subs(f,t,-2*t-4)

subplot(2,1,2),ezplot(y)

结果：

2. 已知两个连续信号

f1(t)??t,0?t?1,f2(t)?sin(2?t)，用MATLAB绘

出下列信号的波形： 1）f3(t)3）

?f1(?t)?f1(t) 2）f4(t)??[f1(?t)?f1(t)]

f5(t)?f2(t)?f3(t) 4）f6(t)?f1(t)?f2(t)

解：代码：syms t f

f1= -t*(heaviside(t)- heaviside(t+1)) f2=sin(2*pi*t) y=subs(f1,t,-t) f3=f1+y

subplot(2,2,1), ezplot(f3) ylabel('f3'); f4=-1*f3

subplot(2,2,2), ezplot(f4) ylabel('f4'); f5=f2*f3

subplot(2,2,3), ezplot(f5) ylabel('f5'); f6=f1*f2

subplot(2,2,4), ezplot(f6) ylabel('f6');

结果：