24小时热门版块排行榜

返回列表

karry2karry

新虫 (初入文坛)

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

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

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

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

回复此楼

» 猜你喜欢

中国科学院物理研究所谌志国研究员团队招收2027年博士研究生已经有6人回复
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究已经有12人回复
物理学I论文润色/翻译怎么收费? 已经有266人回复
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有1人回复
如何从铁电相到顺电相。已经有1人回复
各位大佬，求一份最新Wien2k版本作为学习用途已经有0人回复
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后已经有24人回复
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向已经有22人回复

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);%转化为数值

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