《科学分析工具MATLAB实验》实验指导书实验报告哈工大完整版绝对很完整.docx
《《科学分析工具MATLAB实验》实验指导书实验报告哈工大完整版绝对很完整.docx》由会员分享，可在线阅读，更多相关《《科学分析工具MATLAB实验》实验指导书实验报告哈工大完整版绝对很完整.docx（11页珍藏版）》请在冰点文库上搜索。
《科学分析工具MATLAB实验》实验指导书实验报告哈工大完整版绝对很完整
《科学分析工具MATLAB实验》
实验报告
（18学时系列课程实验）
指导教师：
学生：
班级:
学号:
哈尔滨工业大学机电工程学院机械制造及自动化系
2012年4月
实验一最小二乘法数据处理实验
一、实验目的
通过Matlab上机编程，掌握利用Matlab软件进行数据拟合分析及数据可视化方法
二、实验内容
1.有一组风机叶片的耐磨实验数据，如下表所示，其中X为使用时间，单位为小时h，Y为磨失质量，单位为克g。
要求：
（1）Forpersonaluseonlyinstudyandresearch;notforcommercialuse
（2）
（3）编写一个M文件，对该数据进行合理的最小二乘法数据拟合；
（4）绘出原始数据点和拟合的数据曲线。
表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.对ln（1+x）在[0，1]内的采样数据进行多项式拟合；
3.对指数函数
在[0，1]内的采样数据进行多项式拟合。
三、实验报告内容
1）编写出M文件，对该数据进行合理的最小二乘法数据拟合；
2）绘出原始数据点和拟合的数据曲线。
1）1．clc
clear
xlabel（'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
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（'pleaseinputa='）;
b=input（'pleaseinputb='）;
c=input（'pleaseinputc='）;
q=b^2-4*a*c;
p=sqrt（b^2-4*a*c）;
ifq>0
disp（'两个实根'）
x1=（-b-p）/（2*a）
x2=（-b+p）/（2*a）
elseifq==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
其中x甲1代表在甲机床上加工1号工件，其余类推。
2）
0.5*x甲1+1.1*x乙1+1.0*x丙1<=720
0.5*x乙1+1.2*x乙2+1.3*x乙3<=850
0.6*x丙1+1.3x丙2+0.8*x丙3<=650
x甲1+x乙1+x丙1=340
x乙1+x乙2+x乙3=550
x丙1+x丙2+x丙3=420
3）
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.0000
550.0000
0.0000
54.1355
0.0000
0.0000
55.8645
0.0000
420.0000
fval=
1.2420e+004
所以1号工件在甲机床上加工230件，在乙机床加工54件，丙机床加工56件。
2号工件在甲机床上加工550件。
3号工件在丙机床上加工420件。
总费用为12420。
实验四运动学分析实验
一、实验目的
通过Matlab上机编程，掌握Matlab语言变量、函数、Matlab常用程序结构和流程控制语句、交互控制语句、程序调试运行等。
二、实验内容
凸轮是把一种运动转化为另一种运动的装置，通常凸轮做旋转运动，并转化为从动件的振动、直线运动或者二者的结合。
对基圆为
的凸轮，
为最小圆周半径，与凸轮表面相切，并和凸轮的转轴同心。
则从动件的运动方程为：
对于旋转运动，假定有：
其中，
为从动件的最大位移，
为凸轮转角，
为从动件摆角，并且
。
如果凸轮的角速度
为常量，试求出当
=60°时从动件的位移
、速度
、加速度
、瞬时加速度
的变化曲线。
三、实验报告内容
1）建立M文件，求解从动件的位移、速度、加速度及瞬时加速度；
2）绘制
=60°时从动件的位移
、速度
、加速度
、瞬时加速度
的变化曲线。
clc;
clear;
b=pi/3;
x=linspace（0,2*pi）
fork=1:
1:
100
ifx（k）>=0&x（k）<=b
s（k）=x（k）/b-sin（2*pi*x（k）/b）/（2*pi）
v（k）=（1-cos（2*pi*x（k））/b）/b
a（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/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/b^3*cos（2*pi*（x（k）-b）/b）
升级会员 