实验指导书2013实验报告-哈工大完整版绝对很完整.doc

科学分析工具MATLAB实验实验报告(18学时系列课程实验) 指导教师： 学 生： 班 级: 学 号: 哈尔滨工业大学机电工程学院 机械制造及自动化系2012年4月- 11 -实验一 最小二乘法数据处理实验一、实验目的通过Matlab上机编程，掌握利用Matlab软件进行数据拟合分析及数据可视化方法二、实验内容1.有一组风机叶片的耐磨实验数据，如下表所示，其中X为使用时间，单位为小时h，Y为磨失质量，单位为克g。要求：(1) 编写一个M文件，对该数据进行合理的最小二乘法数据拟合；(2) 绘出原始数据点和拟合的数据曲线。 表1 风机叶片耐磨实验X10000110001200013000140001500016000Y24.026.529.832.434.737.741.1X17000180001900020000210002200023000Y42.844.647.365.887.5137.8174.22.对ln（1+x）在0，1内的采样数据进行多项式拟合；3.对指数函数1-x在0，1内的采样数据进行多项式拟合。三、实验报告内容1）编写出M文件，对该数据进行合理的最小二乘法数据拟合；2）绘出原始数据点和拟合的数据曲线。1）1 clcclearxlabel(h)ylabel(g)x=10000:1000:23000; y=24.0 26.5 29.8 32.4 34.7 37.7 41.1 42.8 44.6 47.3 65.3 87.5 137.8 174.2; p=polyfit(x,y,4) y1=polyval(p,x); plot(x,y,+,x,y1,-r)2 clcclearx=0:0.05:1; y=log10(1+x);p=polyfit(x,y,3) y1=polyval(p,x);plot(x,y,+,x,y1,-r)3clcclearx=0.05:0.05:1; y=1-sqrt(x);p=polyfit(x,y,3) y1=polyval(p,x);plot(x,y,+,x,y1,-r)实验二 Matlab程序设计实验一、实验目的 通过Matlab上机编程，掌握利用Matlab软件进行多项式运算、方程求解的方法，掌握命令型M文件和函数型文件的操作方法及这两种M文件的区别与联系。二、实验内容一元二次方程的求解问题。对于一元二次方程ax2+bx+c=0，针对用户输入的不同a、b、c的值，求出方程的解情况。要求不用roots命令和solve命令，编程完成一元二次方程的求解，给出设计的程序源代码，并针对a、b、c的不同取值情况，给出程序的运行结果。三、实验报告内容要求编写函数型M文件，输入三个系数，运用判别式法判定根的情形，并输出各种类型的根。clccleara=input(please input a=);b=input(please input b=);c=input(please input c=);q=b2-4*a*c;p=sqrt(b2-4*a*c);if q0 disp( 两个实根) x1=(-b-p)/(2*a) x2=(-b+p)/(2*a)else if q=0 disp( 一个实根) x=(-b-p)/(2*a)else disp(无实根)x1=(-b-p)/(2*a)x2=(-b+p)/(2*a)end实验三 生产计划调度优化实验一、实验目的通过Matlab上机编程，掌握利用Matlab软件进行求解优化问题的方法和常用函数、有约束优化问题。二、实验内容某车间有甲、乙、丙三台车床，可用于加工3种工件。这3台机床的可用台时数分别为720h、850h和650h，3种工件的数量分别为340个、550个和420个，已知三种不同机床加工单位数量的不同工件所需的台时数和加工费用如下表所示，问怎样分配机床的加工任务，才能既满足加工工件的要求，又使加工费用最低？ 表2 机床加工情况表机床类型单位工件加工台时数单位工件加工费用可用台时数工件1工件2工件3工件1工件2工件3甲0.51.11.013910720乙0.51.21.314128850丙0.61.30.814137650三、实验报告内容1）建立线性规划问题；2）建立约束条件；3）编写m文件；4）给出运行结果。1） f= 13x甲1+9x甲2+10x甲3+14x乙1+12x乙2+8x乙3+14x丙1+13x丙2+7x丙3其中x甲1代表在甲机床上加工1号工件，其余类推。2）0.5*x甲1+ 1.1*x乙1+1.0* x丙1=7200.5*x乙1+ 1.2*x乙2+1.3* x乙3=8500.6*x丙1+ 1.3x丙2+ 0.8*x丙3=0 & x(k)pi/3 & x(k)2*pi/3 s(k)=2-x(k)/b+sin(2*pi*(x(k)-b)/b)/(2*pi) v(k)=-(1-cos(2*pi*(x(k)-b)/b)/b a(k)=2*pi*sin(2*pi*(x(k)-b)/b)/(b*b) j(k)=4*pi*pi/b3*cos(2*pi*(x(k)-b)/b)else s(k)=0 v(k)=0 a(k)=0 j(k)=0endendsubplot(2,2,1)plot(x,s)subplot(2,2,2)plot(x,v)subplot(2,2,3)plot(x,a)subplot(2,2,4)plot(x,j)实验五 机电系统仿真实验一、 实验目的 通过Matlab上机编程，掌握利用SIMULINK进行机电控制系统仿真，包括仿真的基本概念、机电系统建模基础、利用SIMULINK进行机电控制系统仿真的案例分析。二、实验内容 如下图所示为一标有连接件编号和位置角度的曲柄滑块机构示意图，已知连杆的长度r2=15 mm、r3=55 mm，滑块初始位置r1=70 mm，曲柄初始角度2(0)= 3(0)=0, 曲柄输入角速度2=150 rad/s，试建立系统的运动学模型，并利用SIMULINK进行系统仿真，要求根据输入转速求出滑块的速度和位移曲线。三、实验报告内容1) 建立系统的数学模型；2) 给出利用SIMULINK进行系统仿真的结果。1） 如图是个只有一个自由度的曲柄滑块机构，其输入为2=2，输出分别为3=3、3、v1=r1、r1。设每一连杆（包括固定杆件）均由一位移矢量表示，下图给出了该机构各个杆件之间的矢量关系，则机构的运动学方程可如下导出：曲柄滑块机构的闭环位移矢量方程R2+R3=R1闭环矢量方程的分解r2cos2+r3cos3=r1r2sin2+r3sin3=0曲柄滑块机构的运动学方程对上面的位置方程求时间的导数，即得机构的运动学方程。2）function out=SpeedEqu