function [ basis ]=HubbardBose_Basis(N,L) D=factorial(L+N-1)/factorial(L-1)/factorial(N); %First of all we have to build the basis% basis=zeros(D,L);basis(1,1)=N; k=1; for i=2:D basis(i,:)=basis(i-1,:); if basis(i,k)>0 basis(i,k) = basis(i-1,1)-1; basis(i,k+1) = basis(i,k+1)+1; elseif basis(i,k)==0 ni = find(basis(i,:)>0); ni=ni(1); if (ni==L && basis(i,L)==N) break else basis(i,1)=basis(i,ni)-1; basis(i,ni)=0; basis(i,ni+1)=basis(i,ni+1)+1; end end end end

这个函数是用于构建 Hubbard-Bose 模型的基态基的函数。给定粒子数 N 和格点数 L，该函数返回一个大小为 D × L 的矩阵 basis，其中 D 是基态基的维数。每一行表示一个基态，每一列表示一个格点上的粒子数。

在该函数中，首先计算基态基的维数 D。然后初始化 basis 矩阵，并将第一行的第一个格点设置为 N。接下来，使用循环生成其他基态。在每次循环中，复制上一行的基态到当前行，并根据不同的情况进行更新。

如果当前格点上的粒子数大于 0，则将该格点上的粒子数减1，并将下一个格点上的粒子数加1。如果当前格点上的粒子数等于 0，则找到第一个非空格点，并将其粒子数减1，将当前格点上的粒子数设置为该非空格点减1，将下一个格点上的粒子数加1。

最后，返回生成的基态基矩阵 basis。