lls1.m

来自「一些matlab的经典例程」· M 代码 · 共 26 行

function LLS1()
Ae=3;																			%<2>
Av=5;																			%<3>
[A,b]=LLSDATA(Ae,Av);
x=A\b;
%以下都是为获得较清晰的显示而编写
S='恰定';
if Ae>Av
   S='超定';
elseif Ae<Av
   S='欠定';
end
disp(['外部独立程序求' S '方程 Ax=b 解，其中'])
cs=blanks(Ae)';			%产生一个由方程数决定的空白字符列
ns=fix(Ae/2);				%从零方向最接近空白列长之半的数
As=cs;As(ns)='A';			
bs=cs;bs(ns)='b';
es=cs;es(ns)='=';
disp([As cs es cs num2str(A) cs cs cs cs cs bs cs es cs num2str(b)])
disp('方程的解 ')
nxs=fix(Av/2);
cxs=blanks(Av)';
xs=cxs;xs(nxs)='x';
exs=cxs;exs(nxs)='=';
disp([xs cxs exs cxs num2str(x)])