文档介绍：《科学分析工具MATLAB实验》
实验报告
(18学时系列课程实验)
指导教师:
学生:
班级:
学号:
哈尔滨工业大学机电工程学院机械制造及自动化系
2012年4月
实验一最小二乘法数据处理实验
一、实验目的
通过Matlab上机编程,掌握利用Matlab软件进行数据拟合分析及数据可视化方法
二、实验内容
,如下表所示,其中X为使用时间,单位为小时h,Y为磨失质量,单位为克g。要求:
编写一个M文件,对该数据进行合理的最小二乘法数据拟合;
绘出原始数据点和拟合的数据曲线。
表1 风机叶片耐磨实验
X
10000
11000
12000
13000
14000
15000
16000
Y

X
17000
18000
19000
20000
21000
22000
23000
Y

(1+x)在[0,1]内的采样数据进行多项式拟合;
-x在[0,1]内的采样数据进行多项式拟合。
三、实验报告内容
1)编写出M文件,对该数据进行合理的最小二乘法数据拟合;
2)绘出原始数据点和拟合的数据曲线。
1)1. clc
clear
xlabel('h')
ylabel('g')
x=10000:1000:23000;
y=[ ];
p=polyfit(x,y,4)
y1=polyval(p,x);
plot(x,y,'+',x,y1,'-r')
2
clc
clear
x=0::1;
y=log10(1+x);
p=polyfit(x,y,3)
y1=polyval(p,x);
plot(x,y,'+',x,y1,'-r')
3
clc
clear
x=::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 q>0
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
甲

13
9
10
720
乙

14