78
5 Rotation of a Dipole
sumx2=zeros(N-1,NR);
sumx=zeros(N-1,NR);
sigma2=zeros(N-1,NR);
raizdt=sqrt(dt/tB);
dw=raizdt*randn(N-1,NR);
Trat=zeros(N-1,NR);
Tdip=zeros(N,NR);
frot=kT/Drot;
TBmax=sqrt(2.*frot*kT);
q=1*eo;
p=q*d;
%Tdipo=p*E;
Tdipo=0.;
Trato =p*Eo;
%Trato =0.;
%Tload=(k.e-3)*Trato; % k = 1,2,3,. . . .
xo=2*pi;
%xo=pi/2;
for j=1:NR
t=0.;
x(1,j)=0.;
sumx2(1,j)=0.;
sumx(1,j)=0.;
for i=1:N-1
xij=x(i,j);
%============================================
% DIMENSION
%=============================================
t=t+dt;
phase=omega*t;
Tdip(i,j)=-Tdipo*sin(xij)*sin(phase);
Trat(i,j)= Trato*(cos(2.*pi*xij/xo) - 0.5*cos(4*pi*xij/xo));
%x(i+1,j) = xij+(dt/frot)*Tdip(i,j) + dw(i,j);
%x(i+1,j)=xij + (dt/frot)*Trat(i,j) + dw(i,j);% - (dt/frot)*Tload;
%ratio = dw(i,j)/((dt/frot)*Trat(i,j));
%ratio1 = dw(i,j)/((dt/frot)*Tdip(i,j));
%ratio2 = Trat(i,j)/Tdip(i,j); x(i+1,j)=xij+(dt/frot)*Tdip(i,j)+(dt/frot)*Trat(i,j)+dw(i,j);
% -Tload*(dt/frot);
w(i+1,j) = x(i+1,j)/t;
%w(i+1,j) = (x(i+1,j)-x(i,j))/dt;
%=============================================
% DIMENSIONLESS
%=============================================
%raizdt=sqrt(dt);
%ddw=raizdt*randn(N-1,NR);