- 积分
- 1
- 注册时间
- 2011-8-21
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2011-10-2 10:13:39
|
显示全部楼层
来自 江苏扬州
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) |
评分
-
1
查看全部评分
-
|