信电学院MATLAB第二次试验.doc

实验四 连续时间信号的频域分析

1、（1）编程实现信号的幅度频谱
M文件：
syms t w;
Gt1=sym('Heaviside(t+1)-Heaviside(t-1)');
Gt2=sym('Heaviside(2*t+1)-Heaviside(2*t-1)');
Fw1=fourier(Gt1,t,w)
Fw2=fourier(Gt2,t,w)
FP1=abs(Fw1);
FP2=abs(Fw2);
subplot(2,1,1);ezplot(FP1,[-10*pi,10*pi]);grid;
axis([-10*pi,10*pi,0,2.2]);xlabel('w');ylabel('Fw1');
%title('Gt1的幅频特性曲线');
subplot(2,1,2);ezplot(FP2,[-20*pi,20*pi]);grid;
axis([-20*pi,20*pi,0,1.2]);xlabel('w');ylabel('Fw2');
%title('Gt2的幅频特性曲线');

运行结果：

Fw1 =

2/w*sin(w)

Fw2 =

2/w*sin(1/2*w)
>>
2、(2)利用ifourier()函数，求频谱函数的傅氏反变换。
M文件：
syms t w;
Fw=sym('(-w^2+5*i*w-8)/(-w^2+6*i*w+5)');
ft=ifourier(Fw,w,t)

运行结果：
ft =

dirac(t)+(-3*exp(-t)+2*exp(-5*t))*heaviside(t)

>>

实验五 系统的零极点及频域响应特性

1、（1）由系统函数H（s），求其零极点，并画出零极点图。
M文件：
num=[1 0 1];
den=[1 2 5];
[z,p]=tf2zp(num,den)
zplane(z,p)
运行结果：

z =
0 + 1.0000i
0 - 1.0000i

p =
-1.0000 + 2.0000i
-1.0000 - 2.0000i
>>
2、(1) 已知系统函数H(s)，求其频率特性。
M文件：
num=[2 0];
den=[1 2^0.5 1];
w=logspace(-2,2);
freqs(num,den,w)
运行结果：

实验六 LTI系统的响应

1、（1）已知系统微分方程和激励信号e(t)，用MATLAB绘出系统单位冲激响应和系
统零状态响应的波形。
M文件：
a=[1 4 4];
b=[1 3];
subplot(2,1,1);impulse(b,a,4)
p=0.01;
t=0:p:4;
x=exp(-t);
subplot(2,1,2);lsim(b,a,x,t)

运行结果：

2、（2）用MATLAB求出差分方程所描述的离散系统，在0~20时间范围内的单位函数
相应、阶跃响应和系统零状态响应的数值解，并绘出其波形。
M文件：
a=[1 -0.7 0.1];
b=[7 -2 0];
figure(1);
impz(b,a,-1:20)
hk=impz(b,a,0:20)
k=0:20;
x=ones(size(k));
y=filter(b,a,x)
figure(2);
subplot(2,1,1);dstep(b,a,20)
subplot(2,1,2);stem(k,y)

运行结果：
hk =

7.0000
2.9000
1.3300
0.6410
0.3157
0.1569
0.0783
0.0391
0.0195
0.0098
0.0049
0.0024
0.0012
0.0006
0.0003
0.0002
0.0001
0.0000
0.0000
0.0000
0.0000

y =

Columns 1 through 10

7.0000 9.9000 11.2300 11.8710 12.1867 12.3436 12.4218
12.4609 12.4805 12.4902

Columns 11 through 20

12.4951 12.4976 12.4988 12.4994 12.4997 12.4998 12.4999
12.5000 12.5000 12.5000

Column 21

12.5000

>>