| 查看: 383 | 回复: 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 ] |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有11人回复
-
救命帖
已经有5人回复
-
限项规定
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有35人回复
-
最失望的一年
已经有18人回复
-
求推荐博导
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
疑惑？
已经有5人回复
信彼南山
木虫 (著名写手)
- 应助: 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