没有大错误，将中文标点改下就行了 
新建一个cjxxfun.m文件 function yp=cjxxfun(t,y,flag,q,m) yp=[y(4);y(5);y(6);q*(10*y(3)+1)*y(5)/m-q*((-5*sqrt(y(1).^2+y(2).^2))*cos(pi/4))*y(6)/m;q*((-5*sqrt(y(1).^2+y(2).^2))*sin(pi/4))*y(6)/m;-q*((-5*sqrt(y(1).^2+y(2).^2))*sin(pi/4))*y(5)/m]; 在command窗口输入下面代码结果就有了 clear; q=1.6e-19; m=1.6725e-27; V0=1e6;sita=pi/4; Vy=0;Vz=V0*cos(sita);Vx=V0*sin(sita); c=[0;1;0;Vx;Vy;Vz];%给定初始量 [t,y]=ode45('cjxxfun',[0:1e-7:1.5e-6],c,[],q,m);%解微分方程组 plot3(y(:,1),y(:,2),y(:,3),'linewidth',2);%绘制三维粒子运动轨迹 xlabel('z');ylabel('y');zlabel('x'); grid on;title('磁镜现象') figure subplot(2,1,1); plot(y(:,1),y(:,2),'linewidth',2);%从yoz面看到的轨迹 xlabel('z');ylabel('y'); grid on;title('yoz平面') subplot(2,1,2); plot(y(:,1),y(:,3),'linewidth',2);%从xoz面观察的轨迹 xlabel('z');ylabel('x'); title('xoz平面');grid on; - 追问
-
谢谢你的回答，我不太会用matlab，那个 cjxxfun.m放哪啊，我复制你的程序测试还是不成功
- 回答
-
放在当前工作目录下，current folder目录下，然后直接复制粘贴我贴出来的那些代码
- 追问
-
不好意思啊，当前工作目录是什么目录啊，我用matlab都是直接输程序然后运行的，这个文件夹是从哪来的
- 回答
-
能看到current folder窗口吗？直接把那个m文件copy进去就行了
- 提问者评价
谢谢!
|
