ÕæÕýµÄ¸ßÊÖÊDz»ÉÏÍøµÄ¡¤¡¤¡¤
dingyixiaxiaxiaxia
function =f(n)
ny=n;
nx=3*(ny-1)+1;
A=zeros(nx-1,ny-1);
for i=1: (nx-1)*(ny-1)
p=mod(i,nx-1);
q=fix(i./(nx-1))+1;
if 0<p&&p<ny && p<q&&q<ny
k=q+(p-1)*(nx-1);
A(i,i)=1;
A(i,k)=-1;
end
end
»Ø¸´ _math µÄÌû×Ó
¸ßÊÖѽ£¡£¡£¡{:3_59:}
²é²éÕâ¸öÖµq+(p-1)*(nx-1)£¬³öÏÖÁ˸ºÊý p=mod(i,nx-1);
q=fix(i./(nx-1))+1;ÕâÁ½¾äÇ°ÃæÊDz»ÊÇÒª¼ÓÉÏABS