5.4 Program 5.1, Matlab Code
77
Fig. 5.9 Typical ω stat vs F load curves, We used the same parameters as Fig. 5.4
5.4 Program 5.1, Matlab Code
Solve the Eq. for the rotation of a dipole in an Electric Field E
%help rotor;
% T dip ≡ E
p , T rat ≡
E(rat)
p
% T dipo ≡ E
p (0), T rato ≡
E(rat)
p
(0)
clear;
N=20000; NR=1000; Drot=2.e8; omega=4.e7; d=50.e-10; eo=1.6e-19; kT=4.14e-21;
E=4.e9; Eo=4.e9; pi=3.14159;
tB=1/(2.*Drot);
tE=2*pi/omega;
dt=0.1*tB;
x=zeros(N-1,NR);
w=zeros(N-1,NR);