实验函数式——直接确定性模型.pdfVIP

验 函数式 直接确定性模型 一、 验名称：函数式 直接确定性模型 二、 验目的：掌握一维和二维函数的数据插值 和数据拟合 三、 验内容 1、最小二乘法原理、函数插值和数据拟合的原理。 2、Matlab中数据插值命令是interp1()。 3、Matlab中数据拟合命令是polyfit() 。 4、一元和二元函数的插值和拟合 四、 验步骤 1、最小二乘法、函数插值和数据拟合原理。 (1) 最小二乘法原理(略) (2) 函数插值原理 插值：构造一个相对简单的函数y=g(x)使g 通过全部节点，即使g (x ) = y ，j=0,1, , n j j 用g (x)作为函数f ( x ) 的近似。 (3)、函数(多项式)拟合 多项式拟合: 对给定的数据（x ,y ），j = 0,1, , n，选取 j j 适当阶数的多项式(也可采用其它形式的函 数)，例如二次多项式g(x)=ax^2+bx+c使g(x) 尽可能逼近（拟合）这些数据。 • 例1．世界人口数据拟合问题： 据统计，六十年代世界人口数据如下(单位：亿) 年 1960 1961 1962 1963 1964 1965 1966 1967 1968 人口 29.72 30.61 31.51 32.13 32.34 32.85 33.56 34.20 34.83 根据表中数据，预测公元2000年时的世界人口。 问题分析与数学模型 设人口总数为N(t)，根据人口理论的马尔萨斯模 型，采用指数函数 N(t) = e a + b t 对数据进行拟合。为了计算方便，将上式两边同 取对数，得 ，令 y = ln N 或 N = e y 变换后的拟合函数为 y (t) = a + b t 由人口数据取对数（y = ln N ）计算，得下表 T 1960 1961 1962 1963 1964 1965 1966 1967 1968 y 3.3918 3.4213 3.4503 3.4698 3.4763 3.4920 3.5133 3.5322 3.5505 • 根据表中数据及等式 a + b t k = y k ( k = 1，2， ，9) • 可列出关于两个未知数a 、b的9个方程的超定方程组 （方 程数多于未知数个数的方程组） a + t b = y (j = 1，2， ，9) j j 可用最小二乘法求解 算法与数学模型求解 算法如下： 第一步：输入人口数据，并计算所有人口数据的 对数值； 第二步：建立超定方程组的系数矩阵，并计算对 应的正规方程组的系数矩阵和右端向量； 第三步：求解超定方程组并输出结果：a，b； 第四步：利用数据结果构造指数函数计算2000年 人口近似值N(2000)，结束。 MATLAB程序 t=1960:1968;t0=2000; N= [29.72 30.61 31.51 32.13 32.34 32.85 33.56 34.20 34.83]; y=log(N); A= [ ones(9,1), t ]; d=A\ y ;a=d(1),b=d(2) N0=exp(a+b*t0) x=1960:2001;yy=exp(a+b*x); plot(x,yy,t,N,o,2000,N0,o) 计算结果为 － 33.0383， 0.0186 N(2000) = 63.2336 所以取五位有效数，可得人口数据的指数拟合函数 经计算得2000年人口