MATLAB无法求解的方程

luohui5151

MATLAB无法求解的方程 %要求解的函数方程 x11，x22，x33，x44，x55为未知数，现在已给定x22，x33，x44，x55的初值如何把x11求解出来 G = 1221438094798879/4503599627370496/(exp(2273642015846523/1125899906842624/(exp(6674651505643105/1125899906842624*x22 - 2793337058193389/562949953421312*x11 + 7786462297992035/4503599627370496*x33 + 1653909689282905/281474976710656*x44 + 5433627487937363/562949953421312*x55 - 5001350306684475/2251799813685248) + 1) - 6135921343060489/562949953421312/(exp(937576994128045/281474976710656*x33 - 3467296016277945/562949953421312*x22 - 3240027304186851/562949953421312*x11 - 6577976738527551/562949953421312*x44 + 8239160716855743/4503599627370496*x55 - 3103180927759649/562949953421312) + 1) - 3688538146920391/1125899906842624/(exp(2684082352387689/562949953421312*x55 - 7708168409607693/18014398509481984*x22 - 942897792097483/140737488355328*x33 - 7741229057071969/562949953421312*x44 - 8672930768552633/1125899906842624*x11 - 4840088719479237/1125899906842624) + 1) + 1471696715307291/281474976710656/(exp(3118351633487947/281474976710656*x11 - 7853219751344125/4503599627370496*x22 + 6105326520408029/562949953421312*x33 + 925849721297521/140737488355328*x44 + 9195142606471/2199023255552*x55 + 1405077869555417/140737488355328) + 1) + 7005314404872831/562949953421312/(exp(3035375437702187/4503599627370496*x44 - 7959388699611645/2251799813685248*x22 - 7195125680488049/1125899906842624*x33 - 4662903614550775/562949953421312*x11 - 4825552313583613/4503599627370496*x55 - 5800951974740635/2251799813685248) + 1) - 4880974594433113/1125899906842624/(exp(627703555971279/140737488355328*x22 - 2910382859279367/562949953421312*x11 + 790261749753887/140737488355328*x33 - 5951410103928677/562949953421312*x44 + 2501091896917881/562949953421312*x55 - 3763606444710747/1125899906842624) + 1) - 3027530274576885/562949953421312/(exp(2000499427241703/140737488355328 - 363233467706395/35184372088832*x22 - 2813681892435847/1125899906842624*x33 - 787638263869643/1125899906842624*x44 - 2574835446796537/4503599627370496*x55 - 497214584725331/140737488355328*x11) + 1) + 5857295351200843/1125899906842624/(exp(5791043294061667/562949953421312*x11 - 396618447962835/140737488355328*x22 - 5400721441860269/9007199254740992*x33 - 5127489132397941/1125899906842624*x44 + 6638755810587995/1125899906842624*x55 - 4230081332402245/18014398509481984) + 1) - 6668846988513495/2251799813685248/(exp(7747187959609555/1125899906842624*x11 - 3229808822116109/562949953421312*x22 - 2386131422021977/2251799813685248*x33 + 1899563400329549/4503599627370496*x44 - 4881852760744081/562949953421312*x55 + 7501729668374471/4503599627370496) + 1) + 939560072923647/1125899906842624/(exp(597736820929765/70368744177664*x44 - 3188718215305537/2251799813685248*x22 - 980219022366973/140737488355328*x33 - 337381694789001/70368744177664*x11 - 5327384306613857/4503599627370496*x55 - 5817219189801755/9007199254740992) + 1) + 4884441831596095/2251799813685248/(exp(774502251351843/281474976710656*x11 + 261462308779547/35184372088832*x22 + 6259411426601239/1125899906842624*x33 + 7502111824535615/562949953421312*x44 + 7051888953861725/4503599627370496*x55 - 5262328655115733/2251799813685248) + 1) - 5846558804573373/1125899906842624/(exp(2875255587654425/281474976710656*x22 - 6016778396995649/1125899906842624*x11 + 8738804164252581/4503599627370496*x33 - 6733868477198299/1125899906842624*x44 - 7296685399108217/562949953421312*x55 + 2825245163422549/1125899906842624) + 1) - 2614440078990679/1125899906842624/(exp(7365284222602151/1125899906842624*x22 - 8876855104251123/1125899906842624*x11 - 5522424747779721/1125899906842624*x33 + 7184380823647707/4503599627370496*x44 + 147032754251151/17592186044416*x55 + 8347972414748649/562949953421312) + 1) - 1462547946948677/140737488355328/(exp(2217235135544767/562949953421312*x22 - 6609094369867589/1125899906842624*x11 + 1169326339060419/281474976710656*x33 - 3031721464553685/4503599627370496*x44 + 3802780210235707/562949953421312*x55 - 719869121028201/70368744177664) + 1) + 2051126352039019/281474976710656/(exp(- 7641095530309263/562949953421312*x11 - 4508090997093077/562949953421312*x22 - 365067712644181/35184372088832*x33 - 8398382136173275/4503599627370496*x44 - 3866637561304089/1125899906842624*x55 - 8030501314131339/2251799813685248) + 1) - 6456467689336373/1125899906842624/(exp(4439171828929803/562949953421312*x11 + 6682439266192645/4503599627370496*x22 + 5091859231832373/1125899906842624*x33 - 7862681374852263/2251799813685248*x44 + 2271535358156113/281474976710656*x55 - 4699296319561971/4503599627370496) + 1) + 5034905105637641/562949953421312/(exp(3657763927337439/562949953421312*x22 - 8252160975396503/562949953421312*x11 - 2354412763110275/2251799813685248*x33 - 4686864789195651/1125899906842624*x44 - 1408115645791733/281474976710656*x55 + 3760016435879917/1125899906842624) + 1) + 3790226335040611/1125899906842624/(exp(8093596502563617/1125899906842624*x11 - 259889862927773/35184372088832*x22 - 6718505227583575/1125899906842624*x33 - 643841380235031/281474976710656*x44 - 2537460212885647/562949953421312*x55 + 4028680120065205/4503599627370496) + 1) - 585625909801649/70368744177664/(exp(3375178609602405/281474976710656*x22 - 6391383990304283/2251799813685248*x11 + 4755996295854039/9007199254740992*x33 + 5012470827772507/562949953421312*x44 - 148031837914901/35184372088832*x55 + 3818961718829973/1125899906842624) + 1) - 2931818734747457/562949953421312/(exp(6511572917613039/2251799813685248*x11 - 296543410462969/70368744177664*x22 - 4869711008120353/562949953421312*x33 + 3225806956942819/562949953421312*x44 - 7281569200177939/2251799813685248*x55 + 8611869000116201/1125899906842624) + 1) + 2287331402949521/2251799813685248) + 1) - 2313582895533293/9007199254740992

lin2009

这不算复杂，matlab在非线性方程数值解方面的功能蛮强的啊。
当其它参数均等于1时，x11 = 0.4744。
用fzero很容易计算出来。

lin2009

只要fzero调用的格式符合要求，就能求出来。
可能有多解，与初值选择有关。

feiyuzhen

听说2010a里面有个全局优化工具箱，可以试试

luohui5151

syms x y
f1=sym(x^2+y^2)
x=0.5;y=0.4
pf1=jacobian(f1,[x y]))
为什么我采用jacobian矩阵求导后，程序没有把x，y值代入到jacobian矩阵，
如何才能代入得到数值解了，请大家帮忙下

rocwoods

1. 
   
2. clear
   
3. syms x y
   
4. f1=sym(x^2+y^2);
   
5. pf1=jacobian(f1,[x y]);
   
6. subs(pf1,{x,y},{0.5,0.4})
   
7. 或者
   
8. x = 0.5;y=0.4;
   
9. eval(pf1)
   

复制代码