数学建模社区-数学中国

标题: 一个关于插值的matlab程序问题 [打印本页]

作者: sdccumcm 时间: 2012-2-3 16:59
标题: 一个关于插值的matlab程序问题
例山区地貌：
在某山区测得一些地点的高程如下表。平面区域为
1200<=x<=4000,1200<=y<=3600)
试作出该山区的地貌图和等高线图，并对几种插值方法进行比较。

X
Y 1200 1600 2000 2400 2800 3200 3600 4000
1200 1130 1250 1280 1230 1040 900 500 700
1600 1320 1450 1420 1400 1300 700 900 850
2000 1390 1500 1500 1400 900 1100 1060 950
2400 1500 1200 1100 1350 1450 1200 1150 1010
2800 1500 1200 1100 1550 1600 1550 1380 1070
3200 1500 1550 1600 1550 1600 1600 1600 1550
3600 1480 1500 1550 1510 1430 1300 1200 980

下面是照着书上打下来的代码，运行之后报错，想请高手们帮忙看看：
x=0:400:5600;
y=0:400:4800;
z=[1130 1250 1280 1230 1040 900 500 700
1320 1450 1420 1400 1300 700 900 850
1390 1500 1500 1400 900 1100 1060 950
1500 1200 1100 1350 1450 1200 1150 1010
1500 1200 1100 1550 1600 1550 1380 1070
1500 1550 1600 1550 1600 1600 1600 1550
1480 1500 1550 1510 1430 1300 1200 980];
%用原始数据作山区表面图
figure(1);
meshz(x,y,z)
xlabel('X'),ylabel('Y'),zlabel('Z')
%每隔50加密网络，分别作最近点插值、双线性插值和双三次插值，做出插值后的表面图
xi=0:50:5600;
yi=0:50:4800;

figure(2)
z1i=interp2(x,y,z,xi,yi,'nearest');%最近点插值
surfc(xi,yi,z1i)
xlabel('X'),ylabel('Y'),zlabel('Z')

figure(3)
z2i=interp2(x,y,z,xi,yi,'linear');%线性插值
surfc(xi,yi,z2i)
xlabel('xx'),ylabel('yy'),z;ane;('zz')

figure(4)
z3i=interp2(x,y,z,xi,yi,'cubic')%三次插值
surfc(xi,yi,z3i)
xlabel('**'),ylabel('yyy'),zlabel('zzz')

figure(5)
suplot(1,3,1),contour(xi,yi,z1i,10,'r');
suplot(1,3,2),contour(xi,yi,z2i,10,'r');
suplot(1,3,3),contour(xi,yi,z3i,10,'r');

还有一个问题就是，我在运行这样几个语句时，系统报错了，分析过后可能是我自己电脑的问题，但是又查不出来：
>> x=1:5;
>> yy=1:3;
>> temps=[82 81 80 82 84;79 63 61 65 81;84 84 82 85 86];
>> mesh(x,y,temps)
??? Undefined function or method 'usev6plotapi' for input arguments of type 'double'.

Error in ==> mesh at 53
[v6,args] = usev6plotapi(varargin{:},'-mfilename',mfilename);

谢谢！

作者: 李——建辉 时间: 2012-2-3 19:43
特意来佩服一下···401842064618743749216563207627303297359495429075897925516987348927616539792984

作者: sdccumcm 时间: 2012-2-19 12:15

作者: zmssorry 时间: 2012-2-21 10:09
数据格式不匹配

作者: sdccumcm 时间: 2012-2-21 10:56
喔~~好的谢谢

作者: wyanxia 时间: 2012-2-21 11:08
路过，插值方法都快忘记了

作者: 新生泪 时间: 2012-2-21 15:42
usev6plotapi没有用过啊

作者: 紫藤恋晴天 时间: 2012-8-16 13:32
你仔细看了你的程序没有?细看一下你就会发现你的z的取值不够（你看看你的原题，我没有所以具体不知道），还有就是你的z2i=interp2(x,y,z,xi,yi,'linear）中的yi应该是其转自（字打错了）也就是yi'的'

作者: sdccumcm 时间: 2012-8-16 14:44

紫藤恋晴天发表于 2012-8-16 13:32
你仔细看了你的程序没有?细看一下你就会发现你的z的取值不够（你看看你的原题，我没有所以具体不知道），还 ...

明白了多谢了！

作者: 生于青春 时间: 2012-8-16 15:48

作者: sdccumcm 时间: 2012-8-16 18:58

生于青春发表于 2012-8-16 15:48

~~~~~~~~~~

作者: sdccumcm 时间: 2012-11-19 16:02
今年初发的帖子了…… 回顶个纪念一下下

作者: dabiu 时间: 2012-11-21 20:24
正确的程序如下：
clc
clear
x=0:800:5600; %15
y=0:800:4800;%13
z=[1130       1250       1280       1230       1040       900       500       700
1320       1450       1420       1400       1300       700       900       850
1390       1500       1500       1400       900       1100       1060       950
1500       1200       1100       1350       1450       1200       1150       1010
1500       1200       1100       1550       1600       1550       1380       1070
1500       1550       1600       1550       1600       1600       1600       1550
1480       1500       1550       1510       1430       1300       1200       980];
%用原始数据作山区表面图
[x,y]=meshgrid(x,y); %13*15
figure(1);
mesh(x,y,z)
xlabel('X'),ylabel('Y'),zlabel('Z')
%每隔50加密网络，分别作最近点插值、双线性插值和双三次插值，做出插值后的表面图
xi=0:50:5600;
yi=0:50:4800;
yi=yi';
z1i=griddata(x,y,z,xi,yi,'nearest');%最近点插值
figure(2)
% z1i=interp2(x,y,z,xi,yi,'v4')
surfc(xi,yi,z1i)
xlabel('X'),ylabel('Y'),zlabel('Z')

figure(3)
z2i=interp2(x,y,z,xi,yi,'linear');%线性插值
surfc(xi,yi,z2i)
xlabel('xx'),ylabel('yy'),zlabel('zz')

figure(4)
z3i=interp2(x,y,z,xi,yi,'cubic')%三次插值
surfc(xi,yi,z3i)
xlabel('**'),ylabel('yyy'),zlabel('zzz')

figure(5)
subplot(1,3,1),contour(xi,yi,z1i,10,'r');
subplot(1,3,2),contour(xi,yi,z2i,10,'r');
subplot(1,3,3),contour(xi,yi,z3i,10,'r');

作者: sdccumcm 时间: 2012-11-21 20:52

dabiu 发表于 2012-11-21 20:24
正确的程序如下：
clc
clear

多谢回复啦！！！我好好看下哈！

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)