| 查看: 398 | 回复: 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 ] |
回复此楼» 猜你喜欢
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？
已经有8人回复
-
面上可以超过30页吧？
已经有11人回复
-
版面费该交吗
已经有15人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有18人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有10人回复
-
什么是人一生最重要的？
已经有4人回复
信彼南山
木虫 (著名写手)
- 应助: 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