Научный форум - Форум программистов и сисадминов Киберфорум

Расчитываю матрицу поворота для спина S =1. Поворот выглядит так.

D = 2870;

Q = -5;

A = -2;

Wn = 1;

We = -1400;

gamma_n = -307.7*10^(-6);

gamma_e = 2.8024;



Sx = 0.5^0.5*[0,1,0;1,0,1;0,1,0];

Sz = [1,0,0;0,0,0;0,0,-1];

Sy = (1i/sqrt(2))*[0,-1,0;1,0,-1;0,1,0];

I = eye(3);



S1z = SMGen(Sz,I,n,1);

S2z = SMGen(Sz,I,n,2);



S1x = SMGen(Sx,I,n,1);

S2x = SMGen(Sx,I,n,2);



S1y = SMGen(Sy,I,n,1);

S2y = SMGen(Sy,I,n,2);

I = eye(3^n);



W14 = D+We+A;

W47 = A-D+We;

W25 = D+We;

W58 = -D+We;

W36 = -A+D+We;

W69 = -A-D+We;



W12 = A+Q+Wn;

W23 = A-Q+Wn;

W45 = Q+Wn;

W56 = -Q+Wn;

W78 = -A+Q+Wn;

W89 = -A-Q+Wn;



le = 20000;

ln = 20000;



M1 = Y_rotation(n,8,9,-pi);



t = (2*pi*ln)/Wn;

theta = pi;

h_pulse = abs(theta/(sqrt(2)*t*gamma_n));

phi = 3*(pi/2);

H = A*S1z*S2z+D*S1z*S1z+Q*S2z*S2z+(We-W89)*S1z+(Wn-W89)*S2z;

H = H-gamma_e*h_pulse*(S1x*cos(phi)+S1y*sin(phi));

H = H-gamma_n*h_pulse*(S2x*cos(phi)+S2y*sin(phi));

P1 = expm(-1i*t*H);

Матрицы P1 и M1 при увеличении ln должны совпасть. Но при увеличении ln возникают анамалии. Появляются лишние недиагональные члены. S1 и S2 матрицы спинов 1 и 2 размером 9x9, полученные посредством прямого произведения матриц.

function A = SMGen(Sz,I,n,k)

if k == 1

    A = Sz;

else

    A = I;

end;

for i =1:n

    if i == k

        if i ~=1

            B =FMGen(A,Sz,3^(i-1));

            A = B;

        end;

    else

        B =FMGen(A,I,3^(i-1));

        A = B;

    end;

end;

end

function D = FMGen(Sz,I,n)

for i =1:n

    for j = 1:n

        if j == 1

            A =Sz(i,j)*I;

        else

            B = Sz(i,j)*I;

            C = cat(2,A,B);

            A = C;

            clear C;

        end;

    end;

    if i == 1

        D = A;

    else

        C = cat(1,D,A);

        D = C;

        clear C;

    end;

end;

end

function R = Y_rotation(n,l1,l2,theta)

    phi = pi/2;

    R = eye(3^n);

    R(l1,l1) = cos(theta/2);

    R(l1,l2) = -1i*exp(-1i*phi)*sin(theta/2);

    R(l2,l1) = -1i*exp(1i*phi)*sin(theta/2);

    R(l2,l2) = cos(theta/2);

end

]]>

J.dinAlt Tipi Polirendo	J.dinAlt Внутренняя очистка
Tirr ra	Солнечный бог войны
Pir	Мы (Мы/Они/Множество)
Pir ra	Мы солнце (Мы/Они/Множество солнце)
Sori ma ya	Дарующая мать-космос
Terri sorro puto resto	Земля скорбит об оставленном
Cherri koserro U O	Жизнь (цветение) пронизывает нас
Terri sorro puto resto	Земля скорбит об оставленном
Cherri koserro U E	Жизнь (цветение) Ты, Мы
Tu	Ты
Sa gu tu*	109 бог - Ты
Ma ra tu	Мать солнце - Ты
Tu ki ra tu	Ты их ключ солнца
Tipi Polirendo	Внутренняя очистка
Tipi Ku ra to	Внутри - с солнцем разрушения
Tipi Tolisendo	Внутреннее принятие
Savu Auo	Святое омовение
Tipi Polirendo	Внутренняя очистка
Tipi Ku ra to	Внутри - с солнцем разрушения
Tipi Tolisendo	Внутреннее принятие
Savu Auo	Святое омовение
Ke poro u te sto	Кто/Что тут? Я рядом с тобой
Tu ri ma sa ya	Ты, великий(ая), рождённый(ая) космосом
Pi na ro va tekso	Не тяни огонь в воду
Koli su ra ya	Кто Полубоги Света и Космоса?
Tu	Ты
Sa gu tu*	109 бог - Ты
Ma ra tu	Мать солнце - Ты
Tu ki ra tu	Ты их ключ солнца
Tipi Polirendo	Внутренняя очистка
Tipi Ku ra to	Внутри - с солнцем разрушения
Tipi Tolisendo	Внутреннее принятие
Savu Auo	Святое омовение
Tipi Polirendo	Внутренняя очистка
Tipi Ku ra to	Внутри - с солнцем разрушения
Tipi Tolisendo	Внутреннее принятие
Savu Auo	Святое омовение
Soto fei to	Снаружи - огонь
Soto ra	Снаружи - солнце
Soto hei vo	Снаружи - ветер
Soto va	Снаружи - вода
Tipi Polirendo	Внутренняя очистка
Tipi Ku ra to	Внутри - с солнцем разрушения
Tipi Tolisendo	Внутреннее принятие
Savu Auo	Святое омовение
Tipi Polirendo	Внутренняя очистка
Tipi Ku ra to	Внутри - с солнцем разрушения
Tipi Tolisendo	Внутреннее принятие
Savu Auo	Святое омовение
Ku sa	Жить в гармонии
Su tei ro	Природа будет
Ku sa	Жить в гармонии
Su tei ro	Природа будет
Ku sa	Жить в гармонии
Su tei ro pei du	Природа будет даже тогда, когда не будет нас
Ku sa	Жить в гармонии