- 积分
- 1
- 注册时间
- 2011-8-21
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2011-9-18 10:11:12
|
显示全部楼层
来自 江苏扬州
本帖最后由 stats01 于 2011-9-18 10:11 编辑
clear,clc
fx1=@(b,x)b(1)+b(2)*x.^b(3) + b(4)*x.^3+ b(5)*x.^2 +b(6)*x;
fx2=@(b,x)b(1)+b(7)*x.^b(8) + b(9)*x.^3+ b(10)*x.^2 +b(11)*x;
x=[0 0 0.0772 0.0772 0.1409 0.2439 0.2439 0.3956 0.4106 0.4106 0.4325 0.5772 0.5772 0.6419 0.7439 0.7439 0.9106 0.9106 0.959 1.0772 1.0772 1.2439 1.2439 1.2997 1.4106 1.4106 1.5772 1.5772 1.6507 1.7439 1.7439 1.9106 1.9106 1.9942 2.0479 2.0772 2.0772 2.214 2.2439 2.2439 2.3093 2.3435 2.3811]';
y=[0.8333 0.8333 0.7561 0.9106 1 0.7671 1.103 0.8333 0.8356 1.1454 1.1667 0.8503 1.3031 1.3333 0.8627 1.4354 0.8627 1.4707 1.5 0.8627 1.6182 0.8627 1.6373 1.6667 0.8627 1.7775 0.8627 1.804 1.8333 0.8627 1.9265 0.8627 1.9169 1.8333 1.6667 0.8676 1.6368 1.5 0.9004 1.3987 1.3333 1 1.1667]';
x1=linspace(0,2.438);
b=[0.817731470822958,-0.660171615762068,0.566416207806644,0.159008423164863,-0.676688824702786,1.23360131960106,-560.036919356186,1.99547820437660,-1.16910821815770,559.796257923457,2.12294518543622];
figure(1),clf
plot(x,y,'o','markersize',10,'markerfacecolor','k')
y1=fx1(b,x1);
y2=fx2(b,x1);
hold on
plot(x1,y1,'r-',x1,y2,'r-','linewidth',3)
legend('data','fit','location','best') |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|