Смоделировать потери пучка молекул через прохождение ими через вещество

25.05.2012, 12:12. **Показов** 937. **Ответов** 0

Author24 — интернет-сервис помощи студентам

Суть в том, что я не очень разбираюсь в матлабе, и был бы безмерно благодарен если кто-либо просто скинул готовый физический эксперимент

Пытался смоделировать потери пучка молекул через прохождение ими через вещество (точнее воссоздать этот эксперимент по методичке http://apicentr2.narod.ru/copy... Matlab.pdf ) но как я понимаю успехом моя попытка не увенчалась...

Matlab M

clear;
Na=250; k=(1:Na)';
no=100; Section=0.05;
dx=0.01;
x=zeros(size(k));
sc=ones(size(k));
dW=no*Section*dx;
hl=line(x,k);
set(hl,'Marker','o','MarkerSize',3);
axis([0 1 0 Na+1]);
pause;
while(any(sc)>0);
ra=find(ra-dW<0);
sc(k1)=0;
x=x+sc*dx;
set(hl,'XData',x);
end;
clear;
Na=25;
L=5;k=(1:Na)';no=100;dx=0.01;
SecA=0.1; SecAB=0.15;
x=zeros(size(k));
scA=ones(size(k));
dWa=no*dx*SecA;
dWab=no*dx*SecAB;
while (any(scA)>0 & all(x)<L);
    ra=rand(size(k));
    ka=find(ra-dWa<0);
    scA(ka)=0;
    kb=find(dWa<ra & ra<dWa+dWab);
    scB(kb)=scA(kb);
    scA(kb)=0;
    x=x+scA*dx+scB*dx;
    pause(1);
end;
ra=rand(size(k));
ka=find(ra-Na*Sect*dX<0);
for ii=ka
while 1
eps(ii)=epsmax*rand;
if Fmax*rand<F(E(ii),eps(ii))
break;
end;

	25.05.2012, 12:12