examp4_34.m

来自「薛定宇高等应用数学问题的MATLAB求解例子源码」· M 代码 · 共 19 行

A=[1 2 3 4; 2 2 1 1; 2 4 6 8; 4 4 2 2];  B=[1;3;2;6];
C=[A B]; [rank(A), rank(C)]

Z=null(A,'r')   % 解出规范化的化零空间

x0=pinv(A)*B    % 得出一个特解

a1=randn(1); a2=rand(1);  % 取不同分布的随机数
x=a1*Z(:,1)+a2*Z(:,2)+x0; norm(A*x-B)

Z=null(sym(A))

x0=sym(pinv(A)*B)

a1=randn(1); a2=rand(1);  % 取不同分布的随机数
x=a1*Z(:,1)+a2*Z(:,2)+x0; norm(double(A*x-B))

syms a1 a2;
x=a1*Z(:,1)+a2*Z(:,2)+x0