当前位置: 首页 > 数学 >求助一个函数积分求值问题

求助一个函数积分求值问题

作者 利维坦6021
来源: 小木虫 350 7 举报帖子
+关注

在使用fortran编写UMAT子程序时,遇到了一个比较繁琐的积分问题
 \"求助一个函数积分求值问题\"
要求A的值,其余量均为已知量,积分上下限也已知,积分值已知。
为了节省程序计算量,避免过于繁琐的循环嵌套,该方程能不能解出A的值,即给出A的表达式?
貌似是个超越方程...
 \"求助一个函数积分求值问题-1\"
期待大手子解答! 返回小木虫查看更多

今日热帖
  • 精华评论
  • 利维坦6021

    公式在附件里,第一次发帖没设置好
    求助一个函数积分求值问题-2
    捕获.PNG


    求助一个函数积分求值问题-3
    捕获2.PNG

  • alober

    [latex]Q=-\frac{\varepsilon_0(e^{\frac{pA(\varepsilon_0-D)}{\varepsilon_0}}-1)}{pA} \Rightarrow pAQ+\varepsilon_0(-1+e^{\frac{pA(\varepsilon_0-D)}{\varepsilon_0}})=0,\frac{pA}{\varepsilon_0}=x,\frac{\varepsilon_0-D}{Q}=y \Rightarrow xQ+e^{xyQ}-1=0,xyQ+ye^{xyQ}-y=0,xyQ=\xi \Rightarrow \xi+ye^{\xi}-y=0 \Rightarrow ye^ye^{\xi-y}+(\xi-y)=0 \Rightarrow \xi=y-W(ye^y) \Rightarrow x=\frac{1-\frac{1}{y}W(ye^y)}{Q} \Rightarrow A=\frac{\varepsilon_0(1-\frac{Q}{\varepsilon_0-D}W(\frac{\varepsilon_0-D}{Q}e^{\frac{\varepsilon_0-D}{Q}}))}{pQ}[/latex]

  • 终之太刀—晓

    补充一下,3楼层主的解答过程用了一个Lambert W Function,即所求的参数A并不能用初等函数来表示。
    该函数的一些内容可见网址链接:https://baike.baidu.com/item/%E6 ... /2500676?fr=aladdin

  • 利维坦6021

    引用回帖:
    3楼: Originally posted by alober at 2017-08-19 22:27:49
    Q=-\frac{\varepsilon_0(e^{\frac{pA(\varepsilon_0-D)}{\varepsilon_0}}-1)}{pA} \Rightarrow pAQ+\varepsilon_0(-1+e^{\frac{pA(\varepsilon_0-D)}{\varepsilon_0}})=0,\frac{pA}{\varepsilon_0}=x,\frac{\vareps ...

    谢谢您的热心解答,不胜感激
    我有个疑问是因为这个A的值我是要用fortran语言编写出来的。而朗伯W函数应该无法用程序语言编写的吧.....
    所以我想再咨询您一下有没有通过近似表达的方式把朗伯W函数用较易于使用程序书写的函数表达出来?
    因为并不需要精确的A的值,只需要近似值即可所以再麻烦您一下?
    本来通过循环迭加 A的值是很好求的,但由于求A的程序已经包含在一个循环中,所以使用循环会大大增加运算量,故不得已才向坛友们求助。

  • 利维坦6021

    引用回帖:
    4楼: Originally posted by 终之太刀—晓 at 2017-08-20 15:07:17
    补充一下,3楼层主的解答过程用了一个Lambert W Function,即所求的参数A并不能用初等函数来表示。
    该函数的一些内容可见网址链接:https://baike.baidu.com/item/%E6%9C%97%E4%BC%AFW%E5%87%BD%E6%95%B0/2500676? ...

    谢谢您的热心解答,不胜感激
    我有个疑问是因为这个A的值要用fortran语言编写出来的。而朗伯W函数应该无法用程序语言编写的吧.....
    所以我想再咨询您一下有没有通过近似表达的方式把朗伯W函数用较易于使用程序书写的函数表达出来?
    因为并不需要精确的A的值,只需要近似值即可所以再麻烦您一下?
    本来通过循环迭加 A的值是很好求的,但由于求A的程序已经包含在一个循环中,所以使用循环会大大增加运算量,故不得已才向坛友们求助

  • alober

    引用回帖:
    5楼: Originally posted by 利维坦6021 at 2017-08-20 19:24:51
    谢谢您的热心解答,不胜感激
    我有个疑问是因为这个A的值我是要用fortran语言编写出来的。而朗伯W函数应该无法用程序语言编写的吧.....
    所以我想再咨询您一下有没有通过近似表达的方式把朗伯W函数用较易于使用程序 ...

    python 里有 from scipy.special import lambertw 能用,可以看看相应的源参考。

  • peterflyer

    这个方程明显是个关于A的超越方程,不可能得到A的解析式的。

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