| 查看: 2062 | 回复: 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 |
回复此楼» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有129人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有17人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
此贴已删除
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有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
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
★
感谢参与，应助指数 +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
jerkwin
专家顾问 (正式写手)
-
专家经验: +14 - 计算强帖: 1
- 应助: 454 (硕士)
- 金币: 20699.1
- 散金: 148
- 红花: 81
- 帖子: 813
- 在线: 2648.3小时
- 虫号: 1023452
- 注册: 2010-05-19
- 专业: 理论和计算化学
- 管辖: 分子模拟
4楼2013-04-22 19:57:27
50