exm040208.m

来自「matlab电子书籍」· M 代码 · 共 12 行

A=reshape(1:12,4,3);	
b=(13:16)';	

ra=rank(A)				
rab=rank([A,b])	

xs=A\b;					
xg=null(A);				
c=rand(1);				
ba=A*(xs+c*xg)			
norm(ba-b)