matlab自动控制原理实验一,自动控制原理实验( 实验一).doc

自动控制原理实验( 实验一).doc

实验一 MATLAB的环境熟悉

一、实验目的

1、熟悉MATLAB的运行环境和帮助系统

2、熟悉MATLAB的运算和简单编程

3、掌握MATLAB的矩阵运算和二维绘图

4、掌握熟悉MATLAB的多项式运算

5、掌握熟悉MATLAB的微分方程和差分方程的求解

二、实验内容

1、MATLAB基本命令的例子操作

如：二维绘图plot、hold、text、title、xlabel、ylabel等命令；拉氏变换和拉氏反变换laplace、ilaplace等命令，z变换和z反变换ztrans、iztrans等命令；多项式运算函数root、poly、conv、polyval等函数；微分方程的求解函数dslove和差分方程的求解函数maple。

2、绘制的曲线。其中 1，t [0:0.1:15]，要求同一图画四条线： 0.001(实线)， 0.5(虚线)， 1.001(加号)， 2(星号)； 并要加图形标题和x,y轴标记。提示：用a表示，用b表示，函数表示：y 1- 1/sqrt 1-a^2 *exp -a*b*t .* sin b*sqrt 1-a^2 *t+atan sqrt 1-a^2 /a ;

解：代码：

t [0:0.1:15];

b 1;

a 0.001;

c 0.5;

d 1.001;

e 2;

y 1- 1/sqrt 1-a^2 *exp -a*b*t .* sin b*sqrt 1-a^2 *t+atan sqrt 1-a^2 /a ;

y1 1- 1/sqrt 1-c^2 *exp -c*b*t .* sin b*sqrt 1-c^2 *t+atan sqrt 1-c^2 /c ;

y2 1- 1/sqrt 1-d^2 *exp -d*b*t .* sin b*sqrt 1-d^2 *t+atan sqrt 1-d^2 /d ;

y3 1- 1/sqrt 1-e^2 *exp -e*b*t .* sin b*sqrt 1-e^2 *t+atan sqrt 1-e^2 /e ;

plot t,y,'-',t,y1,'--',t,y2,'+',t,y3,'*' ;

xlabel 't' ;

ylabel 'y' ;

title '曲线Y的绘制' ;

grid;

绘图如图2

图2

3、参考p231例9-15完成二阶系统的响应曲线绘制，t， zeta 的取值同2，并比较。

解：代码

t [0:0.1:15];

num [1];

zeta1 0.001;

den1 [1 2*zeta1 1];

zeta3 0.5;

den3 [1 2*zeta3 1];

zeta5 1.001;

den5 [1 2*zeta5 1];

zeta7 2;

den7 [1 2*zeta7 1];

[y1,x,t] step num,den1,t ;

[y3,x,t] step num,den3,t ;

[y5,x,t] step num,den5,t ;

[y7,x,t] step num,den7,t ;

plot t,y1,t,y3,t,y5,t,y7 ;

grid on;

绘图如图3

图3

4、用微分方程的求解函数dslove求解p37习题2-5(2)和p256习题9-1的解，并分别画出相应的波形 图。其中t [0,100]，取1000个点。提示：求出c的表达式后要重新运行才能画。

解：代码

(1)y dsolve '2*D2y+7*Dy+5*y 2','y 0 0,Dy 0 0'

t [0:0.1:100];

y 4/15*exp -5/2*t -2/3*exp -t +2/5;

plot t,y ;

绘图如图4(1)

如图4(1)

代码：

(2)y dsolve 'D2y+5*Dy+6*y 6','y 0 2,y 0 2' y -4*exp -3*t +5*exp -2*t +1;

plot t,y ;

绘图如图4(2)

如图4(2)

5、求多项式的根。

解：代码

p [1 4 -7 -10];

r roots p

6、求解方程组

提示：求解表达式为：[x1,x2,x3,x4] solve '2*x1-3*x2+x3+2*x4-8', 'x1+3*x2+x4-6',

'x1-x2+x3+8*x4-7', '7*x1+x2-2*x3+2*x4-5', 'x1,x2,x3,x4' ；

解：代码

syms x1 x2 x3 x4;

[x1,x2,x3,x4] solve '2*x1-3*x2+x3+2*x4-8', 'x1+3*x2+x4-6', 'x1-x2+x3+8*x4-7', '7*x1+x2-2*x3+2*x4-5', 'x1,x2,x3,x4'

7、求差分方程：p219页习题8