- 积分
- 0
- 注册时间
- 2009-5-15
- 仿真币
-
- 最后登录
- 1970-1-1
|
各位学友好,我在利用牛顿迭代法回归公式参数时遇到一些麻烦,请帮忙调试一下吧。谢谢了!下面是三个非线性方程,求解k,a,b.
F1=
- (47/5)^a*(886/25)^b*((47/5)^a*(886/25)^b*k - 79/12500) - (47/5)^a*(4903/100)^b*((47/5)^a*(4903/100)^b*k - 47/12500) - (28/5)^a*(182/5)^b*((28/5)^a*(182/5)^b*k - 129/25000) - (58/5)^a*(559/25)^b*((58/5)^a*(559/25)^b*k - 333/25000) - (47/5)^a*(1393/50)^b*((47/5)^a*(1393/50)^b*k - 241/25000) - (64/5)^a*(2731/100)^b*((64/5)^a*(2731/100)^b*k - 301/25000) - (28/5)^a*(4317/100)^b*((28/5)^a*(4317/100)^b*k - 87/25000) - (28/5)^a*(461/10)^b*((28/5)^a*(461/10)^b*k - 209/50000) - (28/5)^a*(1803/100)^b*((28/5)^a*(1803/100)^b*k - 793/50000) - (28/5)^a*(2973/100)^b*((28/5)^a*(2973/100)^b*k - 291/50000) - (58/5)^a*(5099/100)^b*((58/5)^a*(5099/100)^b*k - 193/50000) - (28/5)^a*(559/25)^b*((28/5)^a*(559/25)^b*k - 1067/100000) - (58/5)^a*(883/20)^b*((58/5)^a*(883/20)^b*k - 479/100000) - (58/5)^a*(934/25)^b*((58/5)^a*(934/25)^b*k - 549/100000) - (28/5)^a*(627/25)^b*((28/5)^a*(627/25)^b*k - 1237/100000) - (28/5)^a*(781/50)^b*((28/5)^a*(781/50)^b*k - 1431/100000) - (47/5)^a*(2059/100)^b*((47/5)^a*(2059/100)^b*k - 1433/100000) - (58/5)^a*(2973/100)^b*((58/5)^a*(2973/100)^b*k - 911/100000) - (28/5)^a*(3257/100)^b*((28/5)^a*(3257/100)^b*k - 663/100000) - (28/5)^a*(3929/100)^b*((28/5)^a*(3929/100)^b*k - 601/100000)=0;
F2=
- 5045629660834333/2251799813685248*(47/5)^a*(886/25)^b*((47/5)^a*(886/25)^b*k - 79/12500) - 5045629660834333/2251799813685248*(47/5)^a*(4903/100)^b*((47/5)^a*(4903/100)^b*k - 47/12500) - 7758651007633171/4503599627370496*(28/5)^a*(182/5)^b*((28/5)^a*(182/5)^b*k - 129/25000) - 5519172823270913/2251799813685248*(58/5)^a*(559/25)^b*((58/5)^a*(559/25)^b*k - 333/25000) - 5045629660834333/2251799813685248*(47/5)^a*(1393/50)^b*((47/5)^a*(1393/50)^b*k - 241/25000) - 5740840160890957/2251799813685248*(64/5)^a*(2731/100)^b*((64/5)^a*(2731/100)^b*k - 301/25000) - 7758651007633171/4503599627370496*(28/5)^a*(4317/100)^b*((28/5)^a*(4317/100)^b*k - 87/25000) - 7758651007633171/4503599627370496*(28/5)^a*(461/10)^b*((28/5)^a*(461/10)^b*k - 209/50000) - 7758651007633171/4503599627370496*(28/5)^a*(1803/100)^b*((28/5)^a*(1803/100)^b*k - 793/50000) - 7758651007633171/4503599627370496*(28/5)^a*(2973/100)^b*((28/5)^a*(2973/100)^b*k - 291/50000) - 5519172823270913/2251799813685248*(58/5)^a*(5099/100)^b*((58/5)^a*(5099/100)^b*k - 193/50000) - 7758651007633171/4503599627370496*(28/5)^a*(559/25)^b*((28/5)^a*(559/25)^b*k - 1067/100000) - 5519172823270913/2251799813685248*(58/5)^a*(883/20)^b*((58/5)^a*(883/20)^b*k - 479/100000) - 5519172823270913/2251799813685248*(58/5)^a*(934/25)^b*((58/5)^a*(934/25)^b*k - 549/100000) - 7758651007633171/4503599627370496*(28/5)^a*(627/25)^b*((28/5)^a*(627/25)^b*k - 1237/100000) - 7758651007633171/4503599627370496*(28/5)^a*(781/50)^b*((28/5)^a*(781/50)^b*k - 1431/100000) - 5045629660834333/2251799813685248*(47/5)^a*(2059/100)^b*((47/5)^a*(2059/100)^b*k - 1433/100000) - 5519172823270913/2251799813685248*(58/5)^a*(2973/100)^b*((58/5)^a*(2973/100)^b*k - 911/100000) - 7758651007633171/4503599627370496*(28/5)^a*(3257/100)^b*((28/5)^a*(3257/100)^b*k - 663/100000) - 7758651007633171/4503599627370496*(28/5)^a*(3929/100)^b*((28/5)^a*(3929/100)^b*k - 601/100000)=0;
F3=
- 1004257997774109/281474976710656*(47/5)^a*(886/25)^b*((47/5)^a*(886/25)^b*k - 79/12500) - 4382489226632347/1125899906842624*(47/5)^a*(4903/100)^b*((47/5)^a*(4903/100)^b*k - 47/12500) - 8094249297019231/2251799813685248*(28/5)^a*(182/5)^b*((28/5)^a*(182/5)^b*k - 129/25000) - 6996958222281519/2251799813685248*(58/5)^a*(559/25)^b*((58/5)^a*(559/25)^b*k - 333/25000) - 7492170254429321/2251799813685248*(47/5)^a*(1393/50)^b*((47/5)^a*(1393/50)^b*k - 241/25000) - 3723635771968247/1125899906842624*(64/5)^a*(2731/100)^b*((64/5)^a*(2731/100)^b*k - 301/25000) - 264948582265051/70368744177664*(28/5)^a*(4317/100)^b*((28/5)^a*(4317/100)^b*k - 87/25000) - 539138992942431/140737488355328*(28/5)^a*(461/10)^b*((28/5)^a*(461/10)^b*k - 209/50000) - 3256144230786041/1125899906842624*(28/5)^a*(1803/100)^b*((28/5)^a*(1803/100)^b*k - 793/50000) - 3819228841605155/1125899906842624*(28/5)^a*(2973/100)^b*((28/5)^a*(2973/100)^b*k - 291/50000) - 4426621327271717/1125899906842624*(58/5)^a*(5099/100)^b*((58/5)^a*(5099/100)^b*k - 193/50000) - 6996958222281519/2251799813685248*(28/5)^a*(559/25)^b*((28/5)^a*(559/25)^b*k - 1067/100000) - 8528901047446683/2251799813685248*(58/5)^a*(883/20)^b*((58/5)^a*(883/20)^b*k - 479/100000) - 2038216946648547/562949953421312*(58/5)^a*(934/25)^b*((58/5)^a*(934/25)^b*k - 549/100000) - 3627729118719845/1125899906842624*(28/5)^a*(627/25)^b*((28/5)^a*(627/25)^b*k - 1237/100000) - 3094594603345047/1125899906842624*(28/5)^a*(781/50)^b*((28/5)^a*(781/50)^b*k - 1431/100000) - 3405628254355651/1125899906842624*(47/5)^a*(2059/100)^b*((47/5)^a*(2059/100)^b*k - 1433/100000) - 3819228841605155/1125899906842624*(58/5)^a*(2973/100)^b*((58/5)^a*(2973/100)^b*k - 911/100000) - 7843900599305793/2251799813685248*(28/5)^a*(3257/100)^b*((28/5)^a*(3257/100)^b*k - 663/100000) - 2066572409409471/562949953421312*(28/5)^a*(3929/100)^b*((28/5)^a*(3929/100)^b*k - 601/100000)=0; |
|