24小时热门版块排行榜

>论坛更新日志 (2479)
>导师招生 (254)
>招聘信息布告栏 (226)
>材料综合 (100)
>无机非金属 (90)
>药学 (30)
>虫友互识 (27)
>硕博家园 (22)
>休闲灌水 (22)
>学术会议 (16)
>论文道贺祈福 (15)
>找工作 (14)
>考研 (12)
>论文投稿 (12)
>考博 (11)
>基金申请 (10)

返回列表

karry2karry

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 102.6
散金: 20
帖子: 18
在线: 17.3小时
虫号: 3795233
注册: 2015-04-08
性别: MM
专业: 计算机科学的基础理论

[求助] 请问这个问题，用MATLAB代码怎么实现？已有1人参与

请问这个问题用MATLAB怎么实现

@月只蓝 @beefly 发自小木虫Android客户端

回复此楼

» 猜你喜欢

求助中医院临床中医中标的青年标书学习，谢谢已经有6人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有259人回复
物理学I论文润色/翻译怎么收费? 已经有158人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有47人回复
25届硕士申请，材料计算方向已经有2人回复
课题组招1-2名硕士已经有1人回复
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上已经有0人回复
斯坦福大学与Elsevier世界前2%学者新排名查询已经有1人回复
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有19人回复
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计已经有1人回复
apl投稿状态已经有0人回复

1楼 2016-12-27 13:15:59

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

oldoldstone

木虫 (小有名气)

应助: 2 (幼儿园)
金币: 2113
散金: 466
红花: 4
帖子: 166
在线: 297.2小时
虫号: 414842
注册: 2007-06-28
专业: 生物、医学光子学

【答案】应助回帖

感谢参与，应助指数 +1

我去，matlab数值法最大只能算出170的阶乘，符号法最大也只能算出999的阶乘，应该需要进一步近似的。。。
写了一个符号法的，代码放在同一个文件pp.m
%pp.m
function pp
N=0:999;
A=zeros(length(N),1);
for n=N
A(n+1)=fun(n);
end
plot(N,A)

function result=fun(N)
syms n m lambda c
n=sym(num2str(N));
lambda=sym('3000');
f=lambda^m/factorial(m); %求和函数
S=lambda^n/factorial(n)/symsum(f,m,0,n); %分子的3000我也是用lambda代替了
result=vpa(S,8);%转化为数值

貌似运行需要几分钟，出来几乎是一条直线。