应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 应用数学与力学 » 求助一个函数积分求值问题
【调剂】北京石油化工学院2024年16个专业接受调剂
81/1返回列表
查看: 1087  |  回复: 7
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者利维坦6021将赠送您 5 个金币

利维坦6021

新虫 (初入文坛)

[求助] 求助一个函数积分求值问题已有2人参与

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2017-08-19 20:54:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

利维坦6021

新虫 (初入文坛)
公式在附件里，第一次发帖没设置好
求助一个函数积分求值问题-2
捕获.PNG


求助一个函数积分求值问题-3
捕获2.PNG
赞一下 回复此楼
2楼2017-08-19 20:58:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alober

木虫 (著名写手)
赞一下 回复此楼
3楼2017-08-19 22:27:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

终之太刀—晓

铁杆木虫 (著名写手)

数学爱好者

【答案】应助回帖

感谢参与，应助指数 +1
补充一下，3楼层主的解答过程用了一个Lambert W Function，即所求的参数A并不能用初等函数来表示。
该函数的一些内容可见网址链接：https://baike.baidu.com/item/%E6 ... /2500676?fr=aladdin
赞一下 回复此楼
PreferenceforMathematics
4楼2017-08-20 15:07:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

利维坦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的程序已经包含在一个循环中，所以使用循环会大大增加运算量，故不得已才向坛友们求助。
赞一下 回复此楼
5楼2017-08-20 19:24:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

利维坦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的程序已经包含在一个循环中，所以使用循环会大大增加运算量，故不得已才向坛友们求助。
赞一下 回复此楼
6楼2017-08-20 19:34:33
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alober

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

python 里有 from scipy.special import lambertw 能用，可以看看相应的源参考。
赞一下 回复此楼
7楼2017-08-20 20:16:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer

【答案】应助回帖

感谢参与，应助指数 +1
这个方程明显是个关于A的超越方程，不可能得到A的解析式的。
赞一下 回复此楼
8楼2017-08-20 23:56:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 利维坦6021 的主题更新
81/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 想问一下有没有需要文章但是没时间带学生的研究生导师。 +9 lekinna 2024-04-19 19/950 2024-04-20 05:35 by charles-c
[基金申请] 申请省自然科学基金，研究区能否是省外区域 100+3 喜欢兔兔的我 2024-04-15 13/650 2024-04-19 21:49 by 喜欢兔兔的我
[公派出国] CSC在职申请求问 +7 Amber0919 2024-04-16 7/350 2024-04-19 21:37 by 小狐_小胡
[找工作] 家乡二本高校/沿海传统私企，如何抉择？ 10+6 化学巷 2024-04-15 15/750 2024-04-19 20:43 by jiayouqf7
[硕博家园] 哎 +5 s150535912 2024-04-18 6/300 2024-04-19 20:33 by s150535912
[教师之家] 某系主任被同校教授举报包养情妇、长期嫖娼 +15 zju2000 2024-04-13 17/850 2024-04-19 18:17 by ZHONGWU_U
[考博] 24/25读博求博导 +4 宝23 2024-04-16 4/200 2024-04-19 17:18 by wangzhe_bs
[基金申请] 国家资助博士后BC档出校后资助的概率多大？ +3 卡卡罗特哦 2024-04-16 3/150 2024-04-18 12:58 by wolfgangHugh
[考研] 本科211 277分求调剂 +3 0854通信工程 2024-04-18 4/200 2024-04-18 09:11 by 花花平头哥
[无机/物化] H3PO2在130度附近分解为磷化氢（剧毒），怎样防护，求指点。 1(金币+5) +3 大力2010 2024-04-13 9/450 2024-04-17 21:55 by 大力2010
[考研] 工学322，英语六级求调剂 +3 domax 2024-04-17 10/500 2024-04-17 21:25 by 东方花园小万
[论文投稿] 可以打电话问编辑部是否可以先发录用通知吗 +7 双倍好运锦鲤 2024-04-14 10/500 2024-04-17 13:38 by cjzhu
[考博] 24计算机申博 +4 学无止境er 2024-04-13 6/300 2024-04-16 19:15 by 学无止境er
[考博] 25年申博求助，射频方向 +4 SherlockAH 2024-04-15 10/500 2024-04-16 18:06 by ljysudaee
[考研] 320求调剂 +5 陆志伟 2024-04-15 5/250 2024-04-16 11:11 by 19862091
[考研] 296求调剂 +3 Cclocomotive 2024-04-16 4/200 2024-04-16 10:04 by 19862091
[考研] 273求调剂 +5 Late婉安 2024-04-15 7/350 2024-04-15 13:01 by Late婉安
[考研] 专硕调剂招生 +3 电致发光 2024-04-15 4/200 2024-04-15 07:34 by ashorewmj
[考研] 287求调剂 +6 南沨 2024-04-14 6/300 2024-04-14 23:08 by lincunhui
[考研] 290求调剂 +3 杨yhr 2024-04-14 5/250 2024-04-14 21:50 by coco1981
信息提示
关闭
请填处理意见
关闭