适用于解决的函数问题为:
代码全文见下载
coef= cell(1,5);
for i = 1:5
if ~poly_type
coef{i} = b(i,1) * x(:,1) .^ 2 + b(i,2) * x(:,2) .^ 2 + b(i,3) * x(:,1) .* x(:,2) + ...
b(i,4) * x(:,1) + b(i,5) * x(:,2) + b(i,6); % Quadratic form
else
coef{i} = b(i,1) * x(:,1) .^ 3 + b(i,2) * x(:,1) .^ 2 .* x(:,2) + b(i,3) * x(:,1) .* x(:,2) .^ 2 + ...
b(i,4) * x(:,2) .^ 3 + b(i,5) * x(:,1) .^ 2 + b(i,6) * x(:,1) .* x(:,2) + b(i,7) * x(:,2) .^ 2 + ...
b(i,8) * x(:,1) + b(i,9) * x(:,2) + b(i,10); % Qubic form
end
end