| 查看: 507 | 回复: 2 | ||
kezl铜虫 (正式写手)
|
[求助]
求高人帮我化简优化下面的循环，以减少运行时间，万分感谢！
|
|
|
求高人帮我化简优化下面的循环，以减少运行时间，万分感谢！ 以下Matlab程序以实现对函数yy(x,p,t,y) 中的变量y的积分,最后得到函数aa(x,p,t) clear all % p=-5:0.1:5; % x=-5:0.1:5; % t=0:0.02:1 % y=-5:0.1:5; aa=[]; for n=0:100; p=0.1*(n-50); for nn=0:50 t=0.02*nn for m=0:100; x=0.1*(m-50); gg=0; for mm=0:100; y=0.1*(mm-50); y1=1/sqrt(2)*(1-exp(i*0.5*t)*0.5328*(0.5+tanh(x+y))).*((0.3985*cosh(x+y).*exp(-0.5*0.25*((cosh(x+y)).^2+0.5* (x+y)+0.5*0.5*sinh(2*(x+y)))))); y2=1/sqrt(2)*(1-exp(-i*0.5*t)*0.5328*(0.5+tanh(x+y))).*((0.3985*cosh(x-y).*exp(-0.5*0.25*((cosh(x-y)).^2+0.5*(x-y)+0.5*0.5*sinh(2*(x-y)))))); yy=2*p*exp(2*i*p.*y).*y1.*y2; gg=gg+yy*0.1; end aa(n+1,nn+1,m+1)=gg end end end [ Last edited by kezl on 2013-1-20 at 20:09 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有21人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有21人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
2楼2013-01-20 21:13:41
kezl
铜虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 117.9
- 散金: 10
- 帖子: 450
- 在线: 39.9小时
- 虫号: 1940844
- 注册: 2012-08-16
- 性别: GG
- 专业: 基础物理学
3楼2013-01-22 21:37:11
20