其实用matlab很方便的就求出通解了,用null命令就行
例: 求下面非齐次线性方程的通解:
>> A=[2 4 -1 4 16;-3 -6 2 -6 -23;3 6 -4 6 19;1 2 5 2 19]
A =
2 4 -1 4 16
-3 -6 2 -6 -23
3 6 -4 6 19
1 2 5 2 19
>> b=[-2 7 -23 43]'
b =
-2
7
-23
43
>> B=[A,b]
B =
2 4 -1 4 16 -2
-3 -6 2 -6 -23 7
3 6 -4 6 19 -23
1 2 5 2 19 43
>> rref(B)
ans =
1 2 0 2 9 3
0 0 1 0 2 8
0 0 0 0 0 0
0 0 0 0 0 0
又x=null(B,'r')
x =
-2 -2 -9 -3
1 0 0 0
0 0 -2 -8
0 1 0 0
0 0 1 0
0 0 0 1
故通解为:。。。。。