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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】请高手看看这个小程序问题出在哪里了？
51/1返回列表
查看: 592  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学319求调剂 +4 锦鲤0909 2026-03-17 4/200 2026-03-17 18:21 by 重科小霸王
[考研] 085601材料工程专硕求调剂 +4 慕寒mio 2026-03-16 4/200 2026-03-17 17:17 by ruiyingmiao
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +5 yangfz 2026-03-17 5/250 2026-03-17 17:01 by ruiyingmiao
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[教师之家] 焦虑 +7 水冰月月野兔 2026-03-13 9/450 2026-03-16 10:00 by Quakerbird
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 求调剂（材料与化工327） +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
信息提示
关闭
请填处理意见
关闭