灵鹊做题网matlab求线性齐次方程组!……
来源：学生作业帮 编辑：灵鹊做题网作业帮 分类：数学作业 时间：2024/06/11 15:57:49
matlab求线性齐次方程组!……
现在有一个六阶方阵A,矩阵里含有一个未知数w.我先解出detA=0时的w值,然后我想把w代回矩阵里,求齐次方程Ax=0的解.但是由于精度的问题将解出来的w代回矩阵时,detA不为零,也就解不出齐次方程组的非零解了.请问这种情况该怎么解出对应不同w的方程组的非零解?
用rref返回的是单位阵，还是解不出。有其他办法吗？
现在有一个六阶方阵A,矩阵里含有一个未知数w.我先解出detA=0时的w值,然后我想把w代回矩阵里,求齐次方程Ax=0的解.但是由于精度的问题将解出来的w代回矩阵时,detA不为零,也就解不出齐次方程组的非零解了.请问这种情况该怎么解出对应不同w的方程组的非零解?
用rref返回的是单位阵，还是解不出。有其他办法吗？
把含w的矩阵贴出来看看
再问： A=w*M-C
挺复杂的。
再答： 可以解，可能你的matlab代码没有写对 C=[0.1261 -0.1261 0 0 0 0; -0.1261 0.8763 -0.7502 0 0 0; 0 -0.7502 1.3170 -0.5669 0 0; 0 0 -0.5669 1.3170 -0.7502 0; 0 0 0 -0.7502 1.3170 -0.5669; 0 0 0 0 -0.5669 0.9784]*1E6; M=diag([0.0005 0.0054 0.01 0.0054 0.01 0.0006]); syms w real; y=det(w*M-C); ww=solve(y); for i=1:length(ww) A=double(ww(i))*M-C; null(A) end
再问： A=w*M-C
挺复杂的。再答： 可以解，可能你的matlab代码没有写对 C=[0.1261 -0.1261 0 0 0 0; -0.1261 0.8763 -0.7502 0 0 0; 0 -0.7502 1.3170 -0.5669 0 0; 0 0 -0.5669 1.3170 -0.7502 0; 0 0 0 -0.7502 1.3170 -0.5669; 0 0 0 0 -0.5669 0.9784]*1E6; M=diag([0.0005 0.0054 0.01 0.0054 0.01 0.0006]); syms w real; y=det(w*M-C); ww=solve(y); for i=1:length(ww) A=double(ww(i))*M-C; null(A) end