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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】请高手看看这个小程序问题出在哪里了？
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 603  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

shenailin

木虫 (正式写手)

[交流] 【求助】请高手看看这个小程序问题出在哪里了？ 已有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)
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-04-24 15:51:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主
★
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中只保存了最后一次循环所赋的值。
赞一下(1人) 回复此楼
明德厚学、求是创新
2楼2010-04-24 17:11:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shenailin

木虫 (正式写手)

这就是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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lds588

铁杆木虫 (著名写手)

愚钝
引用回帖:
Originally posted by shenailin at 2010-04-24 19:00:41:
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         - ...

你给出的A一共5列，下面的x0是A的第六列，这怎么可能？
x0=A(:, 6);
y0=A(:, 7);
z0=A(:, 8);
赞一下 回复此楼
[url=http://weibo.com/u/2577306830?s=6uyXnP][img]http://service.t.sina.com.cn/widget/qmd/2577306830/fa2fdb00/1.png[/img][/url]
4楼2010-04-24 21:56:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shenailin

木虫 (正式写手)
对不起，我改动了一点，这里是
x0=A(:, 2);
y0=A(:, 3);
z0=A(:, 4);
其他的没粘过来
赞一下 回复此楼
5楼2010-04-24 22:13:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shenailin 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 305分求调剂 +3 哈_哈_哈_哈_哈 2026-04-04 5/250 2026-04-07 14:49 by 哈_哈_哈_哈_哈
[考研] 求调剂 +7 wwwwabcde 2026-04-07 7/350 2026-04-07 12:32 by dongzh2009
[考研] 347材料专硕求调剂 +10 zj8215216 2026-04-06 10/500 2026-04-07 11:37 by 我是小康
[考研] 362求调剂 +7 我要考大 2026-04-06 11/550 2026-04-07 09:32 by 我要考大
[考研] 0703化学 +9 goldtt 2026-04-02 11/550 2026-04-06 10:35 by 无际的草原
[考研] 0857大类环境工程B区求调剂 +3 龚禹铭 2026-04-05 3/150 2026-04-06 10:22 by 蓝云思雨
[考研] 一志愿南昌大学，085600，344分求调剂 +9 调剂上岸玘 2026-04-05 10/500 2026-04-06 09:30 by dongzh2009
[考研] 求调剂 +10 Hll胡 2026-04-04 10/500 2026-04-05 20:09 by nepu_uu
[考研] 男生，一志愿沪9生物学071000，初试308求调剂 +3 刘墨墨 2026-04-04 3/150 2026-04-05 08:26 by barlinike
[考研] 调剂 +8 熊二想上岸 2026-04-04 8/400 2026-04-05 05:27 by houyaoxu
[考研] 一志愿北京化工大学，初试成绩350求调剂 +9 沿岸?贝壳 2026-04-04 14/700 2026-04-05 01:09 by 沿岸?贝壳
[考研] 调剂 +11 JLLLLLLLLLL 2026-04-03 11/550 2026-04-04 22:21 by hemengdong
[考研] 363求调剂 +7 zh096 2026-04-04 7/350 2026-04-04 17:13 by dongzh2009
[考研] 322求调剂 +6 FZAC123 2026-04-03 6/300 2026-04-03 22:23 by 科研小专家
[考研] 求调剂 +4 15064154688 2026-04-03 5/250 2026-04-03 15:07 by zrongyan
[考研] 313求调剂 +3 ～微微凉～ 2026-04-03 3/150 2026-04-03 11:25 by 啵啵啵0119
[考研] 071000生物学调剂 +8 知昭蔓 2026-04-02 8/400 2026-04-03 10:36 by macy2011
[考研] 调剂 +3 osbbx 2026-04-02 3/150 2026-04-03 07:47 by cc8418
[考研] 293求调剂 +4 珂珂乐 2026-04-02 4/200 2026-04-02 20:10 by 6781022
[考研] 求调剂0703 +5 周嘉尧 2026-03-31 8/400 2026-04-01 20:32 by ltltkkk
信息提示
关闭
请填处理意见
关闭