查看: 76|回复: 4|关注: 0

[已解决] int 求符号积分,无法得到具体的符号表达式,结果还是int表达式

[复制链接]

新手

7 麦片

财富积分


050


1

主题

4

帖子

0

最佳答案
本帖最后由 Abnerss 于 2020-3-26 10:46 编辑

复制代码
EY2是一个函数,x是自变量,tp是符号,求EY2对x在0到tp求符号积分,代码如下
  1. syms x tp  %两个符号变量
  2. 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函数的形式
  1. 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

主题

4

帖子

0

最佳答案
 楼主| 发表于 6 天前 | 显示全部楼层
  1. syms x tp  %两个符号变量
  2. 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的表达式

  3. int(EY2,x,0,tp) %求积分
复制代码

完整代码,不知道为什么发帖的时候贴的代码,老是显示不全

论坛优秀回答者

5

主题

1919

帖子

552

最佳答案
  • 关注者: 159
发表于 6 天前 | 显示全部楼层
没解析解。只能想办法代入具体数值然后求数值解。
提问请:①准确描述问题②提出你的思考(等着抄作业的一律锁帖)③提供代码文本而非截图④及时反馈

新手

7 麦片

财富积分


050


1

主题

4

帖子

0

最佳答案
 楼主| 发表于 6 天前 | 显示全部楼层
TouAkira 发表于 2020-3-26 10:45
没解析解。只能想办法代入具体数值然后求数值解。

我给tp赋一个具体的值,还是只有int的积分形式,得不到具体的答案,是要换成quad函数吗

论坛优秀回答者

5

主题

1919

帖子

552

最佳答案
  • 关注者: 159
发表于 6 天前 | 显示全部楼层 |此回复为最佳答案
Abnerss 发表于 2020-3-25 22:48
我给tp赋一个具体的值,还是只有int的积分形式,得不到具体的答案,是要换成quad函数吗
...

int基本就是拿来算解析解的,帮助文档里面也有无法求得解析解时,用泰勒展式求近似。求数值解用integral就好,quad已经不推荐使用了。
www.mathworks.com/help/matlab/ref/integral.html
提问请:①准确描述问题②提出你的思考(等着抄作业的一律锁帖)③提供代码文本而非截图④及时反馈
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表