开发工具:
文件大小: 1kb
下载次数: 0
上传时间: 2014-12-30
详细说明: t=0:pi/50:2.5*pi; electron_x=30*cos(t); electron_y=30*sin(t); electron_z=0*t; n=numel(t); nucleus=2; electron=-1; K0=8.98755e9; [X,Y,Z]=meshgrid(-52:2:52,-39:2:39,-39:2:39); [sx,sy,sz]=meshgrid(linspace(20,-20,5),linspace(12,-12,4),... linspace(12,-12,4)); [x0,y0,z0]=sphere(40); for k=1:n; clf r1=sqrt(X.^2+Y.^2+Z.^2); r2=sqrt((X-electron_x(k)).^2+(Y-electron_y(k)).^2+(Z-electron_z(k)).^2); F=K0*(nucleus./r1+electron./r2); plot3(electron_x,electron_y,electron_z,'y') axis([-52 52 -39 39 -39 39]) hold on g1=surfl(2*x0, 2*y0,2*z0); set(g1,'facecolor','interp','edgecolor','none'); g2 = surfl(x0+electron_x(k),y0+electron_y(k),z0+electron_z(k)); set(g2,'facecolor','interp','edgecolor','none') colormap(cool) [u,v,w]=gradient(F); h2=streamline(stream3(X,Y,Z,-u,-v,-w,sx,sy,sz)); set(h2,'color','g'); h3=streamline(stream3(X,Y,Z,u,v,w,sx,sy,sz)); set(h3,'color','g'); quiver3(X,Y,Z,-u,-v,-w,'g') view(-18,15) axis equal vis3d off lighting gouraud material shiny zoom(1.5) set(gcf,'color','k') hold off pause(0.0001) end ...展开收缩
(系统自动生成,下载前可以参看下载内容)
下载文件列表
相关说明
- 本站资源为会员上传分享交流与学习,如有侵犯您的权益,请联系我们删除.
- 本站是交换下载平台,提供交流渠道,下载内容来自于网络,除下载问题外,其它问题请自行百度。
- 本站已设置防盗链,请勿用迅雷、QQ旋风等多线程下载软件下载资源,下载后用WinRAR最新版进行解压.
- 如果您发现内容无法下载,请稍后再次尝试;或者到消费记录里找到下载记录反馈给我们.
- 下载后发现下载的内容跟说明不相乎,请到消费记录里找到下载记录反馈给我们,经确认后退回积分.
- 如下载前有疑问,可以通过点击"提供者"的名字,查看对方的联系方式,联系对方咨询.