byes 写道:
小case
似乎解一个非标准的线性方程组会出错？
如 3X1+2X2-X3=1
2X1-X2+X3=0
求指导
[/quote:3biwbr4y]
以上是楼主站内短信我的问题,请以后尽量不要使用站内短信,因为我很少去看站内短信的.现针对问题回答如下:
[code:3biwbr4y]proc iml;
A={3 2 -1,2 -1 1};
B={1,0};
C1=ECHELON(A);
C2=ECHELON(A||B);
r1=0;
do i=1 to nrow(C1);
if any(C1[i,])=1 then r1=r1+1;
end;
r2=0;
do i=1 to nrow(C2);
if any(C2[i,])=1 then r2=r2+1;
end;
if r1^=r2 then C='该线性方程组无解';
else if r1<ncol(A) then C='该线性方程组有无穷个解';
else if r1=ncol(A) then C=solve(A,B);
print C;
quit;[/code:3biwbr4y]作者: shiyiming 时间: 2012-4-13 11:59 标题: Re: SAS能求解线性方程组吗？ 楼主牛B