matlab第6章实例解析.pptVIP

实例解析【例6-1】导线中的电流与时间的函数关系测量如下表所示，已知测量值的精度很高。t00.1250.250.3750.5i06.247.754.850试计算在时间t=0.01k(k=0,1,2,…,50)时的电流值。解：分别利用lagrange插值法，Newton插值法和hermite插值法求解上述问题，编写程序interp1_test.m。运行结果如右图所示。插值多项式的唯一性！【例6-2】利用多种插值方法对

进行插值。解：调用interp1函数求解，编写如下语句：x=-1:0.01:1;xdata=-1:0.2:1;ydata=1./(1+25*xdata.^2);method={nearest,linear,spline,pchip,cubic,v5cubic};%插值方法fork=1:6subplot(3,2,k)%图形分割yi=interp1(xdata,ydata,x,char(method(k)));%插值plot(x,1./(1+25*x.^2),k,x,yi,k:,xdata,ydata,ko)%绘制图形axistight%控制坐标轴title([\fontname{times}\fontsize{12}\it,char(method(k))])%添加标题endplot(x,yi-1./(1+25*x.^2),k)【例6-3】网格节点插值示例。这里利用不同的插值方法对MATLAB自带的演示函数peaks进行二维插值。figure(Position,[100100560630])%设置图形窗口的位置[X,Y]=meshgrid(-3:.5:3);%产生坐标数据矩阵Z=peaks(X,Y);%计算MATLAB自带峰值函数的值[XI,YI]=meshgrid(-3:0.2:3);%加密网格数据method=char(nearest,linear,spline,cubic);%插值方法s(1)=subplot(321),mesh(X,Y,Z),title(原网格图)%绘制原图并添加标题s(2)=subplot(322),mesh(XI,YI,peaks(XI,YI)),title(加密网格图)%绘制加密后的图形fork=3:6s(k)=subplot(3,2,k)%图形分割ZI=interp2(X,Y,Z,XI,YI,method(k-2,:));%二维插值mesh(XI,YI,ZI),title([method(k-2,:),型插值曲面])e{k-2}=ZI-peaks(XI,YI);%误差endaxis(s,[-33-33-1010])%设置坐标轴范围%绘制插值误差曲面图figurefork=1:4h(k)=subplot(2,2,k)%图形分割mesh(XI,YI,e{k})%绘制误差曲面title([method(k,:),型插值误差曲面])%添加标题xlim([-3,3]);ylim([-3,3]);%设置坐标轴范围end【例6-4】散乱节点插值示例。解：以函数为例进行说明，rand(state,0);x=rand(100,1)*4-2;%随机点的x坐标y=rand(100,1)*4-2;%随机点的y坐标z=@(x,y)x.*exp(-x.^2-y.^2);%定义函数ti=-2:.2:2;[XI,YI]=meshgrid(ti,ti);%生成网格数据ZI=griddata(x,y,z(x,y),XI,YI);%随机点插值h(1)=subplot(2,1,1);%图形分割mesh(XI,YI,ZI);title(散乱数据点插值曲面)%绘制图形并添加标题h(2)=subplot(2,1,2);%图形分割mesh(XI,YI,ZI-z(XI,YI));title(散乱数据点插值误差曲面)%绘制误差曲面并添加标题axis(h,tight)实验范例：国土面积的计算解：首先我们利用MATLAB对上下边界分别进行三次样条和分段三次hermite插值得到加密曲线。此时，设上边界函数为f2(x)，下边界函数为f1(x)，则由定积分定义可知曲线所围区域面积为x=[7.010.513.017.53