完善资料让更多小伙伴认识你,还能领取20积分哦, 立即完善>
|
%用matlab数值模拟的方法模拟带电粒子在恒定磁场中的螺旋运动 %function lx global q m B E q=1.6e-2; m=0.02; B=[1;1;0]; E=[1;0;1]; figure strd{1}='Eneq 0,Bneq 0'; strd{2}='E=0,Bneq 0'; strd{3}='Eneq 0, B=0'; for i=1:3 [t,w]=ode23(@ddl***un,[0:0.01:20],[0,0.01,0,6,0,0.01],[],q,m,B(i),E(i)); axes ('unit','normalized','position', [0.045+(i-1)*0.35 0.062 0.2786 0.6583 ]); plot3(w(:,1),w(:,3),w(:,5),'linewidth',2); grid on title(strd{i},'fontsize',12,'fontweight','demi'); xlabel('x'); ylabel('y'); zlabel('z'); view([-51,18]); end function wdot=ddl***un(t,w,q,m,b,e) %该函数实现(3)式的差分迭代 syms q m b e q=1.6e-2; m=0.02; b=[2;1;0]; e=[1;0;1]; wdot=[w(2); q.*b.*w(4)./m;w(4); q.*e./m-q.*b.*w(2)./m; w(6); 0] |
|
相关推荐 |
|
你正在撰写答案
如果你是对答案或其他答案精选点评或询问,请使用“评论”功能。
我用matlab的coder 封装了一个c语言的dll, 但是在用labview调用时,会出现识别不到库的问题,有大神遇到过吗
4710 浏览 0 评论
5517 浏览 0 评论
在matlab中如何计算含有第一类修正的贝塞尔函数的积分算不出的问题?
9321 浏览 0 评论
怎么利用matlab得到95%,80%和70%的置信区间,并生成不同区间下的功率误差贝塔分布?
10761 浏览 0 评论
请问simulink的s-function模块如何添加多输入输出接口
15034 浏览 2 评论
/9
小黑屋| 手机版| Archiver| 电子发烧友 ( 湘ICP备2023018690号 )
GMT+8, 2025-12-12 23:31 , Processed in 0.692453 second(s), Total 73, Slave 55 queries .
Powered by 电子发烧友网
© 2015 bbs.elecfans.com
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 (电路图) 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证:合字B2-20210191

淘帖
3611