24小时热门版块排行榜    

查看: 2210  |  回复: 3

lhf530001541

新虫 (初入文坛)

[求助] matlab求隐函数解析解

函数方程如下:
x=a-b*(exp((c*x+y)/(e*f))-1)-(c*x+y)/d。除x和y,其他的都是常量,想把x用符号表示出来。事实上x和y的关系图如图:纵轴为x,横轴为y。

有尝试着在matlab里用solve计算,得到的结果是:
Warning: Explicit solution could not be found.
> In solve at 83

ans =

[ empty sym ]
如果把方程后面小部分去掉,简化为:‘x-a+b*(exp((c*x+y)/(e*f))-1’,倒是能得到结果:ans =

a + b - (e*f*lambertw(0, (b*c*exp(y/(e*f))*exp((c*(a + b))/(e*f)))/(e*f)))/c。
不会就多了后面那一小点就求不出解析解了吧?~
求各位大牛帮帮忙~谢过了先!

横轴y纵轴x.PNG
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jerkwin

专家顾问 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
csgt0: 金币+1, 谢谢 2013-04-22 17:57:30
你把它化简一下,能得到如下形式的x的方程:
x+b*exp(a*x) = c
然后利用Lambert W function就可以把x表示出来了

» 本帖已获得的红花(最新10朵)

2楼2013-04-21 22:34:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lhf530001541

新虫 (初入文坛)

送红花一朵
引用回帖:
2楼: Originally posted by jerkwin at 2013-04-21 22:34:18
你把它化简一下,能得到如下形式的x的方程:
x+b*exp(a*x) = c
然后利用Lambert W function就可以把x表示出来了

谢谢哈,我试过了,还是不行。
用solve>>x+b*exp(a*x)-c得不到结果,
改成solve('x-exp(c*x+d)/a+b/a',x),也就是"a*x+b=exp(c*x+d)"
得到了:- lambertw(0, -(c*exp(d))/(a*exp((b*c)/a)))/c - b/a
把我的公式改成"a*x+b=exp(c*x+d)"形式后,得到:
x=-(a*d+b*d-y)/(c+d)-e*f/c*lambertw(0,-c*c*b*d/(e*f*(c+d))*exp((d*y+c*a*d+c*b*d)/(e*f*(c+d)))/(e*f)+(a*d+b*d-y)*c/((e*f)*(c+d)));
可是代入数字:y=0:1:40;
a=3.5;
b=0.0000000001;
c=3.5;
d=700;
e=66;
f=0.0257;
得到的x是虚数,这是怎么回事呀?lambertw(0,   )不应该都是实数吗?
3楼2013-04-22 15:53:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jerkwin

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
csgt0: 金币+1, 谢谢 2013-04-25 15:54:38
lhf530001541: 金币+10, 有帮助, 谢谢帮助哈,这个行不通,还是得想别的办法~ 2013-04-26 09:40:21
去读读Lambert W function的性质,你就知道了。Lambert W(x)的定义域为[1/e, inf],不是你随便试试就能对的。
4楼2013-04-22 19:57:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lhf530001541 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金啊基金 +3 longfie172 2026-08-20 3/150 2026-08-21 00:04 by wuchongjun
[基金申请] 让我中一个面上吧! +8 大萍1987 2026-08-20 10/500 2026-08-20 21:48 by 大萍1987
[基金申请] 言之凿凿,明天官网会通知放榜? +15 医学老男孩 2026-08-20 17/850 2026-08-20 20:35 by CaoLG
[基金申请] 时间戳又变了 +12 wuchongjun 2026-08-20 17/850 2026-08-20 20:05 by spzhao
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
信息提示
请填处理意见