Question

用matlab最小二乘法拟合指数函数

本人新手，刚入门，请高手帮帮忙，跪谢。。。
形如y=A*(B^m)*(x^nC)方程中，,A,B,C是常数，x和y也有多组已知的值，在matlab中如何编程用最小二乘法拟合并求得m和n的值，急用，望大神指点，多谢。
帮忙写下程序吧，拜托。。。

Answer 1

myfun。m
function y = myfun(beta,x)
A=....
B=....
C=....

m=beta(1);

n=beta(2);
y=A*(B^m)*(C*x^n)

窗口下执行以下命令
beta0=rand(1，2)
[beta,r,J]=nlinfit(x,y,@myfun,beta0);

Answer 2

呵呵，还需要转换一次啊。及y1=exp(z1)

clear all
x=[1971:1990];
y=[8.5229 8.7177 8.9221 9.0859 9.2420 9.3717 9.4974 9.6259 9.7542 9.8705 10.0072 10.1654 10.3008 10.4357 10.5851 10.7507 10.9300 11.1026 11.2704 11.4333];
fun=inline('exp(a(1)*t+a(2))','a','t')
a=nlinfit(x,y,fun,[0.01 -20])
xx=1970:1990;
yy=exp(a(1)*xx+a(2));
plot(x,y,'o',xx,yy)

z=log(y);
p=polyfit(x,z,1)
z1=polyval(p,x);
y1=exp(z1);
figure
plot(x,y,'*',x,y1)

结果：
a =

0.014631 -26.68

p =

0.01468 -26.777

追问

我的方程指数不是exp的哇

Answer 3

matlab有自带的拟合程序的。