《科学分析工具MATLAB实验》实验指导书2013实验报告-哈工大完整版绝对很完整

-1-《科学分析工具MATLAB实验》实验报告(18学时系列课程实验)指导教师：学生：班级:学号:哈尔滨工业大学机电工程学院机械制造及自动化系2012年4月-1-实验一最小二乘法数据处理实验一、实验目的通过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.对指数函数在[0，1]内的采样数据进行多项式拟合。三、实验报告内容1）编写出M文件，对该数据进行合理的最小二乘法数据拟合；2）绘出原始数据点和拟合的数据曲线。1）1．clcclearxlabel('h')ylabel('g')x=10000:1000:23000;y=[24.026.529.832.434.737.741.142.844.647.365.387.5137.8174.2];p=polyfit(x,y,4)y1=polyval(p,x);plot(x,y,'+',x,y1,'-r')-2-2clcclearx=0:0.05:1;y=log10(1+x);p=polyfit(x,y,3)y1=polyval(p,x);plot(x,y,'+',x,y1,'-r')-3-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')-4-实验二Matlab程序设计实验一、实验目的通过Matlab上机编程，掌握利用Matlab软件进行多项式运算、方程求解的方法，掌握命令型M文件和函数型文件的操作方法及这两种M文件的区别与联系。二、实验内容一元二次方程的求解问题。对于一元二次方程ax2+bx+c=0，针对用户输入的不同a、b、c的值，求出方程的解情况。要求不用roots命令和solve命令，编程完成一元二次方程20axbxc的求解，给出设计的程序源代码，并针对a、b、c的不同取值情况，给出程序的运行结果。三、实验报告内容要求编写函数型M文件，输入三个系数，运用判别式法判定根的情形，并输出各种类型的根。-5-clccleara=input('pleaseinputa=');b=input('pleaseinputb=');c=input('pleaseinputc=');q=b^2-4*a*c;p=sqrt(b^2-4*a*c);ifq0disp('两个实根')x1=(-b-p)/(2*a)x2=(-b+p)/(2*a)elseifq==0disp('一个实根')x=(-b-p)/(2*a)elsedisp('无实根')x1=(-b-p)/(2*a)x2=(-b+p)/(2*a)end实验三生产计划调度优化实验一、实验目的通过Matlab上机编程，掌握利用Matlab软件进行求解优化问题的方法和常用函数、有约束优化问题。二、实验内容某车间有甲、乙、丙三台车床，可用于加工3种工件。这3台机床的可用台时数分别为720h、850h和650h，3种工件的数量分别为340个、550个和420个，已知三种不同机床加工单位数量的不同工件所需的台时数和加工费用如下表所示，问怎样分配机床的加工任务，才能既满足加工工件的要求，又使加工费用最低？表2机床加工情况表-6-机床类型单位工件加工台时数单位工件加工费用可用台时数工件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=650x甲1+x乙1+x丙1=340x乙1+x乙2+x乙3=550x丙1+x丙2+x丙3=4203）clc;clear;f=[13;9;10;14;12;8;14;13;7];A=[0.51.11.0000000;0000.51.21.3000;0000000.61.20.8];b=[720;850;650];Aeq=[100100100;010010010;001001001];beq=[340;550;420];lb=zeros(9,1);[x,fval]=linprog(f,A,b,Aeq,beq,lb)4）x=230.0000550.00000.000054.13550.00000.000055.8645-7-0.0000420.0000fval=1.2420e+004所以1号工件在甲机床上加工230件，在乙机床加工54件，丙机床加工56件。2号工件在甲机床上加工550件。3号工件在丙机床上加工420件。总费用为12420。实验四运动学分析实验一、实验目的通过Matlab上机编程，掌握Matlab语言变量、函数、Matlab常用程序结构和流程控制语句、交互控制语句、程序调试运行等。二、实验内容凸轮是把一种运动转化为另一种运动的装置，通常凸轮做旋转运动，并转化为从动件的振动、直线运动或者二者的结合。对基圆为br的凸轮，br为最小圆周半径，与凸轮表面相切，并和凸轮的转轴同心。则从动件的运动方程为：)()(srLb对于旋转运动，假定有：)2sin(21)(hs0)(2sin21)(hhs20)(s22其中，h为从动件的最大位移，为凸轮转角，为从动件摆角，并且20。如果凸轮的角速度dtd为常量，试求出当=60°时从动件的位移s、速度v、加速度a、瞬时加速度j的变化曲线。三、实验报告内容1)建立M文件，求解从动件的位移、速度、加速度及瞬时加速度；2)绘制=60°时从动件的位移s、速度v、加速度a、瞬时加速度j的变化曲线。clc;clear;b=pi/3;-8-x=linspace(0,2*pi)fork=1:1:100ifx(k)=0&x(k)=bs(k)=x(k)/b-sin(2*pi*x(k)/b)/(2*pi)v(k)=(1-cos(2*pi*x(k))/b)/ba(k)=2*pi*sin(2*pi*x(k)/b)/(b*b)j(k)=4*pi*pi/b^3*cos(2*pi*x(k)/b)elseifx(k)pi/3&x(k)2*pi/3s(k)=2-x(k)/b+sin(2*pi*(x(k)-b)/b)/(2*pi)v(k)=-(1-cos(2*pi*(x(k)-b))/b)/ba(k)=2*pi*sin(2*pi*(x(k)-b)/b)/(b*b)j(k)=4*pi*pi/b^3*cos(2*pi*(x(k)-b)/b)elses(k)=0v(k)=0a(k)=0j(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)-9-实验五机电系统仿真实验一、实验目的通过Matlab上机编程，掌握利用SIMULINK进行机电控制系统仿真，包括仿真的基本概念、机电系统建模基础、利用SIMULINK进行机电控制系统仿真的案例分析。二、实验内容如下图所示为一标有连接件编号和位置角度的曲柄滑块机构示意图，已知连杆的长度r2=15mm、r3=55mm，滑块初始位置r1=70mm，曲柄初始角度θ2(0)=θ3(0)=0º,曲柄输入角速度ω2=150rad/s，试建立系统的运动学模型，并利用SIMULINK进行系统仿真，要求根据输入转速求出滑块的速度和位移曲线。-10-三、实验报告内容1)建立系统的数学模型；2)给出利用SIMULINK进行系统仿真的结果。1）如图是个只有一个自由度的曲柄滑块机构，其输入为，输出分别为、、、。设每一连杆（包括固定杆件）均由一位移矢量表示，下图给出了该机构各个杆件之间的矢量关系，则机构的运动学方程可如下导出：曲柄滑块机构的闭环位移矢量方程闭环矢量方程的分解曲柄滑块机构的运动学方程对上面的位置方程求时间的导数，即得机构的运动学方程。-11-2）functionout=SpeedEquation(theta2,omega1,theta1)l1=15;l2=55;l0=70;b=[-l1*omega1*sin(theta1);l1*omega1*cos(theta1)];A=[l2*sin(theta2)1;-l2*cos(theta2)0];out=A\b;-12-2013年补全实验五