武汉工程大学实验报告数字仿真分析.docVIP

实验报告 专业班级： 姓名 学号： 课程名称：控制系统数字仿真与CAD 实验成绩： 指导老师： 实验名称：Matlab的基本操作与编程 一、实验目的： 1）熟悉MATLAB软件的运行环境和基本操作 2）掌握MATLAB矩阵的输入方式、元素的提取与组合 3）掌握数值运算。 4）掌握MATLAB软件的绘图功能 5）掌握M函数的编写。 二、实验内容： 1）启动MATLAB软件，观察其界面组成及操作方法，了解各部分的功能 2）使用基本的MATLAB命令，并观察记录执行结果 帮助、查询信息类命令：Demo、help、who、whos 显示、记录格式等命令：clc、clear、format 尝试一下其他的命令（dos命令） 3）生成一个5阶魔方矩阵，并提取其第（3、4、5）行，第（2、3、4）列构成的新的矩阵 a=magic(5) a = 17 24 1 8 15 23 5 7 14 16 4 6 13 20 22 10 12 19 21 3 11 18 25 2 9 a(3:5,2:4) ans = 6 13 20 12 19 21 18 25 2 5）用命令行方式求解下式的值 （提示使用syms x 定义一个符号，使用limit函数） syms x limit((cos(x)-exp(x^2/2)/2)/4,x,0) ans = 1/8 6）MATLAB的绘图 二维绘图命令plot：画出在上的图形 x=0:0.01:pi*2; y=sin(x); plot(x,y) 三维绘图命令plot3： 画出三维螺旋线，的图形． mesh命令：绘制，在，区间的曲面 t=0:0.01:pi*4; x=sin(t);y=cos(t); z=t; plot3(x,y,z) x=-5:0.1:5; y=x; [X,Y]=meshgrid(x,y); Z=exp(-(X.^2+Y.^2)/2^2); mesh(X,Y,Z) 7）编写M函数 利用程序流程控制语句编写一个函数myfactorial（n），实现n！（阶乘）。要求使用help命令可以列出相关的帮助信息。 function[x]=myfactorial(n) %the function is to slove n! x=1; for a=1:n; x=x*a; end help myfactorial the function is to slove n! myfactorial(8) ans = 40320三、实验结果及分析 武汉工程大学实验报告 专业班级：楼宇自动化01班 姓名：刘斐 学号：1104170109 课程名称：控制系统数字仿真与CAD 实验成绩： 指导老师：王利恒 实验名称：典型闭环系统的数字仿真及计算机解题 一、实验目的： 1）熟悉典型闭环的仿真过程 2）掌握MATLAB编程实现典型闭环环节仿真 3）利用典型闭环环节仿真程序解题。 4）掌握MATLAB下对控制系统进行时域、频域和根轨迹的分析 二、实验内容： 1）编写典型环节阶跃响应函数 典型环节冲击响应函数function [yout,t] = my_step(num,den,v,t0,tf,h,R,n) 输入参数： num：传递函数的分子系数向量 den：传递函数的分母系数向量 v：反馈比例系数 t0：仿真起始时间 tf：仿真终止时间 h：仿真步长 R:阶跃幅值 n:系统阶次 输出参数： yout：响应输出 t：时间向量 function [yout,t] = my_step(b,a,X0,V,T0,Tf,h,R,n) b=b/a(1);a=a/a(1);A=a(2:n+1); A=[rot90(rot90(eye(n-1,n)));-fliplr(A)]; B=[zeros(1,n-1),1]; m1=length(b); C=[fliplr(b),zeros(1,n-m1)]; Ab=A-B*C*V; X=X0;y=0;t=T0; N=round((Tf-T0)/h); for i=1:N; K1=Ab*X+B*R; K2=Ab*(X+h*K1/2)+B*R; K3=Ab*(X+h*K2/2)+B*R; K4=Ab*(X+h*K