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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】请高手看看这个小程序问题出在哪里了？
51/1返回列表
查看: 695  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于代码变化问题，想知道的进来 +7 且听虎啸 2026-08-07 10/500 2026-08-07 20:43 by 医学老男孩
[基金申请] 基金中了 +13 laoda193707 2026-08-06 13/650 2026-08-07 20:33 by tchen18
[基金申请] 关于filecode +3 布布和一二 2026-08-07 6/300 2026-08-07 17:51 by 医学老男孩
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[论文投稿] 十年后又回来了，论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +23 低垂的野花 2026-07-31 31/1550 2026-08-07 12:22 by speculiar
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的？有收到消息的吗？ +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐，不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
关闭
请填处理意见
关闭