acav.m at c33ae98dafbc4f154b8372938c2493a50f4386a1 - concept-collection/numbl-open-quantum-systems

concept-collection / numbl-open-quantum-systems

12 lines · 263 BBlame History Raw

1function x = acav(i,L,Nphoton,Is,Iatom)

2a = diag(sqrt(1:Nphoton)',1);

3a = kron(Iatom,a);

4Op_total = cell(1,L);

5for site = 1:L

6 Op_total{site} = Is+double(eq(i,site))*(a-Is);

7end

8x = Op_total{1};

9for site = 2:L

10 x = kron(x,Op_total{site});

11end

12end