24小时热门版块排行榜    

CyRhmU.jpeg
查看: 996  |  回复: 5

luohui5151

铁虫 (小有名气)

[交流] 【求助】MATLAB无法求解的方程已有1人参与

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
回复此楼
一辈子做好一件事情
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunxiao

荣誉版主 (著名写手)

小木虫钓鱼协会会长

这也太夸张了吧,你将你的系数简化下,再发给大家看看吧,这一看都晕了,呵呵
西南交大隧道方向招硕、博士
2楼2010-03-27 11:09:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by luohui5151 at 2010-03-26 21:24:36:
MATLAB无法求解的方程
%要求解的函数方程
x11,x22,x33,x44,x55为未知数,现在已给定x22,x33,x44,x55的初值如何把x11求解出来
G =

1221438094798879/4503599627370496/(exp(2273642015846523/112589 ...

不,这是人类无法看懂的方程
3楼2010-03-27 15:55:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
nono2009(金币+2):鼓励应助!欢迎常来仿真模拟版交流。 2010-04-02 12:02
%%%myfun1.m%%%
function f=myfun1(x11)
global x22;
global x33;
global x44;
global x55;
f=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 ;
%%%%%%%%%%%%%%%%%%%%%%%%%


%%%%%%%fsolvemyfun1.m%%%%%%%%%%%
global x22 x33 x44 x55;
x22=2;
x33=3;
x44=4;
x55=5;

x11=fsolve(@myfun1,1,optimset('Display','iter'))

%%%%%%%%求解结果%%%%%%%%%%%%
                                         Norm of      First-order   Trust-region
Iteration  Func-count     f(x)          step         optimality    radius
     0          2    2.42443e-005                     1.98e-005               1
     1          3    2.42443e-005              1      1.98e-005               1
     2          5    3.15259e-006           0.25      4.87e-005            0.25
     3          7    2.10691e-007       0.064729      1.96e-005           0.625
     4          9    2.59214e-010      0.0107716      6.38e-007           0.625
Optimization terminated: first-order optimality is less than options.TolFun.

x =

    0.6960
明德厚学、求是创新
4楼2010-04-02 11:10:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主

help fsolve
看看fsolve的用法,应该可以解。
上面的x22,x33,x44,x55是随便设的,x11的初值也是随便选的。
明德厚学、求是创新
5楼2010-04-02 11:11:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

icqaa

新虫 (小有名气)

eval(表达式就可以了)
6楼2010-11-09 23:29:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 luohui5151 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见