qiuhe.m

来自「这是利用matlab编写的一个拟合曲面的工具箱!,可以很方便的拟合空间曲面了.」· M 代码 · 共 22 行

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