笔趣阁

手机浏览器扫描二维码访问

库代尔玛拉公主选婿（第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？我直接速通！  贫僧法海，万世之尊  七零之我靠纸扎手艺混的风生水起  禁锢之日  御兽：我能增强御兽血脉  我的宗门有点怪！