matlab 拟合求表达式,Matlab离散(x y z)数据,拟合后怎样求拟合曲面表达式?

这是一个比较复杂的拟合问题,合适的模型不好找,花了较大功夫,先前一个较好的结果由于一次误操作被覆盖掉了,只好从头再来。

大概的结果为:

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 上传

442a53943febe9465fc072b4fbe10813.gif

b2a5a3e0dcc7d508e00275fe42fce1b5.gif

edd6622d1e64f7ce0a5de3f06ad6b686.png

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值