matlab中用高斯消元法接方程组程序方程组如下:10 -1 2 0】 6-1 11 -1 3】 25 2 -1 10 -1】 -110 3 -1 8】 15左边的是矩阵A,右边的是b
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/06 14:18:59
![matlab中用高斯消元法接方程组程序方程组如下:10 -1 2 0】 6-1 11 -1 3】 25 2 -1 10 -1】 -110 3 -1 8】 15左边的是矩阵A,右边的是b](/uploads/image/z/5253974-62-4.jpg?t=matlab%E4%B8%AD%E7%94%A8%E9%AB%98%E6%96%AF%E6%B6%88%E5%85%83%E6%B3%95%E6%8E%A5%E6%96%B9%E7%A8%8B%E7%BB%84%E7%A8%8B%E5%BA%8F%E6%96%B9%E7%A8%8B%E7%BB%84%E5%A6%82%E4%B8%8B%EF%BC%9A10+-1+2+0%E3%80%91+6-1+11+-1+3%E3%80%91+25+2+-1+10+-1%E3%80%91+-110+3+-1+8%E3%80%91+15%E5%B7%A6%E8%BE%B9%E7%9A%84%E6%98%AF%E7%9F%A9%E9%98%B5A%2C%E5%8F%B3%E8%BE%B9%E7%9A%84%E6%98%AFb)
matlab中用高斯消元法接方程组程序方程组如下:10 -1 2 0】 6-1 11 -1 3】 25 2 -1 10 -1】 -110 3 -1 8】 15左边的是矩阵A,右边的是b
matlab中用高斯消元法接方程组程序
方程组如下:
10 -1 2 0】 6
-1 11 -1 3】 25
2 -1 10 -1】 -11
0 3 -1 8】 15
左边的是矩阵A,右边的是b
matlab中用高斯消元法接方程组程序方程组如下:10 -1 2 0】 6-1 11 -1 3】 25 2 -1 10 -1】 -110 3 -1 8】 15左边的是矩阵A,右边的是b
function [X]=gaussian(A,B,n)
C=[A B];
for i=2:n
[r,c]=max(abs(C(i-1:n,i-1)));
VEC=C(i-1,:);C(i-1,:)=C(c+i-2,:);C(c+i-2,:)=VEC;
if C(i-1,i-1)==0
return
end
for k=(i-1):(n-1)
C(k+1,:)=C(i-1,:)*-C(k+1,i-1)/C(i-1,i-1)+C(k+1,:);
end
end
D=C;
B2=C(:,n+1);
D(:,n+1)=[];
X=zeros(1,n);
for i=1:n
X(n+1-i)=(B2(n+1-i)-D(n+1-i,:)*X')/C(n+1-i,n+1-i);
end
运行:
>>A=[10,-1,2,0;-1,11,-1,3;2,-1,10,-1;0,3,-1,8];
>>B=[6,25,-11,15]';
>>gaussian(A,B,4)
ans =
1.0000 2.0000 -1.0000 1.0000