计算机仿真作业题目,计算机仿真，第三次作业，第四章16-31

第四章 *

畅 132**224

Contents

例16 计算非周期三角波信号能量

例17 把系统H(s)转变成状态空间形式。

例18 将系统H(z)=2*(z-3)/[(z-2)*(z-1)]转换为状态空间模型[A,B,C,D]

例19 求离散时间系统H(z)=(2+3*z^(-1))/(1+0.4z^(-1)+z^(-2))的零极点向量和增益系数

例20 求两个单输入单输出子系统H1(s)=1/s+1和H2(s)=2/s+2的级联、并联反馈互联后系统的传递函数。

例21 已知FIR数字滤波器的传递函数，求其级联型结构和格型结构。

例22 已知IIR数字滤波器的传递函数，求其级联型结构和并联型结构

例23 已知描述系统的微分方程，求它的传递函数模型，零极点增益模型，极点留数模型和状态空间模型

例24 某LTI系统的单位冲激响应h(t)=exp(-0.1t),输入x(t),初始条件为零，求系统的响应y(t)

例26

有二阶系统H(s)=(2s^2+5s+1)/(s^2+2s+3)求当输入是周期为4s的方波时的输出响应。

例27

有二阶系统，求系统对100点随机噪声的响应曲线。

例28 有二阶系统，求系统的单位冲激响应

例29 有二阶系统，求其单位冲激响应

例30 对例28所示的二阶系统，当初始状态x0=[1;0]时，求系统零输入响应

例31 有二阶系统，当初始状态x0=[1;0]时，求系统零输入响应

例25 已知LTI离散系统的单位冲激响应为：h[n]=0.5^n,求输入信号序列x[n]=1(-5<=n<=4)的系统响应

例16 计算非周期三角波信号能量

dt=0.0001;t=0:dt:1;

x=tripuls(t);

E=sum(abs(x).^2*dt)

E =

0.1667

例17 把系统H(s)转变成状态空间形式。

num=[0,2,3;1,2,1];

den=[1,0.4,1];

[A,B,C,D]=tf2ss(num,den)

A =

-0.4000 -1.0000

1.0000 0

B =

1

0

C =

2.0000 3.0000

1.6000 0

D =

0

1

例18 将系统H(z)=2*(z-3)/[(z-2)*(z-1)]转换为状态空间模型[A,B,C,D]

z=[3];p=[1,2];k=2;

[A,B,C,D]=zp2ss(z,p,k)

A =

3.0000 -1.4142

1.4142 0

B =

1

0

C =

2.0000 -4.2426

D =

0

例19 求离散时间系统H(z)=(2+3*z^(-1))/(1+0.4z^(-1)+z^(-2))的零极点向量和增益系数

num=[2,3];den=[1,0.4,1];

[num,den]=eqtflength(num,den);

[z,p,k]=tf2zp(num,den)

z =

0

-1.5000

p =

-0.2000 + 0.9798i

-0.2000 - 0.9798i

k =

2

例20 求两个单输入单输出子系统H1(s)=1/s+1和H2(s)=2/s+2的级联、并联反馈互联后系统的传递函数。

num1=1;den1=[1,1];

num2=2;den2=[1,2];

[nums,dens]=series(num1,den1,num2,den2)

[nump,denp]=parallel(num1,den1,num2,den2)

[numf,denf]=feedback(num1,den1,num2,den2)

nums =

0 0 2

dens =

1 3 2

nump =

0 3 4

denp =

1 3 2

numf =

0 1 2

denf =

1 3 4

例21 已知FIR数字滤波器的传递函数，求其级联型结构和格型结构。

b=[2,13/12,5/4,2/3];

a=1;

fprintf('级联型结构系数');

[sos,g]=tf2sos(b,a)

fprintf('格型结构系数');

[K]=tf2latc(b,a)

级联型结构系数

sos =

1.0000 0.5360 0 1.0000 0 0

1.0000 0.0057 0.6219 1.0000 0 0

g =

2

格型结构系数

K =

0.2500

0.5000

0.3333

例22 已知IIR数字滤波器的传递函数，求其级联型结构和并联型结构

b=[1,-3,11,-27,18];

a=[16,12,2,-4,-1];

disp('级联型结构系数')

[sos,g]=tf2sos(b,a)

disp('并联型和结构系数')

[R,P,K]=residuez(b,a)

级联型结构系数

sos =

1.0000 -3.0000 2.0000 1.0000 -0.2500 -0.1250

1.0000 0.0000 9.0000 1.0000 1.0000 0.5000

g =

0.0625

并联型和结构系数

R =

-5.0250 - 1.0750i

-5.0250 + 1.0750i

0.9250 + 0.0000i

27.1875 + 0.0000i

P =

-0.5000 + 0.5000i

-0.5000 - 0.5000i

0.5000 + 0.0000i

-0.2500 + 0.0000i

K =

-18

例23 已知描述系统的微分方程，求它的传递函数模型，零极点增益模型，极点留数模型和状态空间模型

num=[2,-5,3];den=[2,3,5,9];

disp('系统传递函数H(s)');

printsys(num,den,'s');

disp('转为零级点增益模型');

[z1,p1,k]=tf2zp(num,den)

disp('转为零级点留数模型');

[r1,p1,h1]=residue(num,den)

disp('转为状态空间模型');

[A,B,C,D]=tf2ss(num,den)

系统传递函数H(s)

num/den =

2 s^2 - 5 s + 3

-----------------------

2 s^3 + 3 s^2 + 5 s + 9

转为零级点增益模型

z1 =

1.5000

1.0000

p1 =

-1.6441 + 0.0000i

0.0721 + 1.6528i

0.0721 - 1.6528i

k =

1

转为零级点留数模型

r1 =

-0.2322 + 0.4716i

-0.2322 - 0.4716i

1.4644 + 0.0000i

p1 =

0.0721 + 1.6528i

0.0721 - 1.6528i

-1.6441 + 0.0000i

h1 =

[]

转为状态空间模型

A =

-1.5000 -2.5000 -4.5000

1.0000 0 0

0 1.0000 0

B =

1

0

0

C =

1.0000 -2.5000 1.5000

D =

0

例24 某LTI系统的单位冲激响应h(t)=exp(-0.1t),输入x(t),初始条件为零，求系统的响应y(t)

dt=0.5;

x=ones(1,fix(10/dt));

h=exp(-0.1*[0:fix(10/dt)]*dt);

y=conv(x,h);

t=dt*([1:length(y)]-1);

plot(t,y),grid;

例26

有二阶系统H(s)=(2s^2+5s+1)/(s^2+2s+3)求当输入是周期为4s的方波时的输出响应。

num=[2,5,1];den=[1,2,3];

t=0:0.1:10;peiod=4;

u=(rem(t,peiod)>=peiod./2);

lsim(num,den,u,t);

title('方波响应')

例27

有二阶系统，求系统对100点随机噪声的响应曲线。

num=[2,-3.4,5.5];

den=[1,-1.2,0.8];

u=randn(1,100);

dlsim(num,den,u);

title('随机噪声响应')

例28 有二阶系统，求系统的单位冲激响应

a=[-0.55,-0.78;0.78,0];b=[1;0];

c=[5.96,6.45];d=[0];

figure,impulse(a,b,c,d);

title('LTI系统的冲激响应');

例29 有二阶系统，求其单位冲激响应

num=[2,-3.5,1.5];

den=[1,-1.7,0.3];

dimpulse(num,den);axis([0,6,-0.5,4.5]);

title('离散LTI系统的冲激响应');

例30 对例28所示的二阶系统，当初始状态x0=[1;0]时，求系统零输入响应

a=[-0.55,-0.78;0.78,0];b=[1;0]

c=[5.96,6.45]; d=[0];

x0=[1;0]; t0=0:0.1:20;

initial(a,b,c,d,x0,t0);

title('LTI系统的零输入响应');

b =

1

0

例31 有二阶系统，当初始状态x0=[1;0]时，求系统零输入响应

a=[-0.55,-0.78;0.78,0];b=[1;0];

c=[5.96,6.45];d=[0];

x0=[1;0]

figure,dinitial(a,b,c,d,x0);

title('离散系统的零输入响应');

x0 =

1

0

例25 已知LTI离散系统的单位冲激响应为：h[n]=0.5^n,求输入信号序列x[n]=1(-5<=n<=4)的系统响应

x=ones(1,10);lx=length(x);

h=0.5.^[0:14];lh=length(h);

lmax=max(lx,lh);

if lx>lh nx=0;nh=lx-lh;

elseif lx

else nx=0;lh=0;

end

lt=lmax;

u=[zeros(1,lt),x,zeros(1,nx),zeros(1,lt)];

t1=(-lt+1:2*lt);

h=[zeros(1,2*lt),h,zeros(1,nh)];

hf=fliplr(h);

y=zeros(1,3*lt);

for k=0:2*lt

p=[zeros(1,k),hf(1:end-k)];

y1=u.*p;

yk=sum(y1);

y(k+lt+1)=yk;

subplot(4,1,1);stem(t1,u)

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

axis([-lt,2*lt,min(u),max(u)]),hold on; ylabel('x[n]');

subplot(4,1,2);stem(t1,p);

axis([-lt,2*lt,min(p),max(p)]); ylabel('h[k-n]');

subplot(4,1,3);stem(t1,y1);

axis([-lt,2*lt,min(y1),max(y1)+eps]); ylabel('s=u.*h[k-n]');

subplot(4,1,4);stem(k,yk);

axis([-lt,2*lt,floor(min(y)+eps),ceil(max(y+eps))]); hold on;

ylabel('y[k]=sum(s)');

if

k==round(0.8*lt) disp('暂停，按任意键继续'),pause

else

pause(1),

end

end

暂停，按任意键继续

Published with MATLAB®

R2014a