6.5 符号代数方程的求解
6.5.1 线性方程组的符号解
【 * 例 6.5.1 -1 】求 
线性方程组的解。本例演示，符号线性方程组的基本解法。
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]);X1=A/b
X1 =
[ 1]
[ 8]
[ 8]
[ 9]
【 * 例 6.5.1 -2 】求解上例前 3 个方程所构成的“欠定”方程组，并解释解的含义。
syms k
A2=A(1:3,:);X2=A2/b(1:3,1) % 求一个特解：最少非零元素的最小二乘解
XX2=X2+k*null(A2) % 构成通解
A2*XX2 % 验算
X2 =
[ 0]
[ 8]
[ 4]
[ 6]
XX2 =
[ k]
[ 8]
[ 4+4*k]
[ 6+3*k]
ans =
[ 0]
[ 10]
[ 0]
6.5.2 一般代数方程组的解
【 * 例 6.5.2 -1 】求方程组 
， 
关于 
的解。
S=solve('u*y^2+v*z+w=0','y+z+w=0','y','z')
disp('S.y'),disp(S.y),disp('S.z'),disp(S.z)
S =
y: [2x1 sym]
z: [2x1 sym]
S.y
[ -1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))-w]
[ -1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))-w]
S.z
[ 1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))]
[ 1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))]
【 * 例 6.5.2 -2 】用 solve 指令重做例 6.5.1-2 。即求 
， 
， 
构成的“欠定”方程组解。
syms d n p q;eq1=d+n/2+p/2-q;eq2=n+d+q-p-10;eq3=q+d-n/4-p;
S=solve(eq1,eq2,eq3,d,n,p,q);S.d,S.n,S.p,S.q
Warning: 3 equations in 4 variables.
> In E:/MAT53/toolbox/symbolic/solve.m at line 110
In E:/MAT53/toolbox/symbolic/@sym/solve.m at line 49
ans =
d
ans =
8
ans =
4*d+4
ans =
3*d+6
【 * 例 6.5.2 -3 】求 
的解。
clear all,syms x;s=solve('(x+2)^x=2','x')
s =
.69829942170241042826920133106081
上一篇： Matlab 调试器的使用
下一篇：Matlab 符号微分方程的求解
