广义布里渊区方程推导过程中一个公式的理解

是对DOI: 10.1103/PhysRevLett.123.066404补充材料公式(S25)的理解

clear;clc;close all
q=2;
N=1;

M=q*N;
syms L

Mat=sym(zeros(2*M,2*M));

for ii=1:M
    Tp=[];
    for jj=1:2*M
       
        %eval(['syms ', 'f',num2str(ii),num2str(jj)]);
        eval(['syms ','f',num2str(ii),'_beta',num2str(jj),'_ES'])
        %eval(['temp=','f',num2str(ii),'_beta',num2str(jj),'ES'])
        eval(['syms ','g',num2str(ii),'_beta',num2str(jj),'_ES'])
        eval(['syms ','beta',num2str(jj)])
        eval(['tp=','beta',num2str(jj),';'])
        Tp=[Tp,tp];
    end
end

Tp=Tp.^L;

for xx=1:2*M%行索引
    for yy=1:2*M%列索引
if xx<=M
eval(['Mat(',num2str(xx),', ',num2str(yy),')=',...
    'f',num2str(xx),'_beta',num2str(yy),'_ES;'])
elseif xx>M
eval(['Mat(',num2str(xx),', ',num2str(yy),')=',...
    'g',num2str(xx-M),'_beta',num2str(yy),'_ES','*','(beta',num2str(yy),')^L;'])

end
    end
end

%latex(Mat)
%collect合并同类项
latex(collect(det(Mat),Tp))

一个例子：

原文地址：https://blog.csdn.net/jack_BOB_LUO/article/details/143667255

免责声明：本站文章内容转载自网络资源，如本站内容侵犯了原著者的合法权益，可联系本站删除。更多内容请关注自学内容网（zxcms.com）！