| 查看: 463 | 回复: 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 ] |
回复此楼» 猜你喜欢
听说今天filecode变了
已经有40人回复
-
filecode
已经有16人回复
-
基金中了
已经有9人回复
-
固定端突然变了，今天
已经有8人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
-
影响面上的因素
已经有11人回复
信彼南山
木虫 (著名写手)
- 应助: 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
50