文档介绍：仅供个人参考
不得用于商业用途
仅供个人参考
不得用于商业用途
《科学分析工具MATLAB?验》
实验报告
（18学时系列课程实验）
指导教师：
学 生：
班 级：
学 号：
哈尔滨工业大学机电工程学院 机械制造及自动化系
2012年4月
实验一最小二乘法数据处理实验
一、实验目的
仅供个人参考
不得用于商业用途
仅供个人参考
不得用于商业用途
通过Matlab上机编程，掌握利用Matlab软件进行数据拟合分析及数据可视化方法
二、实验内容
.有一组风机叶片的耐磨实验数据，如下表所示，其中 X为使用时间，单位为小时h, Y 为磨失质量，单位为克go要求：
For personal use only in study and research; not for commercial use
⑵
(3)编写一个M文件，对该数据进行合理的最小二乘法数据拟合；
(4)绘出原始数据点和拟合的数据曲线。
表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=poly攸x,y,3) y1=polyval(p,x);
plot(x,y,'+',x,y1,'-r')
仅供个人参考
不得用于商业用途
仅供个人参考
不得用于商业用途
3
clc
clear
x=::1;
y=1-sqrt(x);
p=poly攸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命令，编程完成一元二次
2
方程ax +bx+c=0的求解，给出设计的程序源代码，并针对 a、b、c的不同取值情况，
给出程序的运行结果。
三、实验报告内容
要求编写函数型M文件，输入三个系数，运用判别式法判定根的情形，并输出各种类 型的根。
仅供个人参考
不得用于商业用途
仅供个人参考
不得用于商业用途
仅供个人参考
不得用于商业用途
clc
clear
a=input('please input a=');
b=input('please input b=');
c=input('please input c=');
q=bA2-4*a*c;
p=sqrt(bA2-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