版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

oO小柚子Oo

铁虫 (小有名气)

应助: 2 (幼儿园)
金币: 49.8
散金: 10
帖子: 64
在线: 22.8小时
虫号: 3721021
注册: 2015-03-09
专业: 信号理论与信号处理

[求助] 请问怎样改进三重for循环呐？O(∩_∩)O谢谢已有2人参与

for a=1:101
for b=1:101
      for c=1:101
         x=-0.8+0.016*(a-1);
         y=-0.8+0.016*(b-1);
         z=-0.8+0.016*(c-1);
         r=[x y z];
         sume=e0*exp(1i*dot(k0,r))+e1*exp(1i*dot(k1,r))+e2*exp(1i*dot(k2,r))+e3*exp(1i*dot(k3,r))+e4*exp(1i*dot(k4,r));
         intensity(a,b,c)=sume*sume';
这里可以怎么改进一下，然后就快了呢？

回复此楼

» 猜你喜欢

自荐读博已经有8人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复
自然科学基金委宣布启动申请书“瘦身提质”行动已经有4人回复
求个博导看看已经有18人回复

1楼 2016-11-19 21:11:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

mdl9966

铁虫 (初入文坛)

应助: 2 (幼儿园)
金币: 69.1
红花: 1
沙发: 1
帖子: 46
在线: 21.9小时
虫号: 2906415
注册: 2013-12-31
性别: GG
专业: 岩土力学与岩土工程

【答案】应助回帖

感谢参与，应助指数 +1

OpenMP
Cuda
OpenACC

赞一下

回复此楼

2楼2016-11-20 15:14:46

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

oO小柚子Oo

铁虫 (小有名气)

应助: 2 (幼儿园)
金币: 49.8
散金: 10
帖子: 64
在线: 22.8小时
虫号: 3721021
注册: 2015-03-09
专业: 信号理论与信号处理

引用回帖:

2楼: Originally posted by mdl9966 at 2016-11-20 15:14:46
OpenMP
Cuda
OpenACC

请问是什么意思呢？

发自小木虫Android客户端

赞一下

回复此楼

3楼2016-11-20 17:58:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

应助: 591 (博士)
贵宾: 0.03
金币: 18138.7
红花: 97
帖子: 8891
在线: 1485.4小时
虫号: 2400059
注册: 2013-04-04
专业: 计算机应用技术

【答案】应助回帖

感谢参与，应助指数 +1

CODE:

a=[1:101]';

b=[1:101]';

c=[1:101]';

x=-0.8+0.016*(a-1);

y=-0.8+0.016*(b-1);

z=-0.8+0.016*(c-1);

r=[x y z];

sume=e0*exp(1i*dot(k0,r))+e1*exp(1i*dot(k1,r))+e2*exp(1i*dot(k2,r))+e3*exp(1i*dot(k3,r))+e4*exp(1i*dot(k4,r));

intensity(a,b,c)=sume*sume';

赞一下(1人)

回复此楼

4楼2016-11-20 18:20:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

huab1984666

新虫 (著名写手)

应助: 9 (幼儿园)
金币: 241.5
散金: 2094
红花: 12
帖子: 1320
在线: 210.3小时
虫号: 516828
注册: 2008-03-03
性别: GG
专业: 胶体与界面化学

引用回帖:

4楼: Originally posted by FMStation at 2016-11-20 18:20:50
a=';
b=';
c=';
x=-0.8+0.016*(a-1);
y=-0.8+0.016*(b-1);
z=-0.8+0.016*(c-1);
r=;
sume=e0*exp(1i*dot(k0,r))+e1*exp(1i*dot(k1,r))+e2*exp(1i*dot(k2,r))+e3*exp(1i*dot(k3,r))+e4*exp(1i*dot(k4,r));
...

for i=1:1:101
x(:,i)=-0.8+0.016*(a(i)-1);
y(:,i)=-0.8+0.016*(b(i)-1);
z(:,i)=-0.8+0.016*(c(i)-1);
end
r=[x y z];
sume=e0*exp(1i*dot(k0,r))+e1*exp(1i*dot(k1,r))+e2*exp(1i*dot(k2,r))+e3*exp(1i*dot(k3,r))+e4*exp(1i*dot(k4,r));
intensity(a,b,c)=sume*sume';

赞一下

回复此楼

春风又绿江南岸，明月何时照我还。