用matlab求J最小时,系数a1,a2,a3的值,需要程序.二次函数取极值是的系数问题.J=(3*a1+95*a2-a3+4594305155808867/4503599627370496)^2+(4*a1+122*a2-a3+5215933106418897/4503599627370496)^2+(2*a1+109*a2-a3+8822143

Question

用matlab求J最小时,系数a1,a2,a3的值,需要程序.二次函数取极值是的系数问题.
J =(3*a1 + 95*a2 - a3 + 4594305155808867/4503599627370496)^2 + (4*a1 + 122*a2 - a3 + 5215933106418897/4503599627370496)^2 + (2*a1 + 109*a2 - a3 + 8822143484902521/9007199254740992)^2 + (5*a1 + 88*a2 - a3 + 312242324621947/281474976710656)^2 + (4*a1 + 111*a2 - a3 + 1273249975233745/1125899906842624)^2 + (2*a1 + 95*a2 - a3 + 8417790636486269/9007199254740992)^2 + (a1 + 160*a2 - a3 + 4828402264414901/4503599627370496)^2 + (5*a1 + 160*a2 - a3 + 3235986496898921/2251799813685248)^2 + (3*a1 + 122*a2 - a3 + 1218984465723497/1125899906842624)^2,

Answer 1

clear all
clc
J =@(a)(3*a(1) + 95*a(2) - a(3) + 4594305155808867/4503599627370496)^2 + (4*a(1) + 122*a(2) - a(3) + 5215933106418897/4503599627370496)^2 + (2*a(1) + 109*a(2) - a(3) + 8822143484902521/9007199254740992)^2 + (5*a(1) + 88*a(2) - a(3) + 312242324621947/281474976710656)^2 + (4*a(1) + 111*a(2) - a(3) + 1273249975233745/1125899906842624)^2 + (2*a(1) + 95*a(2) - a(3) + 8417790636486269/9007199254740992)^2 + (a(1) + 160*a(2) - a(3) + 4828402264414901/4503599627370496)^2 + (5*a(1) + 160*a(2) - a(3) + 3235986496898921/2251799813685248)^2 + (3*a(1) + 122*a(2) - a(3) + 1218984465723497/1125899906842624)^2;
a=[1,2,3];
[a J]=fminsearch(J,a)
结果：a =
-0.0796 -0.0039 0.3898
J =
0.0034