数学实验(liuliang).pptVIP

实验8 河流流量估计;实验问题 一条100米宽的河道如下图所示 ;clf;clear x=0:5:100; y=[0 2.41 2.96 2.15 2.65 3.12 4.23 5.12 6.21 5.68 4.22 ... 3.91 3.26 2.85 2.35 3.02 3.63 4.12 3.46 2.08 0]; y1=10-y; plot(x,y1,k.-,markersize,15); axis([0 100 0 10]) grid;最佳曲线拟合; 数据插值;多项式插值 事实上，对于已知的n个数据点，总可以唯一地确定一条n-1次多项式曲线 . 因为;function p=lagrange(x,y) L=length(x); A=ones(L); for i=2:L A(:,i)=A(:,i-1).*x; end q=inv(A)*y; for i=1:L p(i)=q(L-i+1); end;例1：已知观测数据 x 1 2 3 4 5 y -1 1.5 2.1 3.6 4.9 求其多项式插值曲线。;例2：已知观测数据 x 0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1 y -.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.3 11.2 求其多项式插值曲线。;河道问题;; 显然，当数据点较多时，这种高次插值多项式曲线在数据点之间产生大的纹波。这时曲线并不能很好地反映数据点的变 化规律，这就是插值问题中典型的“龙格(Runge)现象”。 实验： 对函数y=1/(1+20x^2) 在[-5,5]上以1为步长 划分作插值，观察插 值曲线的变化。; 针对这种情况，在实际插值计算时，常常采用分段插值的方 法，即在相邻两个数据点构成的子区间上分别进行低次多项式插 值，整个区间上的插值函数将是一个分段的多项式函数。 MATLAB提供了一维插值命令： y0=interp1(x,y,x0,method) 其中x,y为已知数据点坐标向量，y0为对应所求点x0处的插值 结果，method插值类型，省略时为分段线性插值。 分段线性插值为每个子区间上低次多项式为一次多项式时， 几何上为相邻两个数据点间用直线连接的情形。 ;河道问题;河道问题;S=trapz(x,y); dx=diff(x);dy=diff(y); L=sqrt(dx.^2+dy.^2); L=sum(L); fprintf(S=%.2f , L=%.2f\n,S,L); 若不采用直线连接数据点，我们可采用某些光滑的、变化平缓 的曲线来拟合数据点。 最常用的方法是用一个3次多项式，来对相继数据点之间的各 段建模，使其满足相邻两个3次多项式在???点处1阶、2阶导数都 相等，这样可以确定内部各段上的3次多项式，并且多项式通过 节点的斜率和曲率是连续的。而第一个和最后一个多项式必须附 加其他约束条件使其确定。这种类型的插值被称为3次样条插值 或样条插值。 ;例3：已知观测数据 x 0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1 y -.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.3 11.2 求其多项式插值曲线。;clear x=0:.1:1; y=[-.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.3 11.2 ]; plot(x,y,k.,markersize,25); axis([-0.1 1.1 -.6 16]); grid;pause ; hold on t=0:.01:1; f=spline(x,y,t); plot(t,f,r-,linewidth,2);;;河道问题;;clf;clear x=0:5:100; y=[0 2.41