实验题目用正交多项式做小二乘曲线拟合.docVIP

实验题目：用正交多项式做小二乘曲线拟合 实验题目： 用正交多项式做最小二乘的曲线拟合 学生组号： 6 完成日期： 2011/11/27 1 实验目的 针对给定数据的煤自燃监测数据中煤温与之间的非线性关系，用正交多项式做最小二乘曲线拟合。 2 实验步骤 2.1 算法原理 设给定n+1个数据点：（），k=0,1,···,n，则根据这些节点作一个m次的最小二乘拟合多项式（x）=+= 　 = 1 \* GB3 ① 其中，mn,一般远小于n.。 若要构造一组次数不超过ｍ的在给定点上正交的多项式函数系｛（ｊ＝０，１，．．．，ｍ）｝，则可以首先利用｛（ｊ＝０，１，．．．，ｍ）｝作为基函数作最小二乘曲线的拟合，即　（x）＝　　　　　 = 2 \* GB3 ② 根据 = 2 \* GB3 ②式，其中的系数（ｊ＝０，１，．．．，ｍ）为 ，ｊ＝０，１，．．．，ｍ　　　　　　　　　　　　　 = 3 \* GB3 ③ 将 = 4 \* GB3 ④代入 = 3 \* GB3 ③后展开就成一般的多项式。 构造给定点上的正交多项式（ｊ＝０，１，．．．，ｍ）的递推公式如下： 　　　　　　 = 4 \* GB3 ④　 其中　　＝，ｊ＝０，１，．．．，ｍ－１　　　　　　　　　　 = 5 \* GB3 ⑤ 　　　　＝，ｊ＝１，２，．．．，ｍ－１　　　　　　　　　　　　　　 = 6 \* GB3 ⑥ ，ｊ＝０，１，．．．，ｍ－１　　　　　　　　　　　　　　 = 7 \* GB3 ⑦ 则实际计算过程中，根据 = 5 \* GB3 ⑤式逐步求出个正交多项式，并用公式 = 4 \* GB3 ④计算出，并将每次计算展开后累加到拟合多项式 = 1 \* GB3 ①中。 2.2 算法步骤 用三个向量Ｂ,Ｔ,Ｓ,存放多项式，，的系数。 构造，设＝，根据 = 4 \* GB3 ④式，得＝１。再根据 = 7 \* GB3 ⑦ = 3 \* GB3 ③ = 5 \* GB3 ⑤式，计算： 　　　＝ｎ＋１ 　　　 　　　 最后将项展开后累加到拟合多项式中，则＝ （２）构造，设＝＋，根据递推公式 = 4 \* GB3 ④，则可知，＝－，＝１。由公式 = 7 \* GB3 ⑦ = 3 \* GB3 ③ = 5 \* GB3 ⑤ = 6 \* GB3 ⑥求得， 最后将展开累加到拟合多项式 = 1 \* GB3 ①中，有 (3)对于j=2,3,……,m,逐步递推（x） 根据递推公式 = 4 \* GB3 ④有 （x）=（x-）（x）-（x） =（x-）(+….+)- (+….+) 假设 （x）=++…..+ + 则可以得到计算（k=0,1,…,j）的公式： k=j-2,….,2,1 然后分别根据　 = 7 \* GB3 ⑦式 = 3 \* GB3 ③式 = 5 \* GB3 ⑤式与 = 6 \* GB3 ⑥式计算下列量： / / 再将项展开后累加到拟合多项式 = 1 \* GB3 ①中，即有 k=j-1,….,1,0 最后，，为了便于循环使用向量B,T,与S,应将向量T传递给B，向量S传递送给T,即 ，k=j-1,…,1,0 ，k=j,…,1,0 在实际计算中，为了防止运算溢出，将用—来代替，其中 = 在这种情况下，拟合多项式的形式为 ++…..+ 2.3 程序流程图 3 实验结果分析 温度数据： 51.1 52.8 54.8 57.2 58.3 62.7 65.2 67.7 70.6 73.5 75.7 78.6 80.8 84.8 87.8 89.5 92.1 96.4 103.1 112.5 120.8 134.7 152.4 159.1 氧气数据： 20.13 19.9 19.61 18.98 19.61 19.32 19.73 19.59 19.38 20.18 19.98 1