24小时热门版块排行榜    

查看: 1070  |  回复: 25

Math露珠

木虫 (知名作家)


[交流] matlab解答计算问题

求助各位大神,用matlab计算时候。
为什么我正常计算psi(x)-log(x)==-0.1302这个表达式得到的结果x是一个复数。
我知道结果,4带入,也同样满足这个表达式呢?不是矛盾吗?
有两个解的话,matlab也会把所有的解都给出啊。谢谢

>> syms x
>> solve(psi(x)-log(x)==-0.1302,x)

ans =

- 226.75334615375301702825277514765 + 0.61433808584738080815547301292808*i

>> psi(4)-log(4)

ans =

   -0.1302

[ Last edited by 月只蓝 on 2017-12-16 at 15:29 ]
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chunguangli

新虫 (小有名气)



Math露珠(金币+1): 谢谢参与
不要用符号运算

发自小木虫IOS客户端

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

16楼2017-12-16 04:01:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Math露珠

木虫 (知名作家)


送红花一朵
引用回帖:
16楼: Originally posted by chunguangli at 2017-12-16 04:01:31
不要用符号运算

大神,那如何计算?

发自小木虫Android客户端
20楼2017-12-16 05:14:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
20楼: Originally posted by Math露珠 at 2017-12-16 05:14:13
大神,那如何计算?
...

用fsolve或fzero函数
23楼2017-12-16 11:50:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫: 金币+0.5, 给个红包,谢谢回帖
CODE:
format long
[x fval]=fsolve(@(x)psi(x)-log(x)+0.1302,1)

计算结果:
x =

   3.999310778077781

残差fval =

   -8.377323362784850e-009

而  psi(4)-log(4)+0.1302 =2.33e-005 远大于 8.378e-009(前者比后者大了上千倍)
可见,该方程的解取3.99931078,比取4合适。

MALTAB对函数psi(x)明确要求x必须是非负的实数:
Y = PSI(X) evaluates the psi function for each element of X.
    X must be real and nonnegative.  SIZE(Y) is the same as SIZE(X).
    The psi function, also know as the digamma function, is the logarithmic
    derivative of the gamma function:

       psi(x) = digamma(x) = d(log(gamma(x)))/dx = (d(gamma(x))/dx)/gamma(x).

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

24楼2017-12-16 12:07:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Math露珠

木虫 (知名作家)


送红花一朵
引用回帖:
24楼: Originally posted by 月只蓝 at 2017-12-16 12:07:33
format long
=fsolve(@(x)psi(x)-log(x)+0.1302,1)


计算结果:
x =

   3.999310778077781

残差fval =

   -8.377323362784850e-009

而  psi(4)-log(4)+0.1302 =2.33e-005 远大于 8.378e-009(前 ...

跪谢

发自小木虫Android客户端
25楼2017-12-16 13:19:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

陈氏帝国

至尊木虫 (职业作家)


26楼2018-01-10 21:49:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
syhorchid2楼
2017-12-16 00:39   回复  
Math露珠(金币+1): 谢谢参与
2017-12-16 00:58   回复  
Math露珠(金币+1): 谢谢参与
发自小木虫IOS客户端
2017-12-16 01:10   回复  
Math露珠(金币+1): 谢谢参与
gepwc955楼
2017-12-16 01:31   回复  
Math露珠(金币+1): 谢谢参与
rakzdw3856楼
2017-12-16 02:05   回复  
Math露珠(金币+1): 谢谢参与
月只蓝: 金币-200, 屏蔽内容, 违规存档, 多次恶意灌水。 2017-12-16 12:26:41
htufuje717楼
2017-12-16 02:19   回复  
Math露珠(金币+1): 谢谢参与
pdmmd078楼
2017-12-16 02:25   回复  
Math露珠(金币+1): 谢谢参与
wttmf509楼
2017-12-16 02:35   回复  
Math露珠(金币+1): 谢谢参与
Jasonluo10楼
2017-12-16 02:38   回复  
Math露珠(金币+1): 谢谢参与
发自小木虫IOS客户端
cjxjhf9211楼
2017-12-16 03:08   回复  
Math露珠(金币+1): 谢谢参与
nnmmjt9812楼
2017-12-16 03:45   回复  
Math露珠(金币+1): 谢谢参与
月只蓝: 金币-200, 屏蔽内容, 违规存档, 多次恶意灌水。 2017-12-16 12:32:14
月只蓝: 金币-200 2017-12-16 12:32:25
月只蓝: 金币-200, 多次恶意灌水。 2017-12-16 12:32:36
月只蓝: 金币-200, 多次恶意灌水。 2017-12-16 12:32:43
月只蓝: 金币-200 2017-12-16 12:32:51
月只蓝: 金币-200 2017-12-16 12:32:56
月只蓝: 金币-200 2017-12-16 12:33:03
月只蓝: 2017-12-16 12:33:10
月只蓝: 金币-200 2017-12-16 12:33:19
月只蓝: 金币-200 2017-12-16 12:33:32
bxhym52513楼
2017-12-16 03:52   回复  
Math露珠(金币+1): 谢谢参与
遗嗜qi14楼
2017-12-16 03:57   回复  
Math露珠(金币+1): 谢谢参与
pvaft37215楼
2017-12-16 03:59   回复  
Math露珠(金币+1): 谢谢参与
泌僮rg17楼
2017-12-16 04:03   回复  
Math露珠(金币+1): 谢谢参与
futkyyb20118楼
2017-12-16 04:17   回复  
Math露珠(金币+1): 谢谢参与
erxcg5619楼
2017-12-16 04:36   回复  
Math露珠(金币+1): 谢谢参与
祝福
ktnka16621楼
2017-12-16 05:23   回复  
Math露珠(金币+1): 谢谢参与
ubjkux8122楼
2017-12-16 05:24   回复  
Math露珠(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 Math露珠 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见