数值计算一维.ppt

钮海 位于…\matlab\toolbox\exlink目录下 前言 内容： 一维问题 高维问题 前言 层次： 插值与拟合理论 Matlab实现 应用实例 前言 问题 测试数据 分类 一个自变量－一维问题 多个自变量－多维问题 方法 1.测试数据为精确数据或者误差可以忽略 方法：插值法 2.测试数据含误差，必须考虑误差 方法：拟合法或回归 一维问题 插值 一维问题-插值 一维问题-插值 插值： 一维问题-插值 首要问题： f的选取？ 一维问题-插值 求解： 一维问题-插值 解： 一维问题-插值 模型分析：Runge现象 测试点越多，多项式次数越高 一维问题-插值 失真解决办法：每次采取数据点比较少 如2个；4个 问题：原先大量数据？ 方法：分组（每组2个；4个－递进进行） 一维问题-插值 再处理！（光滑） 一维问题-插值 数学模型： 一维问题-插值 后注：取自工程中利用弹性竹条（样条）进行插值的方法 实现 编程？ Matlab实现：插值函数（程序） 插值=interp1(自变量数据,因变量试数据,插值点,方法) 实验－问题 实验－脚本文件 %各类一维分段插值比较－－脚本文件（文件名cz.m） x=linspace(-5,5,9);%原始自变量 y=1./(1+x.^2);%原始因变量 x_o=linspace(-5,5,50); y_o=1./(1+x_o.^2);%原始函数 x_in_l=linspace(-5,5,50); y_in_l=interp1(x,y,x_in_l);%线性插值函数 subplot(3,1,1)%子图 plot(x_o,y_o,x_in_l,y_in_l) title(‘linear) x_in_c=x_in_l; y_in_c=interp1(x,y,x_in_c,cubic);%三次插值函数 subplot(3,1,2) plot(x_o,y_o,x_in_c,y_in_c) title(‘cubic) x_in_s=x_in_l; y_in_s=interp1(x,y,x_in_s,spline);%样条插值函数 subplot(3,1,3) plot(x_o,y_o,x_in_s,y_in_s) title(‘spline) 实验－运行 在matlab工作区输入cz 插值方法选择 插值函数 插值=interp1(自变量数据,因变量试数据,插值点,方法) 建模实例 估计水塔的水流量(AMCM92A) 美国某州的各用水管理机构要求各社区提供以每小时多少加仑计的用水率以及每天所用的总水量。但许多社区并没有测量流入或流出当地水塔的水量的设备，他们只能代之以每小时测量水塔中的水位，其精度在5%以内。更为重要的是，无论什么时候，只要水塔中的水位下降到某一最低水位L时，水泵就启动向水塔重新充水直至某一最高水位H，但也无法得到水泵的供水量的测量数据。因此，在水泵正在工作时，人们不容易建立水塔中的水位与水泵工作时的用水量之间的关系。水泵每天向水塔充水两次，每次约二小时。 试估计在任何时刻，甚至包括水泵正在工作的时间内，水从水塔流出的流量f (t)，并估计一天的总用水量和水泵的工作功率。表1给出了某个真实小镇某一天的真实数据。 建模实例 建模实例 表1给出了从第一次测量开始的以秒为单位的时刻，以及该时刻的高度单位为百分之一英尺的水塔中水位的测量值，例如，3316秒后，水塔中的水位达到31.10英尺。水塔是一个垂直圆形柱体，高为40英尺，直径为57英尺。通常当水塔的水位降至约27.00英尺时水泵开始向水塔充水，而当水塔的水位生至约35.50英尺时水泵停止工作。 建模实例 t0=[0,3316,6635,10619,13937,17921,21240,25223,28543,32284,39435,43318,46636,49953,... 53936,57254,60574,64554,68535,71854,75021,85968,89953,93270]; L0=[3175,3110,3054,2994,2947,2892,2850,2797,2752,2697,3550,3445,3350,3260,3167,3087,... 3012,2927,2842,2767,2697,3475,3397,3340]; t=linspace(0,86400,200); L=interp1(t0,L0,t,clubic); plot(t0,L0,*,t,L) 建模实例 一维问题 拟合 一维问题-拟合 一维问题-拟合 拟合： 一维问题-拟合 一维问题-拟合 首要问题： f的选取？ 一维问题-拟合 求解： 二次无约束规划问题（多元函数极值） 实现－多项式拟合 编程？ M