计算机仿真技术第三版第二章,计算机仿真技术第三次作业 （第四章例题16-31）...

专业通信工程

姓名张小芳

学号

1323010214

第四章例题(继续)

Contents

例4-16 非周期三角波信号能量的MATLAB计算。

例4-17 将系统H(s)=[2*s+3;s^2+2*s+1]/(s^2+0.4*s+1)转换成状态空间形式

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

例4-19

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

例4-20

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

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

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

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

例4-24 某LTI系统的单位冲击响应，初始条件为零，求系统的响应。

例4-26

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

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

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

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

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

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

例4-25

某LTI离散系统的单位冲击响应h(n)=0,5^n(n=0,1,2,...14求输入信号序列x[n]=(-5<=n<=4)的系统响应。

例4-16 非周期三角波信号能量的MATLAB计算。

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

x=tripuls(t);

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

E =

0.1667

例4-17 将系统H(s)=[2*s+3;s^2+2*s+1]/(s^2+0.4*s+1)转换成状态空间形式

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

例4-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

例4-19

求离散时间系统H(z)=[2+3*z^(-1)]/[1+0.4*z^(-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

例4-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

例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

例4-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

27.1875

P =

-0.5000 + 0.5000i

-0.5000 - 0.5000i

0.5000

-0.2500

K =

-18

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

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

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

printsys(num,den,'s');

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

[z1,p1,k1]=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.0721 + 1.6528i

0.0721 - 1.6528i

k1 =

1

转为零极点留数模型

r1 =

-0.2322 + 0.4716i

-0.2322 - 0.4716i

1.4644

p1 =

0.0721 + 1.6528i

0.0721 - 1.6528i

-1.6441

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

例4-24 某LTI系统的单位冲击响应，初始条件为零，求系统的响应。

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

title('张小芳 1323010214')

h =

Columns 1 through 8

1.0000

0.9512  0.9048  0.8607  0.8187  0.7788  0.7408  0.7047

Columns 9 through 16

0.6703

0.6376  0.6065  0.5769  0.5488  0.5220  0.4966  0.4724

Columns 17 through 21

0.4493

0.4274  0.4066  0.3867  0.3679

例4-26

有二阶系统H(s)=(2*s^2+5*s+1)/(s^2+2*s+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('方波响应 zhang')

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

num=[2,-3.4,5.5];

den=[1,-1.2,0.8];

u=randn(1,100);

dlsim(num,den,u);

title('随机噪声响应 zhang')

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

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

b=[1;0];

c=[5.96,6.45];

d=[0];

impulse(a,b,c,d);title('LTI系统的冲激响应 zhang')

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

num=[2,-3.5,1.5];

den=[1,-1.7,0.3];

dimpulse(num,den,6);

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

例4-30 对例4-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系统的零输入响应 zhang')

例4-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]

dinitial(a,b,c,d,x0);title('离散系统的零输入响应 zhang');

x0 =

1

0

例4-25

某LTI离散系统的单位冲击响应h(n)=0,5^n(n=0,1,2,...14求输入信号序列x[n]=(-5<=n<=4)的系统响应。

clear

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)');

title('张小芳

1323010214')

if k==round(0.8*lt)

disp('暂停，按任意键继续'),pause

else pause(1),

end

end

暂停，按任意键继续

Published with MATLAB® 7.11