rightmatrix.m

曲线拟合工具

function B=rightmatrix(x,p,y,q,z)
% A*a=B 
% B  为一个列向量 长为p*q
% x y z 为点的坐标 
%p q 为x y指数

if nargin~=5 
    error('error check check! rightmatrix')
end

B=zeros(p*q,1);
for i=1 : p*q
    x_z=quotient(i-1,q);
    y_z=mod(i-1,q);
    B(i,1)=qiuhe(x,x_z,y,y_z,z);
end