24小时热门版块排行榜

返回列表

karry2karry

新虫 (初入文坛)

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

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

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

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

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有90人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有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);%转化为数值

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