浙江科技学院机电系统开放性实验
设计题目： MATLAB控制机电系统实分析
专业：机械设计制造及其自动化班级：机制 122 学生姓名：王朝旭
指导老师：高俊
2014年05月30日
目录
实验一 MATLAB语言基本命令 (3)
1 实验目的 (3)
2 实验器材 (3)
3 实验内容 (3)
4 实验步骤： (3)
5 实验类容： (3)
实验二 MATLAB控制系统工具箱使用 (6)
1 实验目的 (6)
2 实验器材 (6)
3 实验内容 (6)
4 实验步骤及内容............................... 错误！未定义书签。
实验三 MATLAB数字控制器设计初步 (9)
1 实验目的 (9)
2 实验器材 (9)
3 实验内容 (6)
4 实验步骤及内容 (6)
实验四控制系统数学模型的MATLAB实现 (12)
1 实验目的 (12)
2 实验器材 (12)
3 实验步骤及内容 (12)
总结 (14)
实验一 MATLAB语言基本命令
1 实验目的
1. 掌握科学计算的有关方法，熟悉MATLAB语言及其在科学计算中的运用；
2. 掌握MATLAB的命令运行方式和M文件运行方式；
3. 掌握矩阵在MATLAB中的运用。
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.*B
A^2 A.^2 B/A B./A
(2).输入 C=1:2:20，则 C（i）表示什么？其中 i=1,2,3, (10)
(3).查找已创建变量的信息，删除无用的变量；
(4). 试用 help命令理解下面程序各指令的含义：
clear
t =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 实验步骤：
打开MATLAB程序，将实验内容中的题目依次输入MATLAB中，运行得到并记录结果，最后再对所得结果进行验证。
5 实验报告内容
（1）>> A=[7 1 5;2 5 6;3 1 5]
A =
7 1 5
2 5 6
3 1 5
>> B=[1 1 1; 2 2 2;3 3 3]
B =
1 1 1
2 2 2
3 3 3
>> A(2, 3) “取矩阵A中第二行第三列的元素”
ans =
6
>> A(:,2) “取矩阵A中第二列的元素”
ans =
1
5
1
>> A(3,:) “取矩阵A中第三行的元素”
ans =
3 1 5
>> A(:,1:2:3) “使矩阵的第一列到第二列，第二例到第三列，构成新的2*3的矩阵”ans =
7 5
2 6
3 5
>> A(:,3).*B(:,2) “取矩阵A的三列和B的第二列，对应元素一次相乘”
ans =
5
12
15
>> A(:,3)*B(2,:) “取矩阵A的三列和B的第二行，对应元素进行叉乘”
ans =
10 10 10
12 12 12
10 10 10
>> A*B “矩阵进行叉乘”
ans =
24 24 24
30 30 30
20 20 20
>> A.*B “两矩阵对应元素相乘”
ans =
7 1 5
4 10 12
9 3 15
>> A^2 “矩阵A叉乘上矩阵A”
ans =
66 17 66
42 33 70
38 13 46
>> A.^2 “举证A中的元素翻倍”
ans =
49 1 25
4 2
5 36
9 1 25
>> B/A “矩阵B除以矩阵A”
ans =
0.1842 0.2105 -0.2368
0.3684 0.4211 -0.4737
0.5526 0.6316 -0.7105
>> B./A “两矩阵对应元素相除”
ans =
0.1429 1.0000 0.2000
1.0000 0.4000 0.3333
1.0000 3.0000 0.6000
（2）>> C=1:2:20
C =
1 3 5 7 9 11 13 15 17 19 C（i）表示奇数 2*i-1
(3)查找已创建变量的信息，删除无用的变量
clear “清除”
t =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);
(4)polar(t,(cos(2*t)).^0.5)
实验二 MATLAB 控制系统工具箱使用
1 实验目的
（1）了解Matlab 控制系统工具箱中常用函数.
（2）掌握模型转换方法以及奈氏图、伯德图的绘制。
（3）掌握系统脉冲、阶跃等响应的仿真方法
2 实验器材
计算机Win7、Matlab7.0软件
3 实验内容
（1）利用residue 函数求取传递函数的部分分式展开式
（2）利用tf 函数构建系统传递函数
（3）利用impulse 、step 、lsim 等函数仿真系统的单位脉冲、阶跃、速度响应。
（4）画出单位负反馈的开环传递函数的在频率0.01-100rad/s 间的伯德图
（5）画出系统的乃奎斯特图
4 实验步骤及内容
(1) 写出传递函数432
43211395226()10355024
s s s s G s s s s s ++++=++++的部分展开式
a. 列些分子num 、分母den 的系数
b. 利用residue(num,den) 求出零极点和余项
实验运行结果：
r =[ 1.0000 2.5000 -3.0000 0.5000 ] p =[-4.000 -3.0000 -2.000 -1.0000]
k=1
（2）对于系统3237(23)
()(31)(2)(538)s G s s s s s s +=++++
A.首先利用conv 函数展开分母，获得各阶系数，此时num=[14 21];
den=conv(conv(conv([1 0 0], [3 1]),conv([1 2],[1 2])), [5 0 3 8])
b.再利用tf 函数构建传递函数
运行结果: G(S)=100 s 240 s^ 223 s^ 84 s^ 15 s^50
s 20 s^ 223452+++++++
（3）a. 用tf(num,den)函数构建高阶系统传递函数2654322050()1584223240100
s s G s s s s s s ++=+++++ b.构建时间向量t
c 用函数impulse 及step 仿真系统的单位脉冲响应、单位阶跃响应
d. 用lism 函数仿真系统速度响应
e. 程序调试成功后实验结果：
A impulse 单位脉冲信号实验结果
B step 单位阶跃信号仿真结果
C lism 函数仿真系统速度响应
（4）绘制单位负反馈的开环传递函数10(1)()(7)
s G s s s +=+的伯德图。
a. 写出num ，den
b. 用logspace 函数设定频率范围
c 利用bode 函数输入幅值增益、相位
d. 把增益值转换为分贝值
e. 添加网格线，并生成二行一列的图表
f. 用semilogx 在两个图标中分别生成幅值频率和相频特性曲线。
g. 用margin 函数获取开环的幅值裕量和相位裕量。
h. 程序调试后写入实验报告
margin 函数获取开环的幅值裕量和相位裕量。
Gm = Inf Pm =126.0768 wcg =NaN
wcp = 7.2725
实验三 MATLAB 数字控制器设计初步
1 实验目的
（1）掌握离散信号z 变换和z 逆变换的matlab 方法
（2）掌握离散系统的系统函数零极点分布与系统频率特性的matlab 实现方法
（3）掌握matlab 计算差分方程的方法。
2 实验器材
计算机Win7、Matlab7.0软件
3 实验内容
（1）找到matlab 的符号运算工具箱提供的z 变换函数ztrans 与反变换函数iztrans ，通过help 理解其用法，并掌握符号表达式的构建和应用。
F=ztrans(f) %求符号函数f 的z 变换，返回函数的自变量为z ；
F=ztrans(f,w) %求符号函数f 的z 变换，返回函数的自变量为w ；
F=ztrans(f,k,w) %对自变量为k 的符号函数f 求z 变换，返回函数的自变量为w 。
f=iztrans(F) %对自变量为z 的符号函数F 求逆z 变换，返回函数的自变量为n ； f=iztrans(F,k) %对自变量为z 的符号函数F 求逆z 变换，返回函数的自变量为k ； f=iztrans(F,w,k) %对自变量为w 的符号函数F 求逆z 变换，返回函数的自变量为k 。
4 实验步骤及实验内容
（1）.已知序列k k f -=2)(，求其z 变换。
解：在命令窗口中输入如下命令，即可完成f(k)的z 变换
>> syms k
>> f=sym('2^(-k)'); %定义序列k k f -=2)(
>> F=ztrans(f) %求z 变换
F =
z/(z - 1/2)
（2）已知一离散系统的系统函数2
3)(2++=
z z z z H ，求其冲激响应h (k )。 解：运行如下M 文件，
syms k z
H=sym('z/(z^2+3*z+2)');
h=iztrans(H,k) %求逆z 变换
h =
(-1)^k - (-2)^k
(3) MA TLAB 的zplane 函数用于系统函数的零极点图的绘制,调用方式为：
zplane(b,a) 其中， b 、a 分别为系统函数分子、分母多项式的系数向量
例.已知一离散系统的系统函数2
31.07.0)(22+++-=z z z z z H ，试绘制其零极点图。
解：在MATLAB 的命令窗口中输入如下命令，即可得到其零极点图（图7-1）。 >> a=[1 3 2];
>> b=[1 -0.7 0.1];
>> zplane(b,a) %绘制其零极点图
(4) 若描述离散系统的差分方程为)()2(2
1)1(21)(n x n y n y n y =---+，已知激励)(2)(n u n x n =，初始状态y(-1)=1,y(-2)=0，求系统的零输入响应，零状态响应。
【解】将差分方程两边取z 变化的
)()]2()1()([2
1)]1()([21)(121z X y y z z Y z y z Y z z Y =-+-+--++---，进整理，得 )(2
1211121211221)1)(1(21)(21211z X z z z z y y z Y ------++-+-+--=）（，前一部分为零输入响应，后一部分为零状态响应。
将y(-1)=1,y(-2)=0代入第一部分，则零输入响应的表达式为
z
z z z z z z z z z z z z z z z z Y zi 161132)21)(1()1(212121)1(2121211)1(21)(2211-++-=-+--=-+--=-+-=---，因此经过z 的逆变化得，)(])2
1(61)1(32[)(n u n y n n zi +--=
同理，则零状态激励响应响应的表达式为
2
982
1911922212111)(212111)(2121-+--++=-?-+=-+=
----z z z z z z z z z z z X z z z Y zs ， )(])2(9
8)21(91)1(92[)(n u n y n n n zs +--= 如果直接利用iztrans 函数进行编程求解如下
clc
syms z
y=0.5*(1/z-1)/(1+1/2*z-1/2*z^2)
y=iztrans(y)
结果：1/2*charfcn[1](n)-3/4*charfcn[0](n)+1/12*2^n+2/3*(-1)^n
（5）已知22
(1)(3)()s s s s F s +++=，通过部分分式展开法求F(z) 。
F=sym('(s+2)/(s*(s+1)^2*(s+3))');%传递函数F(s)进行符号定义
[numF,denF]=numden(F) ;%提取分子分母
pnumF=sym2poly(numF) ;%将分子转化为一般多项式
pdenF=sym2poly(denF) ;%将分母转化为一般多项式
[R,P,K]=residue(pnumF,pdenF)%部分分式展开
对应部分分式分解结果为：
21111()0.083
0.750.50.6731(1)F s s s s s
=--++++
实验四 控制系统数学模型的MATLAB 实现
（用MATLAB 展开及合并部分分式）
1 实验目的
（1）掌握matlab 展开部分分式的方法
2 实验器材
计算机Win7、Matlab7.0软件
3 实验步骤及实验内容
(1)用matlab 展开部分分式
01m-1m n 01n-1n 11(s)++(s)(s)++m m n B b s b s b s b G A a s a s a s a --++==++……
用num 和den 分别表示G(s)的分子和分母多项式，即
[]01
m num b b b =…
[]01m den a a a =…
MATLAB 提供函数reside 用于实现部分分式展开，其句法为
[r, p, k]=reside(num,den)其中，r, p 分别为展开后的留数及极点构成的列向量；k 为余项多项式行向量。若无重极点，MA TLAB 展开后的一般形式为：
(1)(1)(n)(s)+(s)(1)(2)(n)r r r G K s p s p s p =+++---…
若存在q 个重极点p （j ），则展开式将包括下列各式：
2(j)(j 1)(j q 1)+(j)[(j)][s (j)]q r r r s p s p p ++-++---…
4 实验步骤及实验内容
（1）求下式的部分分式展开
43243211395226()10355024s s s s G s s s s s ++++=++++
解：>> num=[1 11 39 52 26];
>> den=[1 10 35 50 24];
>> [r,p,k]=residue(num,den)
r =
1.0000
2.5000
-3.0000
0.5000
p =
-4.0000
-3.0000
-2.0000
-1.0000
k =
1 所以展开式为：1 2.530.5(s)14331G s s s s -=++++++++
函数residue 也可用于将部分分式合并，其句法为
[num,den]=residue(r,p,k)
>> r=[1 2 3 4]';p=[-1 -2 -3 -4]';k=0
k =
>> [num,den]=residue(r,p,k)
num =
10 70 150 96
den =
1 10 35 50 24
所以合并式为：
32432107015096(s)10355024s s s F s s s s +++=++++
实验报告总结：
通过本次实验，培养了一定的动手能力以及一定的交流能力和团队精神，具备了对本专业、本行业等方面的大系统的理解、掌控能力。初步掌握机械设计制造、计算机与电子科学、电气与自动控制、网络与信息检索、传感与检测等基础知识和技术，具有机械工程产品与系统的设计和应用能力。通过本开放实验，初步掌握MATLAB仿真应用软件的建模，仿真方法与过程。同时也通过本次实验，了解了自身的不足，还需要多多学习。
自动控制理论实验报告 实验一 典型环节的时域响应 院系： 班级： 学号： 姓名：
实验一 典型环节的时域响应 一、 实验目的 1．掌握典型环节模拟电路的构成方法，传递函数及输出时域函数的表达式。 2．熟悉各种典型环节的阶跃响应曲线。 3．了解各项参数变化对典型环节动态特性的影响。 二、 实验设备 PC 机一台，TD-ACC+教学实验系统一套。 三、 实验步骤 1、按图1-2比例环节的模拟电路图将线接好。检查无误后开启设备电源。 注：图中运算放大器的正相输入端已经对地接了100k 电阻。不需再接。 2、将信号源单元的“ST ”端插针与“S ”端插针用“短路块”接好。将信号形式开关设为“方波”档，分别调节调幅和调频电位器，使得“OUT ”端输出的方波幅值为1V ，周期为10s 左右。 3、将方波信号加至比例环节的输入端R(t)， 用示波器的“CH1”和“CH2”表笔分别监测模拟电路的输入R(t)端和输出C(t)端。记录实验波形及结果。 4、用同样的方法分别得出积分环节、比例积分环节、惯性环节对阶跃信号的实际响应曲线。 5、再将各环节实验数据改为如下： 比例环节：；，k R k R 20020010== 积分环节：；，u C k R 22000== 比例环节：；，，u C k R k R 220010010=== 惯性环节：。，u C k R R 220010=== 用同样的步骤方法重复一遍。 四、 实验原理、内容、记录曲线及分析 下面列出了各典型环节的结构框图、传递函数、阶跃响应、模拟电路、记录曲线及理论分析。 1．比例环节 (1) 结构框图： 图1－1 比例环节的结构框图 (2) 传递函数： K S R S C =) () ( K R(S) C(S)
控制工程基础实验指导书 自控原理实验室编印
（内部教材）
实验项目名称: （所属课 程: 院系: 专业班级: 姓名: 学号: 实验日期: 实验地点: 合作者: 指导教师: 本实验项目成绩: 教师签字: 日期: （以下为实验报告正文） 、实验目的 简述本实验要达到的目的。目的要明确，要注明属哪一类实验（验证型、设计型、综合型、创新型）。 二、实验仪器设备 列出本实验要用到的主要仪器、仪表、实验材料等。 三、实验内容 简述要本实验主要内容，包括实验的方案、依据的原理、采用的方法等。 四、实验步骤 简述实验操作的步骤以及操作中特别注意事项。 五、实验结果
给出实验过程中得到的原始实验数据或结果，并根据需要对原始实验数据或结果进行必要的分析、整理或计算，从而得出本实验最后的结论。 六、讨论 分析实验中出现误差、偏差、异常现象甚至实验失败的原因，实验中自己发现了什么问题，产生了哪些疑问或想法，有什么心得或建议等等。 七、参考文献 列举自己在本次准备实验、进行实验和撰写实验报告过程中用到的参考文献资 料。 格式如下 作者，书名（篇名），出版社（期刊名），出版日期（刊期），页码
实验一控制系统典型环节的模拟、实验目的 、掌握比例、积分、实际微分及惯性环节的模拟方法; 、通过实验熟悉各种典型环节的传递函数和动态特性; 、了解典型环节中参数的变化对输出动态特性的影响。 二、实验仪器 、控制理论电子模拟实验箱一台; 、超低频慢扫描数字存储示波器一台; 、数字万用表一只;
、各种长度联接导线。 三、实验原理 运放反馈连接 基于图中点为电位虚地，略去流入运放的电流，则由图 由上式可以求得下列模拟电路组成的典型环节的传递函数及其单位阶跃响应。 、比例环节 实验模拟电路见图所示 U i R i U o 接示波器 以运算放大器为核心元件，由其不同的输入网络和反馈网络组成的各种典型环节，如图所示。图中和为复数阻抗，它们都是构成。 Z2 Z1 Ui ,— U o 接示波器 得:
实验1模拟控制系统在阶跃响应下的特性实验 一、实验目的 根据等效仿真原理，利用线性集成运算放大器及分立元件构成电子模拟器, 以干电池作为输入信号，研究控制系统的阶跃时间响应。 二、实验内容 研究一阶与二阶系统结构参数的改变，对系统阶跃时间响应的影响。 三、实验结果及理论分析 1.一阶系统阶跃响应 a. 电容值1uF，阶跃响应波形： b. 电容值2.2uF，阶跃响应波形:
c. 电容值4.4uF，阶跃响应波形: 2?—阶系统阶跃响应数据表 U r= -2.87V R°=505k? R i=500k? R2=496k 其中
T = R2C U c C：)=「(R/R2)U r 误差原因分析： ①电阻值及电容值测量有误差； ②干电池电压测量有误差； ③在示波器上读数时产生误差； ④元器件引脚或者面包板老化，导致电阻变大； ⑤电池内阻的影响输入电阻大小。 ⑥在C=4.4uF的实验中，受硬件限制，读数误差较大3?二阶系统阶跃响应 a.阻尼比为0.1，阶跃响应波形: b.阻尼比为0.5，阶跃响应波形:
4.二阶系统阶跃响应数据表 E R w ( ?) 峰值时间 U o (t p ) 调整时间 稳态终值 超调（%） 震荡次数 C. d. 阻尼比为0.7，阶跃响应波形: 阻尼比为1.0，阶跃响应波形: CHI 反相 带宽限制 伏/格
四、回答问题 1.为什么要在二阶模拟系统中 设置开关K1和K2 ,而且必须 同时动作？ 答：K1的作用是用来产生阶跃信号，撤除输入信后，K2则是构成了C2的 放电回路。当K1 一旦闭合（有阶跃信号输入），为使C2不被短路所以K2必须断开，否则系统传递函数不是理论计算的二阶系统。而K1断开后，此时要让 C2尽快放电防止烧坏电路，所以K2要立即闭合。 2.为什么要在二阶模拟系统中设置 F3运算放大器？ 答：反相电压跟随器。保证在不影响输入和输出阻抗的情况下将输出电压传递到输入端，作为负反馈。 实验2模拟控制系统的校正实验 一、实验目的 了解校正在控制系统中的作用
实验二二阶系统的瞬态响应分析 一、实验目的 1、熟悉二阶模拟系统的组成。 2、研究二阶系统分别工作在ξ=1，0<ξ<1，和ξ> 1三种状态下的单 位阶跃响应。 3、分析增益K对二阶系统单位阶跃响应的超调量σP、峰值时间tp和调 整时间ts。 4、研究系统在不同K值时对斜坡输入的稳态跟踪误差。 5、学会使用Matlab软件来仿真二阶系统，并观察结果。 二、实验仪器 1、控制理论电子模拟实验箱一台； 2、超低频慢扫描数字存储示波器一台； 3、数字万用表一只； 4、各种长度联接导线。 三、实验原理 图2-1为二阶系统的原理方框图，图2-2为其模拟电路图，它是由惯性环节、积分环节和反号器组成，图中K=R2/R1，T1=R2C1，T2=R3C2。 图2-1 二阶系统原理框图
图2-1 二阶系统的模拟电路 由图2-2求得二阶系统的闭环传递函 12 22 122112 /() (1)()/O i K TT U S K U S TT S T S K S T S K TT ==++++ :而二阶系统标准传递函数为 (1)(2), 对比式和式得 n ωξ== 12 T 0.2 , T 0.5 , n S S ωξ====若令则。调节开环增益K 值，不仅能改变系统无阻尼自然振荡频率ωn 和ξ的值，可以得到过阻尼(ξ>1)、 临界阻尼（ξ=1）和欠阻尼（ξ<1）三种情况下的阶跃响应曲线。 （1）当K >0.625， 0 < ξ < 1，系统处在欠阻尼状态，它的单位阶跃响应表达式为： 图2-3 0 < ξ < 1时的阶跃响应曲线 （2）当K =0.625时，ξ=1，系统处在临界阻尼状态，它的单位阶跃响应表达式为： 如图2-4为二阶系统工作临界阻尼时的单位响应曲线。 (2) +2+=222n n n S S )S (G ωξω ω1 ()1sin( ) (3) 2-3n t o d d u t t tg ξωωωω--=+=式中图为二阶系统在欠阻尼状态下的单位阶跃响应曲线 e t n o n t t u ωω-+-=)1(1)(
一、单项选择题： 1. 某二阶系统阻尼比为0，则系统阶跃响应为 A. 发散振荡 B. 单调衰减 C. 衰减振荡 D. 等幅振荡 2． 一阶系统G(s)=1 +Ts K 的时间常数T 越小，则系统的输出响应达到稳态值的时间 A ．越长 B ．越短 C ．不变 D ．不定 3. 传递函数反映了系统的动态性能，它与下列哪项因素有关？ A.输入信号 B.初始条件 C.系统的结构参数 D.输入信号和初始条件 4．惯性环节的相频特性)(ωθ，当∞→ω时，其相位移)(∞θ为 A ．-270° B ．-180° C ．-90° D ．0° 5．设积分环节的传递函数为G(s)= s 1 ，则其频率特性幅值M(ω)= A. ωK B. 2K ω C. ω1 D. 21ω 6. 有一线性系统，其输入分别为u 1(t)和u 2(t)时，输出分别为y 1(t)和y 2(t)。当输入为a 1u 1(t)+a 2u 2(t)时(a 1,a 2为常数)，输出应为 A. a 1y 1(t)+y 2(t) B. a 1y 1(t)+a 2y 2(t) C. a 1y 1(t)-a 2y 2(t) D. y 1(t)+a 2y 2(t) 7．拉氏变换将时间函数变换成 A ．正弦函数 B ．单位阶跃函数 C ．单位脉冲函数 D ．复变函数 8．二阶系统当0<ζ<1时，如果减小ζ，则输出响应的最大超调量%σ将 A.增加 B.减小 C.不变 D.不定 9．线性定常系统的传递函数，是在零初始条件下 A ．系统输出信号与输入信号之比 B ．系统输入信号与输出信号之比 C ．系统输入信号的拉氏变换与输出信号的拉氏变换之比 D ．系统输出信号的拉氏变换与输入信号的拉氏变换之比 10．余弦函数cos t ω的拉氏变换是 A.ω+s 1 B.2 2s ω+ω C.22s s ω+ D. 2 2s 1ω + 11. 微分环节的频率特性相位移θ(ω)= A. 90° B. -90° C. 0° D. -180° 12. II 型系统开环对数幅频渐近特性的低频段斜率为 A. -40(dB/dec) B. -20(dB/dec) C. 0(dB/dec) D. +20(dB/dec) 13．令线性定常系统传递函数的分母多项式为零，则可得到系统的
MATLAB程序设计语言 实 验 报 告 专业及班级：电子信息工程 姓名：王伟 学号：1107050322 日期 2013年6月20日
实验一 MATLAB 的基本使用 【一】 实验目的 1.了解MATALB 程序设计语言的基本特点，熟悉MATLAB 软件的运行环境； 2.掌握变量、函数等有关概念，掌握M 文件的创建、保存、打开的方法，初步具备将一般数学问题转化为对应计算机模型处理的能力； 3.掌握二维图形绘制的方法，并能用这些方法实现计算结果的可视化。 【二】 MATLAB 的基础知识 通过本课程的学习，应基本掌握以下的基础知识： 一. MATLAB 简介 二. MATLAB 的启动和退出 三. MATLAB 使用界面简介 四. 帮助信息的获取 五. MATLAB 的数值计算功能 六. 程序流程控制 七. M 文件 八. 函数文件 九. MATLAB 的可视化 【三】上机练习 1. 仔细预习第二部分内容，关于MATLAB 的基础知识。 2. 熟悉MATLAB 环境，将第二部分所有的例子在计算机上练习一遍 3. 已知矩阵???? ??????=??????????=123456789,987654321B A 。求A*B ，A .* B ，比较二者结果是否相同。并利用MATLAB 的内部函数求矩阵A 的大小、元素和、长度以 及最大值。 程序代码： >> A=[1 2 3;4 5 6;7 8 9]; >> B=[9 8 7;6 5 4;3 2 1]; >> A*B ans =
30 24 18 84 69 54 138 114 90 >> A.*B ans = 9 16 21 24 25 24 21 16 9 两者结果不同 >> [m,n]=size(A) m = 3 n = 3 >> b=sum(A) b = 12 15 18 >> a=length(A) a = 3 >>max(A) ans =
实验内容 （一）直流电机双环调速系统实验，此时必须松开连轴节！不带动工作台！ 1． 测试电流环特性 ，由于外接霍尔传感器只有一套，有五套PWM 放大器有电流输出（接成跟随器方式，其电流采样输出为25芯D 型插座的17（模拟地），19脚，但模拟地是电流环的模拟地，不是实验箱运算放大器OP07的地！所以，只能用万用表量测。多数同学可用手堵转，给定微小的输入电压（小于±50mV ）加入到电流环输入端,再加大就必须松开手，观察电机转速能否控制？为什么？如果要测试电流环静态特性，必须用台钳夹住电机轴，保证电机堵转。所以此项实验由教师按图22进行，这里只给出以下数据： 图 22 电流环静态特性实验接线图 （1）霍尔传感器的校准 利用直流稳压电源和电流表校准霍尔传感器，该 传感器为LEM-25,当原边为1匝时，量程为25A ，而原边采用5匝时， 量程为5A ；现在按后者的接法实验，M R 约500Ω。 （2）然后利用它来测试PWM 功率放大器的静态传递系数。电流环的静态特性如表2所示。注意电机是堵转的！
1V;得到通频带400Hz. 2.根据给定参数，利用MATLAB设计速度环的校正装置参数,画出校正前后的Bode图调，到实验室自己接线，教师检查无误后，可以通电调试；首先，正确接线保证系统处于负反馈，如果正反馈会产生什么现象？如何通过开环特性判断测速反馈是负反馈？对此有正确定答案后方能够开始实验。 （1）在1 β和β=0.4～0.5时分别调试校正装置的参数，使其单位阶跃输入的 = 响应曲线超调量最小，峰值时间最短，并记录阶跃响应曲线的特征值； 能够用A/D卡把数据采集到计算机中更好！ （2）断开电源，记录最佳的校正装置参数； （3）测试速度环静态特性，为加快测试速度，可直接测试输入电压和测速机电压的关系；在转速低的情况下用手动阻止电机的转动，是否会影响转速？ 为什么？分析速度环的机械特性（转速与负载力矩的关系曲线称为机械特 性），从而说明系统的刚度。 （4）有条件的小组可测试速度环频率特性（只测量幅频特性）。 （二）电压－位置伺服系统实验 开始，也必须脱开电机与工作台的连轴节！直到位置环调试好后，再把连轴节连接好！ 1．断开使能，手动电机转动，检查电子电位计工作的正确性！ 2．让位置环开环，利用调速系统，观察电子电位计在大范围工作的正确性，可利用示波器或万用表测试电位计的输出。 3．位置环要使用实验箱的头2个运算放大器，所以必须注意注意位置反馈的极性；为保证位置反馈是负反馈，必须通过位置系统开环来判断，这时位置调节器只利用比例放大器，如果发现目前的接线是正反馈后，怎么接线？ 4．将位置环的位置反馈正确接到反馈输入端，利用给定指令电位计，移动它，使电机位置按要求转动。正确后，即可把连轴节连接好，连接连轴节时用专用内六角扳手。这时应该断电！ 5．按设计的校正装置连接好，再上电。测试具有比例放大器和近似比例积分调节器时的阶跃响应曲线，并记录之； 6．测试输入电压－位置的传递特性曲线； 7．用手轮加小力矩估计系统的（电弹簧）刚度。 三、实验报告要求 （一）速度环实验 1．对速度环建模，画出速度环方块图，传递函数图 2．画出校正前后的Bode图，设计校正装置及其参数； 3．写出实验原始数据，整理出静态曲线和动态数据； 4．从理论和实际的结合上，分析速度环的特点，并写出实验的收获和改进意见； （二）位置环实验 1．对位置环建模，画出位置环方块图，传递函数图；
中北大学 机械工程与自动化学院 实验指导书 课程名称：《机械工程控制基础》 课程代号：02020102 适用专业：机械设计制造及其自动化 实验时数：4学时 实验室：数字化实验室 实验内容：1.系统时间响应分析 2.系统频率特性分析 机械工程系 2010.12
实验一 系统时间响应分析 实验课时数：2学时 实验性质：设计性实验 实验室名称：数字化实验室 一、实验项目设计内容及要求 1.试验目的 本实验的内容牵涉到教材的第3、4、5章的内容。本实验的主要目的是通过试验，能够使学生进一步理解和掌握系统时间响应分析的相关知识，同时也了解频率响应的特点及系统稳定性的充要条件。 2.试验内容 完成一阶、二阶和三阶系统在单位脉冲和单位阶跃输入信号以及正弦信号作用下的响应，求取二阶系统的性能指标，记录试验结果并对此进行分析。 3.试验要求 学习教材《机械工程控制基础（第5版）》第2、3章有关MA TLAB 的相关内容，要求学生用MA TLAB 软件的相应功能，编程实现一阶、二阶和三阶系统在几种典型输入信号（包括单位脉冲信号、单位阶跃信号、单位斜坡信号和正弦信号）作用下的响应，记录结果并进行分析处理：对一阶和二阶系统，要求用试验结果来分析系统特征参数对系统时间响应的影响；对二阶系统和三阶系统的相同输入信号对应的响应进行比较，得出结论。 4.试验条件 利用机械工程与自动化学院数字化试验室的计算机，根据MA TLAB 软件的功能进行简单的编程来进行试验。 二、具体要求及实验过程 1.系统的传递函数及其MA TLAB 表达 (1)一阶系统 传递函数为：1 )(+= Ts K s G 传递函数的MA TLAB 表达： num=[k]；den=[T,1]；G(s)=tf(num,den) (2)二阶系统 传递函数为：2 2 2 2)(n n n w s w s w s G ++= ξ 传递函数的MA TLAB 表达： num=[2n w ]；den=[1，ξ2wn ，wn^2]；G(s)=tf(num,den) （3）任意的高阶系统 传递函数为：n n n n m m m m a s a s a s a b s b s b s b s G ++++++++= ----11 101110)( 传递函数的MA TLAB 表达： num=[m m b b b b ,,,110- ]；den=[n n a a a a ,,,110- ]；G(s)=tf(num,den)
实验二 MATLAB矩阵及其运算 学号：3121003104 姓名：刘艳琳专业：电子信息工程1班日期：2014.9.20 一实验目的 1、掌握Matlab数据对象的特点以及数据的运算规则。 2、掌握Matlab中建立矩阵的方法以及矩阵处理的方法。 3、掌握Matlab分析的方法。 二实验环境 PC_Windows 7旗舰版、MATLAB 7.10 三实验内容 4、1. （1）新建一个.m文件，验证书本第15页例2-1； （2）用命令方式查看和保存代码中的所有变量；
（3）用命令方式删除所有变量； （4）用命令方式载入变量z。 2. 将x=[4/3 1.2345e-6]在以下格式符下输出：短格式、短格式e方式、长格式、长格式e方式、银行格式、十六进制格式、+格式。 短格式 短格式e 长格式
长格式e方式 银行格式 十六进制格式 3.计算下列表达式的值 （1）w=sqrt(2)*(1+0.34245*10^(-6)) （2）x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a) a=3.5;b=5;c=-9.8; （3）y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a) a=3.32;b=-7.9; （4）z=0.5*exp(2*t)*log(t+sqrt(1+t*t)) t=[2,1-3i;5,-0.65];
4. 已知A=[1 2 3 4 5 ;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20]，对其进行如下操作：（1）输出A在[ 7, 10]范围内的全部元素； （2）取出A的第2，4行和第1，3，5列； （3）对矩阵A变换成向量B，B=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20]； （4）删除A的第2，3，4行元素； (1) (2)
《控制工程基础》实验报告 姓名欧宇涵 914000720206 周竹青 914000720215 学院教育实验学院 指导老师蔡晨晓 南京理工大学自动化学院 2017年1月
实验1：典型环节的模拟研究 一、实验目的与要求： 1、学习构建典型环节的模拟电路； 2、研究阻、容参数对典型环节阶跃响应的影响； 3、学习典型环节阶跃响应的测量方法，并计算其典型环节的传递函数。 二、实验内容： 完成比例环节、积分环节、比例积分环节、惯性环节的电路模拟实验，并研究参数变化对其阶跃响应特性的影响。 三、实验步骤与方法 （1）比例环节 图1-1 比例环节模拟电路图 比例环节的传递函数为：K s U s U i O =)()(，其中1 2R R K =，参数取R 2=200K ，R 1=100K 。 步骤： 1、连接好实验台，按上图接好线。 2、调节阶跃信号幅值(用万用表测)，此处以1V 为例。调节完成后恢复初始。 3、Ui 接阶跃信号、Uo 接IN 采集信号。 4、打开上端软件，设置采集速率为“1800uS”，取消“自动采集”选项。 5、点击上端软件“开始”按键，随后向上拨动阶跃信号开关，采集数据如下图。 图1-2 比例环节阶跃响应
（2）积分环节 图1-3 积分环节模拟电路图 积分环节的传递函数为： S T V V I I O 1 -=，其中T I =RC ，参数取R=100K ，C=0.1μf 。 步骤：同比例环节，采集数据如下图。 图1-4 积分环节阶跃响应 （3）微分环节 图1-5 微分环节模拟电路图 200K R V I Vo C 2C R 1 V I Vo 200K
中北大学机械与动力工程学院 实验报告 专业名称__________________________________ 实验课程名称______________________________ 实验项目数_______________总学时___________ 班级______________________________________ 学号______________________________________ 姓名______________________________________ 指导教师__________________________________ 协助教师__________________________________ 日期______________年________月______日____
实验二二阶系统阶跃响应 一、实验目的 1．研究二阶系统的特征参数如阻尼比ζ和无阻尼自然频率ω n 对系统动态性能 的影响；定量分析ζ和ω n 与最大超调量Mp、调节时间t S 之间的关系。 2．进一步学习实验系统的使用方法。 3．学会根据系统阶跃响应曲线确定传递函数。 二、实验仪器 1．EL-AT-II型自动控制系统实验箱一台 2．PC计算机一台 三、实验原理 1．模拟实验的基本原理： 控制系统模拟实验采用复合网络法来模拟各种典型环节，即利用运算放大器不同的输入网络和反馈网络模拟各种典型环节，然后按照给定系统的结构图将这些模拟环节连接起来，便得到了相应的模拟系统。再将输入信号加到模拟系统的输入端，并利用计算机等测量仪器，测量系统的输出，便可得到系统的动态响应曲线及性能指标。若改变系统的参数，还可进一步分析研究参数对系统性能的影响。 2.时域性能指标的测量方法：超调量% σ： 1)启动计算机，在桌面双击图标 [自动控制实验系统] 运行软件。 2)测试计算机与实验箱的通信是否正常,通信正常继续。如通信不正常查 找原因使通信正常后才可以继续进行实验。 3)连接被测量典型环节的模拟电路。电路的输入U1接A/D、D/A卡的DA1 输出，电路的输出U2接A/D、D/A卡的AD1输入。检查无误后接通电源。 4)在实验课题下拉菜单中选择实验二[二阶系统阶跃响应] 。 5)鼠标双击实验课题弹出实验课题参数窗口。在参数设置窗口中设置相应 的实验参数后鼠标单击确认等待屏幕的显示区显示实验结果。 6)利用软件上的游标测量响应曲线上的最大值和稳态值，带入下式算出超 调量： Y MAX - Y ∞ % σ=——————×100% Y ∞ t P 与t s :利用软件的游标测量水平方向上从零到达最大值与从零到达95%稳 态值所需的时间值，便可得到t P 与t s 。 四、实验内容 典型二阶系统的闭环传递函数为 ω2 n ?（S）= （1） s2＋2ζω n s＋ω2 n
本科实验报告 课程名称：Matlab电子信息应用实验项目：矩阵和数组的操作 实验地点：电机馆跨越机房 专业班级：学号： 学生姓名： 指导教师： 2014年3月26 日
一、实验目的 1．掌握矩阵和数组的一般操作，包括创建、保存、修改和调用等。 2．学习矩阵和数组的加减运算与乘法。 3．掌握对数组中元素的寻访与赋值，会对数组进行一般的操作。 二、预备知识 1．常用的产生特殊矩阵的函数 ?eye(m,n) 单位阵 ?rand（m,n) 随机矩阵 ?randn(m,n) 正态分布的随机矩阵 ?zeros(m,n) 零矩阵 ?ones(m,n) 全部元素都为1的矩阵 ?compan(A) 矩阵A的伴随矩阵 ?bankel(m,n) n维Hankel矩阵 ?invhilb(n) n维逆Hilbert矩阵 ?magic(n) n维Magic矩阵 ?toeplitz(m,n) Toeplitz矩阵 ?wilkinson(n) n维Wilkinson特征值测试矩阵 ?handamard(n) n维Handamard矩阵 ?hilb(n) n维Hilbert矩阵 ?kron(A,B) Kronecker张量积 ?pascal(n) n维Pascal矩阵 ?vander(A) 由矩阵A产生Vandermonde矩阵 2．通过矩阵的结构变换，获得新矩阵 表2 矩阵结构变化产生新矩阵 L=tril(A) L主对角线及以下元素取矩阵A 的元素，其余为0 L=tril(A,k) L及第k条对角线及以下元素取矩阵A的元素，其余为 U=triu(A) U主对角线及以上的元素取矩阵A的元素，其余为0 U=triu(A,k) U第k条对角线及以上的元素取矩阵A的元素，其余为
页眉 实验1 模拟控制系统在阶跃响应下的特性实验一、实验目的 根据等效仿真原理，利用线性集成运算放大器及分立元件构成电子模拟器，以干电池作为输入信号，研究控制系统的阶跃时间响应。 二、实验内容 研究一阶与二阶系统结构参数的改变，对系统阶跃时间响应的影响。 三、实验结果及理论分析 1.一阶系统阶跃响应 a.电容值1uF，阶跃响应波形： b.电容值2.2uF，阶跃响应波形： 页脚 页眉
，阶跃响应波形：电容值c.4.4uF 阶系统阶跃响应数据表2.一稳态终值U（∞）（V）时间常数T(s) 电容值c（uF）理论值实际值实际值理论值0.50 2.87 1.0 0.51 2.90 1.07 2.90 2.2 2.87 1.02 2.06 2.90 2.87 4.4 2.24 元器件实测参数=505kU= -2.87V R? R=496k? =500kR?2o1r其中 T?RC2U(?)??(R/R)U rc21页脚 页眉 误差原因分析： ①电阻值及电容值测量有误差；
②干电池电压测量有误差； ③在示波器上读数时产生误差； ④元器件引脚或者面包板老化，导致电阻变大； ⑤电池内阻的影响输入电阻大小。 ⑥在C=4.4uF的实验中，受硬件限制，读数误差较大。 3.二阶系统阶跃响应 a.阻尼比为0.1，阶跃响应波形： b.阻尼比为0.5，阶跃响应波形： 页脚 页眉 ，阶跃响应波形：0.7c.阻尼比为
，阶跃响应波形：阻尼比为1.0d. 阶系统阶跃响应数据表4.二ξR（?）峰值时间U(t) 调整时间稳态终值超调（%）震荡次数pow M（）t）t（s V（）（s UV）N psps6 62.7 2.8 0.3 0.1 2.95 454k 4.8 1 0.5 0.5 3.3 52.9k 2.95 11.9 0.4 1 0.7 0.3 0.4 24.6k 3.0 2.7 2.92 1.0 1.0 2.98 1.0 2.97k 2.98 页脚 页眉 四、回答问题
《控制工程基础》实验指导书常熟理工学院机械工程学院 2009.9
目录 1.MATLAB时域分析实验 (2) 2.MATLAB频域分析实验 (4) 3.Matlab校正环节仿真实验 (8) 4.附录：Matlab基础知识 (14)
实验1 MATLAB 时域分析实验 一、实验目的 1． 利用MATLAB 进行时域分析和仿真。 要求：(1)计算连续系统的时域响应（单位脉冲输入，单位阶跃输入，任意输入）。 2．掌握Matlab 系统分析函数impulse 、step 、lsim 、roots 、pzmap 的应用。 二、实验内容 1．已知某高阶系统的传递函数为 ()265432 220501584223309240100 s s G s s s s s s s ++=++++++，试求该系统的单位脉冲响应、单位阶跃响应、单位速度响应和单位加速度响应。 MATLAB 计算程序 num=[2 20 50]; den=[1 15 84 223 309 240 100]; t= (0: 0.1: 20); figure (1); impulse (num,den,t); %Impulse Response figure (2); step(num,den,t);%Step Response figure (3); u1=(t); %Ramp.Input hold on; plot(t,u1); lsim(num,den,u1,t); %Ramp. Response gtext(‘t’); figure (4); u2=(t.*t/2);%Acce.Input u2=(0.5*(t.*t)) hold on; plot(t,u2); lsim(num,den,u2,t);%Acce. Response
实验二 MATLAB 程序设计 一、 实验目的 1.掌握利用if 语句实现选择结构的方法。 2.掌握利用switch 语句实现多分支选择结构的方法。 3.掌握利用for 语句实现循环结构的方法。 4.掌握利用while 语句实现循环结构的方法。 5.掌握MATLAB 函数的编写及调试方法。 二、 实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 M 文件的编写： 启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器（Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正 三、 实验内容 1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因 c b a 、、的不同取值而定） ，这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。并输入几组典型值加以检验。 （提示：提示输入使用input 函数） 2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。 要求：（1）用switch 语句实现。 （2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 （提示：注意单元矩阵的用法） 3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 2?1 3?10?5?16?8?4?2?1 6?3?10?5?16?8?4?2?1 运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。 请为关键的Matlab 语句填写上相关注释，说明其含义或功能。 4. 的值，调用该函数后，
南理工控制工程基础实验报告 成绩：《控制工程基础》课程实验报告班级：学号：姓名：南京理工大学2015年12月《控制工程基础》课程仿真实验一、已知某单位负反馈系统的开环传递函数如下G(s)?10 s2?5s?25借助MATLAB和Simulink完成以下要求：(1) 把G(s)转换成零极点形式的传递函数，判断开环系统稳定性。>> num1=[10]; >> den1=[1 5 25]; >> sys1=tf(num1,den1) 零极点形式的传递函数：于极点都在左半平面，所以开环系统稳定。(2) 计算闭环特征根并判别系统的稳定性，并求出闭环系统在0～10秒内的脉冲响应和单位阶跃响应，分别绘出响应曲线。>> num=[10];den=[1,5,35]; >>
sys=tf(num,den); >> t=[0::10]; >> [y,t]=step(sys,t); >> plot(t,y),grid >> xlabel(‘time(s)’) >> ylabel(‘output’) >> hold on; >> [y1,x1,t]=impulse(num,den,t); >> plot(t,y1,’:’),grid (3) 当系统输入r(t)?sin5t时，运用Simulink搭建系统并仿真，用示波器观察系统的输出，绘出响应曲线。曲线：二、某单位负反馈系统的开环传递函数为：6s3?26s2?6s?20G(s)?4频率范围??[,100] s?3s3?4s2?2s?2 绘制频率响应曲线，包括Bode图和幅相曲线。>> num=[6 26 6 20]; >> den=[1 3 4 2 2]; >> sys=tf(num,den); >> bode(sys,{,100}) >> grid on >> clear; >> num=[6 26 6 20]; >> den=[1 3 4 2 2]; >> sys=tf(num,den); >> [z , p , k] = tf2zp(num, den); >> nyquist(sys) 根据Nyquist判据判定系统的稳定性。
机械控制工程基础实验 指导书 Company number【1089WT-1898YT-1W8CB-9UUT-92108】
河南机电高等专科学校《机械控制工程基础》 实验指导书 专业：机械制造与自动化、起重运输机械设计与制造等 机械制造与自动化教研室编 2012年12月
目录
实验任务和要求 一、自动控制理论实验的任务 自动控制理论实验是自动控制理论课程的一部分，它的任务是： 1、通过实验进一步了解和掌握自动控制理论的基本概 念、控制系统的分析方法和设计方法； 2、重点学习如何利用MATLAB工具解决实际工程问题和 计算机实践问题； 3、提高应用计算机的能力及水平。 二、实验设备 1、计算机 2、MATLAB软件 三、对参加实验学生的要求 1、阅读实验指导书，复习与实验有关的理论知识，明确每次实验的目的，了解内容和方法。 2、按实验指导书要求进行操作；在实验中注意观察，记录有关数据和图 像，并由指导教师复查后才能结束实验。 3、实验后关闭电脑，整理实验桌子，恢复到实验前的情况。 4、认真写实验报告，按规定格式做出图表、曲线、并分析实验结果。字迹 要清楚，画曲线要用坐标纸，结论要明确。 5、爱护实验设备，遵守实验室纪律。 实验模块一 MATLAB基础实验 ——MATLAB环境下控制系统数学模型的建立 一、预备知识 的简介
MATLAB为矩阵实验室（Matrix Laboratory）的简称，由美国MathWorks公司出品的商业数学软件。主要用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。 来源：20世纪70年代，美国新墨西哥大学计算机科学系主任Cleve Moler 为了减轻学生编程的负担，用FORTRAN编写了最早的MATLAB。1984年由 Little、Moler、Steve Bangert合作成立了的MathWorks公司正式把MATLAB推向市场。到20世纪90年代，MATLAB已成为国际控制界的标准计算软件。 地位：和Mathematica、Maple并称为三大数学软件，在数学类科技应用软件中，在数值计算方面首屈一指。 功能：矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等。 应用范围：工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 图1-1 MATLAB图形处理示例 的工作环境 启动MATLAB，显示的窗口如下图所示。 MATLAB的工作环境包括菜单栏、工具栏以及命令运行窗口区、工作变量区、历史指令区、当前目录窗口和M文件窗口。 (1)菜单栏用于完成基本的文件输入、编辑、显示、MATLAB工作环境交互性设置等操作。 (2)命令运行窗口“Command Window”是用户与MATLAB交互的主窗口。窗口中的符号“》”表示MATLAB已准备好，正等待用户输入命令。用户可以在“》”提示符后面输入命令，实现计算或绘图功能。 说明：用户只要单击窗口分离键，即可独立打开命令窗口，而选中命令窗口中Desktop菜单的“Dock Command Window”子菜单又可让命令窗口返回桌面（MATLAB桌面的其他窗口也具有同样的操作功能）；在命令窗口中，可使用方向
实验二M A T L A B程序设计含实验报告 集团标准化工作小组 #Q8QGGQT-GX8G08Q8-GNQGJ8-MHHGN#
实验二 MATLAB 程序设计 一、 实验目的 1.掌握利用if 语句实现选择结构的方法。 2.掌握利用switch 语句实现多分支选择结构的方法。 3.掌握利用for 语句实现循环结构的方法。 4.掌握利用while 语句实现循环结构的方法。 5.掌握MATLAB 函数的编写及调试方法。 二、 实验的设备及条件 计算机一台（带有以上的软件环境）。 M 文件的编写： 启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器 （Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正 三、 实验内容 1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因c b a 、、的不同取值而定），这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。并输入几组典型值加以检验。 （提示：提示输入使用input 函数） 2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。 要求：（1）用switch 语句实现。
（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 （提示：注意单元矩阵的用法） 3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 21 3105168421 63105168421 运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。 请为关键的Matlab 语句填写上相关注释，说明其含义或功能。 4. y 5. (Root Mean Square)的计算（1（2）x=rand(1,200)，得到的x 为200个（0，1）之间均匀分布的随机数。 6.根据2 2222 1......3121116n ++++=π，求π的近似值。当n 分别取100、1000、10000时，结果是多少 思考题：
河南机电高等专科学校《机械控制工程基础》
目录 实验任务和要求............................................................................................................................................. 实验模块一MATLAB基础实验............................................................................................................
实验任务和要求 一、自动控制理论实验的任务 自动控制理论实验是自动控制理论课程的一部分，它的任务是： 1、通过实验进一步了解和掌握自动控制理论的基本概念、控制系统的分析方法和设 2、按实验指导书要求进行操作；在实验中注意观察，记录有关数据和图像，并由指 导教师复查后才能结束实验。 3、实验后关闭电脑，整理实验桌子，恢复到实验前的情况。
4、认真写实验报告，按规定格式做出图表、曲线、并分析实验结果。字迹要清楚， 画曲线要用坐标纸，结论要明确。 5、爱护实验设备，遵守实验室纪律。 实验模块一MATLAB基础实验 ——MATLAB环境下控制系统数学模型的 建立 一、预备知识 1.MATLAB的简介 MATLAB为矩阵实验室（Matrix Laboratory）的简称，由美国MathWorks公司出品的商业数学软件。主要用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。 来源：20世纪70年代，美国新墨西哥大学计算机科学系主任Cleve Moler为了减轻学生编程的负担，用FORTRAN编写了最早的MATLAB。1984年由Little、Moler、Steve Bangert合作成立了的MathWorks公司正式把MATLAB推向市场。到20世纪90年代，MATLAB已成为国际控制界的标准计算软件。 地位：和Mathematica、Maple并称为三大数学软件，在数学类科技应用软件中，在数值计算方面首屈一指。 功能：矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言