应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » matlab解答计算问题
261/1返回列表
查看: 1324  |  回复: 25
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 猜你喜欢
1楼 2017-12-16 00:33:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chunguangli

新虫 (小有名气)
★
Math露珠(金币+1): 谢谢参与
不要用符号运算

发自小木虫IOS客户端
回复此楼

» 本帖已获得的红花（最新10朵）

Math露珠
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的回帖

xiegangmai

版主 (职业作家)
★
小木虫: 金币+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朵）

Math露珠
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的回帖

陈氏帝国

至尊木虫 (知名作家)
我来看一下。

发自小木虫IOS客户端
回复此楼
26楼2018-01-10 21:49:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
syhorchid2楼
2017-12-16 00:39   回复  
Math露珠(金币+1): 谢谢参与
186303203343楼
2017-12-16 00:58   回复  
Math露珠(金币+1): 谢谢参与
发自小木虫IOS客户端
xdusgmt2324楼
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露珠 的主题更新
261/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭