24小时热门版块排行榜    

查看: 693  |  回复: 3

dingrui

金虫 (初入文坛)

[交流] 【求助】matlab程序求助,请高人帮助修改,能运行得不到结果~~~ 已有2人参与

[x,y,z]=meshgrid(0:1/100:2,0:1/100:2,0:1/100:2);
A=xlsread('C:\Users\aaa\Desktop\GC坐标,电荷.xls');
x0=A(:, 1);
y0=A(:, 2);
z0=A(:, 3);
q0=A(:, 5);
Vx=0;Vy=0;Vz=0;Fx=0;Fy=0;Fz=0;
for i=1:size(x0,1);
    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).*(x-x0(i))./((x-x0(i)).^3+(y-y0(i)).^3+(z-z0(i)).^3);
    Fy=Fy+q0(i).*(y-y0(i))./((x-x0(i)).^2+(y-y0(i)).^2+(z-z0(i)).^2);
    Fz=Fz+q0(i).*(z-z0(i))./((x-x0(i)).^2+(y-y0(i)).^2+(z-z0(i)).^2);
    V=Vx+Vy+Vz;F=Fx+Fy+Fz;
end
subplot(1,2,1)
slice(x,y,z,V,2,2,2);
subplot(1,2,2)
slice(x,y,z,F,2,2,2);
这是多电荷体系电场的电势能和受力
回复此楼
工作也是学习
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuyalei820

金虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
anyuezhiji(金币+1):给个红包,谢谢回帖交流 2010-04-21 13:24
把文件'C:\Users\aaa\Desktop\GC坐标传上来,不然看不出具体的错误。。。。
2楼2010-04-21 10:34:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuyalei820

金虫 (正式写手)

把文件'C:\Users\aaa\Desktop\GC坐标传上来,不然看不出具体的错误。。。。
3楼2010-04-21 10:34:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingrui

金虫 (初入文坛)

很感谢楼上兄弟 ,问题已经解决了,最近一直不在学校,没有及时给你送金币,真是不好意思啊~~~~
工作也是学习
4楼2010-05-06 09:20:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dingrui 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见