rotop.m

来自「The goal of SPID is to provide the user 」· M 代码 · 共 18 行

function output=rotop(spin_num,xy,angle);
    if xy=='x'
        phi=0;
    end
    if xy=='y'
        phi=90;
    end
    beta=angle;
    c=cos(1/2*beta*pi/180);
    s=sin(1/2*beta*pi/180);
    eplus=exp(+i*phi*pi/180);
    eminus=exp(-i*phi*pi/180);
    op=[c -i*s*eminus ; -i*s*eplus c];
    output=op;
    for n=1:(spin_num-1)
        output=kron(op,output);
    end 
end