手机浏览器扫描二维码访问
别买这章全是（第2页）
%定义时间步长和总时间步数
dt=0.01;
num_steps=100;
%用于存储动画帧的数组
frames=cell（num_steps，1）;
%模拟磁漂移过程并生成动画帧
forstep=1:num_steps
%调整这里以确保矩阵与向量运算维度兼容
vx_drift=B（:，ones（1，num_particles））.*y_particles;
vy_drift=-B（:，ones（1，num_particles））.*x_particles;
%更新粒子位置
x_particles=x_particles+vx_drift*dt;
y_particles=y_particles+vy_drift*dt;
%绘制当前粒子位置和磁场分布（这里用曲面图表示磁场，用散点表示粒子）
figure（1）;
subplot（1，2，1）;
surf（x，y，B）;
shadinginterp;
view（2）;%设置为二维视图以便类似pcolor的展示效果
axis（[0Lx0Lymin（B（:））max（B（:））]）;
title（磁场分布）;
subplot（1，2，2）;
scatter（x_particles，y_particles，filled）;
axis（[0Lx0Ly]）;
title（等离子体粒子位置）;
%获取当前图形作为动画帧
frame=getframe（1）;
frames{step}=frame;
%暂停一下以便观察（可根据需要调整或去掉）
pause（0.01）;
end
%生成动画
video=VideoWriter（plaSma_magnetic_drift.avi）;
video.FrameRate=10;
open（video）;
和腹黑三叔闪婚后真香了 贫僧法海，万世之尊 奸情曝光之后 我见南山 摆摊卖水果 倾月班 我在修仙界养猪 高手达人的奇妙爱情 副本存活率0？我直接速通！ 禁锢之日 李辰赵清澜赵蕊 全村亿万富翁，你管这叫贫困村？ 七零之我靠纸扎手艺混的风生水起 我的宗门有点怪！ 有樗 剑影江湖传说 穿越七零，勾的冷面军少心慌慌 娇娇贵女红了眼，摄政王连夜哄 穿进甜文后成了霸总的恶毒前女友 御兽：我能增强御兽血脉
