小说屋

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

别买这章全是（第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？我直接速通！  御兽：我能增强御兽血脉  星囚禁地  斗罗好孕：我靠生子系统被团宠了  全村亿万富翁，你管这叫贫困村？  射雕之逍遥天下  我的宗门有点怪！  倾月班  娇娇贵女红了眼，摄政王连夜哄  有樗  摆摊卖水果  官道危途  我在修仙界养猪  贫僧法海，万世之尊