如何求齐次线性方程组的非平凡解。我尝试使用solve命令,但它只给出了平凡的解决方案。eigen(A)$vector[,x]只给出了方阵的答案,即偶数确定的系统。
solve
eigen(A)$vector[,x]
wh6knrhe1#
你可以使用奇异值分解,对于ddivya,M <- rbind(c(2,3,4), 1),然后执行,例如:
M <- rbind(c(2,3,4), 1)
s <- La.svd(M, nv=ncol(M)) (v <- s$vt[ncol(M),])
则v的任意倍数都是一个解;您可以验证M %*% v是否为零。您的两个系统的列数都比行数多,并且秩都是满的(s$d中没有元素为零或接近零)。否则,您会在$s$vt$的底部找到更多的基向量,这些基向量跨越了您的解空间。
v
M %*% v
s$d
1条答案
按热度按时间wh6knrhe1#
你可以使用奇异值分解,对于ddivya,
M <- rbind(c(2,3,4), 1)
,然后执行,例如:则
v
的任意倍数都是一个解;您可以验证M %*% v
是否为零。您的两个系统的列数都比行数多,并且秩都是满的(s$d
中没有元素为零或接近零)。否则,您会在$s$vt$的底部找到更多的基向量,这些基向量跨越了您的解空间。