实验六、拟合与插值问题.pptVIP

人口预测线性模型 对于开始提出的实验问题, 代如数据，计算得 syms xx=1949:5:1994;y=[541.67 602.66 672.09 704.99 806.71 908.59 975.42 1034.75 1106.76 1176.74];plot(x,y,r*’)z=polyfit(x,y,1)y1=polyval(z,x)hold onplot(x,y1,b+-) syms x x=0:1:14 y=[39 53 72 96 129 171 232 314 386 502 629 760 920 1065 1232]; plot(x,y,r*) z=polyfit(x,y,2) y1=polyval(z,x) hold on plot(x,y1,b+-) 插值的基本原理 选讲 对表格给出的函数，求出没有给出的函数值 其中x0,y0是已知的节点坐标，是同维向量。y对应于x处的插值。y与x是同维向量。 解决上述问题,我们可分两步: 一 用原始数据绘图作为选用插值方法的参考. 二 确定插值方法进行插值计算 对y=1/(1+x2),-5≤x≤5，用n（=11）个节点（等分）作上述两种插值，用m（=21）个插值点（等分）作图，比较结果。 解：键入并运行如下命令 n=11;m=21;x=-5:10/(m-1):5;y=1./(1+x.^2); xo=-5:10/(n-1):5;yo=1./(1+xo.^2); y1=interp1(xo,yo,x); y2=spline(xo,yo,x); plot(x,y,r,x,y1,b,x,y2,k) 例 3 一水库上游河段降暴雨.，根据预报测算上游流入水库的流量为Q(t) (102立方米/秒) ： t (时) 8 12 16 24 30 44 48 56 60 Q( t ) 36 54 78 92 101 35 25 16 13 利用这个预报值估计14：30 和 20：30 时上游流入水库的流量。 假设: 1 已知数据准确。 2 相邻两个时刻之间的流量没有突然的变化。 t=[8,12,16,24,30,44,48,56,60]; q=[36,54,78,92,101,35,25,16,13]; t1=8:0.5:60; q1=interp1(t,q,t1,linear); plot(t,q,b,t1,q1); hold on; q2=interp1(t,q,t1,spline); plot(t,q,b,t1,q1,t1,q2,r) q1 q2 2、 最佳营销策略 问题：某公司有一批以每桶2元购进的彩漆为了获得较高的利润希望以较高的价格卖出但价格越高,售出量就越少,二者之间的关系由表一给出.于是打算用作广告的办法来促销.而广告费与销售量的关系可由销售增长因子来描述.例如,投入3万元的广告费,销售因子为1.85,意味着做广告后的销售量将是未做广告销售量的1.85倍.根据经验,广告费与销售因子的关系如表2,现请你作出决策：投入多少广告费和售价为多少时所获得的利润最大? 提示 用描点法画出预期销售量——售价；销售增长因子——广告费的关系图. 建立上述两个函数关系. 设x为售价,y为广告费,P为所得利润,建立P关于x和y的函数关系. 用多元函数的极值理论计算并回答此问题. 3、在某处测得海洋不同深度处水温如下： 求深度为500、1000、1500米处的水温。 1. 拉格朗日插值:自编程序,如名为 lagr1.m 的M文件， 第一行为 function y=lagr1(x0,y0,x) 输入:节点x0,y0, 插值点x (均为数组，长度自定义)）； 输出:插值y (与x同长度数组)）。 应用时输入x0,y0,x后,运行 y=lagr1(x0,y0,x) 2. 分段线性插值:已有程序 y=interp1(x0,y0,x) 3. 三次样条插值:已有程序 y=interp1(x0,y0,x, ’spline’) 或 y=spline(x0,y0,x) 用MATLAB作插值计算 例2 表中是待加工零件下轮廓线的一组数据，现需要得到x坐标每改变0.1时所对应的y的坐标. 1.0 14 1.6 1.2 1.8 2.0 2.1 1.7 1.2 0 y 15 13 12 11 7 5 3 0 x 下面是关于插值的两条命令(专门用来解决这类问题)： y=interp1(x0,y0,x)