计控实验报告Matlab系统仿真
- 格式:doc
- 大小:124.00 KB
- 文档页数:3
实验Matlab系统仿真(一)MA TLAB语言编程:
1用MATLAB完成矩阵的各种运算
打开matlab软件,输入矩阵矩阵A=
在“command window(命令窗口)”里输入:
A=[11 12 13 14;21 22 23 24;31 32 33 34;41 42 43 44]%矩阵A A(:1)%取矩阵A的第一列元素
A(2,:)%取矩阵A的第二行素
A(1:2,2:3)%取矩阵A的第一二行第二三列元素
A(:,1:2)%取矩阵A的第一列与第二列元素
2 绘制数学函数的图形
函数y(t)=1-2e-tsin(t) (0=<t<=8)
在command window窗口输入matlab程序如下:
t=0:0.1:8; %建立向量t
y=1-2*exp(-t.*sin(t));
plot(t,y); %利用plot命令绘图
xlabel(‘t’);ylable(‘y=1-2*e^(-t*sin(t))’);
3 理解函数文件与命令文件的区别.并自编函数文件并调用已知非线性方程组{sinx
求系统在初始条件为x(0)=1,x(0)’的数值解。
(1)定义方程组的M函数equ.m, matlab程序如下:
fnction xs=equ(xi)
x=xi(1);y=xi(2);z=xi(3);
xs=zeros(3,1);
xs(1)=sin(x)+y^2+log(z)-7;
xs(2)=3*x+2^y-z^3+1
xs(3)=x+y-5
(2)方程组求解
x0=[1 1 1];
xyz=fsolve(@equ,x0)
(二)数值积分算法与函数调用
若matlab求解初值问题的解,其方程如下:
x’=x-t2 ;x(0)=1, t=[0 3]
1.matlab程序编写如下:
Function xdot=fun21(t,x)
xdot=x-t^2
并以fun21.m存盘
2.写如下程序并以fzsy22.m存盘
t0=0;tf=3;tol=1e-6;
x0=1;trace=1;
[t,x]=ode23(‘fun21’,t0,tf,tol,trace)
plot(t,x)
3.在命令窗口下运行fzsy21即可求出x的解,并画出曲线
(三)simulink建立仿真模型
1、输入simulink或双击命令窗日工具栏右起第三个按钮(NEW SIMULINK MODEL),Library simulink然后指向File菜单下拉菜单open 调出fzsv3 1.mdl文件。
然后在fzsy31.mdl文件的菜单观察并记录有关设置参数,然后指向Start下拉菜单单击一次观察输出波形。
该仿真模型(SIMULUNK系统结构框图)如下:
在Matlab主窗口的Workspace中找到y变量,然后双击,可以看到tout=0:0.1:10对于y 的输出,示波器绘图正是根据tout数组与y数组绘制。
2、以双积分系统为例G(s)=1/s2, 绘制根轨迹
g1=tf([1],[1 0 0])
rlocus(g1);
g2= tf([1 1],[1 5 0 0])
rlocus(g2);
g3= tf([1 5],[1 1 0 0])
rlocus(g3);
3、设Kc={0.1 0.5 1},对这三种补偿器进行仿真,仿真时间可为{10 50 5}
超期补偿器:K(s)=Kc(s+1)/(s+5)。