Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#2 |
高级会员
注册日期: 2007-06-24
年龄: 70
帖子: 188
声望力: 21 ![]() |
![]()
试试下列小程序:
x=[0.0032, 0.0074, 0.0107, 0.0136, 0.0166, 0.0224, 0.0282, 0.0341, 0.0399, 0.0457, 0.0545, 0.0632, 0.072, 0.0807, 0.0903, 0.1102, 0.131, 0.1557, 0.1807, 0.2107, 0.2449, 0.2852, 0.3391, 0.4204, 0.5282, 0.591, 0.671, 0.7899, 0.9693, 1.086, 1.2282, 1.4004, 1.6016, 1.8382, 2.1216, 2.466, 2.8771, 3.3882, 3.9971, 4.7438, 5.6816, 6.8193, 7.9749]'; y=[0.04, 0.07, 0.1, 0.13, 0.16, 0.22, 0.28, 0.34, 0.41, 0.47, 0.56, 0.65, 0.72, 0.82, 0.91, 1.09, 1.27, 1.46, 1.63, 1.82, 2.0, 2.17, 2.35, 2.53, 2.67, 2.73, 2.79, 2.85, 2.91, 2.94, 2.97, 3.0, 3.03, 3.06, 3.09, 3.12, 3.15, 3.18, 3.21, 3.21, 3.27, 3.3, 3.33]'; fx=@(b,x)(b(1)*x+b(2)*x.^2+b(3)*x.^3)./(1+b(4)*x+b(5)*x.^2+b(6)*x.^3); b=[9.7865318, 32.775392, 27.227736, 1.967689, 13.118039, 7.8103828]; x1=linspace(0,8) plot(x,y,'o','markerfacecolor','k','linewidth',3); hold on y1=fx(b,x1); plot(x1,y1,'linewidth',2); |
![]() |
![]() |