matlab连续时间系统复频域分析,连续时间与系统复频域分析的MATLAB实现

MATLA实验报告

实验内容(题目、程序、结果及分析)：

M6-1 已知连续时间信号的s域表示式如下，试用residue求出X(s)的部分分式展开式，并写出x(t)的实数形式表达式。 1

X(s)

41.6667

s 3.7444s 25.7604s 41.6667

3

2

解：num=[41.6667];

den=[1 3.7444 25.760 41.6667] ; [r,p,k]=residue(num,den)

[angle,mag]=cart2pol(real(r),imag(r)) 运行结果：

r =-0.9361 - 0.1895i -0.9361 + 0.1895i 1.8722 p =-0.9361 + 4.6237i -0.9361 - 4.6237i -1.8723 k =[]

angle =-2.9418 2.9418 0

mag =0.9551 0.9551 1.8722 由此可得

X(s)

-0.9361 - 0.1895is 0.9361 4.6237i

0.9551e

-0.9361 0.1895is 0.9361 4.6237i

0.9551e

1.8722 s 1.8723

1.872 2

s 0.936 1 4.6237i s 0.936 1 4.6237is 1.872 3

-j2.9418j2.9418

M6-2 已知某连续时间LTI系统的微分方程为 y''(t)+4y'(t)+3y(t)=2x'(t)+x(t)

x(t)=u(t),y(0)=1,y'(0)=2,试求系统的零输入响应，零状态响应和完全响应，并划出相应的图形。

解：对微分方程两边进行单边laplace变换，得

sY(s) sy(0) y'(0) 4[sY(s) y(0)] 3Y(s) 2s 1 X(s)

2

整理后，得Y(s)

sy(0) y'(0) 4y(0)

s 4s 3

2s 1s 4s 3

22

2s 1s 4s 3

2

X(s)

s 6s 4s 3

2

X(s)

零输入响应的Yzi(s)

2

s域表示式为

s 6s 4s 3

Matlab程序如下：