diff options
Diffstat (limited to 'liouville.m')
| -rw-r--r-- | liouville.m | 31 |
1 files changed, 5 insertions, 26 deletions
diff --git a/liouville.m b/liouville.m index 0673a7e..5987925 100644 --- a/liouville.m +++ b/liouville.m @@ -1,3 +1,7 @@ + +% load useful functions; +useful_functions; + % some physical constants hbar=1; im_one=0+1i; @@ -5,6 +9,7 @@ im_one=0+1i; three_levels; %two_levels; + Nfreq=length(modulation_freq); % now we create Liouville indexes list @@ -35,21 +40,6 @@ L=zeros(N); % NxN matrix Li=zeros(N); % NxN Liouville interactive L0=zeros(N); % NxN Liouville from unperturbed hamiltonian -function ret=decay_total(g_decay,i) - ret=0; - for k=1:size(g_decay)(1) - ret=ret+g_decay(i,k); - endfor -endfunction - -function ret=kron_delta(i,j) - if ((i==j)) - ret=1; - else - ret=0; - endif -endfunction - for p=1:N for s=1:N j=rhoLiouville_r(p); @@ -124,17 +114,6 @@ rhoLiouville_dot(1)=1; %solving for density matrix vector rhoLiouville=L\rhoLiouville_dot; -% this function create from Liouville density vector -% the density matrix with given modulation frequency -function rho=rhoOfFreq(rhoLiouville, freqIndex, Nlevels, Nfreq) - rho=zeros(Nlevels); - for r=1:Nlevels - for c=1:Nlevels - rho(r,c)=rhoLiouville((freqIndex-1)*Nlevels^2+(r-1)*Nlevels+c); - endfor - endfor -endfunction - rho_0=rhoOfFreq(rhoLiouville, 1, Nlevels, Nfreq) rho_1=rhoOfFreq(rhoLiouville, 2, Nlevels, Nfreq) |