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

    function he=qiuhe(x,p,y,q,z)
% he    x^p*y^q 从1－>m的和
% x,y 行向量 长度相同
% p,q 为x,y系数
%x=[2 3]; y=[ 3 4];
%p=2; q= 2;



m=length(x);
if (nargin<4 )&(m~=length(y))            %输入量至少为四，x,y行向量长度必需一样
    error('error check check!');
end
if nargin==4                            %没有 z , 默认为单位行向量
    z=ones(m,1);
end

he=0;
for i=1:m
    he=he+x(i)^p * y(i)^q*z(i);     % 1-->m 求和
end