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

《科学分析工具MATLAB实验》 实验报告 (18学时系列课程实验 指导教师： 学 生： 班 级: 学 号: 哈尔滨工业大学机电项目学院机械制造及自动化系 2018年4月 实验一 最小二乘法数据处理实验 一、实验目的 通过Matlab上机编程，掌握利用Matlab软件进行数据拟合分析及数据可视化方法 二、实验内容 1.有一组风机叶片的耐磨实验数据，如下表所示，其中X为使用时间，单位为小时h，Y为磨失质量，单位为克g。要求： 编写一个M文件，对该数据进行合理的最小二乘法数据拟合； 绘出原始数据点和拟合的数据曲线。 表1 风机叶片耐磨实验 X 10000 11000 12000 13000 14000 15000 16000 Y 24.0 26.5 29.8 32.4 34.7 37.7 41.1 X 17000 18000 19000 20000 21000 22000 23000 Y 42.8 44.6 47.3 65.8 87.5 137.8 174.2 2.对ln1+x）在[0，1]内的采样数据进行多项式拟合； 3.对指数函数1-x在 三、实验报告内容 1）编写出M文件，对该数据进行合理的最小二乘法数据拟合； 2）绘出原始数据点和拟合的数据曲线。 1）1． clc clear xlabel(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 clc clear x=0:0.05:1。 y=log10(1+x。 p=polyfit(x,y,3 y1=polyval(p,x。 plot(x,y,+,x,y1,-r 3 clc clear x=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文件，输入三个系数，运用判别式法判定根的情形，并输出各种类型的根。 clc clear a=input(please input a=。 b=input(please input b=。 c=input(please input c=。 q=b^2-4*a*c。 p=sqrt(b^2-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.5 1.1 1.0 13 9 10 720 乙 0.5 1.2 1.3 14 12 8 850 丙 0.6 1.3 0.8 14 13 7 650 三、实验报告内容 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