exm0651_2.m

来自「学matlab还有比这个更好的吗?我千兴万苦才找到的,这是一个学习matlab的」· M 代码 · 共 12 行

%exm0651_2.m

clear
A=sym([1 1/2 1/2 -1;1 1 -1 1;1 -1/4 -1 1;-8 -1 1 1]);
b=sym([0;10;0;1]);

syms k
A2=A(1:3,:);
X2=A2\b(1:3,1)					%求一个特解：最少非零元素的最小二乘解
XX2=X2+k*null(A2)					%构成通解
A2*XX2								%验算