MATLAB中文论坛

 找回密码
 注册

QQ登录

只需一步,快速开始

扫一扫,访问微社区

查看: 472|回复: 0|关注: 0

[未答复] 曲面拟合结果(招领)

[复制链接]

论坛优秀回答者

15

主题

1579

帖子

91

最佳答案
  • 关注者: 107
发表于 2016-4-21 12:02:51 | 显示全部楼层 |阅读模式
u=[0.5 1 2 3 4 5];
v=[0 5 10 15 20 30 40 50 60 70 80];
y=[0.385 0.420 0.493 0.533 0.569 0.589;
   0.413 0.472 0.546 0.586 0.613 0.631;
   0.531 0.584 0.630 0.660 0.681 0.697;
   0.700 0.7 0.719 0.738 0.749 0.758;
   0.788 0.789 0.793 0.805 0.81 0.814;
   0.892 0.888 0.888 0.891 0.892 0.892;
   0.949 0.948 0.944 0.944 0.944 0.944;
   0.977 0.973 0.967 0.967 0.967 0.967;
   0.989 0.985 0.983 0.983 0.983 0.983;
   0.996 0.992 0.990 0.990 0.990 0.990;
   1 0.998 0.998 0.998 0.998 0.998];

fx=@(b,u,v)b(1)+b(2)*u+b(3)*v+b(4)*(-1+2*u.^2+2*v.^2)+b(5)*(u.^2-v.^2)+b(6)*(2*u.*v)+b(7)*((-2)*u+3*u.*v.^2+3*u.^3)+b(8)*((-2)*v+3*v.^3+3*v.*u.^2)+b(9)*(1-6*u.^2-6*v.^2+6*u.^4+12*u.^2.*v.^2+6*v.^4);
b =[0.3483212387  0.05956257475  0.02197588753  -0.0005987895006  -0.001125911034  -0.0009796514683  5.91302647e-06  -1.291955172e-06  5.355802795e-09]
RSS = 0.0195314293528
MSe = 0.00034266
R^2 = 0.99124
Using following model is simple:
fx=@(b,u,v)(b(1)+b(2)*u+b(3)*v+b(4)*v.^2+b(5)*v.^3+b(6)*u.*v.^2)./(1+b(7)*u+b(8)*v+b(9)*u.*v);
b =[0.2634606304  0.2755042527  0.0285420655  -0.0005987091056  3.215281532e-06  -1.690269766e-05  0.3682891916  -0.008010369816  -0.002494424636]
RSS = 0.0158720518945
MSe = 0.00027846
R^2 =  0.99288
如图:

nh1231_1.jpg
nh1231_2.jpg
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /2 下一条

快速回复 返回顶部 返回列表