求助：matlab函数拟合

wu_580

:dizzy:刚学matlab，完全不知道如何写程序，特来求助！
现有一组数据x=[10 12.5 15 17.5 20 22.5 25 27.5 30 32.5 35 37.5 40 42.5 45 47.5 50];
y=[62.1 77.3 92.5 104 112.9 121.9 125 129.4 134 138.2 142.3 143.2 144.6 147.2 147.8 149.1 150.9];
请高手编程拟合成y=Ax^n形式的函数

wu_580

:'(求助求助。。

wu_580

或者是有一组数据
x=[10 12.5 15 17.5 20 22.5 25 27.5 30 32.5 35 37.5 40 42.5 45 47.5 50];
y=[62.1 77.3 92.5 104 112.9 121.9 125 129.4 134 138.2 142.3 143.2 144.6 147.2 147.8 149.1 150.9];
函数为y=Ax^n，怎么求得参数A和n呢

wu_580

:'(:'(。。。

stats01

用此模型：
fx=@(b,x)b(1)*x.^b(2);
结果为：
b =[ 28.89485503  0.4374601403]
RSS = 870.251379192
MSe =       62.161
R^2 =      0.92392
效果不太好。
若用模型：
fx=@(b,x)b(1)*x.^2./(1+b(2)*x+b(3)*x.^2);
SSy =    11438
b =[  1.084012246  0.0104414964  0.00658753375]
RSS = 11.6518822148
MSe =       0.8963
R^2 =      0.99898
效果好得多(见图)。

wu_580

stats01 发表于 2011-9-25 22:32
用此模型：
fx=@(b,x)b(1)*x.^b(2);
结果为：

非常感谢你的帮忙，激动。。:P

wu_580

stats01 发表于 2011-9-25 22:32
用此模型：
fx=@(b,x)b(1)*x.^b(2);
结果为：

你好，请问你能不能给我这个题目的详细的命令流呢？
我编的程序老是出错，弄了好几天  还没弄好。。

stats01

clear,clc
x=[10 12.5 15 17.5 20 22.5 25 27.5 30 32.5 35 37.5 40 42.5 45 47.5 50];
y=[62.1 77.3 92.5 104 112.9 121.9 125 129.4 134 138.2 142.3 143.2 144.6 147.2 147.8 149.1 150.9];
fx=@(b,x)b(1)*x.^2./(1+b(2)*x+b(3)*x.^2);
b =[1.08  0.01  0.0066]; % initial value
for l=1:5
b=lsqcurvefit(fx,b,x,y)
b=nlinfit(x,y,fx,b)
end
plot(x,y,'o','markerfacecolor','k')
hold on
x1=min(x):range(x)/100:max(x);
y1=fx(b,x1);
plot(x1,y1,'r-','linewidth',2)

wu_580

stats01 发表于 2011-10-2 10:13
clear,clc
x=[10 12.5 15 17.5 20 22.5 25 27.5 30 32.5 35 37.5 40 42.5 45 47.5 50];
y=[62.1 77.3 92.5  ...

谢谢你哦，要不我还不知道要弄多久呢。。