高斯消去、追赶法 matlab
1. 分别用Gauss消去法、列主元Gauss消去法、三角分解方法求解方程组
程序:
(1)Guess消去法:
function x=GaussXQByOrder(A,b)
%Gauss消去法
N = size(A);
n = N(1);
x = zeros(n,1);
for i=1:(n-1)
for j=(i+1):n
if(A(i,i)==0)
disp('对角元不能为0');
return;
end
m = A(j,i)/A(i,i);
A(j,i:n)=A(j,i:n)-m*A(i,i:n);
b(j)=b(j)-m*b(i);
end
end
x(n)=b(n)/A(n,n);
for i=n-1:-1:1
x(i)=(b(i)-sum(A(i,i+1:n)*x(i+1:n)))/A(i,i);
end
命令行输入:
A=[1 -1 2 1;-1 3 0 -3 ;2 0 9 -6;1 -3 -6 19];
b=[1 3 5 7];
b=b';
x=GaussXQByOrder(A,b)
运算结果:
x =
-8.000000000000000
0.333333333333333
3.666666666666667
2.000000000000000
(2)列主元Gauss消去法
程序:
function x=GaussXQLineMain(A,b)
%列主元Gauss消去法
N = size(A);
n = N(1);
x = zeros(n,1);
zz=zeros(1,n);
for i=1:(n-1)
[~,p]=max(abs(A(i:n,i)));
zz=A(i,:);
A(i,:)=A(p+i-1,:);
A(p+i-1,:)=zz;
temp=b(i);
b(i)=b(i+p-1);
b(i+p-1)=temp;
for j=(i+1):n
m = A(j,i)/A(i,i);
A(j,i:n)=A(j,i:n)-m*A(i,i:n);
b(j)=b(j)-m*b(i);
end
end
x(n)=b(n)/A(n,n);
for i=n-1:-1:1
x(i)=(b(i)-sum(A(i,i+1:n)*x(i+1:n)))/A(i,i);
end
命令行:
A=[1 -1 2 1;-1 3 0 -3 ;2 0 9 -6;1 -3 -6 19];
b=[1 3 5 7];
b=b';
x=GaussXQLineMain(A,b)
运行结果:
x =
-8.000000000000005
0.333333333333332
3.666666666666668
2.000000000000000
(3)三角分解方法
程序:
function x = LU(A,b)
%三角分解
N = size(A);
n = N(1);
L = eye(n,n);
U = zeros(n,n);
x = zeros(n,1);
y = zeros(n,1);
U(1,1:n) = A(1,1:n);
L(1:n,1) = A(1:n,1)/U(1,1);
for k=2:n
for i=k:n
U(k,i) = A(k,i)-L(k,1:(k-1))*U(1:(k-1),i);
end
for j=(k+1):n
L(j,k) = (A(j,k)-L(j,1:(k-1))*U(1:(k-1),k))/U(k,k);
end
end
y(1)=b(1)/L(1,1);
for i=2:n
y(i)=b(i)-sum(L(i,1:i-1)*y(1:i-1));
end
x(n)=y(n)/U(n,n);
for i=n-1:-1:1
x(i)=(y(i)-sum(U(i,i+1:n)*x(i+1:n)))/U(i,i);
end
命令行:
A=[1 -1 2 1;-1 3 0 -3 ;2 0 9 -6;1 -3 -6 19];
b=[1 3 5 7];
b=b';
x=LU(A,b)
运行结果:
x =
-8.000000000000000
0.333333333333333
3.666666666666667
2.000000000000000
程序:function [times,wucha]=zhuiganfa(a,b,c,f)
%追赶法:x为所求解,times为所有乘除运算次数(即时间),wucha为误差的2-范数。
n=length(f);
x=zeros(n,1);
y=zeros(n,1);
times=0;
alpha=zeros(1,n);
p=zeros(1,n-1);
alpha(1)=b(1);
for i=2:n
p(i-1)=c(i-1)/alpha(i-1);
alpha(i)=b(i)-a(i-1)*p(i-1);
times=times+1;
end
y(1)=f(1)/b(1);
for i=2:n
y(i)=(f(i)-a(i-1)*y(i-1))/alpha(i);
times=times+1;
end
x(n)=y(n);
for i=n-1:-1:1
x(i)=y(i)-p(i)*x(i+1);
times=times+1;
end
A=zeros(n,n);
A=diag(b,0)+diag(a,-1)+diag(c,1);
wucha=norm((A*x-f'),2);
命令行(n=20):
a=repmat(11,1,19);
b=repmat(-19,1,20);
c=repmat(7,1,19);
f1=repmat(1.1,1,18);
f=[0 f1 1];
[times,wucha]=zhuiganfa(a,b,c,f)
运行结果:
times =
57
wucha =
8.009010697694412e-15
n=50时
命令行:
a=repmat(11,1,49);
b=repmat(-19,1,50);
c=repmat(7,1,49);
f1=repmat(1.1,1,48);
f=[0 f1 1];
[times,wucha]=zhuiganfa(a,b,c,f)
运行结果:
times =
147
wucha =
1.292635294609912e-14
命令行(n=100)
a=repmat(11,1,99);
b=repmat(-19,1,100);
c=repmat(7,1,99);
f1=repmat(1.1,1,98);
f=[0 f1 1];
[times,wucha]=zhuiganfa(a,b,c,f)
结果:
times =
297
wucha =
2.599344850768740e-14
程序:function [count,wucha] = zhouqisanduijaiozhuiganfa(a,b,c,f)
%x为所求解,count为所有乘除运算次数
n=length(f);
x=zeros(n,1);
y=zeros(n,1);
count=0;
l=zeros(1,n-2);
s=zeros(1,n-1);
u=zeros(1,n);
t=zeros(1,n-1);
u(1)=b(1);t(1)=1;
s(1)=1/u(1);y(1)=f(1);
for i=2:n-1
l(i-1)=a(i-1)/u(i-1);
u(i)=b(i)-l(i-1)*c(i-1);
t(i)=-l(i-1)*t(i-1);
s(i)=-s(i-1)*c(i-1)/u(i);
y(i)=f(i)-l(i-1)*y(i-1);
count=count+4;
end
st=0;
for k=1:n-1
st=st+s(k)*t(k);
count=count+1;
end
sy=0;
for k=1:n-2
sy=sy+s(k)*y(k);
count=count+1;
end
u(n)=b(n)-st-s(n-1)*(c(n-1)+t(n-1));
y(n)=f(n)-sy;
x(n)=y(n)/u(n);
for i=n-1:-1:1
x(i)=(y(i)-c(i)*x(i+1)-t(i)*x(n))/u(i);
count=count+1;
end
A=zeros(n,n);
A=diag(b,0)+diag(a,-1)+diag(c,1);
A(n,1)=1;
A(1,n)=1;
wucha=norm((A*x-f'),2);
命令行:
n=10;
a=repmat(11,1,n-1);b=repmat(-19,1,n);
c=repmat(7,1,n-1);f1=repmat(1.1,1,n-2);f=[0 f1 1];
[count,wucha]= zhouqisanduijaiozhuiganfa(a,b,c,f)
运行结果:
count =
58
wucha =
4.525439045433075
n=30
count =
198
wucha =
5.951269557941316
n=100
count =
688
wucha =
5.993271932634396
高斯消去、追赶法 matlab的更多相关文章
- 高斯混合模型(GMM)及MATLAB代码
之前在学习中遇到高斯混合模型,卡了很长一段时间,在这里记下学习中的一些问题以及解决的方法.希望看到这篇文章的同学们对高斯混合模型能有一些基本的概念.全文不废话,直接上重点. 本文将从以下三个问题详解高 ...
- matlab 求解线性方程组之LU分解
线性代数中的一个核心思想就是矩阵分解,既将一个复杂的矩阵分解为更简单的矩阵的乘积.常见的有如下分解: LU分解:A=LU,A是m×n矩阵,L是m×m下三角矩阵,U是m×n阶梯形矩阵 QR分解: 秩分解 ...
- [zz] 混合高斯模型 Gaussian Mixture Model
聚类(1)——混合高斯模型 Gaussian Mixture Model http://blog.csdn.net/jwh_bupt/article/details/7663885 聚类系列: 聚类( ...
- matlab中fft快速傅里叶变换
视频来源:https://www.bilibili.com/video/av51932171?t=628. 博文来源:https://ww2.mathworks.cn/help/matlab/ref/ ...
- 【原创】开源Math.NET基础数学类库使用(16)C#计算矩阵秩
本博客所有文章分类的总目录:[总目录]本博客博文总目录-实时更新 开源Math.NET基础数学类库使用总目录:[目录]开源Math.NET基础数学类库使用总目录 上个月 ...
- C++线性方程求解
介绍 程序SolveLinearEquations解决联立方程.该方案需要一个文本文件,其中包含输入和输出方程解决.这个项目是几年前我写在C#中http://www.codeproject.com/A ...
- Libsvm学习
本篇博客转自 http://www.cppblog.com/guijie/archive/2013/09/05/169034.html 在电脑文件夹E:\other\matlab 20 ...
- 开源Math.NET基础数学类库使用(16)C#计算矩阵秩
原文:[原创]开源Math.NET基础数学类库使用(16)C#计算矩阵秩 本博客所有文章分类的总目录:http://www.cnblogs.com/asxinyu/p/4 ...
- minimize.m:共轭梯度法更新BP算法权值
minimize.m:共轭梯度法更新BP算法权值 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ Carl Edward Rasmussen在高斯机器学 ...
随机推荐
- 【IT笔试面试题整理】二叉搜索树转换为双向链表
[试题描述] 将二叉搜索树转换为双向链表 对于二叉搜索树,可以将其转换为双向链表,其中,节点的左子树指针在链表中指向前一个节点,右子树指针在链表中指向后一个节点. 思路一: 采用递归思想,对于二叉搜索 ...
- linux-openvpn
1.安装openvpn 1)安装需要的依赖,需要用到epel源 #yum -y install epel-release 修改epel.repo文件enabled=1 #yum install ea ...
- [日常] Go语言圣经--复合数据类型,数组习题
go语言圣经-复合数据类型 1.以不同的方式组合基本类型可以构造出来的复合数据类型 2.四种类型——数组.slice.map和结构体 3.数组是由同构的元素组成——每个数组元素都是完全相同的类型——结 ...
- Linux进程间的通信方式和原理
进程的概念 进程是操作系统的概念,每当我们执行一个程序时,对于操作系统来讲就创建了一个进程,在这个过程中,伴随着资源的分配和释放.可以认为进程是一个程序的一次执行过程. 进程通信的概念 进程用户空间是 ...
- 排查CentOS7.0的联网情况
1.ifconfig命令. 查看网络配置是否有问题 在/etc/sysconfig/network-scripts/ifcfg-ens33里面配置好网络,记住onboot=on这个选项一定要设置,不然 ...
- POJ3186(KB12-O DP)
Treats for the Cows Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 5801 Accepted: 30 ...
- Code Signal_练习题_Knapsack Light
You found two items in a treasure chest! The first item weighs weight1 and is worth value1, and the ...
- 高性能JavaScript(字符串和正则表达式)
字符串连接 +/+=操作符连接 str += "one" + "two"; 这是常用的连接字符串的方法,它运行的时候会经历下面四个步骤: 1.在内存中创建一个临 ...
- 语义SLAM的数据关联和语义定位(二)Semantic Localization Via the Matrix Permanent
论文假设和单目标模型 这部分想讲一下Semantic Localization Via the Matrix Permanent这篇文章的一些假设. 待求解的问题可以描述为 假设从姿态\(x\)看到的 ...
- 大数据【一】集群配置及ssh免密认证
八月迷情,这个月会对大数据进行一个快速的了解学习. 一.所需工具简介 首先我是在大数据实验一体机上进行集群管理学习,管理五台实验机,分别为master,slave1,slave2,slave3,cli ...