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

clear;clc;close all
q=2;
N=1;M=q*N;
syms LMat=sym(zeros(2*M,2*M));for ii=1:MTp=[];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
endTp=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;'])endend
end%latex(Mat)
%collect合并同类项
latex(collect(det(Mat),Tp))

一个例子：