| 查看: 624 | 回复: 2 | |||
karry2karry新虫 (初入文坛)
|
[求助]
请问这个问题，用MATLAB代码怎么实现？ 已有1人参与
|
请问这个问题用MATLAB怎么实现 @月只蓝 @beefly 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有232人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有112人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有69人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复
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
