| 查看: 2054 | 回复: 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 |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有94人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有11人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有12人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助 matlab隐函数画图问题
已经有4人回复
- Matlab中如何自己定义一个新函数，然后计算时可以调用？
已经有12人回复
- 求助MATLAB隐函数等
已经有7人回复
- 求matlab画一个函数的图
已经有6人回复
- matlab的函数调用问题
已经有3人回复
- MATLAB中产生高斯白噪声，涉及到awgn和wgn函数
已经有3人回复
- 问个关于matlab里find函数的问题！
已经有4人回复
- 求助matlab隐函数参数拟合
已经有16人回复
- 【分享】Matlab常用函数及常用数学函数【已搜索无重复】
已经有80人回复
- matlab做函数图象的问题
已经有6人回复
- 【求助】matlab求解隐函数时如何设置合适的自变量步长？
已经有7人回复
- 【求助】求解隐函数
已经有6人回复
- 【分享】MATLAB常用函数参考
已经有36人回复
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 计算强帖: 1
- 应助: 454 (硕士)
- 金币: 20699.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.3小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
4楼2013-04-22 19:57:27
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 计算强帖: 1
- 应助: 454 (硕士)
- 金币: 20699.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.3小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
★
感谢参与，应助指数 +1
csgt0: 金币+1, 谢谢 2013-04-22 17:57:30
感谢参与，应助指数 +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
lhf530001541
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 969
- 红花: 1
- 帖子: 29
- 在线: 86.9小时
- 虫号: 1722806
- 注册: 2012-03-28
- 专业: 半导体光电子器件
送红花一朵 |
谢谢哈，我试过了，还是不行。 用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
