哪位高手帮忙解一下详细过程的方程。谢谢!
展开全部
这应该属于线性代数中,求解线性方程组AX=b中的解X(X1,X2……Xn),这个程序基本上可以用matlab中的一句话代替 X=A\b ,只是这个运算在方程组有无穷多解的情况下会报错。 你的程序的详细注释如下: A=[6,-2,2,4;12,-8,4,10;3,-13,3,3;-6,4,2,18]; % 线性方程组系数矩阵 A b=[0;-10;-39;-16]; % AX=b 中的 b B=[A b]; %构造 增广矩阵 B ra=rank(A) % 求解系数矩阵 A 的秩 rb=rank(B) % 求解增广矩阵 B 的秩 n=length(b); % 得到 b 的长度,后面要用到 X=zeros(n,1); % 初始化向量 X 用于存放方程组的解 C=zeros(1,n+1); % 初始化 C 用于下面的两行交换 if ra>rb % 如果增广矩阵的秩大于系数矩阵的秩,则方程组无解,(提示一下,这里程序错了,应该是 ra
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询