| 查看: 537 | 回复: 2 | |||
karry2karry新虫 (初入文坛)
|
[求助]
请问这个问题，用MATLAB代码怎么实现？ 已有1人参与
|
请问这个问题用MATLAB怎么实现 @月只蓝 @beefly 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有229人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
-
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有1人回复
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);%转化为数值 貌似运行需要几分钟，出来几乎是一条直线。 |
2楼2016-12-27 15:17:53
3楼2016-12-27 18:25:06
