这是一个比较复杂的拟合问题,合适的模型不好找,花了较大功夫,先前一个较好的结果由于一次误操作被覆盖掉了,只好从头再来。
大概的结果为:
fxy=@(b,x,y)(b(1)+b(2)*x+b(3)*y+b(4)*y.^2+b(5)*x.*y+b(6)*x.^3+b(7)*y.^3+b(8)*x.*y.^2+b(9)*x.^2.*y+b(10)*y.^4)./(1+b(11)*exp(b(12)*x+b(13)*y+b(14)*x.*y+b(15)*x.^2+b(16)*y.^2+b(17)*x.*y.^2+b(18)*x.^2.*y+b(19)*x.^4+b(20)*y.^4+b(21)*x.*y.^3+b(22)*x.^3.*y+b(23)*x.^2.*y.^3+b(24)*x.^3.*y.^2))+b(25)*exp(-b(26)*(x-b(27)).^2-b(28)*(y-b(29)).^2-b(30)*(x-b(27)).*(y-b(29)))+b(31)*exp(-b(32)*(x-b(33)).^2-b(34)*(y-b(35)).^2-b(36)*(x-b(33)).*(y-b(35)))+b(37)*exp(-b(38)*(x-b(39)).^2-b(40)*(y-b(41)).^2-b(42)*(x-b(39)).*(y-b(41)))+b(43)*exp(-b(44)*(x-b(45)).^2-b(46)*(y-b(47)).^2-b(48)*(x-b(45)).*(y-b(47)));
SSy =
5.0647e+05
b =[-1.906319255 0.01917994757 0.1247831779 -0.0004634608052 0.006252063514 -0.0005492145368 -1.079901965e-06 -3.122455158e-05 0.0001121665551 8.716637171e-09 6.922776767e-17 6.39450919 1.245531974 -0.1824594557 -0.1122153644 -0.006254560698 0.001143391639 0.0005996611783 9.443758225e-05 3.992966502e-08 -2.498338809e-06 -4.247505594e-06 2.293866656e-08 -1.433235995e-07 3.621083059 0.4874165498 3.021838237 0.002923549489 55.49181859 0.07019583991 8.977797772 0.03244844398 33.53256408 0.0008291060394 62.11816028 -0.003804294281 2.852426325 0.003289809204 24.14980188 0.002767079285 174.7631641 0.003190887028 5.771802308 0.5697615602 45.0445104 0.005386779555 37.30806025 -0.08389043595]
RSS = 16204.491819
MSe = 1.3749
R^2 = 0.96801
nh1439.png
(669.65 KB, 下载次数: 6)
2017-12-20 18:46 上传



2268

被折叠的 条评论
为什么被折叠?



