当前位置: 首页 > 计算模拟 >求高手给做一道题,用MATLAB或1stOpt,请把结果直接给出来,我的软件就不行,谢谢

求高手给做一道题,用MATLAB或1stOpt,请把结果直接给出来,我的软件就不行,谢谢

作者 一直在纠结
来源: 小木虫 250 5 举报帖子
+关注

这是我用1stOpt写的一段程序,但是我的软件算出来的结果不唯一,我也不知道到底是什么原因,请哪位大神帮忙给算一下
1stOpt程序:
constant a=1,b=5;
Parameter m,n,c,d,e,f,g,h,i;
Function a=c*exp(31.178-4.606*b)+c^2*exp(36.796-4.606*b)+c^3*exp(42.553-4.606*b)+c^4*exp(47.503-4.606*b)+exp(7.607-2.303*b)+exp(-38.744)+exp(-38.533+2.303*b)+exp(-62.699+4.606*b)+exp(25.721-4.606*b)+exp(-38.264+2.303*b)+exp(-67.656+4.606*b)+d*exp(56.16-6.909*b);
m+n=c+c*exp(21.304-2.303*b)+c*exp(31.178-4.606*b)+2*c^2*exp(36.796-4.606*b)+3*c^3*exp(42.553-4.606*b)+4*c^4*exp(47.503-4.606*b)+2*c^2*d*exp(42.771-4.606*b)+c*d*exp(45.204-4.606*b);
m=d+d*exp(23.762-2.303*b)+d*exp(134.267-4.606*b)+c^2*d*exp(42.771-4.606*b)+c*d*exp(45.204-4.606*b)+d*exp(56.16-6.909*b);
2*a+c*exp(21.304-2.303*b)+exp(-2.303*b)=2*m+exp(-32.271+2.303*b);
e=c*exp(31.178-4.606*b);
f=c^2*exp(36.796-4.606*b);
g=c^3*exp(42.553-4.606*b);
h=c^4*exp(47.503-4.606*b);
i=d*exp(56.16-6.909*b);
上面只是算b=5时,解方程,我现在想让b从5到12变化,间隔0.2,e,f,g,h,i的值,我的软件实在是不行,一样的值算出来不一样的数,麻烦大神补充完整程序,然后算出来这个结果,请直接算出结果,就算给了过程软件不能用也是白搭 返回小木虫查看更多

今日热帖
  • 精华评论
  • 一直在纠结

    上面应该是可以化简的,第一个方程里面就包含有e,f,g,h,i这些变量,刚想到了,辛苦大神再改一下

  • 一直在纠结

    为什么没有人呢

  • dingd

    b=5时似乎都没有精确实数解。公式对否?自己好好检查下。

  • 一直在纠结

    引用回帖:
    4楼: Originally posted by dingd at 2014-04-23 09:07:58
    b=5时似乎都没有精确实数解。公式对否?自己好好检查下。

    呃。。亲爱的版主,您终于出现了,我把ph范围弄错了,应该是从7到14,b从7到14变化,间隔0.2,我想用这些数据作图,不需要十分精确的解,有小数点后三位就可以了,但是这个范围怎么写,我也不会。。。
    方程从头到尾核对了一遍,应该没有错误了
    constant a=1,b=7;
    Parameter m,n,c,d,e,f,g,h,i;
    Function a=e+f+g+h+exp(7.607-2.303*b)+exp(-38.744)+exp(-38.533+2.303*b)+exp(-62.699+4.606*b)+exp(25.721-4.606*b)+exp(-38.264+2.303*b)+exp(-67.656+4.606*b)+i;
    m+n=c+c*exp(21.304-2.303*b)+c*exp(31.178-4.606*b)+2*c^2*exp(36.796-4.606*b)+3*c^3*exp(42.553-4.606*b)+4*c^4*exp(47.503-4.606*b)+2*c^2*d*exp(42.771-4.606*b)+c*d*exp(45.204-4.606*b);
    m=d+d*exp(23.762-2.303*b)+d*exp(134.267-4.606*b)+c^2*d*exp(42.771-4.606*b)+c*d*exp(45.204-4.606*b)+d*exp(56.16-6.909*b);
    2*a+c*exp(21.304-2.303*b)+exp(-2.303*b)=2*m+exp(-32.271+2.303*b);
    e=c*exp(31.178-4.606*b);
    f=c^2*exp(36.796-4.606*b);
    g=c^3*exp(42.553-4.606*b);
    h=c^4*exp(47.503-4.606*b);
    i=d*exp(56.16-6.909*b);
    方程组重写写了,再辛苦您一下

  • dingd

    最后两个点无精确实数解:


    循环常数 b        目标函数值        m        n        c        d        e        f        g        h        i
    7        1.08815041857896E-32        -1.11896652853502        1.1719558682823        -0.0237796683244757        -5.49442917506541E-45        -0.0082057156627447        0.0537265132120331        -0.404229292449488        1.35703565864475        5.55111512312736E-17
    7.2        1.95000406869208E-32        -0.657787425825305        1.56307107190141        -0.0294880845886294        -8.11452547186287E-45        -0.00405028103270097        0.0328850208250031        -0.306816009562855        1.27726868994219        1.06951036724574E-24
    7.4        9.89164390284762E-33        -0.300791696735386        1.86174004373355        -0.0366742714152774        -9.32213699297834E-45        -0.00200506295216517        0.0202467773307324        -0.234936629387645        1.21638180117125        5.55111512312594E-17
    7.6        3.7301727135781E-33        1.94359028462639        3.81898482596624        0.0421670631861849        1.51330488330823E-43        0.000917630406324198        0.0106538738433603        0.142139321123577        0.846145944006841        5.55111512312646E-17
    7.8        6.1815807951956E-67        0.193728865677626        2.26282577206383        -0.0571092952349057        3.78955422963424E-44        -0.000494685482184936        0.00777861658273091        -0.140553769746359        1.13320109190674        -7.86230296378159E-34
    8        1.03359833908599E-70        0.363617506496734        2.39262444752267        -0.0714463068540366        1.7869428033513E-43        -0.000246337240507842        0.00484592017555065        -0.109544196140345        1.10490984353902        -1.016660341896E-35
    8.2        1.14423802950463E-32        0.497074498552954        2.48752189101866        -0.0894954432457326        6.13704668865065E-43        -0.000122822839706104        0.00302654027004954        -0.0856999784052749        1.08277774590605        5.55111512314656E-17
    8.4        1.13958785550941E-71        1.38529414328503        3.42912040481492        0.108674522789926        4.2968680663811E-42        5.93654652787266E-5        0.00177634785983418        0.061078602431067        0.937075359951029        -3.37577755618637E-36
    8.6        4.93038065763132E-32        0.684999091614139        2.59477338781432        -0.140824802498377        5.33792039999609E-42        -3.06205655894182E-5        0.00118729532675495        -0.0529018965056987        1.05173924050175        -5.70585946822379E-25
    8.8        6.30835887238605E-69        1.24455657820896        3.37712851283785        0.173296910207584        2.43651348326087E-41        1.49986717381385E-5        0.000715665396793327        0.0392404673340716        0.960025294629748        -7.94251775885763E-35
    9        5.57191085138282E-70        0.802178119943777        2.61200109708009        -0.222180410331444        3.94545471655265E-41        -7.6541270802095E-6        0.000468239268144061        -0.0329159905621633        1.03245321042704        2.36048954596152E-35
    9.2        1.09607385267155E-32        0.843123790738861        2.58930995647316        -0.279263345418199        1.04181275915997E-40        -3.82940908742938E-6        0.000294450366489993        -0.0260171288520858        1.02572511543282        -7.49944579293137E-25
    9.39999999999999        1.66656023688867E-66        1.12317050999823        3.45044143996149        0.347587243522747        3.48671216236527E-40        1.89718455297027E-6        0.00018156820858485        0.0199681104597551        0.979847491876035        -1.29095322794858E-33
    9.59999999999999        2.26139903220066E-68        0.901242535526899        2.47574133597447        -0.441606278992903        7.02884816676999E-40        -9.59420322698149E-7        0.000116656932566657        -0.0162996817457739        1.01618329205836        1.5037948773065E-34
    9.79999999999999        1.3076597355476E-32        0.921606714200915        2.37891999885285        -0.555514738997521        1.80576066321297E-39        -4.80393567183525E-7        7.34783116551772E-5        -0.0129148119617782        1.01284119948576        -1.07379631641187E-24
    9.99999999999999        2.35623956894689E-68        0.93774914446077        2.24906679588421        -0.698922740175862        4.61607983537739E-39        -2.40579540600274E-7        4.62971504618625E-5        -0.0102380386500319        1.01019129334161        -1.53500474546607E-34
    10.2        5.80633895309345E-33        1.0491056008603        3.91710390135972        0.87592828275067        1.29741201974882E-38        1.20012364718753E-7        2.8944159438632E-5        0.00802163364703053        0.991948347988246        5.55111555449154E-17
    10.4        1.28110721321906E-32        1.03895983036159        4.1360141785191        1.10323306218062        3.22797040876456E-38        6.01662072372545E-8        1.8276216883436E-5        0.00637950385180914        0.993600629318216        -1.1102230275806E-16
    10.6        1.29898949743024E-32        1.03085953002902        4.41554386763406        1.38940435030532        8.04641365278919E-38        3.01607366262953E-8        1.15381717191046E-5        0.00507223054217183        0.994913501625013        -1.11022302492814E-16
    10.8        8.1987619264931E-70        0.974972599116266        1.22637764618195        -1.75322858877193        1.91190959188063E-37        -1.51488532042208E-8        7.31281691242339E-6        -0.00405654760041722        1.00404413976384        -2.86334802729907E-35
    11        9.91753254658644E-33        0.979882950396948        0.777208542197053        -2.20681437659809        4.82748818537367E-37        -7.58988061159153E-9        4.61176474648236E-6        -0.00322007718642992        1.00320519398841        1.87259893322153E-20
    11.2        2.78001555490778E-68        1.01480216380344        5.78802204214813        2.77432182530923        1.25603135106382E-36        3.79799071242519E-9        2.90119644187527E-6        0.00254663853796385        0.997428631865274        1.66733786466128E-34
    11.4        6.6983110616527E-33        1.01115083731299        6.5046115719843        3.49324459196391        3.14417795003961E-36        1.90350476277989E-9        1.83083570427626E-6        0.0020235395682365        0.997926100860203        -6.06779262889514E-21
    11.6        3.81967483794191E-33        1.0078942418821        7.40800567633982        4.39830581719567        7.87368820082206E-36        9.53977949674701E-10        1.15528770528374E-6        0.00160771441869978        0.998279183226935        5.55111514914529E-17
    11.8        5.00147302432804E-32        0.989128134611915        -2.54568838348024        -5.54115648397585        1.94127865952329E-35        -4.78389050367273E-10        7.29874278556426E-7        -0.00127962231855593        1.0010134355543        5.28662967471543E-23
    12        2.01809714998058E-32        0.988944452482783        -3.97797906818926        -6.97507702441908        4.87617519246311E-35        -2.3969438114538E-10        4.60334121996374E-7        -0.00101591063814255        1.00037355204889        -2.14502213706793E-21
    12.2        1.97295157725905E-68        0.987382583442022        -5.78168123343215        -8.77896187323228        1.22310837846732E-34        -1.20082532704463E-10        2.90261477389797E-7        -0.000806242887942828        0.999232946270397        -1.40461794707919E-34
    12.4        2.01392886537046E-33        0.991716039783121        14.0420081840445        11.0421025991835        3.08630535926332E-34        6.01196149498155E-11        1.82782475596419E-7        0.000638586400712479        0.995472183765956        1.4589001950303E-20
    12.6        2.93981607994161E-71        0.977574048059465        -10.9074944552552        -13.8855404831582        7.6431660691803E-34        -3.00922950811415E-11        1.15049437618631E-7        -0.000505452817096776        0.990834444886816        5.42200708222827E-36
    12.8        3.78537816971374E-33        0.971822470578542        20.349527728469        17.4133940451357        1.90889766014442E-33        1.5021138345922E-11        7.20200967838101E-8        0.000396798901251442        0.975464736453316        5.16781885041273E-21
    13        1.53125205090141E-68        0.953379678601214        24.5255979320222        21.7169946051664        4.70471885955135E-33        7.45672441622696E-12        4.45876122442278E-8        0.000306370730022418        0.939300774206198        -1.23743769576538E-34
    13.2        5.92559197439191E-33        0.921551013763185        -24.1877988198999        -26.6581691422972        1.14250798393358E-32        -3.64344242945068E-12        2.67425815970076E-8        -0.000225562542798521        0.848896975780431        -1.86323413166332E-20
    13.4        3.23181912085178E-65        0.879123846243464        32.6349891752993        31.0302780548129        2.73818003543277E-32        1.68806778977553E-12        1.44225465604398E-8        0.000141599248740363        0.620303600254526        5.68490907653921E-33
    13.6        1.46471903666259E-55        0.807126700311523        19.9258707315345        20.5423392338415        6.31576162606172E-32        4.44817675253108E-13        2.5159287167908E-9        1.63524022353659E-5        0.0474230130395801        1.46879194692354E-28
    13.8        0.322731697976617        0.693564198630504        -0.6975100615285        -0.0039457502869636        1.36346112910222E-31        -0.231923729025936        -0.231923728935539        -0.231923729488932        -0.231923731159435        -0.231923730562563
    14        4.17630113542854        0.514291767218135        -0.513724320611137        0.000567446623204058        2.54002353693712E-31        -0.834296227450986        -0.834296227160368        -0.834296215253464        -0.834296225229582        -0.834296216394722,

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓