高斯消元法(matlab程序)
来源:互联网 发布:手机数据可以恢复吗 编辑:程序博客网 时间:2024/06/01 20:36
function x=Untitled2(A,b)
n=size(A,1);
x=zeros(n,1);
for j=1:n-1
T=[];
a=A(j:n,j);
m=find(abs(a)==max(abs(a)));
A([j,m(1)+j-1],:) = A([j+m(1)-1,j],:);
[b(j), b(j+m(1)-1)] = deal(b(m(1)+j-1),b(j));
for i=1:n-j
T(i)=-A(i+j,j)/A(j,j);
A(i+j,:)=A(i+j,:)+A(j,:)*T(i);
b(i+j)= b(i+j)+b(j)*T(i);
end
end
for i=n:-1:1
sum=0;
for j=n:-1:i+1
sum=sum+x(j)*A(i,j);
end
x(i)=(b(i)-sum)/A(i,i);
end
B=[1,1,1,1;2,1,1,1;3,2,1,1;4,3,2,1];
c=[4,5,7,10];
x=Untitled2(B,c);
n=size(A,1);
x=zeros(n,1);
for j=1:n-1
T=[];
a=A(j:n,j);
m=find(abs(a)==max(abs(a)));
A([j,m(1)+j-1],:) = A([j+m(1)-1,j],:);
[b(j), b(j+m(1)-1)] = deal(b(m(1)+j-1),b(j));
for i=1:n-j
T(i)=-A(i+j,j)/A(j,j);
A(i+j,:)=A(i+j,:)+A(j,:)*T(i);
b(i+j)= b(i+j)+b(j)*T(i);
end
end
for i=n:-1:1
sum=0;
for j=n:-1:i+1
sum=sum+x(j)*A(i,j);
end
x(i)=(b(i)-sum)/A(i,i);
end
B=[1,1,1,1;2,1,1,1;3,2,1,1;4,3,2,1];
c=[4,5,7,10];
x=Untitled2(B,c);
阅读全文
0 0
- 高斯消元法(matlab程序)
- matlab好玩小程序(画心形)【matlab】
- Matlab数据处理程序(2)
- matlab练习程序(PSNR)
- Matlab程序学习(一)
- matlab程序
- 矩阵的LU分解(Matlab程序)
- 函数插值计算(Matlab程序)
- matlab练习程序(SUSAN检测)
- matlab练习程序(SUSAN检测)
- 数值方法 回代程序(matlab)
- matlab练习程序(k-means聚类)
- matlab练习程序(radon变换)
- matlab练习程序(k-means聚类)
- matlab练习程序(直方图匹配)
- Matlab 独立运行程序编译(Linux)
- matlab练习程序(模拟退火SA)
- MATLAB频谱分析程序(转)
- Android layout_weight属性使用
- 关于高并发支付、秒杀的一些设计思路
- caffe简介
- 数组与指针区别和联系
- 2017北大信科推免机考+面经
- 高斯消元法(matlab程序)
- 机器学习第一天
- PL/SQL之九 函数
- 3_使用C语言控制LED灯实验
- hdu5934(tarjan算法+缩点)bomb
- 一行代码搞定Parcelable 告别Parcelable繁琐代码!
- 条件变量
- 43. Multiply Strings
- Configure GRUB2 Boot Loader Settings In Ubuntu 16.04