武汉工程大学实验报告
专业班号
组别 01 教师
姓名同组者（个人）
impulse( )绘制
语言程序：
>> num=[0 0 0 1 3 7]； >> den=[1 4 6 4 1 0]；
实验结果分析：
实验结果分析：
ω依次取值1,2,4,6时，系统超调量不变，延从上图可以看出，保持ζ=0.25不变,
n
要求：正文用小四宋体，1.5倍行距，图表题用五号宋体，图题位于图下方，表题位于表上方。
第三章 线性系统的时域分析法 3.1 引言 分析控制系统的第一步是建立模型，数学模型一旦建立，第二步 分析控制性能，分析有多种方法，主要有时域分析法，频域分析法，根轨迹法等。每种方法，各有千秋。均有他们的适用范围和对象。本章先讨论时域法。 实际上，控制系统的输入信号常常是不知的，而是随机的。很难用解析的方法表示。只有在一些特殊的情况下是预先知道的，可以用解析的方法或者曲线表示。例如，切削机床的自动控制的例子。 在分析和设计控制系统时，对各种控制系统性能得有评判、比较的依据。这个依据也许可以通过对这些系统加上各种输入信号比较它们对特定的输入信号的响应来建立。 许多设计准则就建立在这些信号的基础上，或者建立在系统对初始条件变化（无任何试验信号）的基础上，因为系统对典型试验信号的响应特性，与系统对实际输入信号的响应特性之间，存在着一定的关系；所以采用试验信号来评价系统性能是合理的。 3.1.1 典型试验信号 经常采用的试验输入信号： ① 实际系统的输入信号不可知性； ② 典型试验信号的响应与系统的实际响应，存在某种关系； ③ 电压试验信号是时间的简单函数，便于分析。 突然受到恒定输入作用或突然的扰动。如果控制系统的输入量是随时间逐步变化的函数，则斜坡时间函数是比较合适的。 （单位）阶跃函数（Step function ） 0,)(1≥t t 室温调节系统和水位调节系统 （单位）斜坡函数（Ramp function ） 速度 0,≥t t ∝ （单位）加速度函数（Acceleration function ）抛物线 0,2 12 ≥t t （单位）脉冲函数（Impulse function ） 0,)(=t t δ 正弦函数（Simusoidal function ）Asinut ，当输入作用具有周期性变化时。 通常运用阶跃函数作为典型输入作用信号，这样可在一个统一的基础上对各种控制系统的特性进行比较和研究。本章讨论系统非周期信号（Step 、Ramp 、对正弦试验信号相应，将在第五章频域分析法，第六章校正方法中讨论）作用下系统的响应。 3.1.2 动态过程和稳态过程
武汉工程大学实验报告 专业班号 组别01 教师 姓名同组者（个人）
2222)(n n n s s s G ωζωω++= （1）分别绘出)/(2s rad n =ω，ζ分别取0,,,和时的单位阶跃响应曲线，分析参数ζ对系统的影响，并计算ζ=时的时域性能指标ss s p r p e t t t ,,,,σ。 （2）绘制出当ζ=, n ω分别取1,2,4,6时单位阶跃响应曲线，分析参数n ω对系统的影响。 （3）系统的特征方程式为010532234=++++s s s s ，试用二种判稳方式判别该系统的稳定性。 （4）单位负反馈系统的开环模型为 )256)(4)(2()(2++++= s s s s K s G 试分别用劳斯稳定判据和赫尔维茨稳定判据判断系统的稳定性，并求出使得闭环系统稳定的K 值范围。 三、 实验结果及分析 1.可以用两种方法绘制系统的阶跃响应曲线。 （1）用函数step( )绘制 MATLAB 语言程序： >> num=[ 0 0 1 3 7]; >> den=[1 4 6 4 1 ]; >>step(num,den); >> grid; >>xlabel('t/s');ylabel('c(t)');title('step response');
MATLAB运算结果: (2)用函数impulse( )绘制 MATLAB语言程序： >> num=[0 0 0 1 3 7]； >> den=[1 4 6 4 1 0]； >> impulse(num,den)； >> grid； >> xlabel('t/s');ylabel('c(t)');title('step response')；MATLAB运算结果:
《机电系统控制基础》大作业一 基于MATLAB的机电控制系统响应分析 哈尔滨工业大学 2013年11月4日
1 作业题目 1. 用MATLAB 绘制系统2 ()25()() 425 C s s R s s s Φ== ++的单位阶跃响应曲线、单位斜坡响应曲线。 2. 用MATLAB 求系统2 ()25 ()()425 C s s R s s s Φ==++的单位阶跃响应性能指标：上升时间、峰值时间、调节时间和超调量。 3. 数控直线运动工作平台位置控制示意图如下： X i 伺服电机原理图如下： L R (1)假定电动机转子轴上的转动惯量为J 1，减速器输出轴上的转动惯量为J 2，减速器减速比为i ，滚珠丝杠的螺距为P ，试计算折算到电机主轴上的总的转动惯量J ； (2)假定工作台质量m ，给定环节的传递函数为K a ，放大环节的传递函数为K b ，包括检测装置在内的反馈环节传递函数为K c ，电动机的反电势常数为K d ，电动机的电磁力矩常数为K m ，试建立该数控直线工作平台的数学模型，画出其控制系统框图； (3)忽略电感L 时，令参数K a =K c =K d =R=J=1，K m =10，P/i =4π，利用MATLAB 分析kb 的取值对于系统的性能的影响。
2 题目1 单位脉冲响应曲线 单位阶跃响应曲线
源代码 t=[0:0.01:1.6]; %仿真时间区段和输入 nC=[25]; dR=[1,4,25]; fi=tf(nC,dR); %求系统模型 [y1,T]=impulse(fi,t); [y2,T]=step(fi,t); %系统响应 plot(T,y1); xlabel('t(sec)'),ylabel('x(t)'); grid on; plot(T,y2); xlabel('t(sec)'),ylabel('x(t)'); grid on; %生成图形 3 题目2 借助Matlab，可得： ans = 0.4330 0.6860 25.3826 1.0000 即
实验三使用MATLAB解决线性代数问题学院：数计学院班级：1003班姓名：黄晓丹学号：1051020144 实验目的： 学习MATLAB有关线性代数运算的指令，主要学习运用MATLAB解决矩阵除法，线性方程组的通解，矩阵相似 对角化问题，以及解决投入产出分析等应用问题。 实验内容： 矩阵转置：A=[1 2;3 4];B=[4 3;2 1]; >> A',B' ans = 1 3 2 4 ans = 4 3 3 1 矩阵加减：A-B ans= -3 -1 1 3 矩阵乘法：A*B,A.*B(数组乘法)||比较矩阵乘法与数组乘法的区别ans= 8 5 20 13 ans= 4 6 6 4 矩阵除法：A\B,B./A ans=
-6 -5 5 4 ans= 4 1.5 0.6667 0.25 特殊矩阵生成：zeros（m，n）||生成m行n列的矩阵 ones（m，n）||生成m行n列的元素全为一的矩阵 eye（n）||生成n阶单位矩阵 rand（m，n）||生成m行n列[0 ，1]上均匀分布随 机数矩阵 zeros(2,3) ans = 0 0 0 0 0 0 >> ones(3,3) ans = 1 1 1 1 1 1 1 1 1 >> eye(3)
ans = 1 0 0 0 1 0 0 0 1 >> rand(2,4) ans = Columns 1 through 3 0.9501 0.6068 0.8913 0.2311 0.4860 0.7621 Column 4 0.4565 0.0185 矩阵处理：trace（A）||返回矩阵的迹 diag（A）||返回矩阵对角线元素构成的向量 tril（A）||提取矩阵的下三角部分 triu（A）||提取矩阵的上三角部分 flipud（A）||矩阵上下翻转 fliplr（A）||矩阵左右翻转 reshape(A,m,n)||将矩阵的元素重排成m行n列矩阵A=[1 2 3;4 5 6;7 8 9]; >> t=trace(A),d=diag(A),u=triu(A)
第三章线性系统的时域分析方法 教学目的：通过本章学习，熟悉控制系统动态性能指标定义，掌握线性系统稳定的充要条件和劳斯判椐的应用，以及稳态误差计算方法，掌握一阶、 二阶系统的时域分析方法。 教学重点：掌握系统的动态性能指标，能熟练地应用劳斯判椐判断系统稳定性，二阶系统的动态响应特性分析。 教学难点：高阶系统的的动态响应特性分析。 本章知识结构图： 系统结构图闭环传递函数 一阶标准式 二阶标准式 特征方程稳定性、稳定域 代数判据 误差传递函数误差象函数终值定理稳态误差开环传递函数系统型别、开环增益 公式 静态误差系数 第九讲
3.1 系统时间响应的性能指标 一、基本概念 1、时域分析方法：根据系统的数学模型求出系统的时间响应来直接分析和评价系统的方法。 （1）响应函数分析方法：建立数学模型→确定输入信号→求出输出响应→ 根据输出响应→系统分析。 （2）系统测试分析方法：系统加入扰动信号→测试输出变化曲线→系统分析。 系统举例分析：举例：原料气加热炉闭环控制系统 2、分析系统的三大要点 （1）动态性能(快、稳) （2）稳态性能（准） （3）稳定性（稳） 二、动态性能及稳态性能 1、动态过程（过渡过程）：在 典型信号作用下，系统输出从初始状态到最终状态的响应过程。（衰减、发散、等幅振荡） 2、稳态过程：在典型信号作 用下，当t → ∞ 系统输出量表现的方式。表征输出量最终复现输入量的程度。（稳态误差描述） 3、动态稳态性能指标 图3-1温度控制系统原理图 （1）上升时间tr ：从稳态值的10%上升到稳态值的90%所需要的时间。 （2）峰值时间tp ：从零时刻到达第一个峰值h(tp)所用的时间。 （3）超调量δ%：最大峰值与稳态值的差与稳态值之比的百分数。（稳） （3-1） %100)(()(%?∞∞-= h h t h p ） δ
武汉工程大学 实验报告 专业 班号 组别 指导教师 姓名 学号 实验名称 线性系统时域响应分析 一、实验目的 1．熟练掌握step( )函数和impulse( )函数的使用方法，研究线性系统在单位阶跃、单位脉冲及单位斜坡函数作用下的响应。 2．通过响应曲线观测特征参量ζ和n ω对二阶系统性能的影响。 3．熟练掌握系统的稳定性的判断方法。 二、实验内容 1．观察函数step( )和impulse( )的调用格式，假设系统的传递函数模型为 1 4647 3)(2 342++++++=s s s s s s s G 可以用几种方法绘制出系统的阶跃响应曲线试分别绘制。 2．对典型二阶系统 2 22 2)(n n n s s s G ωζωω++= 1）分别绘出)/(2s rad n =ω，ζ分别取0,,,和时的单位阶跃响应曲线，分析参数ζ对系统的影响，并计算ζ=时的时域性能指标ss s p r p e t t t ,,,,σ。 2）绘制出当ζ=, n ω分别取1,2,4,6时单位阶跃响应曲线，分析参数n ω对系统的影响。 3．系统的特征方程式为010532234=++++s s s s ，试用两种判稳方式判别该系统的稳定性。 4．单位负反馈系统的开环模型为 ) 256)(4)(2()(2++++= s s s s K s G
试用劳斯稳定判据判断系统的稳定性，并求出使得闭环系统稳定的K 值范围。 三、实验结果及分析 1．观察函数step( )和impulse( )的调用格式，假设系统的传递函数模型为 14647 3)(2342++++++=s s s s s s s G 可以用几种方法绘制出系统的阶跃响应曲线试分别绘制。 方法一：用step( )函数绘制系统阶跃响应曲线。 程序如下： num=[0 0 1 3 7]; den=[1 4 6 4 1]; t=0::10; step(num,den) grid xlabel('t/s'),ylabel('c(t)') title('Unit-step Response of G(s)=s^2+3s+7/(s^4+4s^3+6s^2+4s+1)') Unit-step Response of G(s)=s 2+3s+7/(s 4+4s 3+6s 2+4s+1) t/s (sec) c (t ) 方法二：用impulse( )函数绘制系统阶跃响应曲线。 程序如下： num=[0 0 0 1 3 7 ]; den=[1 4 6 4 1 0]; t=0::10; impulse(num,den) grid xlabel('t/s'),ylabel('c(t)') title('Unit-impulse Response of G(s)/s=s^2+3s+7/(s^5+4s^4+6s^3+4s^2+s)')
信号时域与频域分析 实验报告 姓名：杨 班级：机械 学号： 213
实验数据中，电机转速为1200r/min，采样频率为1280Hz。Hz3为X位移振幅数据，Hz4为Y位移振幅数据，Hz5为速度振幅数据。 Matlab中信号特征对应函数编程 ma = max(Hz) %最大值 mi = min(Hz) %最小值 me = mean(Hz) %平均值 pk = ma-mi %峰-峰值 va = var(Hz); %方差 st = std(Hz); %标准差 ku = kurtosis(Hz); %峭度 rm = rms(Hz); %均方根 一、X轴位移测量分析 plot(Fs3,Hz3)时域图： ma =52.0261 mi =56.7010 me =1.8200 pk =108.7271 va =1.3870e+03 st =37.2431 ku =1.5462 rm =37.2693 频域图： fs=1280; x=Hz3; N=length(Hz3); df=fs/N; f=0:df:N*df-df; y=fft(x); y=abs(y)*2/N; figure(1); plot(f,y); xlabel('频率/Hz') ylabel('幅值') 频谱幅值取得最大值51.9847um，频率为20Hz，与电机转速对应频率一致，应为电机轴未动平衡所致；二倍频处有较大振幅，可能为轴承间隙过大所致。
二、Y轴位移测量分析 plot(Fs4,Hz4)时域图： ma =61.3987 mi =-74.6488 me =-1.1948 pk =136.0475 av =42.6109 va =2.2428e+03 st =47.3582 ku =1.5135 rm =47.3501 频域图： fs=1280; x=Hz4; N=length(Hz4); df=fs/N; f=0:df:N*df-df; y=fft(x); y=abs(y)*2/N; figure(1); plot(f,y); xlabel('频率/Hz') ylabel('幅值') 频谱幅值取得最大值66.6319um，频率为20Hz，与电机转速对应频率一致，应为电机轴未动平衡所致；二倍频处有较大振幅，可能为轴承间隙过大所致。
实验报告 实验名称 线性系统时域响应分析 一、 实验目的 1．熟练掌握step( )函数和impulse( )函数的使用方法，研究线性系统在单位阶跃、单位脉冲及单位斜坡函数作用下的响应。 2．通过响应曲线观测特征参量ζ和n ω对二阶系统性能的影响。 3．熟练掌握系统的稳定性的判断方法。 二、 实验内容 1．观察函数step( )和impulse( )的调用格式，假设系统的传递函数模型为 1 4647 3)(2 342++++++=s s s s s s s G 可以用几种方法绘制出系统的阶跃响应曲线？试分别绘制。 2．对典型二阶系统 2 22 2)(n n n s s s G ωζωω++= 1）分别绘出)/(2s rad n =ω，ζ分别取0,0.25,0.5,1.0和2.0时的单位阶跃响应曲线，分析参数ζ对系统的影响，并计算ζ=0.25时的时域性能指标 ss s p r p e t t t ,,,,σ。 2）绘制出当ζ=0.25, n ω分别取1,2,4,6时单位阶跃响应曲线，分析参数n ω对系统的影响。 3．系统的特征方程式为010532234=++++s s s s ，试用两种判稳方式判别该系统的稳定性。 4．单位负反馈系统的开环模型为 ) 256)(4)(2()(2 ++++= s s s s K s G 试用劳斯稳定判据判断系统的稳定性，并求出使得闭环系统稳定的K 值范围。
三、 实验结果及分析 1.观察函数step( )和impulse( )的调用格式，假设系统的传递函数模型为 1 4647 3)(2342++++++=s s s s s s s G 可以用几种方法绘制出系统的阶跃响应曲线？试分别绘制。 方法一： num=[1 3 7]; den=[1 4 6 4 1]; step(num,den) grid xlabel('t/s'),ylabel('c(t)') title('Unit-step Respinse of G(s)=(s^2+3s+7)/(s^4+4s^3+6s^2+4s+1)') 方法二： num=[1 3 7]; den=[1 4 6 4 1 0]; impulse(num,den) grid xlabel('t/s'),ylabel('c(t)') title('Unit-impulse Respinse of G(s)/s=(s^2+3s+7)/(s^5+4s^4+6s^3+4s^2+s)')
实验2 MATLAB数值计算、符号运算功能 一、实验目的 1、掌握建立矩阵、矩阵分析与处理的方法。 2、掌握线性方程组的求解方法。 3、掌握数据统计和分析方法、多项式的常用运算。 4、掌握求数值导数和数值积分、常微分方程数值求解、非线性代数方程数值求解的方法。 5、掌握定义符号对象的方法、符号表达式的运算法则及符号矩阵运算、符号函数极限及导数、符号函数定积分和不定积分的方法。 二、预习要求 （1）复习4、5、6章所讲内容； （2）熟悉MATLAB中的数值计算和符号运算的实现方法和主要函数。 三、实验内容 1、已知 29618 20512 885 A -?? ?? =?? ?? - ?? ，求A的特征值及特征向量，并分析其数学意义。 >> A=[-29,6,18;20,5,12;-8,8,5]; >> [V,D]=eig(A) V = 0.7130 0.2803 0.2733 -0.6084 -0.7867 0.8725 0.3487 0.5501 0.4050 D = -25.3169 0 0 0 -10.5182 0 0 0 16.8351 V为A的特征向量，D为A的特征值，3个特征值是-25.3169、10.5182和16.8351。 >> A*V ans = -18.0503 -2.9487 4.6007 15.4017 8.2743 14.6886 -8.8273 -5.7857 6.8190 >> V*D
ans = -18.0503 -2.9487 4.6007 15.4017 8.2743 14.6886 -8.8273 -5.7857 6.8190 经过计算，A*V=V*D 。 2、 不用rot90函数，实现方阵左旋90°或右旋90°的功能。例如，原矩阵为A ，A 左旋后得到B ，右旋后得到C 。 147102581136912A ????=??????，101112789456123B ??????=??????，321654987121110B ??????=?????? 提示：先将A 转置，再作上下翻转，则完成左旋90°；如将A 转置后作左右翻转，则完成右旋转90°，可用flipud 、fliplr 函数。 >> a=[1 4 7 10;2 5 8 11;3 6 9 12] a= 1 4 7 10 2 5 8 11 3 6 9 12 >> B=rot90(a) B = 10 11 12 7 8 9 4 5 6 1 2 3 >>C= rot90(s,3) C= 3 2 1 6 5 4 9 8 7 12 11 10
信号的频域 在电子学、控制系统及统计学中，频域是指在对函数或信号进行分析时，分析其和频率有关部份，而不是和时间有关的部份，和时域一词相对。函数或信号可以透过一对数学的运算子在时域及频域之间转换。例如傅里叶变换可以将一个时域信号转换成在不同频率下对应的振幅及相位，其频谱就是时域信号在频域下的表现，而反傅里叶变换可以将频谱再转换回时域的信号。 以信号为例，信号在时域下的图形可以显示信号如何随着时间变化，而信号在频域下的图形（一般称为频谱）可以显示信号分布在哪些频率及其比例。频域的表示法除了有各个频率下的大小外，也会有各个频率的相位，利用大小及相位的资讯可以将各频率的弦波给予不同的大小及相位，相加以后可以还原成原始的信号。在频域的分析中，常会用频谱分析仪来将实际的信号转换为频域下的频谱。 频域，尤其在射频和通信系统中运用较多，在高速数字应用中也会遇到频域。频域最重要的性质是：它不是真实的，而是一个数学构造。时域是惟一客观存在的域，而频域是一个遵循特定规则的数学范畴。 正弦波是频域中唯一存在的波形，这是频域中最重要的规则，即正弦波是对频域的描述，因为时域中的任何波形都可用正弦波合成。这是正弦波的一个非常重要的性质。然而，它并不是正弦波的独有特性，还有许多其他的波形也有这样的性质。正弦波有四个性质使它可以有效地描述其他任一波形： （1）时域中的任何波形都可以由正弦波的组合完全且惟一地描述。 （2）任何两个频率不同的正弦波都是正交的。如果将两个正弦波相乘并在整个时间轴上求积分，则积分值为零。这说明可以将不同的频率分量相互分离开。
（3）正弦波有精确的数学定义。 （4）正弦波及其微分值处处存在，没有上下边界。 使用正弦波作为频域中的函数形式有它特别的地方。若使用正弦波，则与互连线的电气效应相关的一些问题将变得更容易理解和解决。如果变换到频域并使用正弦波描述，有时会比仅仅在时域中能更快地得到答案。 而在实际中，首先建立包含电阻，电感和电容的电路，并输入任意波形。一般情况下，就会得到一个类似正弦波的波形。而且，用几个正弦波的组合就能很容易地描述这些波形。 许多物理元件的特性会随着输入讯号的频率而改变，例如电容在低频时阻抗变大，高频时阻抗变小，而电感恰好相反，高频时阻抗变大，低频时阻抗变小。一个线性非时变系统的特性也会随频率而变化，因此也有其频域下的特性，频率响应的图形即为其代表。频率响应可以视为是一个系统在输入信号振幅相同、频率不同时，其输出信号振幅的变化，可以看出系统在哪些频率的输出较大。有些系统的定义就是以频域为主，例如低通滤波器只允许低于一定频率的讯号通过。 不论是进行拉普拉斯转换、Z转换或是傅立叶变换，其产生的频谱都是一个频率的复变函数，表示一个信号（或是系统的响应）的振幅及其相位。不过在许多的应用中相位的资讯并不重要，若不考虑相位的资讯，都可以将频谱的资讯只以不同频率下的振幅（或是功率密度）来表示。 功率谱密度是一种常应用在许多非周期性也不满足平方可积性（square-integrable）讯号的频域表示法。只要一个讯号是符合广 义平稳随机过程的输出，就可以计算其对应的功率谱密度。 时域分析与频域分析是对模拟信号的两个观察面。时域分析是以时间轴为坐标表示动态信号的关系；频域分析是把信号变为以频率轴为坐标表示出来。一般来说，时域的表示较为形象与直观，频域分析
Matlab线性代数实验指导书 理学院线性代数课程组 二零零七年十月
目录 一、基础知识 (1) 1.1、常见数学函数 (1) 1.2、系统在线帮助 (1) 1.3、常量与变量 (2) 1.4、数组（矩阵）的点运算 (3) 1.5、矩阵的运算 (3) 二、编程 (4) 2.1、无条件循环 (4) 2.2、条件循环 (5) 2.3、分支结构 (5) 2.4、建立M文件 (6) 2.5、建立函数文件 (6) 三、矩阵及其运算 (7) 3.1、矩阵的创建 (7) 3.2、符号矩阵的运算 (11) 四、秩与线性相关性 (14) 4.1、矩阵和向量组的秩以及向量组的线性相关性 (14) 4.2、向量组的最大无关组 (14) 五、线性方程的组的求解 (16) 5.1、求线性方程组的唯一解或特解（第一类问题） (16) 5.2、求线性齐次方程组的通解 (18) 5.3、求非齐次线性方程组的通解 (19) 六、特征值与二次型 (22) 6.1、方阵的特征值特征向量 (22) 6.2、正交矩阵及二次型 (23)
一、基础知识 1.1常见数学函数 函数数学计算功能函数数学计算功能 abs(x) 实数的绝对值或复数的幅值floor(x) 对x朝-∞方向取整acos(x) 反余弦arcsinx gcd(m,n) 求正整数m和n的最大公约数acosh(x) 反双曲余弦arccoshx imag(x) 求复数x的虚部angle(x) 在四象限内求复数x的相角lcm(m,n)求正整数m和n的最小公倍 自然对数(以e为底数) asin(x) 反正弦arcsinx log(x) 常用对数(以 10 为底数) asinh(x) 反双曲正弦arcsinhx log10(x) atan(x) 反正切arctanx real(x) 求复数 x 的实部atan2(x,y) 在四象限内求反正切rem(m,n) 求正整数m和n的m/n之余数atanh(x) 反双曲正切arctanhx round(x) 对x四舍五入到最接近的整数 符号函数：求出 x 的符号ceil(x) 对x朝+∞方向取整 sign(x) conj(x) 求复数x的共轭复数 sin(x) 正弦sinx 反双曲正弦sinhx cos(x) 余弦cosx sinh(x) cosh(x) 双曲余弦coshx sqrt(x) 求实数x的平方根exp(x) 指数函数e x tan(x) 正切tanx fix(x) 对 x 朝原点方向取整 tanh(x) 双曲正切tanhx 如：输入 x=[-4.85 -2.3 -0.2 1.3 4.56 6.75]，则： ceil(x)= -4 -2 0 2 5 7 fix(x) = -4 -2 0 1 4 6 floor(x) =-5 -3 -1 1 4 6 round(x) = -5 -2 0 1 5 7 1.2 系统的在线帮助 1.2.1 help 命令： 1.当不知系统有何帮助内容时，可直接输入 help以寻求帮助: >> help(回车) 2.当想了解某一主题的内容时，如输入： >> help syntax (了解Matlab的语法规定) 3.当想了解某一具体的函数或命令的帮助信息时，如输入： >> help sqrt (了解函数sqrt的相关信息) 1.2.2 lookfor 命令 现需要完成某一具体操作，不知有何命令或函数可以完成，如输入： >> lookfor line (查找与直线、线性问题有关的函数) 1.3 常量与变量
肌电信号的时域和频域分析 摘要：肌电信号是产生肌肉力的电信号根源，它是肌肉中很多运动单元动作电位在时间和空间上的叠加，反映了神经，肌肉的功能状态，在基础医学研究、临床诊断和康复工程中有广泛的应用。 其种类重要有两种：一，临床肌电图检查多采用针电极插入肌肉检测肌 电图，其优点是干扰小，定位性好，易识别，但由于它是一种有创伤的检测 方法，其应用收到了一定的限制。二，表面肌电则是从人体皮肤表面通过电 极记录下来的神经肌肉活动时发放的生物电信号，属于无创伤性，操作简单，病人易接受，有着广泛的应用前景。 本次设计基于matlab用小波变换对肌电信号进行消噪处理，分别选用20N 的肌电信号数据和50N的肌电数据进行对比，最后在GUI界面上完成相应的功能处理。 关键字：肌电信号 Matlab 小波去噪 GUI 第一章绪论 肌电信号是产生肌肉力的电信号根源，它是肌肉中很多运动单元动作电位在时间和空间上的叠加，反映了神经，肌肉的功能状态，在基础医学研究、临床诊断和康复工程中有广泛的应用。 其种类重要有两种：一，临床肌电图检查多采用针电极插入肌肉检测肌电图，其优点是干扰小，定位性好，易识别，但由于它是一种有创伤的检测方法，其应用收到了一定的限制。二，表面肌电则是从人体皮肤表面通过电极记录下来的神经肌肉活动时发放的生物电信号，属于无创伤性，操作简单，病人易接受，有着广泛的应用前景。 肌电信号本身是一种较微弱的电信号。检测和记录表面肌电信号，需要考虑的主要问题是尽量消除噪声和干扰的影响, 提高信号的保真度[1]。
第二章肌电信号的时域分析 2.1 肌电信号时域图的显示及比较 肌电信号采用两个不同的数据进行比较，通过比较时域图及其特性来进行分析[2]。其图像如下所示： 如上图所示：肌电数据分别是同一个体在20N的力和50N的力所反映的图像。可以看出在不同作用力时，其图像的差别很大。 2.2 时域参数 2.2.1 均值 对于一个随机变量来说，均值是一个很重要的数值特征。粗略的说，就是来描述一个群体的平均水平。其严格的数学定义非常的简单，就是一个随机变量关于概率测度的积分。这样的积分在测度轮或者实分析里是没有什么直观的解释的。而在概率论里却成为了一个群体的主要指标。在此处，均值表示肌电信号的平均水平。 2.2.2 标准差 标准差（Standard Deviation），也称均方差（mean square error），是各数据偏离平均数的距离的平均数，它是离均差平方和平均后的方根，用σ表示。
第三章 线性系统的时域分析与校正习题及答案 3-1 已知系统脉冲响应t 25.1e 0125.0)t (k -=，试求系统闭环传递函数)s (Φ。 解 [])25.1s /(0125.0)t (k L )s (+==Φ 3-2 设某高阶系统可用下列一阶微分方程)t (r )t (r )t (c )t (c T +τ=+? ? 近似描述，其中，1)T (0<τ-<。试求系统的动态性能指标s r d t ,t ,t 。 解 设单位阶跃输入s s R 1)(= 当初始条件为0时有： 1 Ts 1 s )s (R )s (C ++τ= 1Ts T s 1s 11Ts 1s )s (C +τ--=?++τ= ∴ T /t e T T 1)t (h )t (c -τ--== T )0(h τ=，1)(h =∞，20T T )]0(h )(h [05.0τ -=-∞=? 1) 当 d t t = 时 2T T e T T 1)]0(h )(h [5.0)0(h )t (h t /t d τ += τ--=-∞+=- T /t d e 2 1 -= ; 693T .0t d = 2) 求r t （即)t (c 从1.0)(h ∞到9.0)(h ∞所需时间) 当T /t 2e T T 1)0(h )]0(h )(h [9.0)t (h -τ-- =+-∞=; 当T /t 1e T T 1)0(h )]0(h )(h [1.0)t (h -τ--=+-∞=; )T 1(.0T ln T t 2τ+τ-=， τ +τ -=)T 9(.0T ln T t 1 则 2T .29ln T t t t 12r ==-= 3) 求 s t T /t s s e T T 1)0(h )]0(h )(h [95.0)t (h -τ-- =+-∞= 3T 05.ln0T t s ==∴ 3-3 一阶系统结构如图所示。要求系统闭环增益2k =Φ，调节时间4.0t s ≤s ，试确定参数21k ,k 的值。 解 由结构图写出闭环系统传递函数 1k k s k 1k k s k s k k 1s k )s (212211211 +=+=+ =Φ
第二章：信号时域与频域分析实验指导书 一.实验目的 本实验结合《机械系统故障诊断》课程第二章“信号时域与频域分析”的课堂教学内容，通过实验进一步了解振动信号的获取过程与时域、频域分析方法，加深对所学的理论知识的掌握与理解。 二.教学基本要求 要求学生学习并掌握信号调理、采集与时域分析、频域分析方法。搭建由振动传感器、数据采集箱、计算机组成的信号调理与采集系统，测量故障模拟试验系统的振动信号，用Matlab软件编写信号的时域分析和频域分析程序，学会数字信号的获取与分析方法，掌握振动信号的测量系统搭建的基本方法。 三.实验内容 搭建用于信号调理、采集测试系统，测量振动信号，用Matlab软件编写信号的时域分析和频域分析软件，并对所测信号进行时域分析和频域分析，撰写实验报告。具体要求如下： 1.利用实验室现有设备搭建由振动传感器、信号调理箱、A/D板、计算机组成的振动数据采集系统； 2.采集转子试验台的电涡流信号、速度传感器信号、加速度传感器信号，对比不同传感器信号区别，总结不同传感器适用场合。 3.找出电涡流信号/速度传感器/加速度传感器信号的时域信号特征（波形、峰值、脉冲、峭度等）、频域信号特征。 4.对加速度信号/速度传感器信号做自相关、互相关分析。 5. 在Matlab软件中编写时域、频谱、自相关、互相关分析软件； 6. 利用自编软件分析所测数据并编写实验报告。 四. 使用的主要仪器 电涡流式位移、速度、加速度传感器、信号采集箱、计算机。 五.实验报告要求 1．实验报告内容包括计算分析的图、表或数值结果，以及对结果的简要分析、自编软件； 2．实验报告应独立完成；六.实验注意事项 1．开启电源前检查传感器安装、电源线、信号线连接是否正确； 2．实验完成后，关闭仪器的电源、清洁好实验台。
实验二 线性系统时域响应分析 一、实验目的 1．熟练掌握step( )函数和impulse( )函数的使用方法，研究线性系统在单位阶跃、单位脉冲及单位斜坡函数作用下的响应。 2．通过响应曲线观测特征参量ζ和n ω对二阶系统性能的影响。 二、基础知识及MATLAB 函数 （一）基础知识 时域分析法直接在时间域中对系统进行分析，可以提供系统时间响应的全部信息，具有直观、准确的特点。为了研究控制系统的时域特性，经常采用瞬态响应（如阶跃响应、脉冲响应和斜坡响应）。本次实验从分析系统的性能指标出发，给出了在MATLAB 环境下获取系统时域响应和分析系统的动态性能和稳态性能的方法。 用MATLAB 求系统的瞬态响应时，将传递函数的分子、分母多项式的系数分别以s 的降幂排列写为两个数组num 、den 。由于控制系统分子的阶次m 一般小于其分母的阶次n ，所以num 中的数组元素与分子多项式系数之间自右向左逐次对齐，不足部分用零补齐，缺项系数也用零补上。 1．用MATLAB 求控制系统的瞬态响应 1）阶跃响应 求系统阶跃响应的指令有： step(num,den) 时间向量t 的范围由软件自动设定，阶跃响应曲线 随即绘出 step(num,den,t) 时间向量t 的范围可以由人工给定（例如 t=0:0.1:10） [y ，x]=step(num,den) 返回变量y 为输出向量，x 为状态向量 在MATLAB 程序中，先定义num,den 数组，并调用上述指令，即可生成单位阶跃输入信号下的阶跃响应曲线图。 考虑下列系统： 25 425 )()(2++=s s s R s C 该系统可以表示为两个数组，每一个数组由相应的多项式系数组成，并且以s
第三章 线性系统的时域分析与校正 习题及答案 3-1 已知系统脉冲响应 t e t k 25.10125.0)(-= 试求系统闭环传递函数)(s Φ。 解 Φ()()./(.)s L k t s ==+00125125 3-2 设某高阶系统可用下列一阶微分方程 T c t c t r t r t ?? +=+()()()()τ 近似描述，其中，1)(0<-<τT 。试证系统的动态性能指标为 T T T t d ?? ? ?????? ??-+=τln 693.0 t T r =22. T T T t s ?? ??? ? -+=)ln( 3τ 解 设单位阶跃输入s s R 1)(= 当初始条件为0时有： 1 1 )()(++=Ts s s R s C τ 1 11 11)(+--= ? ++= ∴ Ts T s s Ts s s C ττ C t h t T T e t T ()()/==---1τ 1) 当 t t d = 时 h t T T e t t d ()./==---051τ 12=--T T e t T d τ/ ; T t T T d -??? ??-=-τln 2ln ????? ???? ??-+=∴ T T T t d τln 2ln
2) 求t r （即)(t c 从1.0到9.0所需时间) 当 T t e T T t h /219.0)(--- ==τ; t T T T 201=--[ln()ln .]τ 当 T t e T T t h /111.0)(---==τ; t T T T 109=--[ln()ln .]τ 则 t t t T T r =-==21 09 01 22ln ... 3) 求 t s T t s s e T T t h /195.0)(---==τ ]ln 3[]20ln [ln ]05.0ln [ln T T T T T T T T T t s τ ττ-+=+-=--=∴ 3-3 一阶系统结构图如图3-45所示。要求系统闭环增益2=ΦK ，调节时间4.0≤s t s ，试确定参数21,K K 的值。 解 由结构图写出闭环系统传递函数 111)(212211211 +=+=+ =ΦK K s K K K s K s K K s K s 令闭环增益21 2 == ΦK K ， 得：5.02=K 令调节时间4.03 32 1≤= =K K T t s ，得：151≥K 。 3-4 在许多化学过程中，反应槽内的温度要保持恒定， 图3-46（a ）和（b ）分别为开环和闭环温度控制系统结构图，两种系统正常的K 值为1。 （1） 若)(1)(t t r =，0)(=t n 两种系统从响应开始达到稳态温度值的63.2％各需多长时间？ （2） 当有阶跃扰动1.0)(=t n 时，求扰动对两种系统的温度的影响。
第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos 1000 x mx y e =，求''y syms x diff(exp(x)*cos(902*x/1000),2) ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算 22 11 00 x y e dxdy +?? dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.1394 1.4 计算4 2 2 4x dx m x +? syms x int(x^4/(902^2+4*x^2)) ans = (91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求 syms x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x) 1.6 0x =的泰勒展式(最高次幂为4).
第3章 线性系统的时域分析 学习要点 1控制系统时域响应的基本概念，典型输入信号及意义； 2控制系统稳定性的概念、代数稳定判据及应用； 3控制系统的时域指标，一阶二阶系统的阶跃响应特性与时域指标计算； 4高阶系统时域分析中主导极点和主导极点法； 5 控制系统稳态误差概念、计算方法与误差系数，减小稳态误差的方法。 思考与习题祥解 题 思考与总结下述问题。 （1）画出二阶系统特征根在复平面上分布的几种情况，归纳ξ值对二阶系统特征根的影响规律。 （2）总结ξ和n ω对二阶系统阶跃响应特性的影响规律。 （3）总结增加一个零点对二阶系统阶跃响应特性的影响规律。 （4）分析增加一个极点可能对二阶系统阶跃响应特性有何影响 （5）系统误差与哪些因素有关试归纳减小或消除系统稳态误差的措施与方法。 （6）为减小或消除系统扰动误差，可采取在系统开环传递函数中增加积分环节的措施。请问，该积分环节应在系统结构图中如何配置，抗扰效果是否与扰动点相关 答：（1）二阶系统特征根在复平面上分布情况如图所示。 图 二阶系统特征根在复平面上的分布 当0ξ=，二阶系统特征根是一对共轭纯虚根，如图中情况①。 当01ξ<<，二阶系统特征根是一对具有负实部的共轭复数根，变化轨迹是 以n ω为半径的圆弧，如图中情况②。 当1ξ=，二阶系统特征根是一对相同的负实根，如图中情况③。 当1ξ>，二阶系统特征根是一对不等的负实根，如图中情况④。
（2）ξ和n ω是二阶系统的两个特征参量。 ξ是系统阻尼比，描述了系统的平稳性。 当0ξ=，二阶系统特征根是一对共轭纯虚根，二阶系统阶跃响应为等幅振荡特性，系统临界稳定。 当01ξ<<，二阶系统特征根是一对具有负实部的共轭复数根，二阶系统阶跃响应为衰减振荡特性，系统稳定。ξ越小，二阶系统振荡性越强，平稳性越差； ξ越大，二阶系统振荡性越弱，平稳性越好。因此，二阶系统的时域性能指标超 调量由ξ值唯一确定，即001_ 100%2 ?=-π ξξ σe 。在工程设计中，对于恒值控制系 统，一般取 ξ＝～；对于随动控制系统ξ＝～。 n ω是系统无阻尼自然振荡频率，反映系统的快速性。当ξ一定，二阶系统的 时域性能指标调节时间与n ω值成反比，即34 s n t ξω≈:。 （3）二阶系统增加一个零点后，增加了系统的振荡性，将使系统阶跃响应的超调量增大，上升时间和峰值时间减小。 所增加的零点越靠近虚轴，则上述影响就越大；反之，若零点距离虚轴越远，则其影响越小。 （4）二阶系统增加一个极点后，减弱了系统的振荡性，将使系统阶跃响应的超调量减小，上升时间和峰值时间减小； 所增加的极点越靠近虚轴，则上述影响就越大；反之，若极点距离虚轴越远，则其影响越小。 （5）系统误差与系统的误差度（开环传递函数所含纯积分环节的个数或系统型别）、开环放大系数，以及作用于系统的外部输入信号有关。如果是扰动误差还与扰动作用点有关。 因此，减小或消除系统稳态误差的措施与方法有：增大开环放大系数，增加系统开环传递函数中的积分环节，引入按给定或按扰动补偿的复合控制结构。 无论采用何种措施与方法减小或消除系统稳态误差，都要注意系统须满足稳定的条件。 （6）采取在系统开环传递函数中增加积分环节的措施来减小或消除系统扰动误差时，所增加的积分环节须加在扰动作用点之前。若所增加的积分环节加在扰动作用点之后，则该积分环节无改善抗扰效果作用。这一点可以通过误差表达式分析得到。 题系统特征方程如下，试判断其稳定性。 （a ）0203.002.023=+++s s s ； （b ）014844122345=+++++s s s s s ； （c ）025266.225.11.0234=++++s s s s 解：（a ）稳定； （b ）稳定； （c ）不稳定。 题 系统结构如题图所示。控制器)1 1()(s T K s G i p c + =，为使该系统稳定，控制器参数p K 、i T 应满足什么关系