MATLAB编程基础第7讲插值拟合与初值常微分方程的求解.ppt

* * （2）load filename variables 将以前用save命令保存的变量variables从磁盘文件中调入MATLAB工作空间。 用load 命令调入的变量，其名称为用save命令保存时的名称，取值也一样。 Variables所表示的变量列表中，不能用逗号，各个不同的变量之间只能用空格来分隔。 未列出variables时，表示将磁盘文件中的所有变量都调入工作空间。 3、清屏和退出 clc quit 或 exit * * ?load —— ?load data —— ?load data a b —— mat文件是标准的二进制文件，还可以ASCII码形式保存。 即可恢复保存过的所有变量 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 * * * MATLAB编程基础之 插值、拟合与初值常微分方程的求解 第七讲 梁丙臣 * * 3.8.3多项式的拟合 需求：实验数据总结为规律曲线等等 最小二乘法：所用曲线限定为多项式，在数据点上拟合值和函数值有最小误差平方和。 P=polyfit(x,y,n)：把自变量x和函数值y拟合成n阶多项式，向量为p。 * * 例3-35 x = (0:0.1:1); y = [-0.4 1.9 3.2 6.2 7.1 7.3 7.7 9.6 9.5 9.3 12]; % 给出一组11个点数据 y2 = polyfit(x,y,2) % 计算2阶拟合的多项式向量 x1= 0:0.01:1; f2 = polyval(y2,x1); % 2阶拟合曲线在各点的函数值 y10 = polyfit(x,y,10) % 计算10阶拟合的多项式向量 f10 = polyval(y10,x1); % 10阶拟合曲线在各点的函数值 plot(x,y,o,x1,f2,:,x1,f10,k) * * 例3-36 x = (0: 0.1: 2.5); % 给出一组数据，为误差函数的一个区间 y = erf(x); p = polyfit(x,y,6) % 计算该区间内6阶拟合多项式的向量 x1 = (0: 0.1: 5); % 将区间增长一倍 y1 = erf(x1); % 计算误差函数在新区间内的函数值 f = polyval(p,x1); % 计算6阶拟合曲线在新区间内的取值 plot(x1,y1,o,x1,f,-) % 绘图，比较前后区间内的曲线拟合效果，如图3-7 axis([0 5 0 2]) * * 例：拟合以下数据 x 0.5 1.0 1.5 2.0 2.5 3.0 y 1.75 2.45 3.81 4.80 8.00 8.60 x=[0.5 1.0 1.5 2.0 2.5 3.0]; y=[1.75 2.45 3.81 4.80 8.00 8.60]; a=polyfit(x,y,2) a = 0.4900 1.2501 0.8560 x1=[0.5:0.05:3.0]; y1=a(3)+a(2)*x1+a(1)*x1.^2; plot(x,y,*) hold on plot(x1,y1,-r) * * 例：根据经验公式y=a+bx2，拟合如下数据： xi 19 25 31 38 44 yi 19.0 32.3 49.0 73.3 98.8 x=[19 25 31 38 44]; y=[19.0 32.3 49.0 73.3 98.8]; x1=x.^2; x1=[ones(5,1),x1]; ab=x1\y ab = 0.5937 0.0506 x0=[19:0.2:44]; y0=ab(1)+ab(2)*x0.^2; plot(x,y,o,x0,y0,-r) * * 3.8.4多项式的插值 插值的定义——是对某些集合给定的数据点之间函数的估值方法。 当不能很快地求出所需中间点的函数时，插值是一个非常有价值的工具。 Matlab提供了一维、二维、 三次样条等许多插值选择 * * 1.一维插值 yi=interp1(x,Y,xi,method)：求已知同维数据x和Y，运用method指定的方法计算插值点xi处的数值yi。 method有如下四种： nearest 最近点差值，取与最近的已知数据点的值 linear 线性差值，直线连接数据点，插值点值位于 直线上 spline 样条插