function qiuni =INV_GET(a)
N=length(a);
M=eye(N);
%得到上三角矩?
for i=1:N
max=a(i,i);
A=i;
for j=i+1:N
if(abs(a(j,i))>abs(max))%找最大值
max=a(j,i);
A=j;
end
end
for m=1:N
temp1=a(i,m);%交换最大值值所在的行和当前行
a(i,m)=a(A,m);
a(A,m)=temp1;
temp2=M(i,m);
M(i,m)=M(A,m);
M(A,m)=temp2;
end
for k=i+1:N
jiaquan=a(k,i)/a(i,i);%其他行加权得到上三角
for n=1:N
a(k,n)=a(k,n)-jiaquan*a(i,n);
M(k,n)=M(k,n)-jiaquan*M(i,n);
end
end

end
temp3=a;
%得到单位矩阵
for s=N:-1:1
xishu1 =a(s,s);
for p=s:N
a(p,s)=a(p,s)/xishu1;
end
for q=1:N
M(s,q)=M(s,q)/xishu1;
end
for w=s-1:-1:1
xishu=a(w,s);
a(w,s)=0;
for t=1:N
M(w,t)=M(w,t)-xishu*M(s,t);
end
end
end
temp2=a;
qiuni=M;

参考的c语言程序和注释

#include<stdio.h>

#include<math.h>

long int const N=1000;     //定义矩阵最大为1000阶

int n;              //n表示矩阵的行数和列数。

double juzhen[N][N]; //定义一个1000阶矩阵

double danwei[N][N];   //定义一个单位矩阵

bool zhaozuidazhi(int s)//定义一个布尔型从s行到n行选择最大的元素作为主元的函数。

{

    int i,j,A;

    double mas,temp;   //temp为中间变量,实现行交换

    mas=fabs(juzhen[s][s]);    //数学函数:fabs  功能:求浮点数juzhen[s][s]的绝对值

    //计算juzhen[s][s], 当juzhen[s][s]不为负时返回juzhen[s][s],否则返回-juzhen[s][s]

    A=s;

    for(i=s+1;i<n;i++)

    {

       if(mas<fabs(juzhen[i][s]))

       {

           mas=fabs(juzhen[i][s]);

           A=i;

       }

    }

    if(mas==0)

       return false;

    //交换两行

    for(j=0;j<n;j++)

    {

       temp=juzhen[s][j];

       juzhen[s][j]=juzhen[A][j];

       juzhen[A][j]=temp;

      

       temp=danwei[s][j];

       danwei[s][j]=danwei[A][j];

       danwei[A][j]=temp;

    }

    return true;

}

void jisuan(int s)  //消元计算

{

    int i,j;

    double mas=juzhen[s][s],r;

    for(i=s+1;i<n;i++)

    {

       r=juzhen[i][s]/mas;

       for(j=0;j<n;j++)

       {    //利用消元计算方阵使之成为上三角矩阵,最后使主对角线上的元素相乘就是最终结果 

           juzhen[i][j]=juzhen[i][j]-juzhen[s][j]*r;

           danwei[i][j]=danwei[i][j]-danwei[s][j]*r;

       }

    }

}

void HH(int s)

{

    int i,j;

    double mas;

    mas=juzhen[s][s];

    for(i=s;i<n;i++)

       juzhen[s][i]=juzhen[s][i]/mas;

    for(i=0;i<n;i++)

       danwei[s][i]=danwei[s][i]/mas;

    for(i=s-1;i>=0;i--)

    {

       mas=juzhen[i][s];

       juzhen[i][s]=0;

        for(j=0;j<n;j++)

           danwei[i][j]=danwei[i][j]-mas*danwei[s][j];

    }

}

int main()

{

    int i,j;

zl: printf("第一行输入矩阵的阶数,然后回车输入原始矩阵:\n");

    scanf("%d",&n);

    for(i=0;i<n;i++)

    {

       for(j=0;j<n;j++)

       {

           if(i==j) danwei[i][j]=1;

           else danwei[i][j]=0;

           scanf("%lf",&juzhen[i][j]);//C语言的函数参数是传值而不是传引用的,

           //通常函数无法修改和操作参数.scanf可以修改参数的原因是传给scanf的

           // 那个参数是一个指针,scanf通过传过来的指针来修改指针指向的内容。所以,

           //把x取地址再传过去。&是取juzhen[i][j]的地址,%lf是说juzhen[i][j]是一个双精度浮点型.

       }

    }

   

    //将原始矩阵转成上三角矩阵

    for(i=0;i<n-1;i++)

    {

       //每一次选一个绝对值最大的值作为a[i]即主元

       if(zhaozuidazhi(i)) //函数调用

       {

           jisuan(i);   //函数调用

 

       }

       else

       {

           printf("该矩阵不可逆");//高斯消去法不能将原始矩阵化为三角形的格式,

                                 //那就代表原始矩阵是一个不可逆的矩阵,

           return 0;

       }

    }

    //将原始矩阵转化成单位矩阵

    for(i=n-1;i>=0;i--)

    {

       HH(i);

    }

         printf("\n原始矩阵的逆矩阵为:\n");//输出逆矩阵语句

         for(i=0;i<n;i++)

         {    

             for(j=0;j<n;j++)

                printf("%.2lf ",danwei[i][j]);//%.2lf在printf()里表示按下面格式

             // 输出一个long double型数:整数部分输出全部输出小数部分输出2位,

             //没有小数的输出两个0,不足两位的后面补0,大于两位的截短到两位。

              //测试结果会出现两位小数

     

             printf("\n");

         }

         goto zl;

}

matlab中求逆矩阵的高斯消元法实现的代码的更多相关文章

  1. MATLAB中求矩阵非零元的坐标

    MATLAB中求矩阵非零元的坐标: 方法1: index=find(a); [i,j]=ind2sub(size(a),index); disp([i,j]) 方法2: [i,j]=find(a> ...

  2. MATLAB中的多项式运算

    作者:长沙理工大学 交通运输工程学院 王航臣 1.多项式求根 在MATLAB中求取多项式的根用roots函数. 函数:roots 功能:一元高次方程求解. 语法:roots(c) 说明:返回一个列向量 ...

  3. matlab中求解线性方程组的rref函数

    摘自:http://www.maybe520.net/blog/987/ matlab中怎么求解线性方程组呢? matlab中求解线性方程组可应用克拉默法则(Cramer's Rule)即通过det( ...

  4. matlab中矩阵的表示与简单操作

    原文地址为:matlab矩阵的表示和简单操作 一.矩阵的表示在MATLAB中创建矩阵有以下规则: a.矩阵元素必须在”[ ]”内: b.矩阵的同行元素之间用空格(或”,”)隔开: c.矩阵的行与行之间 ...

  5. MATLAB中FFT的使用方法

    MATLAB中FFT的使用方法 说明:以下资源来源于<数字信号处理的MATLAB实现>万永革主编 一.调用方法X=FFT(x):X=FFT(x,N):x=IFFT(X);x=IFFT(X, ...

  6. Matlab中函数定义方法

    Matlab自定义函数的六种方法 n1.函数文件+调用函数(命令)文件:需单独定义一个自定义函数的M文件: n2.函数文件+子函数:定义一个具有多个自定义函数的M文件: n3.Inline:无需M文件 ...

  7. MATLAB中白噪声的WGN和AWGN函数的使用

    MATLAB中白噪声的WGN和AWGN函数的使用如下: MATLAB中产生高斯白噪声非常方便,可以直接应用两个函数,一个是WGN,另一个是AWGN.WGN用于产生高斯白噪声,AWGN则用于在某一 信号 ...

  8. matlab中的xcorr 自相关函数

    转载自 http://blog.163.com/to_be_myself/blog/static/176060227201101762159227/ Matlab中用于计算自相关函数的指令是xcorr ...

  9. MATLAB中plot()画图的颜色线型和希腊字母参数设置

    y         黄色           ·             点线      m         粉红           ○             圈线      c          ...

随机推荐

  1. 【Flask】Flask上下文

    # 上下文: ### Local对象:在`Flask`中,类似于`request`的对象,其实是绑定到了一个`werkzeug.local.Local`对象上.这样,即使是同一个对象,那么在多个线程中 ...

  2. Java四则运算总结

    Java结对编程四则运算第二周总结 需求分析 要随机产生规定数量的式子. 要对分式进行计算. 要支持多语言. 要自动判断对错. 设计思路 确定要选用的语言: 要产生随机数: 在产生随机数的时候就要产生 ...

  3. Oracle大总结

    maven的常见两个指令说明 mvn install 是将你打好的jar包安装到你的本地库中,一般没有设置过是在 用户目录下的 .m2\下面.mvn package 只是将你的代码打包到输出目录,一般 ...

  4. Android Studio Gradle Could not reserve enough space for object heap

    Studio 创建第一个工程报错 Error:Unable to start the daemon process.This problem might be caused by incorrect ...

  5. DPDK l3fwd

    l3fwd负责三层转发,比l2fwd要复杂点. /*- * BSD LICENSE * * Copyright(c) 2010-2014 Intel Corporation. All rights r ...

  6. response.setHeader()用法

    response.setHeader()下载中文文件名乱码问题 收藏  1. HTTP消息头 (1)通用信息头 即能用于请求消息中,也能用于响应信息中,但与被传输的实体内容没有关系的信息头,如Data ...

  7. 使用SpringMVC时报错HTTP Status 405 - Request method 'GET' not supported

    GET方法不支持.我出错的原因在于,在JSP中我希望超链接a以post方式提交,但是这里写js代码时出错. <script type="text/javascript"> ...

  8. mongodb笔记(一) 分片 &&文档连接

    版本:mongodb3.4 ; 分片: 工作顺序:router=>config=>shards 一,配置config: 3.4中config必须为replSet.下面配置两个config. ...

  9. .SourceInsight添加.S文件

    在Option->Document Option添加配置.S然后再去添加文件

  10. C#学习历程(二)[基础知识]

    c#中类型的转换 1.Convert.ToInt32(string s) 这个方法的返回值是int类型,要用int类型的变量接收 如: string strNum=Console.ReadLine() ...