为什么一直出错呀。。。。那位高手帮帮忙。。

gh.sha

function [A,nx,ny]=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

水木年华zzu

查查这个值q+(p-1)*(nx-1)，出现了负数  p=mod(i,nx-1);
    q=fix(i./(nx-1))+1;这两句前面是不是要加上ABS

gh.sha

回复 _math 的帖子
高手呀！！！

   

zsy312

function [A,nx,ny]=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

柯雷

dingyixiaxiaxiaxia

yujianyuanhaha

真正的高手是不上网的···

伤感的幸福

应该有问题啊
实现什么功能啊

_math

function [A,nx,ny]=aaaa(n)
ny=n;
nx=3*(ny-1)+1;
A=zeros((nx-1)*(ny-1),(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
楼主貌似if语句那里写错了。。。
我的可以运行的。。
aaaa(21)
运行后矩阵中有190个1和-1

linmatsas

0<p<ny那个程序可以用这样的符号的………………肯定分开写呀………………

gh.sha

走过路过的高手帮帮忙喽。。。。