基于MATLAB的系统仿真实验指导书

基于MATLAB的系统仿真实验实验指导书说明MATLAB是MATrix LABoratory的缩写,是一种基于矩阵的数学与工程计算系统,可以用作动态系统的建模与仿真。

研究系统的结构和参数的变化对系统性能的影响可利用MATLAB强大的计算和作图功能,因此本实验采用MATLAB仿真研究连续控制系统和离散控制系统的性能分析过程。

通过该实验提高学生对控制系统的分析与设计能力,加深对《自动控制原理》课程内容的理解。

该实验教学大纲中只安排了两个学时,本实验的内容比较多,要求同学们参考有关MATLAB应用的书籍利用课后时间完成。

同时可以将MATLAB用于习题的验证和课程的学习中。

仿真实验基于MATLAB 连续系统的仿真一、 实验目的在研究系统的结构和参数的变化对系统性能的影响时,采用解析和作图的方法比较麻烦,而且误差也大,用MATLAB 仿真实现则简单方便,精度高。

本实验采用MATLAB 实现控制系统的数学描述、控制系统的时域分析及根轨迹和频率特性分析。

通过该实验,加深学生对系统阶次,型号,参数与系统性能的关系的理解。

二、实验环境在计算机Windows 环境下安装好MATLAB6.3以上版本后,双击MATLAB 图标或成“开始”菜单打开MATLAB ,即可进入MATLAB 集成环境。

三、实验原理三、MATLAB 应用实例1.拉氏变换和反变换例 求22)(2++=t t t f 的拉氏变换 解键入 syms s t;ft=t^2+2*t+2; st=laplace(ft,t,s)运行结果为 st=2/s^3+2/s^2+2/s例 求)2)(34(6)(2++++=s s s s s F 的拉氏反变换 解键入 syms s t;Fs=(s+6)/(s^2+4*s+3)/(s+2); ft=ilaplace(Fs,s,t)运行结果为 ft=3/2*exp(-3*t)+5/2*exp(-t)-4*exp(-2*t)2。

求根运算例 求多项式 43)(23++=s s s p 的根,再由根建多项式。

解键入 p=[1 3 0 4]; r=root(p) 运行结果为 r=-3.35530.1777+1.0773i 1.7777-1.0773i键入 p=poly(r)运行结果为 p=1.0000 3.0000 0.0000 40000例 实现多项式相乘:)4)(123(2+++s s s ,并求5-=s 时的值。

解键入 p=[3 2 1]; q=[1 4]; n=conv(p,q) 运行结果为 n=3 14 9 4 键入 vlaue=polyval(n,-5) 运行结果为 value=-663.微分方程求解例 解下列微分方程:1)(2)(3)(322=++t y dt t dy dtt y d 初始条件 0)0()0('==y y 解键入 y=dsolve('3*D2y+3*Dy+2*y=1','y(0)=0,Dy(0)=0') 运行结果为 y=1/2-1/2*exp(-1/2*t)*cos(1/6*15^(1/2)*t)-1/10*15^(1/2)*exp(-1/2*t)*sin(1/6*15^(1/2)*t)4.传递函数例 求如图所示系统的传递函数。

解键入numg=[1];deng=[500 0 0]; numh=[1 1];denh=[1 2];[num,den]=feedback(numg,deng,numh,denh,-1); printsys(num,den) 运行结果为 num/den=s +2------------------------------500s^3 + 1000 s^2 + s+1其中, 由函数 printsys(num,den) 打印出传递函数。

5。

系统输出响应及性能分析例 二阶系统闭环传递函数的标准形式为 2222)()()(nn n s s s R s C s ωξωω++==Φ 若n ω确定,系统的瞬态响应和ξ的取值有关。

下面用MATLAB 分析在不同的ξ值时,系统的单位阶跃响应。

所用的MA TLAB 程序如下:%n ω=1t=0:0.1:12;num=[1];zeta1=0;den1=[1 2*zeta1 1]; zeta3=0.3;den3=[1 2*zeta3 1]; zeta5=0.5;den5=[1 2*zeta5 1]; zeta7=0.7;den7=[1 2*zeta7 1]; zeta9=1.0;den9=[1 2*zeta9 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); [y9,x,t]=step(num,den9,t); plot(t,y1,t,y3,t,y5,t,y7,t,y9,) grid on; 运行结果见图。

例 已知AAK s s K s 55.345)(2++=Φ,分别计算K A =200时,系统的性能指标t p ,t s ,%σ。

解 MA TLAB 程序如下:t=0:0.01:2; num=[1000];den=[1 34.5 1000]; [y,x,t]=step(num,den,t); plot(t,y); %求超调量 maxy=max(y); yss=y(length(t));pos=100*(maxy-yss)/yss; %求峰值时间 for i=1:1:201if y(i)==maxy,n=i;end endtp=(n-1)*0.01; %求调节时间 for i=1:1:201if (y(i)<1.05 & y(i)>0.95),m=i;end break; endts=(m-1)*0.01;计算结果为t p =0.12s ,t s =0.17s ,%σ=12.93%例 系统闭环特征方程为 0242)(23=+++=s s s s q 用MA TLAB 判系统稳定性。

解键入 den=[1 1 2 24]; roots(den) 运行结果为 ans=-3.00001. 000+2.6458i 1.0000-2.6458i可见,系统有两个右根,故系统不稳定。

例 求出下列系统11)(1+=s s G ,)1(1)(2+=s s s G ,)1(14)(23++=s s s s G 单位阶跃 响应及稳态误差。

解MATLAB 程序如下:t=0:0.1:20;[num1,den1]=cloop([1],[1 1]); [num2,den2]=cloop([1],[1 1 0]); [num3,den3]=cloop([4 1],[1 1 0 0]); y1=step(num1,den1,t); y2=step(num2,den2,t); y3=step(num3,den3,t); subplot(311); plot(t,y1);subplot(312);plot(t,y2); subplot(313);plot(t,y3); er1=y1(length(t))-1; er2=y2(length(t))-1; er3=y3(length(t))-1;运行结果见图。

在命令窗口可得er1 = -0.5000 % 0型系统的稳态误差 er2 = 2.4294e-005 % Ⅰ型系统的稳态误差 er3 = -3.3525e-004 % Ⅱ型系统的稳态误差例 用SIMULINK 对第七节中单闭环无静差调速系统进行仿真。

解 在SIMULINK 下建立系统仿真结构图如图3-46所示。

仿真的结果见图3-47。

6。

绘制系统的根轨迹图例 设一单位反馈系统如图4-31所示,试绘制该系统的根轨迹。

解键入num=[1];den=[1 3 2 0]; rlocus(num,den) 运行结果见图4-32。

例 已知系统前向通道传递函数为5)(+=s K s G c ,)8(1)(0++=s s s s G ,反馈通道传递函数为21)(+=s s H ,试绘出系统闭环根轨迹图。

解键入 Gc=tf(1,[1 5]);Go=tf([1 1],[1 8 0]); H=tf(1,[1 2]); rlocus(Gc*Go*H);v=[-10 10; –10 13];axis(v); grid on运行结果见图4-33。

图4-33 例4-18系统根轨迹7.绘制系统的频率特性图例5 振荡环节传递函数的标准形式为2222)()()(nn n s s s R s C s ωζωω++==Φ 当n ω确定时,系统的频率特性与ζ的取值有关。

用MATLAB 绘出1=n ω时,在不同ζ取值下,系统的Bode 图和Nyquist 图。

解 1)作Bode 图程序: num=[1];zeta1=0.1;den1=[1 2*zeta1 1]; zeta3=0.3;den3=[1 2*zeta3 1];zeta5=0.5;den5=[1 2*zeta5 1]; zeta7=0.7;den7=[1 2*zeta7 1]; zeta9=1.0;den9=[1 2*zeta9 1];[mag1,phase1,w1]=bode(num,den1); [mag3,phase3,w3]=bode(num,den3); [mag5,phase5,w5]=bode(num,den5); [mag7,phase7,w7]=bode(num,den7); [mag9,phase9,w9]=bode(num,den9); subplot(211);semilogx(w1,20*log10(mag1),w3,20*log10(mag3),…w5,20*log10(mag5),w7,20*log10(mag7), w9,20*log10(mag9)); subplot(212);semilogx(w1,phase1,w3,phase3,w5,phase5,w7,phase7,w9,phase9); 运行结果见图。

2)作Nyquist 图程序: num=[1];zeta1=0.4;den1=[1 2*zeta1 1]; zeta6=0.6;den6=[1 2*zeta6 1]; zeta8=0.8;den8=[1 2*zeta8 1]; [re1,im1]=nyquist(num,den1); [re2,im2]=nyquist(num,den6);[re3,im3]=nyquist(num,den8);plot(re1,im1,re2,im2,re3,im3); 运行结果见图。

合集下载

现代控制理论-基于MATLAB的实验指导书课程设计指导书

现代控制理论-基于MATLAB的实验指导书课程设计指导书

现代控制理论基于MATLAB的实验指导书第一部分实验要求1.实验前做好预习。

2.严格按照要求操作实验仪器,用毕恢复原状。

3.实验完成后,由指导教师检查实验记录、验收仪器后,方可离开。

4.实验报告应包括以下内容:1)实验目的;2)实验原理图;3)实验内容、步骤;4)仿真实验结果(保留仿真实验波形,读取关键参数);5)仿真实验结果分析。

第二部分MATLAB平台介绍实际生产过程中,大部分的系统是比较复杂的,并且要考虑安全性、经济性以及进行实验研究的可能性等,这在现场实验中往往不易做到,甚至根本不允许这样做。

这时,就需要把实际系统建立成物理模型或数学模型进行研究,然后把对模型实验研究的结果应用到实际系统中去,这种方法就叫做模拟仿真研究,简称仿真。

到目前为止,已形成了许多各具特色的仿真语言。

其中美国Mathworks软件公司的动态仿真集成软件Simulink与该公司著名的MATLAB软件集成在一起,成为当今最具影响力的控制系统应用软件。

国内MA TLAB软件的著名论坛为“MATLAB中文论坛”,网址为:https:///forum.php,建议同学们注册并参与论坛相关内容的讨论。

图1 MA TLAB仿真环境第三部分 实验实验一线性系统的时域分析实验目的熟悉MATLAB 环境,掌握用MATLAB 控制系统工具箱进行线性定常系统的时域分析、能控性与能观性分析、稳定性分析的方法。

实验要求完成指导书规定的实验内容,记录并分析实验结果,写出实验报告。

实验内容1.已知系统的状态模型,求系统在单位阶跃输入下的各状态变量、输出响应曲线。

例:[]⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡+⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡--=⎥⎦⎤⎢⎣⎡2121214493.69691.1,0107814.07814.05572.0x x y u x x x x 。

键入:a = [-0.5572, -0.7814; 0.7814,0]; b = [1; 0]; c = [1.9691,6.4493]; d = 0;[y, x, t]=step(a, b, c, d); plot(t, y); grid (回车,显示输出响应曲线。

《Matlab与控制系统仿真》实验指导书

《Matlab与控制系统仿真》实验指导书

机械与汽车工程学院《Matlab控制系统仿真》实验指导书学院班级姓名学号浙江科技学院机械与汽车工程学院制实验一 MATLAB语言基本命令1 实验目的1. 掌握科学计算的有关方法,熟悉MA TLAB语言及其在科学计算中的运用;2. 掌握MATLAB的命令运行方式和M文件运行方式;3. 掌握矩阵在MA TLAB中的运用。

2 实验器材计算机WinXP、Matlab7.0软件3 实验内容(1). 输入A=[7 1 5;2 5 6;3 1 5],B=[1 1 1; 2 2 2;3 3 3],在命令窗口中执行下列表达式,掌握其含义:A(2, 3) A(:,2) A(3,:) A(:,1:2:3)A(:,3).*B(:,2) A(:,3)*B(2,:) A*B A.*BA^2 A.^2 B/A B./A(2).输入C=1:2:20,则C(i)表示什么?其中i=1,2,3, (10)(3).查找已创建变量的信息,删除无用的变量;(4). 试用help命令理解下面程序各指令的含义:cleart =0:0.001:2*pi;subplot(2,2,1);polar(t, 1+cos(t))subplot(2,2,2);plot(cos(t).^3,sin(t).^3)subplot(2,2,3);polar(t,abs(sin(t).*cos(t)))subplot(2,2,4);polar(t,(cos(2*t)).^0.5)4 实验步骤:打开MA TLAB程序,将实验内容中的题目依次输入MATLAB中,运行得到并记录结果,最后再对所得结果进行验证。

5 实验报告要求记录实验数据,理解其含义实验二 MATLAB语言程序设计1 实验目的(1)掌握Matlab程序的编制环境和运行环境。

(2)掌握Matlab程序的编写方法。

(3)能编写基本的数据处理Matlab程序。

(4)能编写基本的数据可视化Matlab程序。

2 实验器材计算机WinXP、Matlab7.0软件3 实验内容(1) Matlab脚本文件编写和执行(2) Matlab 函数文件的编写和调用(3) nargm和nargout函数使用方法(4) 局部变量与全局变量使用4 实验步骤1、Matlab命令文件编写(1) 建立自己工作目录,如/Mywork。

《MATLAB语言与控制系统仿真》实验指导书01

《MATLAB语言与控制系统仿真》实验指导书01

《MATLAB语言与控制系统仿真》实验指导书实验一、MATLAB语言环境与基本运算一、实验目的及要求1.学习了解MATLAB语言环境2.练习MATLAB命令的基本操作3.练习MATLAB数值运算相关内容4.练习MATLAB符号运算相关内容5.撰写实验报告二、实验内容1.熟悉Matlab语言环境2.1).学习了解MATLAB语言环境MATLAB语言操作界面(主界面的各个窗口)变量查询命令who, whos 目录与目录结构搜索路径联机帮助2).MATLAB基本操作命令demos,clc,clf,clear,contro-c(^c),diary3.Matlab数值运算与符号运算1).MATLAB数值运算相关内容MATLAB变量及变量赋值初等矩阵函数ones, zeros, eye, rand, randn, size矩阵的基本运算矩阵的特征运算det, eig, rank, svd 矩阵的分解运算基本数学函数矩阵的运算与点运算多项式的运算2).MATLAB符号运算相关内容符号变量的定义基本符号运算符号运算的扩展实验二、MATLAB编程与绘图一、实验目的及要求1.练习MATLAB编程2.练习MATLAB二维绘图与三维绘图3.撰写实验报告二、实验内容1.MATLAB编程MATLAB关系运算MATLAB逻辑运算switch分支选择语句if条件语句for循环语句while循环语句2.MATLAB二维绘图与三维绘图MATLAB二维绘图plot, subplot, figureMATLAB三维绘图plot3, mesh, surf z=rand(10,10);>> mesh(z)[X,Y,Z]=cylinder(r,30);surf(X,Y,Z)实验三、控制系统时间响应与频率响应分析一、实验目的及要求1.练习控制系统仿真绘图2.练习控制系统的数学模型的描述3.SIMULINK系统仿真4.撰写实验报告二、实验内容1.控制系统数学模型与时域分析时间响应绘图impulse, step, lsim控制系统的数学模型tf, zp, ss; printsys, conv, feedback,控制系统的数学模型之间的相互转换2.控制系统频域分析与稳定性分析频率特性绘图bode, nichols, nyquist根轨迹图rlocus3.SIMULINK系统仿真SIMULINK的功能模块利用SIMULINK进行控制系统仿真。

2021级《系统仿真与MATLAB语言》实验指导书

2021级《系统仿真与MATLAB语言》实验指导书

2021级《系统仿真与MATLAB语言》实验指导书系统仿真与matlab语言实验指导书韶关学院2021级自动化专业系统仿真与matlab语言实验说明1、所有同学必须出席实验课,按照定出顺序上机。

2、每次上机请登记在实验室记录本上,固定机位,一人一机。

3、本课程7次上机实验(其中前6次实验要提交纸质实验报告,第7次实验递交电子版实验报告)。

4、在机房内务必保持安静,不得大声喧哗,不得从事与实验无关事项,如打游戏、看电影等。

5、写好预习报告:针对实验内容建议用自己的语言详细的载明实验目的、实验原理,并针对实验内容预先搞好程序的撰写。

所有实验项目必须提早编写程序不写下复习报告严禁出席实验6、数据记录记录实验的原始数据,实验期间当场提交。

拒绝抄袭。

实验过程中要能回答我的提问。

7、实验内容总结所搞实验的内容提问有关思考题2021、9实验一matlab语言工作环境和基本操作一、实验目的:熟悉matlab的工作环境,学会使用matlab进行一些简单的运算。

二、实验内容:matlab的启动和退出,熟悉matlab的桌面(desktop),包括菜单空间(workspace)等;顺利完成一些基本的矩阵操作方式;自学采用在线协助系统。

三、实验步骤:1、启动matlab,熟识matlab的桌面。

2、在命令窗口执行命令完成以下运算,观察workspace的变化,记录运算结果。

(1)(365-52?2-70)?3(2)area=pi*2.5^2(3)未知x=3,y=4,在matlab中求z:x2y3z?2?x?y?(4)输入c=1:2:20,则c(i)表示什么?其中i=1,2,3,…,10;(5)掌握matlab常用命令>>who%列举工作空间中变量>>whos%列出工作空间中变量,同时包括变量详细信息>>savetest%将工作空间中变量存储到test.mat文件中>>loadtest%从test.mat文件中读取变量到工作空间中>>clear%清除工作空间中变量>>help函数名%对选好函数的功能、调用格式及有关函数得出表明>>lookfor%搜寻具备某种功能的函数但却不晓得该函数的精确名称例如:lookforlyapunov可以列举与lyapunov有关的所有函数。

matlab控制系统计算机仿真实验-完整版

matlab控制系统计算机仿真实验-完整版

MALTAB 仿真实验指导书实验一实验题目:欧拉法&梯形法的MATLAB 实现实验目的:1.熟练掌握MATLAB 的使用方法2.牢记欧拉法、梯形法的计算过程3.熟悉欧拉法、梯形法以及实现二阶动态响应的程序编写 实验内容:已知被控对象的系数矩阵分别为A=[-5 -2 -1 -0.5;4 0 0 0;0 2 0 0;0 0 1 0 ]B=[1;0;0;0];C=[0 0 0.25 0.5];D=0;根据欧拉法、梯形法的递推公式,应用MATLAB 语言编写相应的仿真程实验要求:1.取计算步长65.0=h ,初值均为零,输入为阶跃信号,取25=u ,研究系统25秒的动态过程。

2.取计算步长01.0=h ,初值均为零,输入为阶跃信号,取25=u ,研究系统25秒的动态过程。

实验算法:欧拉法递推公式:),(1k k k k y t hf y y +=+梯形法的递推公式: )],(),([2),(011101++++++=+=k k k k k k k k k k y t f y t f h y y y t hf y y实验方法:利用所学过数值积分方法(欧拉法、梯形法),通过MATLAB 语言对给定的系统进行仿真实验步骤:1.了解并掌握基本数值积分的方法,即欧拉法、梯形法,并做比较,了解它们之间的联系与区别和优缺点,其中重点掌握梯形法。

2.通过给定的系统,利用欧拉法、梯形法编写相应MATLAB 语言,实现仿真,得出相应的仿真曲线。

3.比较仿真实验结果,并得出结论。

4.撰写实验报告。

实验程序:1.欧拉法A=[-5 -2 -1 -0.5;4 0 0 0;0 2 0 0;0 0 1 0];B=[1;0;0;0];C=[0 0 0.25 0.5];D=0;x0=[0;0;0;0];% x0为状态变量的初值,此处以列向量表示;u=25;% u为输入向量;t0=0;% t0为仿真时间的起始时刻;tf=15;% tf为仿真时间的结束时刻;h=0.65;% h=0.01 h为仿真时所取的仿真步长;m=(tf-t0)/h;[r,c]=size(A);for i=1:mfor j=1:rx(j)=x0(j)+h*(A(j,:)*x0+B(j,:)*u);endy(i)=C*x';x0=x';t(i)=i*h;endplot(t,y)grid ontitle('useEuler')2.梯形法A=[-5 -2 -1 -0.5;4 0 0 0;0 2 0 0;0 0 1 0];B=[1;0;0;0];C=[0 0 0.25 0.5];D=0;x0=[0;0;0;0];% x0为状态变量的初值,此处以列向量表示;u=25;% u为输入向量;t0=0;% t0为仿真时间的起始时刻;tf=15;% tf为仿真时间的结束时刻;h=0.65;% h=0.01 h为仿真时所取的仿真步长;m=(tf-t0)/h;[r,c]=size(A);for i=1:mfor j=1:rx(j)=x0(j)+h*(A(j,:)*x0+B(j,:)*u);endx1=x';for k=1:rxx(k)=x0(k)+0.5*h*((A(k,:)*x0+B(k,:)*u)+(A(k,:)*x1+B(k,:)*u)); endy(i)=C*xx';x0=xx';t(i)=i*h;endplot(t,y)grid ontitle('useLadder')实验报告要求:1.书写实验报告,其中包括实验题目,实验目的,实验内容,实验要求,实验思路,实验方法,实验步骤,实验程序等。

MATLAB实验指导书(共5篇)

MATLAB实验指导书(共5篇)

MATLAB实验指导书(共5篇)第一篇:MATLAB实验指导书MATLAB 实验指导书皖西学院信息工程学院实验一 MATLAB编程环境及简单命令的执行一、实验目的1.熟悉MATLAB编程环境二、实验环境1.计算机2.MATLAB7.0集成环境三、实验说明1.首先应熟悉MATLAB7.0运行环境,正确操作2.实验学时:2学时四、实验内容和步骤1.实验内容(1)命令窗口的使用。

(2)工作空间窗口的使用。

(3)工作目录、搜索路径的设置。

(4)命令历史记录窗口的使用。

(5)帮助系统的使用。

(6)了解各菜单的功能。

2.实验步骤(1)启动MATLAB,熟悉MATLAB的桌面。

(2)进入MATLAB7.0集成环境。

(3)在命令窗口执行命令完成以下运算,观察workspace的变化,记录运算结果。

1)(365-52⨯2-70)÷3 2)>>area=pi*2.5^2 3)已知x=3,y=4,在MATLAB中求z:x2y3 z=2(x-y)4)将下面的矩阵赋值给变量m1,在workspace中察看m1在内存中占用的字节数。

⎡162313⎤⎢511108⎥⎥m1=⎢⎢97612⎥⎢⎥414151⎣⎦执行以下命令>>m1(2 , 3)>>m1(11)>>m1(: , 3)>>m1(2 : 3 , 1 : 3)>>m1(1 ,4)+ m1(2 ,3)+ m1(3 ,2)+ m1(4 ,1)5)执行命令>>helpabs 查看函数abs的用法及用途,计算abs(3 + 4i)6)执行命令>>x=0:0.1:6*pi;>>y=5*sin(x);>>plot(x,y)7)运行MATLAB的演示程序,>>demo,以便对MATLAB有一个总体了解。

五、思考题1、以下变量名是否合法?为什么?(1)x2(2)3col(3)_row (4)for2、求以下变量的值,并在MATLAB中验证。

控制系统MATLAB仿真实验指导书.

1实验5. 控制理论仿真实验1 控制系统的建模一、实验目的1.学习在MATLAB 命令窗口建立系统模型的方法;2.学习如何在三种模型之间相互转换;3.学习如何用SIMULINK 仿真工具建模。

二、相关知识1.传递函数模型设连续系统的传递函数为:nn n n m m m m a s a s a s a b s b s b s b s den s num s G ++++++++==----11101110)()()( 设离散系统的传递函数为:n n n n m m m m a z a z a z a b z b z b z b z den z num z G ++++++++==----11101110)()()( 则在MATLAB 中,都可直接用分子/分母多项式系数构成的两个向量num 与den 构成的矢量组[num ,den ]表示系统,即num =],,,[10m b b b den =],,,[10n a a a建立控制系统的传递函数模型(对象)的函数为tf (),调用格式为:sys=tf (num ,den )sys=tf (num ,den ,Ts)sys=tf(othersys)sys=tf (num ,den )返回的变量sys 为连续系统的传递函数模型。

sys=tf (num ,den ,Ts)返回的变量sys 为离散系统的传递函数模型,Ts 为采样周期,当Ts=-1或Ts=[]时,系统的采样周期未定义。

sys=tf(othersys)将任意的控制系统对象转换成传递函数模型。

离散系统的传递函数的表达式还有一种表示为1-z 的形式(即DSP 形式),转换为DSP 形式的函数命令为filt(),调用格式为:sys=filt(num ,den )sys=filt(num ,den ,Ts)sys=filt(num ,den )函数用来建立一个采样时间未指定的DSP 形式传递函数。

sys=filt(num ,den ,Ts)函数用来建立一个采样时间为Ts 的DSP 形式传递函数。

控制系统MATLAB仿真实验指导书

则因输入命令axis(V),式中V是一个四元向量。axis(V)把坐标轴定在规定的范围内。
2.4 MATLAB
下面的实验中,我们将较为详细地学习使用MATLAB的语言基础和控制系统仿真,但是MATLAB的内容及其丰富,在学习和使用中,可以利用MATLAB的联机帮助功能。
MATLAB的联机帮助既可以由help命令来直接获得,又可以由MATLAB图形界面下的Help菜单来查询,还可以用lookfor命令查询有关的关键词:
r
绿色
g
蓝色
b
白色
w
无色
i
参见plot的帮助文件。
7.自动绘图算法
在MATLAB中,图形是自动定标的。在另一幅图形画出之前,这幅图形作为现行图将保持不变,但是在另一幅图形画出后,原图形将被删除,坐标轴自动地重新定标。
8.手工坐标定标
如果需要在下列语句指定的范围内绘制曲线:
V=[x-min x-max y-min y-max]
>>help plot
>>lookfor plot
3.1
一、实验目的
1.学习了解MATLAB语言环境;
2.练习MATLAB命令的基本操作;
3.练习m文件的基本操作。
二、实验内容
1. 学习了解MATLAB语言环境
(1)MATLAB语言操作界面
开机双击相应图标即可进入MATLAB命令窗口,如图1.1所示。在命令提示符位置键入命令,完成下述练习。
图1.1MATLAB窗口
(2)练习DOS相似命令
MATLAB语言与DOS操作系统有如下相似的命令,在操作界面上练习这些命令。
dirdir('c:\windows')
typetypeexpm.m

控制系统仿真与工具(matlab)实验指导书

控制系统仿真与工具实验指导书目录实验一熟悉MATLAB语言工作环境和特点 (1)实验二图形绘制与修饰 (4)实验三系统的时间响应分析 (8)实验四系统的时间响应分析 (12)实验五SIMULINK仿真基础 (14)实验一熟悉MATLAB语言工作环境和特点一、实验目的通过实验使学生熟悉MA TLAB语言的工作环境,并了解MATLAB语言的特点,掌握其基本语法。

二、实验设备PC机MATLAB应用软件三、实验内容本实验从入门开始,使学生熟悉MA TLAB的工作环境,包括命令窗、图形窗和文字编辑器、工作空间的使用等。

1、命令窗(1)数据的输入打开MATLAB后进入的是MA TLAB的命令窗,命令窗是用户与MATLAB做人机对话的主要环境。

其操作提示符为“》”。

在此提示下可输入各种命令并显示出相应的结果,如键入:x1=sqrt(5),x2=1.35,y=3/x2显示结果为:x1=2.2361x2=1.3500y=2.2222上命令行中两式之间用逗号表示显示结果,若用分号,则只运行而不用显示运行结果。

如键入A=[1,2,3;4,5,6;7,8,9],则显示为说明:●直接输入矩阵时,矩阵元素用空格或逗号分隔,矩阵行用分号相隔,整个矩阵放在方括号中。

注意:标点符号一定要在英文状态下输入。

●在MA TLAB中,不必事先对矩阵维数做任何说明,存储时自动配置。

●指令执行后,A被保存在工作空间中,以备后用。

除非用户用clear指令清除它,或对它重新赋值。

●MATLAB对大小写敏感。

(2)数据的显示在MA TLAB工作空间中显示数值结果时,遵循一定的规则,在缺省的情况下,当结果是整数,MATLAB将它作为整数显示;当结果是实数,MATLAB以小数点后4位的精度近似显示。

如果结果中的有效数字超出了这一范围,MATLAB以类似于计算器的计算方法来显示结果。

也可通过键入适当的MA TLAB命令来选择数值格式来取代缺省格式。

如键入format bank命令(数据格式显示命令,数据小数部分以两个十进制小数表示),若在前面键入的三个表达式前键入上述格式命令,则显示结果为:x1=2.24x2=1.35y=2.22前面显示的结果为默认显示结果,它的小数部分为四位数,还有format compact 它以紧凑格式显示结果,默认显示方式为稀疏格式。

MATLAB与系统仿真实验指导书

实验指导书(标准格式)《MATLAB与系统仿真》实验指导书课程编号:ME2121025课程名称:MATLAB与系统仿真适应专业:电气工程及其自动化课程类别:限选实验教学种类:上机课程总学时:30+6实验学时:6执笔人:勾燕洁西安电子科技大学机电工程学院2006 年9 月目录实验一MA TLAB的基本使用及矩阵操作------------------- 3 实验二MA TLAB编程与图形处理---------------------------- 9 实验三MA TLAB中的数值运算与系统仿真---------------- 15实验一MA TLAB的基本使用及矩阵操作一、实验目的和要求熟悉MA TLAB的界面和基本操作,掌握矩阵的建立方法及各种运算。

1.熟悉MA TLAB软件的界面和帮助系统。

2.掌握MA TLAB软件中关于矩阵建立、矩阵初等变换以及矩阵算术、关系、逻辑运算的各种命令。

3.掌握MA TLAB软件中M函数和M文件的编写,以及程序结构与控制,学会编写一般程序。

二、实验内容1.启动与退出2.数、数组、矩阵的输入3.MA TLAB的基本命令4.矩阵大小的测试5.矩阵元素的操作6.特殊矩阵的产生7.矩阵的算术运算8.矩阵的关系运算9.矩阵的逻辑运算10.常用函数三、实验仪器、设备(软、硬件)及仪器使用说明PC机一台,注意正确开、关机及打开软件。

四、实验原理无五、实验方法与步骤1.启动与退出双击MA TLAB图标,进入MA TLAB命令窗口,即可输入命令,开始运算;观察各个窗口以及菜单;单击File菜单中的Exit,或使用MATLAB的Exit命令退出。

2.数、数组、矩阵的输入(1)数的输入>>a=25>>b=3-9i问题1.1:输入“>>a=25;”,结果有什么区别?(2)数组的输入>>c=[1,2,3,4;5,6,7,8]>>d=0:1:10>>e=linspace(1,11,6)问题1.2:体会以上输入方法有什么区别和联系?若A为在0~2 之间均匀分布的22个数据,B=(1.3,2.5,7.6,2,-3),C=(23,20,17,14,11,8,5,2),各用何种方法输入较简单?(3)矩阵的输入>>A=[1,2,3;4,5,6;7,8,9]问题1.3:输入“A=[1 2 3;4 5 6;7 8 9]”,结果相同吗?MA TLAB中对变量名有什么规定?3.MA TLAB的基本命令(1)变量的存储与调用>>save data a b c>>load data(2)MA TLAB命令编辑使用↑键和↓键查看或执行以前执行过的命令(3)显示格式设定>>d=1223.58>>format short e>>d>>format rat>>d(4)工作空间管理>>who a b>>whos>>workspace>>path>>clear %注意工作空间的变化(5)窗口清屏命令>>clc问题1.4:执行命令“clc”和“clear”结果有何不同?(6)联机求助>>help eig>>demo4.矩阵大小的测试(1)>>f=length(A)(2)>>[n,m]=size(A)(3)>>g=nidms(A)问题1.5:以上三个命令有什么区别?5.矩阵元素的操作>>A (1,:)>>A ([1,3],:)>>A (2:3,1:2)问题1.6:输入“A (2,3)”,结果任何,输入“A (6)”,结果相同吗?>>A ([1,3],:)=A ([3,1],:)>>A (2,:)=4>>A (find(A==4))=0>>A (3,:)=[]>>reshape (A,3,2)问题1.7:如何将A的第2列和第3列互换?输入“reshape (2,4)”能否得到正确结果?>>A (4,5)=2>>[A(1:3,2:4),A(2:4,1:3);A,A(:,2)]>>diag(A,1)>>tril(A,2)>>triu(A,2)>>flipud(A)>>fliplr(A)>>rot90(A)问题1.8:输入“rot90(A,2)”和“rot90(A,-2)”结果有区别吗?6.特殊矩阵的产生>>B=eye(5)>>C=ones(2,3)>>D=zeros(3,2)>>E=rand(2,5)>>randn(3,4)>>vander([2 3 9])问题1.9:产生一个在[10,20]内均匀分布的4阶随机矩阵。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
相关文档
最新文档