文档章节

行列式求解,列主元高斯消元法

quanwei9958
 quanwei9958
发布于 2015/09/05 20:19
字数 240
阅读 717
收藏 0

  高斯消元法的弊端就是针对系数矩阵A,当遇到A(n,n)=0的情况时边无法处理(出现除0错误),以及有效减少计算机在处理浮点运算时出现舍入误差。 

 列主消元法代码:

function  [x]=ext_gauss(A,b)
  n=size(A,1);
  for k=1:n-1
    [value,index]=min(abs(A(k:n,k)));
    index+=k-1;
    if value==0 
      break
    endif
    A([k,index],:)=A([index,k],:);
    b([k,index])=b([index,k]);
    %%gause
    m=A(k+1:n,k)/A(k,k);
    A(k+1:n,k+1:n)-=m*A(k,k+1:n);
    A(k+1:n,k)=zeros(n-k,1);
    b(k+1:n)-=m*b(k);
    %%end gauss
  endfor
  
  x=zeros(n,1);
  x(n)=b(n)/A(n,n);
  for k=n-1:-1:1
    x(k)=(b(k)-A(k,k+1:n)*x(k+1:n))/A(k,k);
  endfor
endfunction

  解决了高斯消元法中存在不能消元的现象。 本程序只能处理满秩系数矩阵。



© 著作权归作者所有

quanwei9958

quanwei9958

粉丝 29
博文 79
码字总数 36250
作品 0
程序员
私信 提问
行列式求解,高斯消元法

程序部分 gauss.m function [x]=gauss(A,b)n=size(A,1);for k=1:n-1 endforx=zeros(n,1);x(n)=b(n)/A(n,n);for k=n-1:-1:1 endforendfunction 参数: A 输入系数矩阵 b 右端项 程序原理就是利......

quanwei9958
2015/09/05
301
0
数论常用内容——高斯消元

高斯消元法 数学上,高斯消元法,是线性代数中的一个算法,可用来为线性方程组求解 高斯消元法求解线性方程组时,首先需要根据方程,列出增广矩阵。然后再利用初等行变换把增广矩阵转换为行阶...

tick_tock97
2017/05/06
0
0
Java 矩阵行列式算法(非高斯消元)

最近由于项目任务较少,手上有不少空闲的时间,所以抽空研究了一下矩阵行列式的算法。 先来说说行列式,以下摘自百度百科: 行列式在数学中,是由解线性方程组产生的一种算式。行列式的特性可...

Acce1erator
2015/12/25
424
0
Gauss Elimination算法分析与实现

高斯消去法分为两个过程:第一步是前向消元(forward elimination),也就是将系数矩阵转化成上三角矩阵的过程;第二步是回代(back substitution)过程,自底向上求解方程组的过程。 选择主...

SVD
2015/11/03
1K
0
克莱姆法则、逆矩阵、体积[MIT线代第二十课]

0、前言 MIT线性代数课程精细笔记[第十八九课]笔记见行列式介绍[MIT线代第十八九课],该笔记是连载笔记,本文由坤博所写,希望对大家有帮助。 一、知识概要 上一节中我们介绍了行列式的求法,...

忆臻
2017/12/15
0
0

没有更多内容

加载失败,请刷新页面

加载更多

jquery底层函数-jQuery.ajaxSetup()

jQuery.ajaxSetup( options )返回: 描述: 为以后要用到的Ajax请求设置默认的值 添加的版本: 1.1jQuery.ajaxSetup( options ) options 类型: PlainObject 一个用来配置Ajax请求的"{键:值}"对,......

前端老手
29分钟前
5
0
Qt编写自定义控件56-波浪曲线

一、前言 波浪曲线控件,其实是之前一个水波进度条控件的一个核心,其实就是利用正弦曲线来生成对应的坐标进行绘制,把这个功能单独提取出来,是为了更详细的研究各种正弦余弦等拓展效果,当...

飞扬青云
35分钟前
5
0
深入理解计算机系统cp1:存储单位与编码

摘要: 理解计算机是如何存储数据的。 原文:深入理解计算机系统cp1:存储单位与编码 作者:Chor Fundebug经授权转载,版权归原作者所有。 1. 存储单位 位:即 bit,表示二进制位,要么是 0 ...

Fundebug
38分钟前
5
0
Python查询比特币实时价格

在本文中我们将学习如何使用使用coinmarketcap提供的比特币行情API,编写Python程序来获取像比特币、莱特币或以太币之类的区块链数字货币的实时行情/实时价格。 1、比特币行情API调用代码 我...

汇智网教程
今天
8
0
爬虫平台Crawlab核心原理--自动提取字段算法

⚠注意: 可配置爬虫现在仅在Python版本(v0.2.1-v0.2.4)可用,在最新版本Golang版本(v0.3.0)还暂时不可用,后续会加上,请关注近期更新 背景 实际的大型爬虫开发项目中,爬虫工程师会被要...

tikazyq
今天
6
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部