incorrect result in solving system using MATLAB ?
2 次查看(过去 30 天)
显示 更早的评论
I am trynig to solve the following system (mat=coefficient matrix,B=output matrix)
mat =
1 2 6
2 5 14
5 7 24
B=
0
0
0
the expected result is not zeros for sure because this system is linear dependent . And I used mat/B and linsolve(mat,B).However , I have got the answer
ans=
0
0
0
and this warning
Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND =
1.009294e-18.
could any one help me why I got this result and what I should to get the correct answer?
0 个评论
采纳的回答
Torsten
2017-11-15
Try
Z = null(mat)
It should give you a vector Z different from the zero vector that solves
mat*Z = 0
Best wishes
Torsten.
8 个评论
Steven Lord
2017-11-16
Was the result of null(mat, 'r') closer to what you found when you did it manually? From the help for the null function:
"The orthonormal basis is preferable numerically, while the rational basis may be preferable pedagogically."
The rational basis is what you get when you add the 'r' flag.
更多回答(2 个)
Steven Lord
2017-11-14
That answer is one solution to your system of equations. It is not the only solution to that system.
If you add x and any linear combination of columns from the nullspace of your matrix that will be another solution to your system. Take a look at the null function for more information.
0 个评论
Walter Roberson
2017-11-14
编辑:Walter Roberson
2017-11-14
rank(mat) is 2 because the third column is twice the sum of the other two columns. You cannot use the \ operator with singular matrix.
If the b matrix was not all zero then pinv(mat)*b might work, but since it is all 0 the result is going to be all 0.
With the all-zero b, you are effectively looking null spaces, for which you should look at null(mat)
2 个评论
Walter Roberson
2017-11-14
11 * [1 2 6] - 3 * [2 5 14] = [5 7 24] . Therefore the last row does not add any new information to the matrix and the rank is 2.
"and if I get result rather than zeros this means it is Linear dependent"
No, it would mean that they were linearly independent
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!