| 查看: 742 | 回复: 3 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
dingzemin银虫 (正式写手)
|
[交流]
【求助】matlab求解方程，方程里面含有psi函数（ploygamma函数），不好解啊
|
||
|
求解函数： y=1-x/5-1.44*{0.25/x-0.05[-psi(1)+0.5*psi(1-0.1*x)+ 0.5*psi(1+0.1*x)]} 其中psi是ploygamma函数，在matlab里面可以直接调用。y是给定的值，其变化范围是[-0.5, 1]， x是要求出来的变量。对应每一个y，存在两个x值。我不知道怎么求解这两个x的值。调用solve函数也解不出来。 我试着按照常规思路求解了一下：先把y取成一个定值0.1，调用solve求解，如下： y=0.1 f=solve('1-x/5-1.44*(0.25/x-0.05*(-psi(1)+1/2*psi(1-x/10)+1/2*psi(1+x/10)))-y', 'x') 结果显示是：Warning: Explicit solution could not be found. 看到网上有人建议将这个问题看成是一个非线性最优化问题，即求得一个x，使得方程左端函数值与右端常数值的平方和最小。函数如下： y=0.1 f=@(x)norm(1-x/5-1.44*(0.25/x-0.05*(-psi(1)+1/2*Psi(1-x/10)+1/2*Psi(1+x/10)))-y)^2; [x,err]=fminsearch(f,1e-5) 但是这种方法得出的x值只有一个值，实际上x存在两个值。改怎么解决呢？ 请大家帮忙啊！！ [ Last edited by dingzemin on 2010-1-17 at 20:11 ] |
回复此楼» 猜你喜欢
关于代码变化问题，想知道的进来
已经有20人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有11人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★ ★ ★ ★ ★
sunxiao(金币+2,VIP+0):谢谢详细解释，欢迎常来 1-18 04:35
dingzemin(金币+3,VIP+0): 1-18 13:06
sunxiao(金币+2,VIP+0):谢谢详细解释，欢迎常来 1-18 04:35
dingzemin(金币+3,VIP+0): 1-18 13:06
|
通过图片可以知道这个方程x=[0 10],y=[-inf 0.4616].最高点为(1.3 0.4616)(放大图片可以看出来),当y<0.4614时,x有两个值,分别在[0 1.3]及[1.3 10]区间上。 当y=0.1时,可以采用以下求出两个解 >> f=fsolve('1-x/5-1.44*(0.25/x-0.05*(-psi(1)+1/2*psi(1-x/10)+1/2*psi(1+x/10)))-0.1', 0.5) Optimization terminated: first-order optimality is less than options.TolFun. f = 0.443865675644518 >> f=fsolve('1-x/5-1.44*(0.25/x-0.05*(-psi(1)+1/2*psi(1-x/10)+1/2*psi(1+x/10)))-0.1', 3) Optimization terminated: first-order optimality is less than options.TolFun. f = 3.967250801226048 |
3楼2010-01-17 22:40:51
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
2楼2010-01-17 22:36:11
dingzemin
银虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 482.8
- 散金: 897
- 帖子: 373
- 在线: 69.5小时
- 虫号: 793813
- 注册: 2009-06-13
- 专业: 冶金物理化学与冶金原理
4楼2010-01-18 13:05:07
100