matlab用自编函数最小二乘法拟合

在 MATLAB 中，使用自编函数最小二乘法拟合可以分为以下几个步骤：

1. 定义自编函数

function [p,fit]=myfit(x,y,n)
% x,y为样本数据，n为拟合多项式的次数
% 返回拟合参数p和拟合结果fit

% 构造拟合矩阵
X=zeros(length(x),n+1);
for i=0:n
    X(:,i+1)=x.^i;
end

% 计算拟合参数
p=inv(X'*X)*X'*y;

% 计算拟合结果
fit=X*p;
end

2. 调用自编函数进行拟合

% 样本数据
x=[1 2 3 4 5];
y=[1.2 1.9 3.2 4.1 5.2];

% 调用自编函数进行拟合
n=2; % 拟合二次曲线
[p,fit]=myfit(x,y,n);

% 绘制拟合结果
plot(x,y,'o',x,fit)
legend('样本数据','拟合结果')

在上述代码中，我们使用了自编函数 `myfit` 进行拟合。该函数首先构造拟合矩阵，然后使用最小二乘法计算拟合参数，最后计算拟合结果。在调用自编函数时，我们需要提供样本数据和拟合多项式的次数。最后，我们使用 `plot` 函数绘制样本数据和拟合结果。