计算机硬件技术基础实验报告模板
东北大学
计算机硬件技术基础实验报告
专业班级：成型1002 姓名：廖力学号：20101691
实验题目：顺序程序设计日期：2012年10月9日
一、实验目的
(1) 掌握顺序程序设计方法；
(2) 学习数据传送及算术和逻辑运算指令的用法；
(3) 熟悉在PC机上建立、汇编、连接、调试和运行8086汇编语言程序的过程。
二、实验内容
设X和Y均为8位无符号数，写一个求表达式Z=X+Y值的程序。
三、算法描述流程图
开始
设置数据段地址
首地址取X
加法X+Y
和保存到Z地址
结束
四、实验步骤
1. DSEG SEGMENT ;数据段
X DW 1234H
Y DW 5678H
Z DW ?
DSEG ENDS
CSEG SEGMENT ;代码段
ASSUME CS:CSEG,DS:DSEG
START: MOV AX,DSEG ;段寄存器初值
MOV DS,AX
MOV AX,X
ADD AX,Y ;求x+y和
MOV Z,AX ;保存
MOV AH,4CH
INT 21H ;程序结束退出
CSEG ENDS
END START
2. 保存源文件，执行masm汇编源程序，生成obj目标文件，结果无错误。
3. 执行link连接程序，生成exe可执行文件，结果无错误。
4. 执行debug调试程序，使用t参数单步执行程序，在AL寄存器中得到计算结
果46H，结果正确，执行d参数命令察看Z所在的地址单元中的数据为0002字节单元中存放46H。
五、实验结果
1执行masm编译源程序得到如下界面：
2执行link连接程序，生产exe可执行文件得到如下界面：
3执行debug调试程序，使用t参数单步执行程序，在AL寄存器中得到计算机结果68ACH，结果正确，得到如下界面：
4执行d参数命令查看Z所在的地址单元中的数据为0002字节单元中存放68ACH，得到如下界面：
六、思考题
本程序没有考虑进位的情况。如果考虑可能产生的进位，应该如何修改程序？
答：如果考虑进位的情况，应该定义Z为DW类型，在程序中执行完ADD AL,Y后，执行ADDC命令。
七、扩展题
设X和Y均为三字节无符号数，编制程序求其和，并将结果放置在Z开始的存储单元。
具体程序如下：
DSEG SEGMENT ;数据段
X DB 12H,34H,56H
Y DB 34H,56H,78H
Z DB 3 DUP(?)
DSEG ENDS
CSEG SEGMENT ;代码段
ASSUME CS:CSEG,DS:DSEG
START: MOV AX,DSEG ;段寄存器初值
MOV DS,AX
MOV AL,X
ADD AL,Y ;求x+y和
MOV Z,AL ;保存
MOV AL,X+1
ADC AL,Y+1; 求x+y和
MOV Z+1，AL; 保存
MOV AL,X+2;
ADC AL,Y+2; 求x+y和
MOV Z+2，AL; 保存
MOV AH,4CH
INT 21H; 程序结束退出
CSEG ENDS
执行debug调试程序，执行d参数命令，在Z对应的储存单元中得到计算结果468ACEH,结果正确。
八、实验心得
通过本次试验，我初步理解了汇编语言中编程的基本步骤，掌握了通过地址对数据赋值的方法。掌握汇编语言编制基本顺序结构的算法，包括地址取值和加法的算法步骤等。但是在本次试验中，我更加认识到自己对于汇编语言知识的掌握不够，这就要求我加倍努力，多读课本，多向老师和同学学习，尽快弥补自身不足。
现代控制理论实验报告
实验报告 ( 2016-2017年度第二学期) 名称：《现代控制理论基础》 题目：状态空间模型分析 院系：控制科学与工程学院 班级： ___ 学号： __ 学生姓名： ______ 指导教师： _______ 成绩： 日期： 2017年 4月 15日
线控实验报告 一、实验目的： l.加强对现代控制理论相关知识的理解； 2.掌握用 matlab 进行系统李雅普诺夫稳定性分析、能控能观性分析； 二、实验内容 1 第一题：已知某系统的传递函数为G (s) S23S2 求解下列问题： （1）用 matlab 表示系统传递函数 num=[1]; den=[1 3 2]; sys=tf(num,den); sys1=zpk([],[-1 -2],1); 结果： sys = 1 ------------- s^2 + 3 s + 2 sys1 = 1 ----------- (s+1) (s+2) （2）求该系统状态空间表达式： [A1,B1,C1,D1]=tf2ss(num,den); A = -3-2 10 B = 1 C = 0 1
第二题：已知某系统的状态空间表达式为： 321 A ,B,C 01：10 求解下列问题： （1）求该系统的传递函数矩阵： （2）该系统的能观性和能空性： （3）求该系统的对角标准型： （4）求该系统能控标准型： （5）求该系统能观标准型： （6）求该系统的单位阶跃状态响应以及零输入响应：解题过程： 程序： A=[-3 -2;1 0];B=[1 0]';C=[0 1];D=0; [num,den]=ss2tf(A,B,C,D); co=ctrb(A,B); t1=rank(co); ob=obsv(A,C); t2=rank(ob); [At,Bt,Ct,Dt,T]=canon(A,B,C,D, 'modal' ); [Ac,Bc,Cc,Dc,Tc]=canon(A,B,C,D, 'companion' ); Ao=Ac'; Bo=Cc'; Co=Bc'; 结果： （1） num = 0 01 den = 1 32 （2）能控判别矩阵为： co = 1-3 0 1 能控判别矩阵的秩为： t1 = 2 故系统能控。 （3）能观判别矩阵为： ob = 0 1
计算方法上机实验报告
. / 《计算方法》上机实验报告 班级：XXXXXX 小组成员：XXXXXXX XXXXXXX XXXXXXX XXXXXXX 任课教师：XXX 二〇一八年五月二十五日
前言 通过进行多次的上机实验，我们结合课本上的内容以及老师对我们的指导，能够较为熟练地掌握Newton 迭代法、Jacobi 迭代法、Gauss-Seidel 迭代法、Newton 插值法、Lagrange 插值法和Gauss 求积公式等六种算法的原理和使用方法，并参考课本例题进行了MATLAB 程序的编写。 以下为本次上机实验报告，按照实验内容共分为六部分。 实验一： 一、实验名称及题目： Newton 迭代法 例2.7(P38):应用Newton 迭代法求在附近的数 值解，并使其满足. 二、解题思路： 设'x 是0)(=x f 的根,选取0x 作为'x 初始近似值,过点())(,00x f x 做曲线)(x f y =的切线L ,L 的方程为))((')(000x x x f x f y -+=,求出L 与x 轴交
点的横坐标) (') (0001x f x f x x - =,称1x 为'x 的一次近似值,过点))(,(11x f x 做曲线)(x f y =的切线,求该切线与x 轴的横坐标) (') (1112x f x f x x - =称2x 为'x 的二次近似值,重复以上过程,得'x 的近似值序列{}n x ,把) (') (1n n n n x f x f x x - =+称为'x 的1+n 次近似值，这种求解方法就是牛顿迭代法。 三、Matlab 程序代码： function newton_iteration(x0,tol) syms z %定义自变量 format long %定义精度 f=z*z*z-z-1; f1=diff(f);%求导 y=subs(f,z,x0); y1=subs(f1,z,x0);%向函数中代值 x1=x0-y/y1; k=1; while abs(x1-x0)>=tol x0=x1; y=subs(f,z,x0); y1=subs(f1,z,x0); x1=x0-y/y1;k=k+1; end x=double(x1) K 四、运行结果：
力学实验报告
力学实验报告 篇一：工程力学实验(全) 工程力学实验学生姓名：学号：专业班级：南昌大学工程力学实验中心目录实验一金属材料的拉伸及弹性模量测定试验实验二金属材料的压缩试验实验三复合材料拉伸实验实验四金属扭转破坏实验、剪切弹性模量测定实验五电阻应变片的粘贴技术及测试桥路变换实验实验六弯曲正应力电测实验实验七叠（组）合梁弯曲的应力分析实验实验八弯扭组合变形的主应力测定实验九偏心拉伸实验实验十偏心压缩实验实验十二金属轴件的高低周拉、扭疲劳演示实验实验十三冲击实验实验十四压杆稳定实验实验十五组合压杆的稳定性分析实验实验十六光弹性实验实验十七单转子动力学实验实验十八单自由度系统固有频率和阻尼比实验 1 2 6 9 12 16 19 23 32 37 41 45 47 49 53 59 62 65实验一金属材料的拉伸及弹性模量测定试验实验时间：设备编号：温度：湿度：一、实验目的二、实验设备和仪器三、实验数据及处理引伸仪标距l =mm 实验前 2低碳钢弹性模量测定 E? 实验后 ?F?l = (?l)?A 屈服载荷和强度极限载荷 3载荷―变形曲线(F―Δl曲线)及结果四、问题讨论（1）比较低碳钢与铸铁在拉伸时的力学性能；（2）试从不同的断口特征说明金属的两种基本破坏形式。 4篇二：工程力学实验报告工程力学实验报告自动化12级实验班 1-1 金属材料的拉伸实验一、试验目的 1．测定低碳钢（Q235 钢）的强度性能指标：上屈服强度ReH，下屈服强度ReL和抗拉强度Rm 。 2．测定低碳钢（Q235 钢）的塑性性能指标：断后伸长率A和断面收缩率Z。 3．测定铸铁的抗拉强度Rm。 4．观察、比较低碳钢（Q235 钢）和铸铁的拉伸过程及破坏现象，并比较其机械性能。 5．学习试验机的使用方法。二、设备和仪器 1．试验机（见附录）。 2．电子引伸计。 3．游标卡尺。三、试样 (a) (b) 图1-1 试样拉伸实验是材料力学性能实验中最基本的实验。为使实验结果可以相互比较，必须对试样、试验机及实验方法做出明确具体的规定。我国国标GB/T228－2002 “金属材料室温拉伸试验方法”中规定对金属拉伸试样通常采用圆形和板状两种试样，如图（1-1）所示。它们均由夹持、过渡和平行三部分组成。夹持部分应适合于试验机夹头的夹持。过渡部分的圆孤应与平行部分光滑地联接，以保证试样
现代控制理论实验
华北电力大学 实验报告| | 实验名称状态空间模型分析 课程名称现代控制理论 | | 专业班级：自动化1201 学生姓名：马铭远 学号：2 成绩： 指导教师：刘鑫屏实验日期：4月25日
状态空间模型分析 一、实验目的 1.加强对现代控制理论相关知识的理解； 2.掌握用 matlab 进行系统李雅普诺夫稳定性分析、能控能观性分析； 二、实验仪器与软件 1. MATLAB7.6 环境 三、实验内容 1 、模型转换 图 1、模型转换示意图及所用命令 传递函数一般形式： MATLAB 表示为： G=tf(num,den)，，其中 num,den 分别是上式中分子，分母系数矩阵。 零极点形式： MATLAB 表示为：G=zpk(Z,P,K) ，其中 Z，P ，K 分别表示上式中的零点矩阵，极点矩阵和增益。 传递函数向状态空间转换：[A,B,C,D] = TF2SS(NUM,DEN)； 状态空间转换向传递函数：[NUM,DEN] = SS2TF(A,B,C,D,iu)---iu 表示对系统的第 iu 个输入量求传递函数；对单输入 iu 为 1。
例1：已知系统的传递函数为G（S）= 2 2 3 24 11611 s s s s s ++ +++ ,利用matlab将传递函数 和状态空间相互转换。 解：1.传递函数转换为状态空间模型： NUM=[1 2 4];DEN=[1 11 6 11]; [A,B,C,D] = tf2ss(NUM,DEN) 2.状态空间模型转换为传递函数： A=[-11 -6 -11;1 0 0;0 1 0];B=[1;0;0];C=[1 2 4];D=[0];iu=1; [NUM,DEN] = ss2tf(A,B,C,D,iu); G=tf(NUM,DEN) 2 、状态方程状态解和输出解 单位阶跃输入作用下的状态响应： G=ss(A,B,C,D);[y,t,x]=step(G);plot(t,x). 零输入响应 [y,t,x]=initial(G,x0)其中，x0 为状态初值。
数值计算实验报告
(此文档为word格式，下载后您可任意编辑修改！) 2012级6班###(学号)计算机数值方法 实验报告成绩册 姓名:宋元台 学号： 成绩：
数值计算方法与算法实验报告 学期： 2014 至 2015 第 1 学期 2014年 12月1日课程名称: 数值计算方法与算法专业:信息与计算科学班级 12级5班 实验编号： 1实验项目Neton插值多项式指导教师：孙峪怀 姓名：宋元台学号：实验成绩： 一、实验目的及要求 实验目的： 掌握Newton插值多项式的算法，理解Newton插值多项式构造过程中基函数的继承特点，掌握差商表的计算特点。 实验要求： 1. 给出Newton插值算法 2. 用C语言实现算法 二、实验内容 三、实验步骤(该部分不够填写.请填写附页)
1.算法分析： 下面用伪码描述Newton插值多项式的算法： Step1 输入插值节点数n，插值点序列{x(i),f(i)},i=1,2,……,n,要计算的插值点x. Step2 形成差商表 for i=0 to n for j=n to i f(j)=((f(j)-f(j-1)(x(j)-x(j-1-i)); Step3 置初始值temp=1,newton=f(0) Step4 for i=1 to n temp=(x-x(i-1))*temp*由temp(k)=(x-x(k-1))*temp(k-1)形成 (x-x(0).....(x-x(i-1)* Newton=newton+temp*f(i); Step5 输出f（x）的近似数值newton(x)=newton. 2.用C语言实现算法的程序代码 #include
现代控制理论实验报告
现代控制理论实验报告
实验一系统能控性与能观性分析 一、实验目的 1.理解系统的能控和可观性。 二、实验设备 1.THBCC-1型信号与系统·控制理论及计算机控制技术实验平台； 三、实验容 二阶系统能控性和能观性的分析 四、实验原理 系统的能控性是指输入信号u对各状态变量x的控制能力，如果对于系统任意的初始状态，可以找到一个容许的输入量，在有限的时间把系统所有的状态引向状态空间的坐标原点，则称系统是能控的。 对于图21-1所示的电路系统，设iL和uc分别为系统的两个状态变量，如果电桥中 则输入电压ur能控制iL和uc状态变量的变化，此时，状态是能控的。反之，当 时，电桥中的A点和B点的电位始终相等，因而uc不受输入ur的控制，ur只能改变iL的大小，故系统不能控。 系统的能观性是指由系统的输出量确定所有初始状态的能力，如果在有限的时间根据系统的输出能唯一地确定系统的初始状态，则称系统能观。为了说明图21-1所示电路的能观性，分别列出电桥不平衡和平衡时的状态空间表达式： 平衡时：
由式（2）可知，状态变量iL和uc没有耦合关系，外施信号u只能控制iL的变化，不会改变uc的大小，所以uc不能控。基于输出是uc，而uc与iL无关连，即输出uc中不含有iL的信息，因此对uc的检测不能确定iL。反之式（1）中iL与uc有耦合关系，即ur的改变将同时控制iL和uc的大小。由于iL与uc的耦合关系，因而输出uc的检测，能得到iL 的信息，即根据uc的观测能确定iL（ω） 五、实验步骤 1．用2号导线将该单元中的一端接到阶跃信号发生器中输出2上，另一端接到地上。将阶跃信号发生器选择负输出。 2．将短路帽接到2K处，调节RP2，将Uab和Ucd的数据填在下面的表格中。然后将阶跃信号发生器选择正输出使调节RP1，记录Uab和Ucd。此时为非能控系统，Uab和Ucd没有关系（Ucd始终为0）。 3．将短路帽分别接到1K、3K处，重复上面的实验。 六、实验结果 表20－1Uab与Ucd的关系 Uab Ucd
计算方法实验报告
计算方法实验报告（四） 方程和方程组的迭代解法 一、实验问题 利用简单迭代法，两种加速技术，牛顿法，改进牛顿法，弦割法求解习题5-1，5-2，5-3中的一题，并尽可能准确。 选取5-3：求在x=1.5附近的根。 二、问题的分析（描述算法的步骤等） （1）简单迭代法算法： 给定初始近似值,求的解。 Step 1 令i=0； Step 2 令(计算)； Step 3 如果,则迭代终止，否则重复Step 2。 （2）Aitken加速法算法 Step 1 令k=0，利用简单迭代算法得到迭代序列； Step 2 令-(计算得到一个新的序列，其中k=0,1,2…)；Step 3 如果,则迭代终止，否则重复Step 2。 (3)插值加速法算法 Step 1 令k=0，利用简单迭代算法得到迭代序列； Step 2 令+(计算得到一个新的序列，其中k=1,2,3…)； Step 3 如果,则迭代终止，否则重复Step 2。 (4)牛顿法算法
Step 1给定初始近似值； Step 2令,其中k计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 （5）改进牛顿法的算法 Step 1给定初始近似值； Step 2令,其中k迭代计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 （6）弦割法算法（双点弦割法） Step 1给定初始近似值，； Step 2令其中k计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 三、程序设计 (1)简单迭代法 利用迭代公式进行迭代运算。 #include
计算方法实验报告格式
计算方法实验报告格式 小组名称: 组长姓名(班号)： 小组成员姓名(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一个完整的实验，应包括数据准备、理论基础、实验内容及方法，最终对实验结果进行分析，以达到对理论知识的感性认识，进一步加深对相关算法的理解，数值实验以实验报告形式完成，实验报告格式如下： 一、实验名称 实验者可根据报告形式需要适当写出． 二、实验目的及要求 首先要求做实验者明确，为什么要做某个实验，实验目的是什么，做完该实验应达到什么结果，在实验过程中的注意事项，实验方法对结果的影响也可以以实验目的的形式列出． 三、算法描述（实验原理与基础理论） 数值实验本身就是为了加深对基础理论及方法的理解而设置的，所以要求将实验涉及到的理论基础，算法原理详尽列出． 四、实验内容 实验内容主要包括实验的实施方案、步骤、实验数据准备、实验的算法以及可能用到的仪器设备． 五、程序流程图 画出程序实现过程的流程图，以便更好的对程序执行的过程有清楚的认识，在程序调试过程中更容易发现问题． 六、实验结果 实验结果应包括实验的原始数据、中间结果及实验的最终结果，复杂的结果可以用表格
形式列出，较为简单的结果可以与实验结果分析合并出现． 七、实验结果分析 实验结果分析包括对对算法的理解与分析、改进与建议． 数值实验报告范例 为了更好地做好数值实验并写出规范的数值实验报告，下面给出一简单范例供读者参考． 数值实验报告 小组名称: 小组成员(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一、实验名称 误差传播与算法稳定性． 二、实验目的 1．理解数值计算稳定性的概念． 2．了解数值计算方法的必要性． 3．体会数值计算的收敛性与收敛速度． 三、实验内容 计算dx x x I n n ? += 1 10 ，1,2,,10n = ． 四、算法描述 由 dx x x I n n ? += 1 10 ，知 dx x x I n n ?+=--101110，则
土力学实验报告
园林学院 土力学实验报告 学生姓名 学号2009041001 专业班级土木工程091 指导教师李西斌 组别第三组 成绩
实验目录 前言 (1) 实验一含水量试验 (2) 实验二密度实验 (5) 实验三液限和塑限试验 (7) 实验四固结试验 (13) 实验五直接剪切试验 (18)
前言 土是矿物颗粒所组成的松散颗粒集合体，其物理力学性质与其他材料不同；土力学是利用力学的基本原理和土工试验技术来研究土的强度和变形及其规律性的一门应用学科。 土的天然含水率、击实性、压缩性、抗剪强度是水利工程中的四大问题，他们的好坏与否直接关系到水利工程的经济效益与安全问题，因此在工程中作好土料的指标实验，确定出相应标对水利工程具有十分重要的意义。
实验一 含水量试验 一、概述 土的含水率 是指土在温度105~110℃下烘干至恒量时所失去的水质量与达 到恒量后干土质量的比值，以百分数表示。 含水率是土的基本物理性质指标之一，它反映了土的干、湿状态。含水率的变化将使土物理力学性质发生一系列变化，它可使土变成半固态、可塑状态或流动状态，可使土变成稍湿状态、很湿状态或饱和状态，也可造成土在压缩性和稳定性上的差异。含水率还是计算土的干密度、孔隙比、饱和度、液性指数等不可缺少的依据，也是建筑物地基、路堤、土坝等施工质量控制的重要指标。 二、实验原理 土样在在105℃～110℃温度下加热，土中自由水会变成气体挥发，土恒重后， 即可认为是干土质量s m ，挥发掉的水分质量为w s m m m =-。 三、实验目的 测定土的含水量，供计算土的孔隙比、液性指数、饱和度等不可缺少的一个基本指标。并查表可确定地基土的允许承载力 四、实验方法 含水率实验方法有烘干法、酒精燃烧法、比重法、碳化钙气压法、炒干法等，其中以烘干法为室内实验的标准方法。在此仅用烘干法来测定。 烘 烘干法是将实样放在温度能保持105~110℃的烘箱中烘至恒量的方法，是室内测定含水率的标准方法。 （一）仪器设备 （1）保持温度为105~110℃的自动控制电热恒温烘箱； （2）称量200g 、最小分度值0.01g 的天平； （3）玻璃干燥缸；
现代控制理论实验报告河南工业大学
河南工业大学 现代控制理论实验报告姓名：朱建勇 班级：自动1306 学号：201323020601
现代控制理论 实验报告 专业: 自动化 班级: 自动1306 姓名: 朱建勇 学号: 201323020601 成绩评定: 一、实验题目： 线性系统状态空间表达式的建立以及线性变换 二、实验目的 1. 掌握线性定常系统的状态空间表达式。学会在MATLAB 中建立状态空间模型的方法。 2. 掌握传递函数与状态空间表达式之间相互转换的方法。学会用MATLAB 实现不同模型之 间的相互转换。 3. 熟悉系统的连接。学会用MATLAB 确定整个系统的状态空间表达式和传递函数。 4. 掌握状态空间表达式的相似变换。掌握将状态空间表达式转换为对角标准型、约当标准 型、能控标准型和能观测标准型的方法。学会用MATLAB 进行线性变换。 三、实验仪器 个人笔记本电脑 Matlab R2014a 软件 四、实验内容 1. 已知系统的传递函数 (a) ) 3()1(4)(2++=s s s s G
(b) 3486)(22++++=s s s s s G
(c) 6 1161)(232+++++=z z z z z z G （1）建立系统的TF 或ZPK 模型。 （2）将给定传递函数用函数ss( )转换为状态空间表达式。再将得到的状态空间表达式用函 数tf( )转换为传递函数，并与原传递函数进行比较。 （3）将给定传递函数用函数jordants( )转换为对角标准型或约当标准型。再将得到的对角 标准型或约当标准型用函数tf( )转换为传递函数，并与原传递函数进行比较。 （4）将给定传递函数用函数ctrlts( )转换为能控标准型和能观测标准型。再将得到的能控标 准型和能观测标准型用函数tf( )转换为传递函数，并与原传递函数进行比较。
数学计算方法实验报告
数学计算方法实验报告 习题二 2.估计用二分法求方程f(x)=x3+4x2-10=0在区间[1,2]内根的近似值，为使方程不超过10时所需的二分次数。f(x k) 程序过程： function two (tolerance) a=1;b=2;counter=0; while (abs(b-a)>tolerance) c=(a+b)/2; fa=a^3+4*a^2-10;
fb=b^3+4*b^2-10; fc=c^3+4*c^2-10; if ((fa==0|fb==0)) disp(counter); elseif (fa*fc<0) b=c;counter=counter+1; elseif (fb*fc<0) a=c;counter=counter+1; elseif (fb==0) disp(counter); end end solution=(a+b)/2; disp(solution); disp(counter); 实验结果： 6.取x0=1.5,用牛顿迭代法求第三中的方程根.f(x)=x3+4x2-10=0的近似值（精确到||x k+1-x k|≦10-5,并将迭代次数与3题比较。 程序过程： function six (g) a=1.5; fa=a^3+4*a^2-10;
ga=3*a^2+8*a; b=a-fa/ga; k=1; while(abs(b-a)>g) a=b; fa=a^3+4*a^2-10; ga=3*a^2+8*a; b=a-fa/ga; k=k+1; end format long; disp(a); disp(k); 实验结果：程序结果计算结果 8.用弦割法求方程f(x)=x3-3x2-x+9=0在区间[-2,-1]内的一个实根近似值x k,|f(x k)|≦10-5. 程序过程： function eight (t) a=-2; b=-1; fa=a^3-3*a^2-a+9; fb=b^3-3*b^2-b+9; c=b-fb*(b-a)/(fb-fa); k=1; while(abs(c-b)>t) a=b; b=c; fa=a^3-3*a^2-a+9; fb=b^3-3*b^2-b+9; c=b-fb*(b-a)/(fb-fa); k=k+1; end
现代控制理论课程报告
现代控制理论课程总结 学习心得 从经典控制论发展到现代控制论,是人类对控制技术认识上的一次飞跃。现代控制论是用状态空间方法表示,概念抽象,不易掌握。对于《现代控制理论》这门课程，在刚拿到课本的时候，没上张老师的课之前，咋一看，会认为开课的内容会是上学期学的控制理论基础的累赘或者简单的重复，更甚至我还以为是线性代数的复现呢！根本没有和现代控制论联系到一起。但后面随着老师讲课的风格的深入浅出，循循善诱，发现和自己想象的恰恰相反，张老师以她特有的讲课风格，精心准备的ppt 课件，向我们展示了现代控制理论发展过程，以及该掌握内容的方方面面，个人觉得，我们不仅掌握了现代控制理论的理论知识，更重要的是学会了掌握这门知识的严谨的逻辑思维和科学的学习方法，对以后学习其他知识及在工作上的需要大有裨益，总之学习了这门课让我受益匪浅。 由于我们学习这门课的课时不是很多，并结合我们学生学习的需求及所要掌握的课程深入程度，张老师根据我们教学安排需要，我们这学期学习的内容主要有：1.绪论；2.控制系统的状态表达式；3.控制系统状态表达式的解；4.线性系统的能空性和能观性；5.线性定常系统的综合。而状态变量和状态空间表达式、状态转移矩阵、系统的能控性与能观性以及线性定常系统的综合是本门课程的主要学习内容。当然学习的内容还包括老师根据多年教学经验及对该学科的研究的一些深入见解。 在现代科学技术飞速发展中,伴随着学科的高度分化和高度综合,各学科之间相互交叉、相互渗透,出现了横向科学。作为跨接于自然科学和社会科学的具有横向科学特点的现代控制理论已成为我国理工科大学高年级的必修课。 经典控制理论的特点 经典控制理论以拉氏变换为数学工具，以单输入－单输出的线性定常系统为主要的研究对象。将描述系统的微分方程或差分方程变换到复数域中，得到系统的传递函数，并以此作为基础在频率域中对系统进行分析和设计，确定控制器的结构和参数。通常是采用反馈控制，构成所谓闭环控制系统。经典控制理论具有明显的局限性，突出的是难以有效地应用于时变系统、多变量系统，也难以揭示系统更为深刻的特性。当把这种理论推广到更为复杂的系统时，经典控制理论就显得无能为力了，这是因为它的以下几个特点所决定。 1．经典控制理论只限于研究线性定常系统，即使对最简单的非线性系统也是无法处理的；这就从本质上忽略了系统结构的内在特性，也不能处理输入和输出皆大于1的系统。实际上，大多数工程对象都是多输入－多输出系统，尽管人们做了很多尝试，但是，用经典控制理论设计这类系统都没有得到满意的结果；2．经典控制理论采用试探法设计系统。即根据经验选用合适的、简单的、工程上易于实现的控制器，然后对系统进行分析，直至找到满意的结果为止。虽然这种设计方法具有实用等很多完整，从而促使现代控制理论的发展：对经典理论的精确化、数学化及理论化。优点，但是，在推理上却是不能令人满意的，效果也
计算方法实验报告
实验报告 一、求方程f(x)=x^3-sinx-12x+1的全部根, ε=1e -6 1、 用一般迭代法; 2、 用牛顿迭代法; 并比较两种迭代的收敛速度。 一、首先，由题可求得：12cos 3)(2 ' --=x x x f . 其次，分析得到其根所在的区间。 ① 令()0=x f ，可得到x x x sin 1123 =+-. ② 用一阶导数分析得到1123 +-x x 和x sin 两个函数的增减区间；再用二阶导数分析得到 两个函数的拐点以及凹凸区间. ③ 在直角坐标轴上描摹出01123 =+-x x 和0sin =x 的图，在图上可以看到他们的交点，然后估计交点所在的区间，即是所要求的根的区间。经过估计，得到根所在的区间为 []3,4--，[]1,0和[]4,3. 1、 一般迭代法 （1）算法步骤： 设ε为给定的允许精度，迭代法的计算步骤为： ① 选定初值0x .由()0=x f 确定函数()x g ，得等价形式()x g x =. ② 计算()0x g .由迭代公式得()01x g x =. ③ 如果ε≤-01x x ，则迭代结束，取1x 为解的近似值；否则，用1x 代替0x ，重复步骤②和步骤③. （2）程序代码： ① 在区间[]3,4--内， 代码： clc
x0=-3.5; %初值0x iter_max=100; %迭代的最大次数 ep=1e-6; %允许精度 ε k=0; while k<=iter_max %k 从0开始到iter_max 循环 x1=(sin(x0)+12*x0-1).^(1/3); %代入0x ，算出1x 的值 if abs(x1-x0) 试验目的： 1. 测定低碳钢（塑性材料）的弹性摸量E；屈服极限σs 等机械性能。 2．测定灰铸铁（脆性材料）的强度极限σb 3．了解塑性材料和脆性材料压缩时的力学性能。 材料拉伸与压缩实验指导书 低碳钢拉伸试验 拉伸试验的意义: 单向拉伸试验是在常温下以缓慢均匀的速度对专门制备的试件施加轴向载荷，在试件加载过程中观测载荷与变形的关系，从而决定材料有关力学性能。通过拉伸试验可以测定材料在单向拉应力作用下的弹性模量及屈服强度、抗拉强度、延伸率、截面收缩率等指标。其试验方法简单且易于得到较可靠的试验数据，所以是研究材料力学性能最基本、应用最广泛的试验。 操作步骤： 1．试验设备：WDW－3050电子万能试验机 2．试件准备：用游标卡尺测量试件试验段长度l0和截面直径d0，并作记录。 3．打开试验机主机及计算机等相关设备。 4．试件安装（详见WDW3050电子万能试验机使用与操作三．拉伸试件的安装）。 5．引伸计安装（用于测量E, 详见WDW3050电子万能试验机使用与操作四．引伸计安装）。 6．测量参数的设定: 7．再认真检查一遍试件安装等试验准备工作。 8．负荷清零，轴向变形清零，位移清零。 9．开始进行试验，点击试验开始。 10．根据提示摘除引伸计。 11．进入强化阶段以后，进行冷作硬化试验，按主机控制面板停止，再按▼，先卸载到10kN，再加载，按▲，接下来计算机控制，一直到试件断裂（此过程中计算机一直工作，注意观察负荷位移曲线所显示的冷作硬化现象．）． 12．断裂以后记录力峰值。 13．点击试验结束（不要点击停止）。 14．材料刚度特征值中的弹性模量E的测定 试验结束后，在试验程序界面选定本试验的试验编号，并选择应力─应变曲线。在曲线上较均匀地选择若干点，记录各点的值，分别为及 (如i =0,1,2,3,4),并计算出相应的 计算E i的平均值，得到该材料的弹性模量E的值。 15．材料强度特征值屈服极限和强度极限的测定 试验结束后，在试验程序界面选定本试验的试验编号，并选择负荷─位移曲线，找到的曲线屈服阶段的下屈服点，即为屈服载荷F s, 找到的曲线上最大载荷值，即为极限载荷P b. 计算屈服极限：；计算强度极限：； 16．材料的塑性特征值延伸率及截面收缩率的测定 试件拉断后，取下试件，沿断裂面拼合，用游标卡尺测定试验段长度，和颈缩断裂处截面直径。 计算材料延伸率 计算截面收缩率 低碳钢拉伸试验报告 试验目的： 1. 掌握电子万能试验机操作； 2. 理解塑性材料拉伸时的力学性能； 3. 观察低碳钢拉伸时的变形特点； 4. 观察低碳钢材料的冷作硬化现象； 5. 测定低碳钢材料弹性模量E ； 6. 测定材料屈服极限和强度极限； 7. 测定材料伸长率δ和截面收缩率Ψ 试验设备： 实验三 利用MATLAB 导出连续状态空间模型的离散化模型 实验目的： 1、基于对象的一个连续时间状态空间模型，导出其相应的离散化状态空间模型； 2、通过编程、上机调试，掌握离散系统运动分析方法。 实验原理： 给定一个连续时间系统的状态空间模型： ()()()()()() x t Ax t Bu t y t Cx t Du t =+=+ （3.1） 状态空间模型（3.1）的输入信号()u t 具有以下特性： ()(),u t u kT kT t kT T =≤≤+ （3.2） 已知第k 个采样时刻的状态()x kT 和第k 个采样时刻到第1k +个采样时刻间的输入()()u t u kT =，可得第1k +个采样时刻(1)k T +处的状态 (1)((1))((1))()((1))()k T kT x k T k T kT x kT k T Bu d τττ++=Φ+-+Φ+-? （3.3） 其中： ((1))((1))A k T kT AT k T kT e e +-Φ+-== ((1))((1))A k T k T e ττ+-Φ+-= 由于输入信号在两个采样时刻之间都取常值，故对式（3.3）中的积分式进行一个时间变量替换(1)k T στ=+-后，可得 0((1))()()()AT A x k T e x kT e d Bu kT τ σσ+=+? （3.4） 另一方面，以周期T 对输出方程进行采样，得到 ()()()y kT Cx kT Du kT =+ 在周期采样的情况下，用k 来表示第k 个采样时刻kT 。因此，连续时间状态空间模型 （3.1）的离散化方程可以写成 (1)()()()()()()() x k G T x k H T u k y k Cx k Du k +=+=+ （3.5） 其中： 0()()()AT A G T e H T e d B τσσ==? （3.6） 已知系统的连续时间状态空间模型，MATLAB 提供了计算离散化状态空间模型中状态矩阵和输入矩阵的函数： [G ,H]=c2d(A,B,T) 其中的T 是离散化模型的采样周期。 实验步骤 1、导出连续状态空间模型的离散化模型，采用MA TLAB 的m-文件编程； 2、在MA TLAB 界面下调试程序，并检查是否运行正确。 例3.1 已知一个连续系统的状态方程是 010()()()2541x t x t u t ????=+????--???? 若取采样周期0.05T =秒，试求相应的离散化状态空间模型。 编写和执行以下的m-文件： A=[0 1;-25 –4]; B=[0;1]; [G ,H]=c2d(A,B,0.05) 得到 G= 0.9709 0.0448 -1.1212 0.7915 H= 0.0012 0.0448 因此，所求的离散化状态空间模型是 0.97090.04480.0012(1)()()1.12120.79150.0448x k x k u k ????+=+????-???? 材料力学实验报告答案 Prepared on 22 November 2020 材料力学实验报告 评分标准 拉伸实验报告 一、实验目的（1分） 1. 测定低碳钢的强度指标（σs、σb）和塑性指标（δ、ψ）。 2. 测定铸铁的强度极限σb。 3. 观察拉伸实验过程中的各种现象，绘制拉伸曲线（P－ΔL曲线）。 4. 比较低碳钢与铸铁的力学特性。 二、实验设备（1分） 机器型号名称电子万能试验机 测量尺寸的量具名称游标卡尺精度0.02 mm 三、实验数据（2分） 四、实验结果处理 （4分） 0A P s s = σ ＝300MPa 左右 0 A P b b = σ ＝420MPa 左右 ％10000 1?-= L L L δ ＝20～30％左右 ％＝ 1000 1 0?-A A A ψ ＝60～75％左右 五、回答下列问题（2分，每题分） 1、画出（两种材料）试件破坏后的简图。 略 2、画出拉伸曲线图。 3、试比较低碳钢和铸铁拉伸时的力学性质。 低碳钢在拉伸时有明显的弹性阶段、屈服阶段、强化阶段和局部变形阶段，而铸铁没有明显的这四个阶段。 4、材料和直径相同而长短不同的试件，其延伸率是否相同为什么 相同 延伸率是衡量材料塑性的指标，与构件的尺寸无关。 压缩实验报告 一、实验目的（1分） 1. 测定压缩时铸铁的强度极限σb 。 2. 观察铸铁在压缩时的变形和破坏现象，并分析原因。 二、实验设备 （1分） 机器型号名称电子万能试验机 （分） 测量尺寸的量具名称 游标卡尺 精度 0.02 mm （分） 三、实验数据（1分） 四、实验结果处理 （2分） A P b b = σ =740MPa 左右 五、回答下列思考题（3分） 1.画出（两种材料）实验前后的试件形状。 略 2. 绘出两种材料的压缩曲线。 略 3. 为什么在压缩实验时要加球形承垫 实验六利用MATLAB设计状态观测器 ******* 学号 1121***** 实验目的： 1、学习观测器设计算法； 2、通过编程、上机调试，掌握基于观测器的输出反馈控制系统设计方法。 实验原理： 1、全阶观测器模型： () ()x Ax Bu L y Cx A LC x Bu Ly =++-=-++ 由极点配置和观测器设计问题的对偶关系，也可以应用MATLAB 中极点配置的函数来确定所需要的观测器增益矩阵。例如，对于单输入单输出系统，观测器的增益矩阵可以由函数 L=(acker(A ’,C ’,V))’ 得到。其中的V 是由期望的观测器极点所构成的向量。类似的，也可以用 L=(place(A ’,C ’,V))’ 来确定一般系统的观测器矩阵，但这里要求V 不包含相同的极点。 2、降阶观测器模型： ???w Aw By Fu =++ b x w Ly =+ 基于降阶观测器的输出反馈控制器是： ????()[()]()b a b b a b w A FK w B F K K L y u K w K K L y =-+-+=--+ 对于降阶观测器的设计，使用MATLAB 软件中的函数 L=(acker(Abb’,Aab’,V))’ 或 L=(place(Abb’,Aab’,V))’ 可以得到观测器的增益矩阵L 。其中的V 是由降阶观测器的期望极点所组成的向量。 实验要求 1．在运行以上例程序的基础上，考虑图6.3所示的调节器系统，试针对被控对象设计基于全阶观测器和降 阶观测器的输出反馈控制器。设极点配置部分希望的闭环极点是1,22j λ=-± （a ） 对于全阶观测器，1 8μ=-和 28μ=-； （b ） 对于降阶观测器，8μ=-。 比较系统对下列指定初始条件的响应： （a ） 对于全阶观测器： 1212(0)1,(0)0,(0)1,(0)0x x e e ==== （b ） 对于降阶观测器： 121(0)1,(0)0,(0)1x x e === 进一步比较两个系统的带宽。 《数据分析方法与技术》上机实验——实验1描述性统计方法 学号： 姓名： 日期： 实验项目（一）：描述性统计方法 一、实验内容 1．实验目的 掌握常用的描述性图表展示方法的原理及操作，包括：频数分布表、分组频数表、列联表、茎叶图、箱线图、误差图、散点图等； 掌握常用的描述性统计方法的原理及操作，包括：算术平均值、中位数、众数、四分位数、极差、平均差、方差、标准差、标准分数、离散系数等。 2. 实验内容和要求 实验内容：基于标准数据集，属性描述性图表展示方法（数分布表、分组频数表、列联表、茎叶图、箱线图、误差图、散点图等），对统计指标（算术平均值、中位数、众数、极差、平均差、方差、标准差、标准分数、离散系数、偏态峰态）进行计算。 实验要求：掌握各种描述性统计指标的计算思路及其在SPSS或EXCEL环境下的操作方法，掌握输出结果的解释。 二、实验过程 1、数据集介绍 1.数据库标题:鲍鱼数据 2.该数据库共计4177行数据 3.该数据有八个属性（包含性别共有九项） 4.以下是关于属性的描述，包括属性的名称，数据类型，测量单元和一个简短的描述： Name Data TypeMeas.Description ---- --------- ----- ----------- Sex nominal M, F, and I (infant)鲍鱼宝宝 Length continuousmm Longest shell measurement最长壳 Diameter continuousmm perpendicular to length垂直长度 Height continuousmm with meat in shell有肉的壳高度 Whole weightcontinuousgramswhole abalone整个鲍鱼 Shucked weightcontinuousgramsweight of meat肉的重量 Viscera weightcontinuousgramsgut weight (after bleeding)放血后内脏重 Shell weightcontinuousgramsafter being dried弄干后重量 Rings integer +1.5 gives the age in years +1.5=年龄 5.数据的值域 . 现代控制理论实验报告 组员： 院系：信息工程学院 专业： 指导老师： 年月日 实验1 系统的传递函数阵和状态空间表达式的转换 [实验要求] 应用MATLAB 对系统仿照[例1.2]编程，求系统的A 、B 、C 、阵；然后再仿照[例1.3]进行验证。并写出实验报告。 [实验目的] 1、学习多变量系统状态空间表达式的建立方法、了解系统状态空间表达式与传递函数相互转换的方法； 2、通过编程、上机调试，掌握多变量系统状态空间表达式与传递函数相互转换方法。 [实验内容] 1 设系统的模型如式(1.1)示。 p m n R y R u R x D Cx y Bu Ax x ∈∈∈?? ?+=+= (1.1) 其中A 为n ×n 维系数矩阵、B 为n ×m 维输入矩阵 C 为p ×n 维输出矩阵，D 为传递阵，一般情况下为0，只有n 和m 维数相同时，D=1。系统的传递函数阵和状态空间表达式之间的关系如式(1.2)示。 D B A SI C s den s num s G +-== -1)() () (()( (1.2) 式(1.2)中，)(s num 表示传递函数阵的分子阵，其维数是p ×m ；)(s den 表示传递函数阵的按s 降幂排列的分母。 2 实验步骤 ① 根据所给系统的传递函数或（A 、B 、C 阵），依据系统的传递函数阵和状态空间表达式之间的关系如式(1.2)，采用MATLA 的file.m 编程。注意：ss2tf 和tf2ss 是互为逆转换的指令； ② 在MATLA 界面下调试程序，并检查是否运行正确。 ③ [1.1] 已知SISO 系统的状态空间表达式为(1.3)，求系统的传递函数。工程力学拉伸实验报告
现代控制理论实验报告3
材料力学实验报告答案
利用MATLAB设计状态观测器—现代控制理论实验报告
数据分析方法与技术- 实验报告模板
《现代控制理论》实验报告