曲面拟合算法程序

该程序利用matlab开发环境编制而得

function A=leftmatrix（xpyq）
% A*a=B  a 即为系数列矩阵  
% A为左边（p-1）（q-1） 乘 （p-1）（q-1） 的矩阵
% xy 为长度一样的列矩阵 也就是给定离散点的xy坐标
% pq为拟合的函数中xy的指数

m=length（x）;
if （nargin~=4） & （m~=length（y））
    error（‘error check check!‘）;
end

A_length=p*q;                        % A 为p*q阶的方阵
A=zeros（A_lengthA_length）;          % 赋值0
for i=1 : p*q
    for j= 1 : p*q
        x_z=quotient（j-1q）+quotient（i-1q）;     % x 的指数   quotient为求商
        y_z=mod（j-1q）+mod（i-1q）;            % y 的指数
        A（ij）=qiuhe（xx_zyy_z）;            
    end
end
    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        614  2004-04-20 13:39  Matlab  拟合\3d_fitting\leftmatrix.m

     文件        830  2004-04-20 16:44  Matlab  拟合\3d_fitting\main.m

     文件        439  2004-04-20 13:38  Matlab  拟合\3d_fitting\qiuhe.m

     文件         65  2004-04-20 13:39  Matlab  拟合\3d_fitting\quotient.m

     文件        296  2004-04-20 13:39  Matlab  拟合\3d_fitting\rightmatrix.m

    ..A..H.       162  2008-07-24 14:56  Matlab  拟合\3d_fitting\~$readme.doc

    ..A..H.     59392  2004-04-20 22:01  Matlab  拟合\3d_fitting\~WRL0977.tmp

    ..A..H.     59392  2008-07-24 15:34  Matlab  拟合\3d_fitting\~WRL2290.tmp

     文件      59392  2004-04-20 22:01  Matlab  拟合\3d_fitting\原理.doc

     目录          0  2008-07-24 15:47  Matlab  拟合\3d_fitting

     目录          0  2008-07-24 15:47  Matlab  拟合

----------- ---------  ---------- -----  ----

               180800                    12