三维数据拟合问题

wgwhl

已知测试数据：x，y，z 分别是 5 X 9 矩阵，现要在 x，y 平面内绘制 z 的等值线。我的拟合方法是按照下面的方法进行的，不知道是否正确，请高手指点。

1. 
   
2. for index = 1 : 5
   
3.     xi = x(: , index);
   
4.     yi = y(: , index);
   
5.     zi = y(: , index);
   
6. 
   
7.     % 1、用 x，y 拟合曲线多项式
   
8.     FX = polyfit(xi,  yi,  4);
   
9. 
   
10.     % 2、根据拟合曲线多项式取1000个点的值用于画曲线图
    
11.     Xi = linspace(min(xi),  max(yi),  1000);    % 得到 1000 个 x 值
    
12.     Yi = polyval(FX, Xi);  % 得到 1000 个 y 值
    
13.    
    
14.     % 3、用 Xi，Yi 拟合Zi
    
15.     F = [ones(size(zi)), xi, yi];
    
16.     a = F \ zi;
    
17.     Zi = a(1) + a(2) * Xi + a(3) * Yi;   %  得到 1000 个 z 值
    
18. 
    
19.     % 4、将每一行拟合后的数据加入各自的拟合数据矩阵
    
20.     X = [X;  Xi];   % 5 X 1000 矩阵数据
    
21.     Y = [Y;  Yi];   % 5 X 1000 矩阵数据
    
22.     Z = [Z;  Zi];   % 5 X 1000 矩阵数据
    
23. end
    
24. % 后面的代码就是用  plot 画 X，Y的 平面图，用 contour 画 Z 在X，y平面的等值线图
    

复制代码

因为数学基础不好，不知道这样拟合出来的数据是否合理？另外：如何做，使结果是 1000 X 1000 的矩阵数据呢？请指点。
谢谢！

另外：如果 X，Y平面是对数坐标，用 contour 画出来的等值线是正确的吗？

[ 本帖最后由 wgwhl 于 2010-8-6 16:14 编辑 ]

ChaChing

没LZ的原始数据, 不好试LZ的程序!:loveliness:
Xi = linspace(min(xi),  max(yi),  1000);是否有误?
还有试试griddata可否满足LZ需求?

wgwhl

是的，这个地方敲错了，因为源程序变量不是这个，在发帖的时候改成x,y,z了，
Xi = linspace(min(xi),  max(yi),  1000);  应为  Xi = linspace(min(xi),  max(xi),  1000);

我已经想到了其他办法，先试试再来这里问吧