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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » matlab解答计算问题
261/1返回列表
查看: 1292  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 05:38 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-17 05:14 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-17 04:14 by k06BKNblrGRH
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-17 03:52 by k06BKNblrGRH
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 23jxep3nCNZb 2026-09-14 5/250 2026-09-17 03:50 by k06BKNblrGRH
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-17 02:38 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-13 6/300 2026-09-17 02:14 by Tql5LQhh5rLK
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-13 8/400 2026-09-17 02:04 by Tql5LQhh5rLK
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-17 01:26 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-16 17:27 by i3KBRD19ERkG
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 15:48 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-16 06:51 by LH5NkK5Ud4bh
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:36 by LH5NkK5Ud4bh
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:15 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:50 by DgNGHc3h5tPl
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 02:14 by DgNGHc3h5tPl
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
信息提示
关闭
请填处理意见
关闭