24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1146  |  回复: 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露珠 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 河北省自然科学基金 +5 Peterchao 2026-05-18 8/400 2026-05-24 11:58 by 晓晓爱翠翠
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[教师之家] 论文撤稿了 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:06 by Equinoxhua
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路! +4 zju2000 2026-05-21 5/250 2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 08:01 by 9ps9vgkqva
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 pmo95bazuy 2026-05-23 7/350 2026-05-24 06:35 by fpo5ljpv91
[基金申请] 揭秘青基评审内幕:几个A才能顺利中标 +3 国自然国社科中 2026-05-23 4/200 2026-05-23 15:37 by 2000zf36392
[基金申请] 青B发送上会通知了吗 +5 chemBioBro 2026-05-22 7/350 2026-05-23 12:35 by zhuifengzhy
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[论文投稿] 投稿求助,期刊 +4 希冀,有书读 2026-05-20 8/400 2026-05-22 10:16 by 希冀,有书读
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +14 1234567wang 2026-05-17 16/800 2026-05-21 17:58 by 脆脆的饼干
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[基金申请] 评审有感 +15 popular289 2026-05-18 26/1300 2026-05-21 10:35 by 西葫芦炒鸡蛋
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
信息提示
请填处理意见