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