24小时热门版块排行榜    

查看: 2034  |  回复: 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的回帖

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的回帖
查看全部 4 个回答

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂院校信息 +4 CX 330 2026-03-21 4/200 2026-03-21 23:48 by ms629
[考研] 一志愿北京化工大学070300 学硕336求调剂 +3 vv迷 2026-03-21 4/200 2026-03-21 23:05 by f19980501
[考研] 材料工程专硕 348分求调剂 +3 冬辞. 2026-03-17 5/250 2026-03-21 18:47 by 学员8dgXkO
[考研] 生物学一志愿985,分数349求调剂 +3 zxts12 2026-03-21 3/150 2026-03-21 16:34 by 33来了真来了
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
信息提示
请填处理意见