MATLAB一维插值和二维插值 比较

插值问题描述:已知一个函数上的若干点,但函数具体表达式未知,现在要利用已知的若干点求在其他点处的函数值,这个过程就是插值的过程.

 

1.一维插值

一维插值就是给出y=f(x)上的点(x1,y1),(x2,y2),…,(xn,yn),由此求出y=f(x)在点xa处的值ya的值.

实现一维插值使用interp1命令,使用参数为interp1(x,y,xa,’method’) ,其中x和y是已知点对应横纵坐标,xa为代求值的横坐标,method参数代表插值类型,参数可以选择的选项如下表,若缺省则为linear.

method                含义                     特   点

linear               线性插值             快,精度不高

cubic           三次多项式插值      较慢,精度高,平滑

spline          三次样条插值          最慢,精度高,最平滑

 

例如,求y=xln(x)的插值问题,代码如下

    %为了说明插值问题,根据已知函数构造一组插值点

    %而实际问题中,只有这些点,函数是未知的

    x = 0.2:0.4:2;

    y = x.*log(x);%构造一组插值点

    plot(x,y,’o’);

    xa =0.2:0.1:2;%d代求y值得x值点

    ya = interp1(x,y,xa,’cubic’);

    hold on;

    plot(xa,ya,’b-‘);

    box on;

    hold on;

    fplot(@(x)x*log(x),[0.2,2]);

    legend(‘插值点’,’三次多项式插值’,’原图’);%图例

 

2.二维插值

    二维插值就是给出z=f(x,y)上的点(x1,y1,z1),…,(xn,yn,zn),由此求出在(xa,ya)处求出za的值.实现一维插值使用interp2命令,使用参数为interp1(x,y,z,xa,ya,’method’),参数含义与interp1类似.

    例如,求z = exp(-x^2-y^2)的空间插值问题,代码如下

    x = linspace(-2,2,10);

    y = linspace(-2,2,10);

    [X,Y] = meshgrid(x,y);

    Z = exp(-X.^2 – Y.^2);

    subplot(1,2,1);

    surf(X,Y,Z);

    title(‘插值点曲图’);

    hold on;

    xa = linspace(-2,2,50);

    ya = linspace(-2,2,50);

    [Xa,Ya] = meshgrid(xa,ya);

    Za = interp2(X,Y,Z,Xa,Ya,’cubic’);

    subplot(1,2,2);

    surf(Xa,Ya,Za);

    title(‘三次插值曲面’);

 

转载自 https://blog.csdn.net/ten_sory/article/details/54429482

    原文作者:这是啥?
    原文地址: https://blog.csdn.net/qq_41764205/article/details/82667330
    本文转自网络文章,转载此文章仅为分享知识,如有侵权,请联系博主进行删除。
点赞