Ja1 =
2*10^b*(10^b*a - 54/25) + 2*20^b*(20^b*a - 86/25) + 2*30^b*(30^b*a - 83/20) + 2*15^b*(15^b*a - 143/50) + 2*5^b*(5^b*a - 127/100) + 2*55^b*(55^b*a - 116/25) + 2*45^b*(45^b*a - 229/50) + 2*50^b*(50^b*a - 231/50) + 2*25^b*(25^b*a - 387/100) + 2*35^b*(35^b*a - 437/100) + 2*40^b*(40^b*a - 451/100) + 2*0^(2*b)*a
Jb1 =
10^b*a*log(100)*(10^b*a - 54/25) + 5^b*a*log(25)*(5^b*a - 127/100) + 15^b*a*log(225)*(15^b*a - 143/50) + 20^b*a*log(400)*(20^b*a - 86/25) + 30^b*a*log(900)*(30^b*a - 83/20) + 25^b*a*log(625)*(25^b*a - 387/100) + 35^b*a*log(1225)*(35^b*a - 437/100) + 40^b*a*log(1600)*(40^b*a - 451/100) + 45^b*a*log(2025)*(45^b*a - 229/50) + 50^b*a*log(2500)*(50^b*a - 231/50) + 55^b*a*log(3025)*(55^b*a - 116/25)
对Ja1=0 和Jb1=0 联立的方程组 求解 最好用牛顿迭代。
2*10^b*(10^b*a - 54/25) + 2*20^b*(20^b*a - 86/25) + 2*30^b*(30^b*a - 83/20) + 2*15^b*(15^b*a - 143/50) + 2*5^b*(5^b*a - 127/100) + 2*55^b*(55^b*a - 116/25) + 2*45^b*(45^b*a - 229/50) + 2*50^b*(50^b*a - 231/50) + 2*25^b*(25^b*a - 387/100) + 2*35^b*(35^b*a - 437/100) + 2*40^b*(40^b*a - 451/100) + 2*0^(2*b)*a
Jb1 =
10^b*a*log(100)*(10^b*a - 54/25) + 5^b*a*log(25)*(5^b*a - 127/100) + 15^b*a*log(225)*(15^b*a - 143/50) + 20^b*a*log(400)*(20^b*a - 86/25) + 30^b*a*log(900)*(30^b*a - 83/20) + 25^b*a*log(625)*(25^b*a - 387/100) + 35^b*a*log(1225)*(35^b*a - 437/100) + 40^b*a*log(1600)*(40^b*a - 451/100) + 45^b*a*log(2025)*(45^b*a - 229/50) + 50^b*a*log(2500)*(50^b*a - 231/50) + 55^b*a*log(3025)*(55^b*a - 116/25)
对Ja1=0 和Jb1=0 联立的方程组 求解 最好用牛顿迭代。
