nihe.m

来自「matlab的60个小程序」· M 代码 · 共 21 行

x=1:8;
y=[15.3,20.5,27.4,36.6,49.1,65.6,87.8,117.6];
y1=log(y);
n=1;%表示用一次函数拟和
m=length(x);
A=zeros(m,n+1);
for i=1:2
    A(:,i)=x.^(i-1)';
end
B=A'*A;
C=A'*y1';
re=B\C;
a=exp(re(1));
b=re(2);
Y=a*exp(b*x);
err=norm(Y-y,2)
plot(x,y,'r+',x,Y,'g');
xlabel('x');
ylabel('y');
legend('原始数据','拟和曲线');