关于二维插值函数interp2的疑问

Pig_U4zfj

如题，请问一下数据如果是散乱的，并且不能保证是递增的，怎么利用interp2进行差值，我理解的是先用meshgrid处理数据，但是数据不能保证x，y都是递增的，比如x=[400 300 200 500 350 150 450 360 200],y=[1000 2000 3000  1000 2000 3000 1000 2000 3000],z=[3 2 3 4 5 6 3 4 1],请问怎么操作才能画出meshgrid进而用二维插值函数interp2

Pig_U4zfj

x=[399.8        354.1        318.5        278.1        236.2        203.6        185.3        157.2        117.2        80.8        409.1        365.7        328.3        284.1        243.7        203.2        164.3        123.9        83.5        39.7        408.3        368.3        328.3        289        244.4        208.8        167.7        132.1        89.5        46.1        425.6        380.3        332.7        290.9        244.4        205.1        160.2        114.5        68.8        30.7        420.7        379.6        334.6        291.6        244.4        202.8        157.5        116        74.1        37.8        404.6        360.5        322.7        283        243.3        205.5        162.1        124.7        86.8        52.4        378        344.7        310.3        264.3        226.1        186.8        154.2        115.3        76.3        34.1        315.6        275.5        242.5        210.3        178.5        145.6        118.6        72.6        52.8        22.4]';
y=[1400        1400        1400        1400        1400        1400        1400        1400        1400        1400        1600        1600        1600        1600        1600        1600        1600        1600        1600        1600        1800        1800        1800        1800        1800        1800        1800        1800        1800        1800        2000        2000        2000        2000        2000        2000        2000        2000        2000        2000        2200        2200        2200        2200        2200        2200        2200        2200        2200        2200        2400        2400        2400        2400        2400        2400        2400        2400        2400        2400        2600        2600        2600        2600        2600        2600        2600        2600        2600        2600        2800        2800        2800        2800        2800        2800        2800        2800        2800        2800]';
z=[222.8        220.4        232.4        228.5        227.8        232.6        248.5        245.9        272.4        329.7        222        221.7        235.4        226.5        230.5        236.8        249.1        276.1        407.9        487        226        225.3        226.4        233.9        242.1        283.3        253.9        271.4        323.5        468.6        206.5        231.1        231.1        233        242        244.9        265        299.8        398        596.8        234.7        259.8        235.5        237.6        242.8        292.3        277.9        308.7        396.2        605.9        174.2        242.2        252.1        287.4        253.6        264.6        290.6        316.8        378        518.8        256.9        253.7        253.5        260        303.8        280.7        300.6        346.6        435.6        812.9        257.9        295.3        282.4        288.7        301.9        329.7        357        475.4        580.3        1080.1]';
我想算的是这个数据，想用x,y,z的数据然后进行二维插值，

Pig_U4zfj

我试过先将z按x一维插值使得z等于特定数值x坐标下的值，但是这样误差增大了，不知道还有其他办法吗

TouAkira

1. N = 80;
   
2. [X,Y] = meshgrid(min(x):range(x)/(N-1):max(x), min(y):range(y)/(N-1):max(y));
   
3. Z = griddata(x,y,z,X,Y,'cubic');
   
4. 
   
5. surf(X,Y,Z);hold on;
   
6. plot3(x,y,z,'o');

复制代码

登录/注册后可看大图

Pig_U4zfj

TouAkira 发表于 2019-6-24 18:12

学习了一下，是interp2只能用在网格插值是吗，griddata用于散乱数据插值，