实验八 连续系统的复频域分析

a=-0.5:0.08:0.5;

b=-2:0.08:2;

[a,b]=meshgrid(a,b);

d=ones(size(a));

c=a+i*b;

c=c.*c;

c=c+d;

c=1./c

c=abs(c);

surf(a,b,c);

axis=([-0.5,0.5,-2,2,0.15]); title('单边正弦信号拉氏变换图'); colormap(hsv);

2.

a=0:0.5:5;

b=-20:0.1:20;

[a,b]=meshgrid(a,b);

c=a+i*b;

c=(1-exp(-2*c)./c);

c=abs(c)

mesh(a,b,c);

sufr(a,b,c);

view(-10,20);

axis([-0.5:-20,20,0.2]);

title('拉氏变换S域像函数');

w=-20:0.1:20;

Fw=(2*sin(w).*exp(i*w));

plot(w,abs(Fw))

title=('傅里叶变换');

xlabel('频率w');

3.

a=[1,2,-3,2,1];

b=[1,4];

sjdt(a,b);

a=[1,5,16,30];

b=[5,20,25,0];

sjdt(a,b);

4.

a=[8,2,3,15];

b=[1,3,2];

[p,q]=sjdt(a,b);

5.

a=[1,0];

b=[1];

impulse(b,a)

a=[1,2];

b=[1];

impulse(b,a)

a=[1,-2];

b=[1];

impulse(b,a)

a=[1,1,16,25];

b=[1];

impulse(b,a)

a=[1,0,16];

b=[1];

impulse(b,a)

a=[1,-1,16,25];

b=[1];

impulse(b,a)

6.

q=[0,0];

p=[-100,100];

f1=0;

f2=1000;

k=0.1;

splxy(f1,f2,k,p,q) q=[0,0];

f1=0;

f2=1000;

k=0.1;

p=[-500,-1000]; splxy(f1,f2,k,p,q) q=[0,0];

f1=0;

f2=1000;

k=0.1;

p=[-2000,-4000]; splxy(f1,f2,k,p,q) 7.

8

a=[1,3];

b=[1,0,4,0];

[v,p,k]=residue(b,a)

a=[1,5,4];

b=[1,5,6,0];

[v,p,k]=residue(b,a)

9.

b=[1,4];

a=[1,5,6,0];

sys=tf(b,a);

pzmap(sys);

syms st;

H=sym('(s+4)/(s^3+2*s^2+2*s)'); h=ilaplace(H)

ezplot(h,[-8,8]);

axis([0,8,-1,3]);