| 查看: 557 | 回复: 4 | |||
[交流]
【求助】请高手看看这个小程序问题出在哪里了？ 已有3人参与
|
|
A=xlsread('D:\proton trasfer\1T.xls'); x0=A(:, 6); y0=A(:, 7); z0=A(:, 8); q0=A(:, 3); Vx=0;Vy=0;Vz=0;Fx=0;Fy=0;Fz=0;V1x=0;V1y=0;V1z=0;V2x=0;V2y=0;V2z=0;V3x=0;V3y=0;V3z=0; for i=1:size(x0,1);j=26;x=x0(j);y=y0(j);z=z0(j);q=q0(j); Vx=Vx+q0(i).*(x-x0(i))./((x-x0(i)).^2+(y-y0(i)).^2+(z-z0(i)).^2); Vy=Vy+q0(i).*(y-y0(i))./((x-x0(i)).^2+(y-y0(i)).^2+(z-z0(i)).^2); Vz=Vz+q0(i).*(z-z0(i))./((x-x0(i)).^2+(y-y0(i)).^2+(z-z0(i)).^2); Fx=Fx+q0(i).*q.*(x-x0(i))./((x-x0(i)).^2+(y-y0(i)).^2+(z-z0(i)).^2); Fy=Fy+q0(i).*q.*(y-y0(i))./((x-x0(i)).^2+(y-y0(i)).^2+(z-z0(i)).^2); Fz=Fz+q0(i).*q.*(z-z0(i))./((x-x0(i)).^2+(y-y0(i)).^2+(z-z0(i)).^2); V(j)=(Vx.^2+Vy.^2+Vz.^2).^(1/2); F(j)=(Fx.^2+Fy.^2+Fz.^2).^(1/2); end V(j) F(j) |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
★
anyuezhiji(金币+1):谢谢回帖交流 2010-04-24 20:40
shenailin(金币+4): 2010-04-24 22:13
anyuezhiji(金币+1):谢谢回帖交流 2010-04-24 20:40
shenailin(金币+4): 2010-04-24 22:13
|
没你的xls文件。 1、程序中j=26;x=x0(j);y=y0(j);z=z0(j);q=q0(j);该放在循环前，不用在每个循环中赋值计算。 2、Vx=Vx+q0(i).*(x-x0(i))./((x-x0(i)).^2+(y-y0(i)).^2+(z-z0(i)).^2);中，不存在数组之间的对应元素运算，所以运算符前的“.”号没必要。 3、for循环中， V(j)=(Vx.^2+Vy.^2+Vz.^2).^(1/2); F(j)=(Fx.^2+Fy.^2+Fz.^2).^(1/2); 每次循环都给V(j),F(j)赋值，而j是不变的，循环结束后V，F中只保存了最后一次循环所赋的值。 |
2楼2010-04-24 17:11:15
这就是Excel文件的内容，麻烦你在给看看还有什么问题！
|
1 0.4397 -1.74856 -2.80728 0 2 -0.5663 -2.14786 -4.14676 0 3 0.4745 0.99366 -4.50379 0 4 0.7356 -1.85968 -0.57935 0 5 -0.6270 -2.54142 -1.74274 0 6 0.7081 0.38903 -1.60605 0 7 0.2450 -1.06427 -6.02435 0 8 0.4328 -3.53884 0.52009 0 9 0.4517 -2.04553 1.49752 0 10 -0.0062 -0.35467 -2.83915 0 11 0.4693 -3.11245 -4.46919 0 12 0.4607 0 0.43235 0 13 -0.6776 1.62009 -1.49525 0 14 0.3414 -1.05544 -4.94205 0 15 -0.7934 -2.52735 0.57338 0 16 -0.6823 -0.4764 -0.51248 0 17 -0.5175 0.03228 -4.17066 0 18 -0.7884 3.00583 1.08143 0 19 -0.3505 2.81789 3.49099 0 20 0.1284 1.96486 4.55372 0 21 0.5274 2.22501 2.17612 0 22 0.8946 0.05087 3.06935 0 23 -0.6403 0.6175 4.34858 0 24 -0.7239 -1.18897 2.96936 0 25 -0.7167 0.89245 1.99694 0 26 0.4431 2.58741 0.14528 0 27 0.4175 4.0125 1.17373 0 28 0.2395 3.89369 3.63269 0 29 0.2320 2.30262 5.5872 0 30 0.4487 -0.03881 5.12388 0 |
3楼2010-04-24 19:00:41
lds588
铁杆木虫 (著名写手)
愚钝
- 应助: 15 (小学生)
- 金币: 10949.6
- 散金: 2018
- 红花: 5
- 帖子: 1444
- 在线: 206.3小时
- 虫号: 890720
- 注册: 2009-11-01
- 专业: 凝聚态物性I:结构、力学和
4楼2010-04-24 21:56:02
5楼2010-04-24 22:13:32