系统工程实验报告 学院：管工学院 班级:工业工程102班 姓名：管华同 学号：109094042
实验一：解释结构模型 一、实验目的： 熟悉EXCEL，掌握解释结构模型规范方法。 二、实验内容： 1.已知可达矩阵如下表1 12345678 111010000 201000000 311110000 401010000 501011000 601011111 701011011 800000001 2. EXCEL中对错误！未找到引用源。中的可达矩阵用实用方法建立其递阶结构模型。（1）对可达矩阵进行缩减，得到缩减矩阵 12345678 111010000 201000000 311110000 401010000 501011000 601011111 701011011 800000001 （2）按小到大给每行排序 1 2 3 4 5 6 7 8 每行的和 2 0 1 0 0 0 0 0 0 1 8 0 0 0 0 0 0 0 1 1 4 0 1 0 1 0 0 0 0 2 1 1 1 0 1 0 0 0 0 3 5 0 1 0 1 1 0 0 0 3 3 1 1 1 1 0 0 0 0 4 7 0 1 0 1 1 0 1 1 5 6 0 1 0 1 1 1 1 1 6
（3）调整行列构成对角单位矩阵 2 8 4 1 5 3 7 6 每行的和 2 1 0 0 0 0 0 0 0 1 8 0 1 0 0 0 0 0 0 1 4 1 0 1 0 0 0 0 0 2 1 1 0 1 1 0 0 0 0 3 5 1 0 1 0 1 0 0 0 3 3 1 0 1 1 0 1 0 0 4 7 1 1 1 0 1 0 1 0 5 6 1 1 1 0 1 0 1 1 6 （4）画出递阶结构有向图 28 4 15 37 6（4）递阶结构模型完成。第一级第五级第二级 第三级第四级
实验3 典型闭环系统的阶跃响应的仿真 一、实验目的 1. 了解MATLAB 在仿真中的具体应用 2. 熟悉MATLAB 语言环境 3. 掌握M 文件的应用 二、实验步骤 1. 对控制系统的典型结构形式二次模型化,经一定方式把数学模型转化为便于在计算 机上运行的表达形式。 2. 讨论采用数值积分法求解系统响应的仿真程序实现,绘制仿真框图。 x 0t f t j k 1 1,++k k y x f t t =图3-1 典型闭环系统的仿真程序框图
3.编写MATLAB程序语句,实现对典型闭环系统的阶跃响应的仿真 Filename:sa.m 1)输入数据 a=[a0,a1,…,an]; %% n+1维分母系数向量 b=[b0,b1,…,bm]; %% m+1维分子系数向量 X0=[x0,x1,…,xn]; %% 状态向量初值 V=V0; %% 反馈系数 n=n0; %% 系统阶次 T0=t0; %% 起始时间 Tf=tf; %% 终止时间 h=h0; %% 计算步长 R=r; %% 阶跃输入函数的幅值 2)形成开、闭环系数阵 b=b/a(1); a=a/a(1); A=a(2:n+1); %% 首一化处理 A=[rot90(rot90(eye(n-1,n)));-fliplr(A)]; %%形成能控标准型A阵 B=[zeros(1,n-1),1]'; %%形成输入阵B m1=length(b); %%分子向量维数M+1 C=[fliplr(b),zeros(1,n-m1)]; %%形成输出阵C Ab=A-B*C*V; %%形成闭环系数阵 X=X0'; y=0; t=T0; %%设初值，准备开始递推运算 3）运算求解 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*K3)+B*R; %%求各次斜率K X=X+h*(K1+2*K2+2*K3+K4)/6; %%求状态 y=[y,C*X]; %%求输出并以向量形式保存 t=[t,t(i)+h]; end
《控制系统仿真与CAD》 实验课程报告
一、实验教学目标与基本要求 上机实验是本课程重要的实践教学环节。实验的目的不仅仅是验证理论知识，更重要的是通过上机加强学生的实验手段与实践技能，掌握应用 MATLAB/Simulink 求解控制问题的方法，培养学生分析问题、解决问题、应用知识的能力和创新精神，全面提高学生的综合素质。 通过对MATLAB/Simulink进行求解，基本掌握常见控制问题的求解方法与命令调用，更深入地认识和了解MATLAB语言的强大的计算功能与其在控制领域的应用优势。 上机实验最终以书面报告的形式提交，作为期末成绩的考核内容。 二、题目及解答 第一部分：MATLAB 必备基础知识、控制系统模型与转换、线性控制系统的计算机辅助分析 1. >>f=inline('[-x(2)-x(3);x(1)+a*x(2);b+(x(1)-c)*x(3)]','t','x','flag','a','b','c');[t,x]=ode45( f,[0,100],[0;0;0],[],0.2,0.2,5.7);plot3(x(:,1),x(:,2),x(:,3)),grid,figure,plot(x(:,1),x(:,2)), grid
2. >>y=@(x)x(1)^2-2*x(1)+x(2);ff=optimset;https://www.doczj.com/doc/ba4345377.html,rgeScale='off';ff.TolFun=1e-30;ff.Tol X=1e-15;ff.TolCon=1e-20;x0=[1;1;1];xm=[0;0;0];xM=[];A=[];B=[];Aeq=[];Beq=[];[ x,f,c,d]=fmincon(y,x0,A,B,Aeq,Beq,xm,xM,@wzhfc1,ff) Warning: Options LargeScale = 'off' and Algorithm = 'trust-region-reflective' conflict. Ignoring Algorithm and running active-set algorithm. To run trust-region-reflective, set LargeScale = 'on'. To run active-set without this warning, use Algorithm = 'active-set'. > In fmincon at 456 Local minimum possible. Constraints satisfied. fmincon stopped because the size of the current search direction is less than twice the selected value of the step size tolerance and constraints are satisfied to within the selected value of the constraint tolerance.
信号与系统仿真实验报告1．实验目的 了解MATLAB的基本使用方法和编程技术，以及Simulink平台的建模与动态仿真方法，进一步加深对课程内容的理解。 2．实验项目 信号的分解与合成，观察Gibbs现象。 信号与系统的时域分析，即卷积分、卷积和的运算与仿真。 信号的频谱分析，观察信号的频谱波形。 系统函数的形式转换。 用Simulink平台对系统进行建模和动态仿真。 3．实验内容及结果 3.1以周期为T，脉冲宽度为2T1的周期性矩形脉冲为例研究Gibbs现象。 已知周期方波信号的相关参数为：x(t)=∑ak*exp(jkω),ω=2*π/T,a0=2*T1/T,ak=sin(kωT1)/kπ。画出x(t)的波形图（分别取m=1,3,7,19,79,T=4T1）,观察Gibbs现象。 m=1; T1=4; T=4*T1;k=-m:m; w0=2*pi/T; a0=2*T1/T; ak=sin(k*w0*T1)./(k*pi); ak(m+1)=a0; t=0:0.1:40; x=ak*exp(j*k'*w0*t); plot(t,real(x)); 3.2求卷积并画图 (1)已知：x1(t)=u(t-1)-u(t-2), x2(t)=u(t-2)-u(t-3)求：y(t)=x1(t)*x2(t)并画出其波形。 t1=1:0.01:2; f1=ones(size(t1)); f1(1)=0; f1(101)=0; t2=2:0.01:3; f2=ones(size(t2)); f2(1)=0; f2(101)=0; c=conv(f1,f2)/100;
t3=3:0.01:5; subplot(311); plot(t1,f1);axis([0 6 0 2]); subplot(312); plot(t2,f2);axis([0 6 0 2]); subplot(313); plot(t3,c);axis([0 6 0 2]); (2）已知某离散系统的输入和冲击响应分别为：x[n]=[1,4,3,5,1,2,3,5], h[n]=[4,2,4,0,4,2].求系 统的零状态响应，并绘制系统的响应图。 x=[1 4 3 5 1 2 3 5]; nx=-4:3; h=[4 2 4 0 4 2]; nh=-3:2; y=conv(x,h); ny1=nx(1)+nh(1); ny2=nx(length(nx))+nh(length(nh)); ny=[ny1:ny2]; subplot(311); stem(nx,x); axis([-5 4 0 6]); ylabel('输入') subplot(312); stem(nh,h); axis([-4 3 0 5]); ylabel('冲击效应') subplot(313); stem(ny,y); axis([-9 7 0 70]); ylabel('输出'); xlabel('n'); 3.3 求频谱并画图 (1) 门函数脉冲信号x1(t)=u(t+0.5)-u(t-0.5) N=128;T=1; t=linspace(-T,T,N); x=(t>=-0.5)-(t>=0.5); dt=t(2)-t(1); f=1/dt; X=fft(x); F=X(1:N/2+1); f=f*(0:N/2)/N; plot(f,F)
过程控制系统Matlab/Simulink 仿真实验 实验一 过程控制系统建模 ............................................................................................................. 1 实验二 PID 控制 ............................................................................................................................. 2 实验三 串级控制 ............................................................................................................................. 6 实验四 比值控制 ........................................................................................................................... 13 实验五 解耦控制系统 . (19) 实验一 过程控制系统建模 指导内容：（略） 作业题目一： 常见的工业过程动态特性的类型有哪几种？通常的模型都有哪些？在Simulink 中建立相应模型，并求单位阶跃响应曲线。 作业题目二： 某二阶系统的模型为2 () 22 2n G s s s n n ?ζ??= ++，二阶系统的性能主要取决于ζ，n ?两个参数。试利用Simulink 仿真两个参数的变化对二阶系统输出响应的影响，加深对二阶 系统的理解，分别进行下列仿真： （1）2n ?=不变时，ζ分别为0.1, 0.8, 1.0, 2.0时的单位阶跃响应曲线； （2）0.8ζ=不变时，n ?分别为2, 5, 8, 10时的单位阶跃响应曲线。
实验一 Matlab使用方法和程序设计 一、实验目的 1、掌握Matlab软件使用的基本方法； 2、熟悉Matlab的数据表示、基本运算和程序控制语句 3、熟悉Matlab绘图命令及基本绘图控制 4、熟悉Matlab程序设计的基本方法 二、实验内容 1、帮助命令 使用help命令，查找sqrt（开方）函数的使用方法； 2、矩阵运算 （1）矩阵的乘法 已知A=[1 2;3 4]; B=[5 5;7 8]; 求A^2*B （2）矩阵除法 已知A=[1 2 3;4 5 6;7 8 9]; B=[1 0 0;0 2 0;0 0 3]; A\B,A/B （3）矩阵的转置及共轭转置 已知A=[5+i,2-i,1;6*i,4,9-i]; 求A.', A' （4）使用冒号选出指定元素 已知：A=[1 2 3;4 5 6;7 8 9]; 求A中第3列前2个元素；A中所有列第2，3行的元素； （5）方括号[] 用magic函数生成一个4阶魔术矩阵，删除该矩阵的第四列 3、多项式 （1）求多项式p(x) = x3 - 2x - 4的根 （2）已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ， 求矩阵A的特征多项式; 求特征多项式中未知数为20时的值； 4、基本绘图命令 （1）绘制余弦曲线y=cos(t)，t∈[0，2π] （2）在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t∈[0，2π] 5、基本绘图控制 绘制[0，4π]区间上的x1=10sint曲线，并要求： （1）线形为点划线、颜色为红色、数据点标记为加号； （2）坐标轴控制：显示范围、刻度线、比例、网络线 （3）标注控制：坐标轴名称、标题、相应文本； 6、基本程序设计 （1）编写命令文件：计算1+2+?+n<2000时的最大n值； （2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。 三、预习要求 利用所学知识，编写实验内容中2到6的相应程序，并写在预习报告上。
Matlab通信原理仿真 学号： 2142402 姓名：圣斌
实验一Matlab 基本语法与信号系统分析 一、实验目的： 1、掌握MATLAB的基本绘图方法； 2、实现绘制复指数信号的时域波形。 二、实验设备与软件环境： 1、实验设备：计算机 2、软件环境：MATLAB R2009a 三、实验内容： 1、MATLAB为用户提供了结果可视化功能，只要在命令行窗口输入相应的命令，结果就会用图形直接表示出来。 MATLAB程序如下： x = -pi::pi; y1 = sin(x); y2 = cos(x); %准备绘图数据 figure(1); %打开图形窗口 subplot(2,1,1); %确定第一幅图绘图窗口 plot(x,y1); %以x，y1绘图 title('plot(x,y1)'); %为第一幅图取名为’plot(x,y1)’ grid on; %为第一幅图绘制网格线 subplot(2,1,2) %确定第二幅图绘图窗口 plot(x,y2); %以x，y2绘图 xlabel('time'),ylabel('y') %第二幅图横坐标为’time’,纵坐标为’y’运行结果如下图： 2、上例中的图形使用的是默认的颜色和线型，MATLAB中提供了多种颜色和线型，并且可以绘制出脉冲图、误差条形图等多种形式图： MATLAB程序如下： x=-pi:.1:pi; y1=sin (x); y2=cos (x); figure (1); %subplot (2,1,1); plot (x,y1); title ('plot (x,y1)'); grid on %subplot (2,1,2); plot (x,y2);
电子科技大学学院学生实验报告 学院：机电工程学院专业:17自动化课程名称：自动控制原理实验与仿真
JET 性能指标 Pole Dampi ng (rad/sec on ds) Freque ncy (sec on ds) Time Con sta nt 1.50e+00 + 9.89e+00i 1.50e-01 1.00e+01 6.67e-01 ?1.50e+00 - 9.89e+00i 1.50e-01 1.00e+01 6.67e-01 2当=0, 0.25 , 0.5 , 0.75 , 1 , 1.25时，对应系统的闭环极点和自然振荡频率n见表，编程求取对应系统的阶跃响应曲线， 并分析n—定时，变化对系统性能的影响。。 曲线:
.w 结论：可见当n/(rad/s) —定时，系统随着阻尼比E增大，闭环极点的实部在S左半平面的位置更加远离原点，虚部减小到0，超调量减小，调节时间缩短，稳定性更好。 3. 0.25, n 10,30,50时，对应点的单位阶跃响应曲线并分析不变，n对系统性能的影响 曲线： 结论：可见当E—定时，随着n/(rad/s)增大，系统响应加速，振荡频率增大，系统调整时间缩短，但是超
调量没有变化。 3.试做出以下系统的单位阶跃响应'并与原系统G(s)= s2—2s—I。的阶跃响应曲线比较，做出实验结果分析? 1)系统分别增加零点z -5z -2, G(s) 22(S 5)G(s) 25(S 2) s2 2s 10 s2 2s 10 代码及曲线： 代码： sys=tf(10,[1 2 10]); step(sys) hold on sysc=tf([2,10],[1 2 10]); step(sysc) hold on sysx=tf([5,10],[1 2 10]); step(sysx) hold off title('单位阶跃系统增加零点比较’)； lab仁'增加零点-2';text(1,1.8,lab1) Iab2='增加零点-5';text(0.25,1.1,lab2) lab3='原系统：text(1.5,1.3,lab3) 曲线:
中南大学系统仿真实验报告 指导老师胡杨 实验者 学号 专业班级 实验日期 2014.6.4 学院信息科学与工程学院
目录 实验一MATLAB中矩阵与多项式的基本运算 (3) 实验二MATLAB绘图命令 (7) 实验三MATLAB程序设计 (9) 实验四MATLAB的符号计算与SIMULINK的使用 (13) 实验五MATLAB在控制系统分析中的应用 (17) 实验六连续系统数字仿真的基本算法 (30)
实验一MATLAB中矩阵与多项式的基本运算 一、实验任务 1．了解MATLAB命令窗口和程序文件的调用。 2．熟悉如下MATLAB的基本运算： ①矩阵的产生、数据的输入、相关元素的显示； ②矩阵的加法、乘法、左除、右除； ③特殊矩阵：单位矩阵、“1”矩阵、“0”矩阵、对角阵、随机矩阵的产生和运算； ④多项式的运算：多项式求根、多项式之间的乘除。 二、基本命令训练 1．eye(m) m=3; eye(m) ans = 1 0 0 0 1 0 0 0 1 2．ones(n)、ones(m,n) n=1;m=2; ones(n) ones(m,n) ans = 1 ans = 1 1
3．zeros(m,n) m=1,n=2; zeros(m,n) m = 1 ans = 0 0 4．rand(m,n) m=1;n=2; rand(m,n) ans = 0.8147 0.9058 5．diag(v) v=[1 2 3]; diag(v) ans = 1 0 0 0 2 0 0 0 3 6．A\B 、A/B、inv(A)*B 、B*inv(A) A=[1 2;3 4];B=[5 6;7 8]; a=A\B b=A/B c=inv(A)*B d=B*inv(A) a = -3 -4 4 5 b = 3.0000 -2.0000 2.0000 -1.0000
系统工程仿真实验报告 姓名：_蒋智颖_ 学号：_110061047_ 成绩：___________ 实验一：基于VENSIM的系统动力学仿真 一、实验目的 VENSIM是一个建模工具，可以建立动态系统的概念化的，文档化的仿真、分析和优化模型。PLE（个人学习版）是VENSIM的缩减版，主要用来简单化学习动态系统，提供了一种简单富有弹性的方法从常规的循环或储存过程和流程图建立模型。本实验就是运用VENSIM进行系统动力学仿真，进一步加深对系统动力学仿真的理解。 二、实验软件 VENSIM PLE 三、原理 1、在VENSIM中建立系统动力学流图； 2、写出相应的DYNAMO方程； 3、仿真出系统中水准变量随时间的响应趋势； 四、实验内容及要求 某城市国营和集体服务网点的规模可用SD来研究。现给出描述该问题的DYNAMO方程及其变量说明。 L S·K=S·J+DT*NS·JK N S=90 R NS·KL=SD·K*P·K/(LENGTH-TIME·K) A SD·K=SE-SP·K C SE=2 A SP·K=SR·K/P·K A SR·K=SX+S·K C SX=60 L P·K=P·J+DT*NP·JK N P=100 R NP·KL=I*P·K C I=0.02 其中：LENGTH为仿真终止时间、TIME为当前仿真时刻，均为仿真控制变量；S为个体服务网点数(个)、NS为年新增个体服务网点数(个/年)、SD为实际千人均服务网点与期望差(个/千人)、SE为期望的千人均网点数、SP为的千人均网点数(个/千人)、SX为非个体服务网点数(个)、SR为该城市实际拥有的服务网点数(个)、P为城市人口数(千人)、NP为年新
实验一 控制系统的阶跃响应 一、实验目的 1. 掌握控制系统多项式模型和零极点模型的建立方法及它们之间的相互转换。 2．观察学习控制系统的单位阶跃响应。 3．记录单位阶跃响应曲线。 4．掌握时间响应分析的一般方法。 5．分析系统阶跃响应曲线与传递函数参数的对应关系。 二、实验设备 PC 机，MATLAB 仿真软件。 三、实验内容 1．作以下二阶系统的单位阶跃响应曲线 10 210)(2++=s s s G 2．分别改变该系统的ζ和n ω，观察阶跃响应曲线的变化。 3．作该系统的脉冲响应曲线。 四、实验原理 1. 建立系统模型 在MATLAB 下，系统数学模型有三种描述方式，在本实验中只用到多项式模型和零极点模型。 （1）多项式模型 num 表示分子多项式的系数，den 表示分母多项式的系数，以行向量的方式 输入。例如，程序为 num=[0 1 3]； %分子多项式系数 den=[1 2 2 1]； %分母多项式系数 printsys (num, den) %构造传递函数并显示 （2）零极点模型 z 表示零点，p 表示极点，以行向量的方式输入，k 表示增益。例如，程序为 k=2； %赋增益值，标量 z=[1]； %赋零点值，向量 p=[-1 2 -3]； %赋极点值，向量 [num, den]=zp2tf(z, p, k)； %零极点模型转换成多项式模型 printsys(num, den) %构造传递函数并显示 （3）相关MATLAB 函数 函数tf(num, den) 用来建立控制系统的多项式模型； 函数zpk(z, p, k)用来建立控制系统的零极点模型； [num, den]=zp2tf (z, p, k) %零极点模型转换成多项式模型 [z, p, k]=tf2zp (num, den) %多项式模型转换成零极点模型 [num, den]=ord2(ωn , ξ) %用来建立二阶系统标准模型
哈尔滨理工大学实验报告 控制系统仿真 专业：自动化12-1 学号：1230130101 姓名：
一．分析系统性能 课程名称控制系统仿真实验名称分析系统性能时间8.29 地点3# 姓名蔡庆刚学号1230130101 班级自动化12-1 一．实验目的及内容： 1. 熟悉MATLAB软件的操作过程； 2. 熟悉闭环系统稳定性的判断方法； 3. 熟悉闭环系统阶跃响应性能指标的求取。 二．实验用设备仪器及材料： PC, Matlab 软件平台 三、实验步骤 1. 编写MATLAB程序代码； 2. 在MATLAT中输入程序代码，运行程序； 3.分析结果。 四．实验结果分析： 1.程序截图
得到阶跃响应曲线 得到响应指标截图如下
2.求取零极点程序截图 得到零极点分布图 3.分析系统稳定性 根据稳定的充分必要条件判别线性系统的稳定性最简单的方法是求出系统所有极点,并观察是否含有实部大于0的极点，如果有系统不稳定。有零极点分布图可知系统稳定。
二．单容过程的阶跃响应 一、实验目的 1. 熟悉MATLAB软件的操作过程 2. 了解自衡单容过程的阶跃响应过程 3. 得出自衡单容过程的单位阶跃响应曲线 二、实验内容 已知两个单容过程的模型分别为 1 () 0.5 G s s =和5 1 () 51 s G s e s - = + ，试在 Simulink中建立模型，并求单位阶跃响应曲线。 三、实验步骤 1. 在Simulink中建立模型，得出实验原理图。 2. 运行模型后，双击Scope，得到的单位阶跃响应曲线。 四、实验结果 1．建立系统Simulink仿真模型图，其仿真模型为
一、系统描述 1.1.系统背景 本系统将基于下面的卫星屏幕快照创建一个模型。当前道路网区域的两条道路均为双向，每个运动方向包含一条车道。Tapiolavagen路边有一个巴士站，Menninkaisentie路边有一个带五个停车位的小型停车场。 1.2.系统描述 （1）仿真十字路口以及三个方向的道路，巴士站，停车点；添加小汽车、公交车的三维动画，添加红绿灯以及道路网络描述符； （2）创建仿真模型的汽车流程图，三个方向产生小汽车，仿真十字路口交通运行情况。添加滑条对仿真系统中的红绿灯时间进行实时调节。添加分析函数，统计系统内汽车滞留时间，用直方图进行实时展示。 二、仿真目标 1、timeInSystem值：在流程图的结尾模块用函数统计每辆汽车从产生到丢弃的，在系统中留存的时间。 2、p_SN为十字路口SN方向道路的绿灯时间，p_EW为十字路口EW方向道路的绿灯时间。 3、Arrival rate：各方向道路出现车辆的速率（peer hour）。
三、系统仿真概念分析 此交通仿真系统为低抽象层级的物理层模型，采用离散事件建模方法进行建模，利用过程流图构建离散事件模型。 此十字路口交通仿真系统中，实体为小汽车和公交车，可以源源不断地产生；资源为道路网络、红绿灯时间、停车点停车位和巴士站，需要实施分配。系统中小汽车（car）与公共汽车（bus）均为智能体，可设置其产生频率参数，行驶速度，停车点停留时间等。 四、建立系统流程 4.1.绘制道路 使用Road Traffic Library中的Road模块在卫星云图上勾画出所有的道路，绘制交叉口，并在交叉口处确保道路连通。 4.2.建立智能体对象 使用Road Traffic Library中的Car type模快建立小汽车（car）以及公共汽车（bus）的智能体对象。 4.3.建立逻辑 使用Road Traffic Library中的Car source、Car Move To、Car Dispose、
《工程系统建模与仿真》实验报告 姓名XXXXXXX 学号XXXXXXX 班级XXXXXXX 专业XXXXXXX 报告提交日期XXXXXXX
实验一 扭摆法测定物体的转动惯量 一、 实验名称 扭摆法测定物体的转动惯量 二、 同组成员 学号 姓名 XXXXXX XXX XXXXXX XXX XXXXXX XXX XXXXXX XXX XXXXXX XXX XXXXXX XXX XXXXXX XXX XXXXXX XXX XXXXXX XXX 三、 实验器材 1) 转动惯量测试仪 2) 数字式电子台秤 3) 游标卡尺 4) 扭摆及几种有规则的待测转动惯量的物体：金属载物圆盘、塑料圆柱体、 木球、验证转动惯量平行轴定理用的金属细杆，杆上有两块可以自由移动的金属滑块。 四、 实验原理 转动惯量的测量，一般都是使刚体以一定形式运动，通过表征这种运动特征的物理量与转动惯量的关系，进行转换测量。本实验使物体作扭转摆动，由于摆动周期及其它参数的测定计算出物体的转动惯量。 扭摆的构造如图 1-1所示，在垂直轴1上装有一根薄片状的螺旋弹簧2，用以产生恢复力矩。在轴的上方可以装上各种待测物体。垂直轴与支座间装有轴承，以降低摩擦力矩。3为水平仪，用来调整系统平衡。 将物体在水平面内转过一定角度θ后，在弹簧的恢复力矩作用下物体就开始绕垂直轴作周期往返扭转运动。 根据虎克定律，弹簧受扭转而产生的恢复力矩M 与所转过的角度θ成正 比，即：M=-Kθ (1) 上式中，K 为弹簧的扭转常数。 由转动定律M ＝Iβ得：β=M /I (2) 令ω2=K /I ，忽略轴承的摩擦阻力矩，由式(1)、(2)得： 2 22 d K dt I θβθωθ= =-=- 图 1-1 上述方程表示扭摆运动具有角简谐振动的特性，角加速度与角位移成正比， 且方向相反。此方程的解为：θ=Acos (ωt +?)。 式中，A 为谐振动的角振幅，φ为初相位角，ω为角速度，此谐振动的周期
一、摘要 2 二、课程设计任务 3 1．问题描述 3 2．设计要求 3 三、课程设计内容 4 1、系统的模型表示 4 2、利用Matlab进行仿真设计 4 3、利用Simulink进行仿真设计 9 总结与体会 10 参考文献 10
本课题以汽车运动控制系统的设计为应用背景，利用MATLAB语言对其进行设计与仿真.首先对汽车的运动原理进行分析，建立控制系统模型，确定期望的静态指标稳态误差和动态指标搬调量和上升时间，最终应用MATLAB环境下的.m文件来实现汽车运动控制系统的设计。其中.m文件用step函数语句来绘制阶跃响应曲线，根据曲线中指标的变化进行P、PI、PID校正；同时对其控制系统建立Simulink进行仿真且进行PID参数整定。仿真结果表明，参数PID控制能使系统达到满意的控制效果，对进一步应用研究具有参考价值，是汽车运动控制系统设计的优秀手段之一。 关键词：运动控制系统 PID仿真稳态误差最大超调量
一、课程设计任务 1. 问题描述 如下图所示的汽车运动控制系统，设该系统中汽车车轮的转动惯量可以忽略不计，并且假定汽车受到的摩擦阻力大小与汽车的运动速度成正比，摩擦阻力的方向与汽车运动的方向相反，这样，该汽车运动控制系统可简化为一个简单的质量阻尼系统。 根据牛顿运动定律，质量阻尼系统的动态数学模型可表示为： ???==+v y u bv v m 系统的参数设定为：汽车质量m =1000kg ， 比例系数b =50 N ·s/m ， 汽车的驱动力u =500 N 。 根据控制系统的设计要求，当汽车的驱动力为500N 时，汽车将在5秒内达到10m/s 的最大速度。由于该系统为简单的运动控制系统，因此将系统设计成10％的最大超调量和2％的稳态误差。这样，该汽车运动控制系统的性能指标可以设定为： 上升时间：t r <5s ； 最大超调量：σ％<10％； 稳态误差：e ssp <2％。 2.设计要求 1．写出控制系统的数学模型。 2．求系统的开环阶跃响应。 3．PID 控制器的设计 （1）比例（P ）控制器的设计 （2）比例积分（PI ）控制器的设计 （3）比例积分微分（PID ）控制器的设计 利用Simulink 进行仿真设计。 二、课程设计内容 1.系统的模型表示
昆明理工大学电力工程学院学生实验报告 实验课程名称：控制系统仿真实验 开课实验室：年月日
实验一 电路的建模与仿真 一、实验目的 1、了解KCL 、KVL 原理； 2、掌握建立矩阵并编写M 文件； 3、调试M 文件，验证KCL 、KVL ； 4、掌握用simulink 模块搭建电路并且进行仿真。 二、实验内容 电路如图1所示，该电路是一个分压电路，已知13R =Ω，27R =Ω，20S V V =。试求恒压源的电流I 和电压1V 、2V 。 I V S V 1 V 2 图1 三、列写电路方程 （1）用欧姆定律求出电流和电压 （2）通过KCL 和KVL 求解电流和电压
四、编写M文件进行电路求解（1）M文件源程序 （2）M文件求解结果 五、用simulink进行仿真建模（1）给出simulink下的电路建模图（2）给出simulink仿真的波形和数值
六、结果比较与分析
实验二数值算法编程实现 一、实验目的 掌握各种计算方法的基本原理，在计算机上利用MATLAB完成算法程序的编写拉格朗日插值算法程序，利用编写的算法程序进行实例的运算。 二、实验说明 1．给出拉格朗日插值法计算数据表； 2．利用拉格朗日插值公式，编写编程算法流程，画出程序框图，作为下述编程的依据； 3．根据MATLAB软件特点和算法流程框图，利用MATLAB软件进行上机编程； 4．调试和完善MATLAB程序； 5．由编写的程序根据实验要求得到实验计算的结果。 三、实验原始数据 上机编写拉格朗日插值算法的程序，并以下面给出的函数表为数据基础，在整个插值区间上采用拉格朗日插值法计算(0.6) f，写出程序源代码，输出计算结果： 四、拉格朗日插值算法公式及流程框图
《嵌入式操作系统》实验报告 班级计算机 学号 姓名 指导教师庄旭菲
内蒙古工业大学信息工程学院计算机系 2018年6月 实验一 Linux内核移植与编译实验 1. 实验目的 了解 Linux 内核相关知识与内核结构 了解 Linux 内核在 ARM 设备上移植的基本步骤和方法 掌握 Linux 内核裁剪与定制的基本方法 2. 实验内容 分析 Linux 内核的基本结构，了解 Linux 内核在 ARM 设备上移植的一些基本步骤及常识。 学习 Linux 内核裁剪定制的基本配置方法，利用 UP-Magic210 型设备配套 Linux 内核进行自定义功能(如helloworld 显示)的添加，并重新编译内核源码，生成内核压缩文件 zImage，下载到 UP-Magic210 型设备中测试。 3. 实验步骤 实验目录：/UP-Magic210/SRC/kernel/编译内核：在宿主机端为UP-Magic210 设备的Linux 内核编写简单的测试驱动（内核）程序并修改内核目录中相关文件，添加对测试驱动程序的支持。 (1)、使用 vim 编辑器手动编写实验代码
内如如下： #include
控制系统数字仿真实验报告 班级：机械1304 姓名：俞文龙 学号： 0801130801
实验一数字仿真方法验证1 一、实验目的 1.掌握基于数值积分法的系统仿真、了解各仿真参数的影响； 2.掌握基于离散相似法的系统仿真、了解各仿真参数的影响； 3.熟悉MATLAB语言及应用环境。 二、实验环境 网络计算机系统（新校区机电大楼D520），MATLAB语言环境 三实验内容 （一）试将示例1的问题改为调用ode45函数求解，并比较结果。 实验程序如下； function dy = vdp(t,y) dy=[y-2*t/y]; end [t,y]=ode45('vdp',[0 1],1); plot(t,y); xlabel('t'); ylabel('y');
（二）试用四阶RK 法编程求解下列微分方程初值问题。仿真时间2s ，取步长h=0.1。 ?????=-=1 )0(2y t y dt dy 实验程序如下: clear t0=0; y0=1; h=0.1; n=2/h; y(1)=1; t(1)=0; for i=0:n-1 k1=y0-t0^2; k2=(y0+h*k1/2)-(t0+h/2)^2; k3=(y0+h*k2/2)-(t0+h/2)^2;
k4=(y0+h*k3)-(t0+h)^2; y1=y0+h*(k1+2*k2+2*k3+k4)/6; t1=t0+h; y0=y1; t0=t1; y(i+2)=y1; t(i+2)=t1; end y1 t1 figure(1) plot(t,y,'r'); xlabel('t'); ylabel('y'); （三）试求示例3分别在周期为5s的方波信号和脉冲信号下的响应，仿真时间20s，采样周期Ts=0.1。
电子科技大学中山学院学生实验报告 学院: 机电工程学院 专业: 17自动化 课程名称:自动控制原理实验与仿真 班级: 姓名: 学号: 组别: 实验名称:基于MATLAB 的控制系统单位阶跃响应分析 实验时间: 成 绩: 批改时间: 一、 实验目的 (1)学会使用MATLAB 编程绘制控制系统的单位阶跃响应曲线。 (2)研究二阶控制系统中ζ ,ωn 对系统阶跃响应的影响。 (3)掌握准确读取动态特性指标的方法。 (4)分析二阶系统闭环极点与闭环零点对系统动态性能的影响。 二、 实验条件 实验设备:每人一台计算机奔腾系列以上计算机,配置硬盘≥2G,内存≥64M 。 实验软件:WINDOWS 操作系统(WINDOWS XP 或WINDOWS 2000),并安装MATLAB 语言编程环境。 三、实验内容 21001.(),3G s s s = +已知系统开环传递函数为试绘制单位负反馈闭环系统的阶跃响应曲线，并测出动态性能指标。 代码、曲线及性能指标: 代码 sys=tf(100,[1 3 0]); sysc=feedback(sys,1); damp(sysc) step(sysc) 曲线 性能指标 Pole Damping Frequency Time Constant (rad/seconds) (seconds) -1、50e+00 + 9、89e+00i 1、50e-01 1、00e+01 6、67e-01 -1、50e+00 - 9、89e+00i 1、50e-01 1、00e+01 6、67e-01
2.=n n ζωωζ当0，0.25，0.5，0.75，1，1.25时，对应系统的闭环极点 和自然振荡频率见表，编程求取对应系统的阶跃响应曲线，并分析一定时，变化对系统性能的影响。。 ζ 闭环极点 /(/)n rad s ω 阶跃响应曲线 =0ζ j ± 10 等幅振荡 =0.25ζ 2.59.68j -± 10 衰减振动 =0.5ζ 58.66j -± 10 衰减振动 =0.75ζ 7.5 6.61j -± 10 衰减振动 =1ζ 两实重根-10 10 单调上升 =1.25ζ 两不等实根 -5与-20 5,20 单调上升 曲线: 结论:可见当/(/)n rad s ω一定时,系统随着阻尼比ξ增大,闭环极点的实部在S 左半平面的位置更加远离原点,虚部减小到0,超调量减小,调节时间缩短,稳定性更好。 0.25,10,30,50n n ζωζ==3.时，对应点的单位阶跃响应曲线并分析不变，对系统性能的影响。 曲线:
《MATLAB与控制系统仿真》 实验报告 班级： 学号： 姓名： 时间：2013 年 6 月
目录实验一MATLAB环境的熟悉与基本运算（一）实验二MATLAB环境的熟悉与基本运算（二）实验三MATLAB语言的程序设计 实验四MATLAB的图形绘制 实验五基于SIMULINK的系统仿真 实验六控制系统的频域与时域分析 实验七控制系统PID校正器设计法 实验八线性方程组求解及函数求极值
实验一MATLAB环境的熟悉与基本运算（一） 一、实验目的 1．熟悉MATLAB开发环境 2．掌握矩阵、变量、表达式的各种基本运算 二、实验基本原理 1.熟悉MATLAB环境: MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器、文件和搜索路径浏览器。 2.掌握MATLAB常用命令 表1 MATLAB常用命令 变量与运算符 3．1变量命名规则 3．2 MATLAB的各种常用运算符 表3 MATLAB关系运算符 表4 MATLAB逻辑运算符
| Or 逻辑或 ~ Not 逻辑非 Xor逻辑异或 符号功能说明示例符号功能说明示例 ：1:1:4;1:2:11 . ；分隔行.. ，分隔列… （）% 注释 [] 构成向量、矩阵！调用操作系统命令 {} 构成单元数组= 用于赋值 的一维、二维数组的寻访 表6 子数组访问与赋值常用的相关指令格式 三、主要仪器设备及耗材 计算机 四．实验程序及结果 1、新建一个文件夹（自己的名字命名，在机器的最后一个盘符） 2、启动MATLAB，将该文件夹添加到MATLAB路径管理器中。 3、学习使用help命令。