25.07.2013 Views

Ph.D. Thesis - Physics

Ph.D. Thesis - Physics

Ph.D. Thesis - Physics

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

J = J0;<br />

Jt(1) = J;<br />

Bt(1) = B;<br />

for kk = 2:length(t)<br />

Fr = 1-A*grad*sin(Omega*t(kk)+pi/2);<br />

Jt(kk) = J*Fr^2*(1+A*sin(Omega*t(kk)+pi/2))^(-3);<br />

Bt(kk) = B;<br />

Ht = Jt(kk)*kron(sz,sz)+ B*(kron(eye([2 2]), sx) + kron(sx,eye([2 2])));<br />

curr = expm(-i*Ht*tstep)*curr;<br />

expsz2 = [expsz2, curr’*sztot1*curr];<br />

end<br />

Jm = mean(Jt);<br />

theerr = expsz2-expsz1;<br />

err_result = mean(theerr);<br />

233

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!