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