文档介绍：计算方法大作业
第一部分、计算方法学****回顾
:
误差分析,插值法与拟合,数值积分,数值微分,线性方程组的直接
解法和迭代解法,非线性方程求根,矩阵特征值问题计算、常微分方程初值问题数值解法
:
(1)最小二乘拟合方法的实际应用
曲线拟合就是拟合测量数据曲线。所选择的曲线有时通过数据点,但在其他点上,曲线接近它们而不必通过它们13,41~在大多数情况下,选择曲线使得数据点的平方误差和最小。这种选择就是最小二乘曲线拟合。下面介绍一下最小二乘法拟合的基本原理。设已知个数据点)(i=0,1,„, 一1),求(m一1) 次最
小二乘拟合多项式:其中
设拟合多项式为各正交多项式:
的线性组合:
则继续往向下推导得:
继续推导最后可得一般形式的m一1次多项式:
即为最小二乘拟合多项式其拟合精度由下式来评定:
应用实例: 某建筑物176 d水平位移测量数据如下表所示,在程序编制过程中,为了防止运算溢出,xi用来代
替,其中x=(i=0n-1)xi/n。
此时,拟合多项式的形式为:
此时,拟合多项式的形式为:
运用最小二乘多项式拟合时,拟合多项式的次数越高,其拟合精度未必越高。以拟合最高次数l9次为例,拟合系数如表2,拟合的精度评定见表3。
根据水平位移的观测数据,实现了累计观测时间与水平位移的曲线拟合,在有限的测量数据条件下,表述了时间与该建筑物水平位移之间的函数关系。曲线拟合的最小二乘法在解决这类问题的数据处理和误差分析中应用非常广泛,提高了数据处理的效率和精确度,最d"-乘曲线拟合实现方法简明、适用,可应用于类似的测量数据处理和实验研究。
举例:在万能拉拨机中有一个园柱形凸轮,其底园半径R=300mm,凸轮的上端面不在同一平面上,而要根据动杆位移变化的需要进行设计制造。按设计要求,将底园周18等分,旋转一周。第个分点对应柱高,数据见下表。为了数控加工,需要计算出园周上任一点的柱高。凸轮高度的数据(单位:mm)
分点
0和18
1
2
3
4
5
柱高
分点
6
7
8
9
10
11
柱高
分点
12
13
14
15
16
17
柱高
我们将园周展开,借助MATLAB软件画出对应的柱高曲线散点图(左下图)。
clear;close;
x=linspace(0,2*pi*300,19);
y=[ ,,,,,,,,,,,,,,,,,,];
plot(x,y,’o’);axis([0,2000,0,550]);
可见,可以用三次多项式插值,下面给出借助MATLAB软件画出的柱高插值曲线图(右上图)。
xi=0:2*pi*300;
yi=interp1(x,y,xi,’cubic’);
plot(xi,yi);
(2)数据插值方法及应用
多项式插值
设有次多项式
通过所有个点,那么就有
可以证明当且时,这样的多项式存在且唯一。若要求得到函数表达式,可直接解上面方程组。若只要求得函数在插值点处数值,可用下列Lagrange插值公式
多项式插值光滑但不具有收敛性,一般不宜采用高次多项式(如)插值。
第二部分、水槽的液位控制
1问题的提出
,对耦合水槽建模如下[
]:
(3-39)
其中,为1号和2号水槽的截面积,,分别是两个水槽的液面高度,,分别是两个水泵向水槽中注水的体积流速(),,,分别是与,,对应的比例常数,通过升高或者降低挡板可以控制两个水槽之间的流速,关闭2号水泵,通过1号水泵向1号水槽中注水以控制2号水槽液面高度
2用 MATLAB语言编程的方法实现水槽液位模糊控制系统的设计与仿真
(1) 建立模糊推理系统(FIS)
要想设计整个模糊控制系统,首先需要建立模糊推理系统,这也是建模的主要工作。采用MATLAB语言编程的方法
建立模糊推理系统的过程如下所述。
1) 建立新的 FIS结构
设所建的 FIS结构为 FC. fis ,则命令行为:FC = newfis(’ FC’) ;这时便建立了一个新的名叫 FC的 mamdani 型系统。此时它的输入、输出及规则都还没有定义。
2) 添加输入输出变量
由于本