信号与系统实训报告

实验一信号的可视化及时域运算

－、实验目的

1. 掌握用M ATLAB 实现连续与离散信号可视化的方法。

2. 掌握用M ATLAB 实现信号的时域运算及变换的方法。

二、实验内容

1、7.1 (1) (2)

2、7.2 (1) (3)

3、7.5 (2)

三、实验原理

主要写：

1、本实验内容与理论课的那些知识点相对应，表述这些知识点

2、在编写Matlab程序时用到的函数及主要功能

P为时间采样间隔，pi为MATLAB的预定义变量π，plot可以绘出信号的时域波形。

四、实验记录（程序、波形等）

7.1（1）

% dm07101

p=0.001;

t=-pi:p:pi

f=4*sin(2*pi*t-pi/4)

plot(t,f)

title('f(t)=4sin(2*pi*t-pi/4')

xlabel('t')

axis([-pi,pi,-5,5])

波形

-3

-2

-1

01

2

3

-5-4-3-2-101234

5f(t)=4sin(2*pi*t-pi/4

t

7.1(2)

%dm07101 p=0.001; t=-pi:p:pi

f=(1-exp(-2*t)).*heaviside 波形

7.2（1） syms t

f=sym('4*sin(2*pi*t-pi/4)')

ezplot(f,[-pi,pi])

set(gcf,'color','w')

图形

4 sin(2 π t-π/4)

t

7.2（3）

syms t

f=sym('Heaviside(t+2)-Heaviside(t-2)') ezplot(f,[-10,10])

set(gcf,'color','w')

Heaviside(t+2)-Heaviside(t-2)

-10-8-6-4-20246810

t

7.5（2）

function fexp(d,w,t1,t2,a)

t=t1:0.01:t2;

f=a*exp((d+i*w)*t);

fr=real(f);

fi=imag(f);

fa=abs(f);

fn=angle(f);

subplot(2,2,1)

plot(t,fr)

axis([t1,t2,-(max(fa)+0.5),max(fa)+0.5]),

title('ʵ²¿');

subplot(2,2,2)

plot(t,fi)

axis([t1,t2,-(max(fa)+0.5),max(fa)+0.5]),

title('Ð鲿');

subplot(2,2,3)

plot(t,fa)

axis([t1,t2,0,max(fa)+1]),

title('Ä£');

subplot(2,2,4)

plot(t,fn)

axis([t1,t2,-(max(fn)+1),max(fn)+1]),

title('Ïà½Ç');

fexp(-2,2,0,5,1)

波形

五、实验分析及实验小结

主要写：

1、编程时遇到的问题

2、参数对实验结果的影响

3、实验结果如何印证理论课知识点

编出了程序显示不出波形，

参数设定的不合适会使波形显示不全或者太小。

实验二 LTI 系统时域分析及MATLAB 实现

－、实验目的

1. 利用 M ATLAB 求连续系统的冲激响应与阶跃响应、求离散系统的单位响应；

二、实验内容 1、8.1 (2) 代码：

a=[1 4 2]; b=[1 0 3]; sys=tf(b,a); p=0.01; t=0:p:10; f=exp(-t);

lsim(sys,f,t);

运行结果：

0.10.20.30.40.50.60.70.80.91Linear Simulation Results

Time (seconds)

A m p l i t u d e

8.1 (4) 代码：

a=[1 5 6]; b=[6]; sys=tf(b,a); p=0.01; t=0:p:10; f=10*cos(2*t); lsim(sys,f,t);

运行结果：

012345

678910

-10

-8-6-4-202468

10Linear Simulation Results

Time (seconds)

A m p l i t u d e

2、8.3 (1) 代码：

a=[1 3 2]; b=[1];

subplot(2,2,1) impulse(b,a) subplot(2,2,2) step(b,a,0:0.01:10)

运行结果：

02460

0.10.20.3

0.4Impulse Response

Time (seconds)

A m p l i t u d e

05

10

0.20.40.6

0.8Step Response

Time (seconds)

A m p l i t u d e

8.3(3) 代码：

a=[1 4 5]; b=[1 0]; subplot(2,2,1) impulse(b,a) subplot(2,2,2) step(b,a,0:0.01:10)

运行结果：