资源描述：
《系统仿真实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、中南大学系统仿真实验报告指导老师成新明实验者冯虎学号专业班级电气1304班实验日期2015.4.20学院信息科学与工程学院30目录实验一MATLAB中矩阵与多项式的基本运算.................................................3实验二MATLAB绘图命令..................................................7实验三MATLAB程序设计..................................................9实验四MATLAB的符号计算与SIMU
2、LINK的使用.................................................13实验五MATLAB在控制系统分析中的应用.................................................17实验六连续系统数字仿真的基本算法.................................................3030实验一MATLAB中矩阵与多项式的基本运算一、实验任务1．了解MATLAB命令窗口和程序文件的调用。2．熟悉如下MATLAB的基本运算：①矩阵的产生、数据的输入、相关元素
3、的显示；②矩阵的加法、乘法、左除、右除；③特殊矩阵：单位矩阵、“1”矩阵、“0”矩阵、对角阵、随机矩阵的产生和运算；④多项式的运算：多项式求根、多项式之间的乘除。二、基本命令训练1．eye(m)m=3;eye(m)ans=1000100012．ones(n)、ones(m,n)n=1;m=2;ones(n)ones(m,n)ans=1ans=113．zeros(m,n)m=1,n=2;zeros(m,n)30m=1ans=004．rand(m,n)m=1;n=2;rand(m,n)ans=0.81470.90585．diag(v)v=[123];diag(v)
4、ans=1000200036．AB、A/B、inv(A)*B、B*inv(A)A=[12;34];B=[56;78];a=ABb=A/Bc=inv(A)*Bd=B*inv(A)a=-3-445b=3.0000-2.00002.0000-1.0000c=-3.0000-4.00004.00005.000030d=-1.00002.0000-2.00003.00007．roots(p)symsx;a=3*x.^3+2*x+1;p=[3,0,2,1];roots(p)ans=0.2012+0.8877i0.2012-0.8877i-0.40238．polyA=[
5、12;34];poly(A)ans=1.0000-5.0000-2.00009．conv、deconvA=[12];B=[56];a=conv(A,B)b=deconv(A,B)a=51612b=0.200010．A*B与A.*B的区别A=[12];B=[56]';a=A*BA=[12];B=[56];b=A.*Ba=3017b=51211．who与whos的使用A=[12;34];whowhosYourvariablesare:ANameSizeBytesClassAttributesA2x232double12．disp、size(a)、length(a)
6、的使用a='ABCDEF';disp(a)a=[1234];B=size(a)C=length(a)ABCDEFB=14C=4实验二MATLAB绘图命令一、实验任务熟悉MATLAB基本绘图命令，掌握如下绘图方法：1．坐标系的选择、图形的绘制；2．图形注解（题目、标号、说明、分格线）的加入；3．图形线型、符号、颜色的选取。二、基本命令训练301．plot2．loglog3．semilogx4．semilogy5．polar6．title7．xlabel8．ylabel9．text10．grid11．bar12．stairs13．contour三、实验举例1．t=
7、[0:pi/360:2*pi*22/3];x=93*cos(t)+36*cos(t*4.15);y=93*sin(t)+36*sin(t*4.15);plot(y,x),grid;%绘制二维坐标网格图2．t=0:0.05:100;x=t;y=2*t;z=sin(2*t);plot3(x,y,z,'r-.')%绘制三维坐标图3．t=0:pi/20:2*pi;y=sin(x);stairs(x,y)%绘制阶梯图304．th=[pi/200:pi/200:2*pi]';r=cos(2*th);polar(th,r),grid%在网格里画极坐标图5．th=[0:pi/
8、10:2*pi];x=exp(j*th
