# 7 麦片

050

1主题 0最佳答案
 本帖最后由 Abnerss 于 2020-3-26 10:46 编辑 复制代码EY2是一个函数，x是自变量，tp是符号，求EY2对x在0到tp求符号积分，代码如下 syms x tp  %两个符号变量 EY2=exp((55071200000000*(-1)^(27/50)*10^(19/25)*3129^(27/50)*gamma(23/50))/(10211638563*(45*log(1850) - 45*log(1829) - 52*log(41113) + 52*log(41300))^(123/50)) - (206517*311656599332291100290357921385515998191072476426947805281733115522809489835979210700000000000000000000000^(1/50)*x^(73/50)*igamma(23/50, x*((27*log(1829))/20000 - (27*log(1850))/20000 + (39*log(41113))/25000 - (39*log(41300))/25000)))/(1418283133750*(-x)^(73/50)*((27*log(1850))/20000 - (27*log(1829))/20000 - (39*log(41113))/25000 + (39*log(41300))/25000)^(123/50)) - (8979*(-1)^(23/50)*2^((447*x)/100000)*7^((39*x)/25000)*25^((291*x)/100000)/31^((27*x)/20000)*37^((27*x)/20000)*59^((21*x)/100000)/41113^((39*x)/25000)*311656599332291100290357921385515998191072476426947805281733115522809489835979210700000000000000000000000^(1/50)*x^(48/25))/(28365662675*(-x)^(73/50)*((27*log(1829))/20000 - (27*log(1850))/20000 + (39*log(41113))/25000 - (39*log(41300))/25000)^2) - (246*(-1)^(23/50)*2^((447*x)/100000)*7^((39*x)/25000)*25^((291*x)/100000)/31^((27*x)/20000)*37^((27*x)/20000)*59^((21*x)/100000)/41113^((39*x)/25000)*311656599332291100290357921385515998191072476426947805281733115522809489835979210700000000000000000000000^(1/50)*x^(73/25))/(1134626507*(-x)^(73/50)*((27*log(1829))/20000 - (27*log(1850))/20000 + (39*log(41113))/25000 - (39*log(41300))/25000)));复制代码 结果如下，得不到具体的表达式，还是int函数的形式 int(exp(- (2245954264241636518445140297121792*x^(73/50)*igamma(23/50, -(3241358370611*x)/144115188075855872))/(461633230275677700738238125*(-x)^(73/50)) + (2^((447*x)/100000)*7^((39*x)/25000)*25^((291*x)/100000)/31^((27*x)/20000)*37^((27*x)/20000)*59^((21*x)/100000)/41113^((39*x)/25000)*x^(48/25)*(- 152985168682033793032572425729/15860216652252021509609660 - 1513752883175318685211812869569i/19825270815315026887012075))/(-x)^(73/50) + (2^((447*x)/100000)*7^((39*x)/25000)*25^((291*x)/100000)/31^((27*x)/20000)*37^((27*x)/20000)*59^((21*x)/100000)/41113^((39*x)/25000)*x^(73/25)*(2682479670041140082555354332031/18052268468082400806439221198848 + 5308503261546323204431977179967i/4513067117020600201609805299712))/(-x)^(73/50) - (2521052312205759/2147483648 - 623630747855207i/67108864)), x, 0, tp) 复制代码因为我是想得到一个符号表达式，然后进行使用自带的optimtool中的ga遗传算法求得到最小值，但是现在无法得到具体的符号表达式，无法进行下一步的遗传算法

# 7 麦片

050

1主题 0最佳答案
楼主| 发表于 6 天前 | 显示全部楼层
 syms x tp  %两个符号变量 EY2=exp((55071200000000*(-1)^(27/50)*10^(19/25)*3129^(27/50)*gamma(23/50))/(10211638563*(45*log(1850) - 45*log(1829) - 52*log(41113) + 52*log(41300))^(123/50)) - (206517*311656599332291100290357921385515998191072476426947805281733115522809489835979210700000000000000000000000^(1/50)*x^(73/50)*igamma(23/50, x*((27*log(1829))/20000 - (27*log(1850))/20000 + (39*log(41113))/25000 - (39*log(41300))/25000)))/(1418283133750*(-x)^(73/50)*((27*log(1850))/20000 - (27*log(1829))/20000 - (39*log(41113))/25000 + (39*log(41300))/25000)^(123/50)) - (8979*(-1)^(23/50)*2^((447*x)/100000)*7^((39*x)/25000)*25^((291*x)/100000)/31^((27*x)/20000)*37^((27*x)/20000)*59^((21*x)/100000)/41113^((39*x)/25000)*311656599332291100290357921385515998191072476426947805281733115522809489835979210700000000000000000000000^(1/50)*x^(48/25))/(28365662675*(-x)^(73/50)*((27*log(1829))/20000 - (27*log(1850))/20000 + (39*log(41113))/25000 - (39*log(41300))/25000)^2) - (246*(-1)^(23/50)*2^((447*x)/100000)*7^((39*x)/25000)*25^((291*x)/100000)/31^((27*x)/20000)*37^((27*x)/20000)*59^((21*x)/100000)/41113^((39*x)/25000)*311656599332291100290357921385515998191072476426947805281733115522809489835979210700000000000000000000000^(1/50)*x^(73/25))/(1134626507*(-x)^(73/50)*((27*log(1829))/20000 - (27*log(1850))/20000 + (39*log(41113))/25000 - (39*log(41300))/25000))); %函数EY2，关于x的表达式 int(EY2,x,0,tp) %求积分复制代码 完整代码，不知道为什么发帖的时候贴的代码，老是显示不全

# 论坛优秀回答者

5主题 552最佳答案
• 关注者： 159

 没解析解。只能想办法代入具体数值然后求数值解。

# 7 麦片

050

1主题 0最佳答案
楼主| 发表于 6 天前 | 显示全部楼层
 TouAkira 发表于 2020-3-26 10:45 没解析解。只能想办法代入具体数值然后求数值解。 我给tp赋一个具体的值，还是只有int的积分形式，得不到具体的答案，是要换成quad函数吗

# 论坛优秀回答者

5主题 552最佳答案
• 关注者： 159