第6章实例解析.pptVIP

【例6-1】导线中的电流与时间的函数关系测量如下表所示，已知测量值的精度很高。 【例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}; % 插值方法 for k=1:6 subplot(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) % 绘制图形 axis tight % 控制坐标轴 title([\fontname{times}\fontsize{12}\it,char(method(k))]) % 添加标题 end 【例6-3】网格节点插值示例。 这里利用不同的插值方法对MATLAB自带的演示函数peaks进行二维插值。 figure(Position,[100 100 560 630]) % 设置图形窗口的位置 [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(加密网格图) % 绘制加密后的图形 for k=3:6 s(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); % 误差 end axis(s,[-3 3 -3 3 -10 10]) % 设置坐标轴范围 % 绘制插值误差曲面图 figure for k=1:4 h(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.0 10.5 13.0 17.5 34 40.5 44.5 48 56 61 68.5 76.5 80.5... 91 96 101 104 106 111.5 118 123.5 136.5 142 146 150 157 158]; % x轴坐标 y1=[44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 ... 32 65 55 54 52 50 66