www.gusucode.com > matlab编写的一个拟合曲面的工具箱源码程序 > 3d_fitting/leftmatrix.m

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

m=length(x);
if (nargin~=4) & (m~=length(y))
    error('error check check!');
end

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