gh.sha ·¢±íÓÚ 2010-12-17 11:11

Ϊʲôһֱ³ö´íѽ¡£¡£¡£¡£ÄÇλ¸ßÊÖ°ï°ïæ¡£¡£

function =aaaa(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<ny && p<q<ny
        k=q+(p-1)*(nx-1);
       A(i,i)=1;
       A(i,k)=-1;
    end
end

gh.sha ·¢±íÓÚ 2010-12-17 11:12

¸ßÊÖÔÚÄÄÀïѽ£¿£¿£¿

gh.sha ·¢±íÓÚ 2010-12-17 11:12

:sleepy::sleepy:

gh.sha ·¢±íÓÚ 2010-12-17 11:13

ÄÇλ¸ßÊÖ°ï°ïæඡ£¡£¡£¡£:)

madio ·¢±íÓÚ 2010-12-17 11:15

ÌáʾÄÄÀïÓÐ´í£¿

gh.sha ·¢±íÓÚ 2010-12-17 11:16

??? Attempted to access A(60,-58); index must be a positive integer or logical.

Error in ==> aaaa at 10
       A(i,q+(p-1)*(nx-1))=-1;

gh.sha ·¢±íÓÚ 2010-12-17 11:16

ÔËÐоͳöÎÊÌâ¡£¡£¡£¡£:(

gh.sha ·¢±íÓÚ 2010-12-17 11:21

¸ßÊÖÔÚÄÄÀïѽ£¿¸ßÊÖÔÚÄÄÀ

gh.sha ·¢±íÓÚ 2010-12-17 11:22

¸ßÊÖÔÚÄÄÀïѽ£¿¸ßÊÖÔÚÄÄÀ

gh.sha ·¢±íÓÚ 2010-12-17 11:23

¸ßÊÖÔÚÄÄÀïѽ£¿¸ßÊÖÔÚÄÄÀ
Ò³: [1] 2 3
²é¿´ÍêÕû°æ±¾: Ϊʲôһֱ³ö´íѽ¡£¡£¡£¡£ÄÇλ¸ßÊÖ°ï°ïæ¡£¡£