我正在尝试在细胞壁的范围内生成随机游走模型。我有以下代码用于单元格内的离散步进:
x(1000)=0;
t(1000)=0;
dt=1;
kB=1.381e-23;
T=310;
R=4.89e-9;
visc=0.007;
friction=6*pi*visc*R;
D=kB*T/friction;
dxdiff=(2*D*dt)^0.5;
Force_applied=1e-15;
for i=(2:1000)
vx=Force_applied/friction;
dx=vx*dt;
x(i)=x(i-1)+dxdiff*(2*(round(rand)-0.5))+dx;
t(i)=(i-1)*dt;
end
哪个有效,只有它没有绑定。我如何将它绑定以使其保持在具有给定直径值的球形单元内?