| 查看: 407 | 回复: 1 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】Matlab关于级数求和问题
|
|||
|
要计算一个关于K,L,kama,alpha从负无穷到正无穷的求和(原公式见附件)，并且分别独立变化，由于这四个同时都在贝塞尔函数中不能分离，按照以下程序运算耗时较长，还请各位帮忙。以下为求和的部分程序！ alpha=-10:10; for m=-10:1:10 for n=-10:1:10 for k=-10:10 for kama=-10:10 for l=-10:10 fai(m+11,n+11)=fai(m+11,n+11)+a0.*(sum((i.^(m+n+m1+n1+k+l).*(-1).^(alpha.*miu+kama.*you+k).*besselj(k,2.*c.*Z).*besselj(l,2.*c.*Z))... .*((-1).^(-m1-k-l).*besselj(m-m1-k-l+2.*alpha.*(miu+2),2.*a.*Z)+besselj(m+m1+k+l+2+2.*alpha.*(miu+2),2.*a.*Z))... .*((-1).^(-n1-k+l).*besselj(n-n1-k+l+2.*kama.*(you+2),2.*b.*Z)+besselj(n+n1+k-l+2+2.*kama.*(you+2),2.*b.*Z)))); end end end end end [ Last edited by lxy6522 on 2009-10-27 at 10:09 ] |
回复此楼» 猜你喜欢
好奇怪的filecode
已经有5人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有21人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有10人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有3人回复
-
国自然结果
已经有5人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有12人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有8人回复
-
据悉今年马上要出结果了
已经有8人回复
-
关于代码变化问题，想知道的进来
已经有23人回复
-
FileCode能看出啥？
已经有6人回复
2楼2009-10-26 20:31:19
20