南邮matlab考试,南邮MATLAB实验三

-1/6

11/36

3.已知有理分式R(x) N(x)3，其中N(x) (3xD(x) x)(x3 0.5)，D(x) (x2 2x 2)(5x3 2x2 1)。求该分式的商多项式和余多项式。

format rat

p1=conv([3,0,1,0],[1,0,0,0.5]);

p2=conv([1,2,-2],[5,2,0,1]);

[q,r]=deconv(p1,p2);

cq='商多项式';cr='余多项式';

disp([cq,poly2str(q,'s')]),disp([cr,poly2str(r,'s')])

商多项式 0.6 s - 1.44

余多项式 -4.4409e-016 s^6 + 8.8818e-016 s^5 +21.88 s^4 - 5.34 s^3 -

5.52 s^2 + 4.58 s - 2.88

qp2=conv(q,p2);

pp1=qp2+r;

pp1==p1

ans =

1 1 1 1 1 1 1

4.求解Lorenz模型的状态方程，并图示.

初值在(0,1)中取值

function xdot=DxDt(t,x)

a=16;b=4;c=45;

xdot=[a*(x(2)-x(1);c*x(1)-x(2)-x(1)*x(3);x(1)*x(2)-b*x(3)];

tspan=[0,30];

x0=[0.5;0.4;0.6];

[tt,xx]=ode45(@DxDt,tspan,x0);

plot(tt,xx(:,1))

xlabel('t'),title('x1(t)')