高斯消元法求逆矩阵 matlab,高斯消元法求逆矩阵
有多组测试数据。每组测试数据先输入一个整数n,表示方阵的阶。然后下面输入n阶方阵。输出其逆矩阵。若无逆矩阵,则输出No inverse matrix。
#include
#include
#include
using namespace std;
const double eps = 1e-6;
bool is_zero( const double num )
{
return fabs(num) < eps;
}
void create( double ** & matrix, const int n )
{
matrix = new double* [n];
for ( int i = 0; i < n; ++i )
matrix[i] = new double[n];
}
void input ( double ** matrix, const int n )
{
for ( int i = 0; i < n; ++i )
{
for ( int j = 0; j < n; ++ j )
cin >> matrix[i][j];
}
}
bool inverse ( double ** matrix1, double ** matrix2, const int n )
{
int i, j;
for ( i = 0; i < n; ++ i )
{
for ( j = 0; j < n; ++ j )
{
if ( i == j )
matrix2[i][j] = 1;
else
matrix2[i][j] = 0;
}
}
for ( i = 0; i < n; ++i )
{
int rowmaxpos = i;
for ( j = i + 1; j < n; ++j )
{
if ( matrix1[i][j] > matrix1[i][rowmaxpos] )
rowmaxpos = j;
}
for ( j = i; j < n; ++ j )
{
swap( matrix1[j][rowmaxpos], matrix1[j][i]);
swap( matrix2[j][rowmaxpos], matrix2[j][i]);
}
if ( !is_zero(matrix1[i][i]) )
{
int divisor = matrix1[i][i];
for ( j = i; j < n; ++ j )
{
matrix1[i][j] /= divisor;
matrix2[i][j] /= divisor;
}
for ( j = i + 1; j < n; ++ j )
{
int multiple = matrix1[j][i];
for ( int k = i; k < n; ++ k )
{
matrix1[i][j] -= matrix1[i][k] * multiple;
matrix2[i][j] -= matrix2[i][k] * multiple;
}
}
}
else
return false;
}
return true;
}
void output( double ** matrix, const int n )
{
for ( int i = 0; i < n; ++i )
{
for ( int j = 0; j < n; ++ j )
cout << matrix[i][j] << ' ';
cout<
}
}
void destroy( double ** matrix, const int n )
{
for ( int i = 0; i < n; ++ i )
delete [] matrix[i];
delete [] matrix;
}
int main()
{
int n;
double ** matrix1;
double ** matrix2;
while ( cin >> n )
{
create( matrix1, n );
create( matrix2, n );
input( matrix1, n);
if ( inverse(matrix1, matrix2, n) )
output( matrix2, n );
else
cout << "No inverse matrix" << endl;
destroy( matrix1, n );
destroy( matrix2, n );
}
return 0;
}
高斯消元法求逆矩阵 matlab,高斯消元法求逆矩阵相关推荐
- matlab求曲线极值程序,matlab函数求极值matlab函数求极值.ppt
matlab函数求极值matlab函数求极值 * * 函数的极值 1.一元函数的极值 函数命令:fminbnd 调用格式:[x,feval,exitflag,output]=fminbnd(fun,x ...
- matlab求两点间距离,matlab如何求一个N*2的矩阵的任意两点间的距离?
matlab如何求一个N*2的矩阵的任意两点间的距离? mip版 关注:64 答案:3 悬赏:0 解决时间 2021-02-23 20:55 已解决 2021-02-23 16:37 假设有个矩 ...
- matlab用lu求逆矩阵,Matlab实现——求矩阵的逆(LU分解)
Program ( :Solve By Factorization with Pivoting ) 思路及原理: 就得到: 程序: function X=Ni(A) %Input - A is an ...
- matlab用lu求逆矩阵,Matlab实现——求矩阵的逆(LU分解) | 学步园
Program ( :Solve By Factorization with Pivoting ) 思路及原理: 就得到: 程序: function X=Ni(A) %Input - A is an ...
- matlab中求立方根,matlab实现求复数立方根的函数代码怎么写
z2] a=[a1; a2] 来表示z 和 a=a1 + i a2 然后我们用两个实数z1和z2来写f(z) 所以 f(z: function [z, z2, 我们来求解根, niter]=ComCu ...
- Matlab三元隐函数求极值,matlab用三重循环求一个三元函数的最大值所对应的x1,x2,x3...
用MATLAB实现for循环 t=2;whileS(t)>Pstrong&&t 求一个MATLAB循环语句表示这个矩阵200分 这样的,i和j是内部虚数变量,避免轻易使用.cle ...
- java 多阶矩阵求逆_求N阶矩阵的逆矩阵 java
求N阶矩阵的逆矩阵 java 别妄想泡我 911次浏览 2020年08月01日 16:49 最佳经验 本文由作者推荐 1毫升等于多少升-岂止的意思 import r; public class 求逆矩 ...
- 求n阶矩阵的逆矩阵(C语言实现)
求n阶矩阵的逆矩阵 求行列式的值部分:https://blog.csdn.net/qq_30291317/article/details/114364893?ops_request_misc=& ...
- matlab对多项式求导,matlab中多项式求导
1 0.5 0 -0.5 -1 -1.5 -2 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 4.对比用多项式函数的 polyder 函数及符号函数中的 diff 函数,求导 x2+2x ...
最新文章
- python详细安装教程环境配置-python环境配置详细步骤
- Latex 数学符号表
- 人类大脑每日24小时工作节奏表
- [转]OpenGL库介绍
- 推荐几个干货超多助你成长的前端大佬
- SQL 语句中 where 条件后 写上1=1 是什么意思
- commons-lang常用方法
- android wifi 静态地址设置_安卓手机连接wifi故障的6种常见解决方法
- 数据可视化demo_火出圈的大屏你真的会做吗?这才是老板最爱的可视化大屏
- python 提示框如何顶层显示_python tkinter之顶层菜单、弹出菜单实例
- tomcat体系结构
- python信噪比signaltonoise, SNR
- ubuntu16.04查看opencv版本
- 做管理,“重人情”还是“重制度”?
- python求梯形面积_Python - 教你用Python进行简单的面积计算
- 解决SharePoint 2003的爬网性能问题- 之四
- python语言编程中的感叹号_python 感叹号
- 大觉山漂流,男人的欢笑女人的尖叫
- 【Markdown语法】5分钟快速入门保姆级教程(建议收藏...)
- 数学笔记28——不定式和洛必达法则
热门文章
- 手机订货系统优势是什么?订货软件优势有啥?
- 等待用户输入名字、地点、爱好、根据用户的名字和爱好进行任意显示 如:敬爱可亲的xxxx,最喜欢在xxxx地方干xxxx
- php 解析word文件,php解析word文档
- 《谁说菜鸟不会数据分析》之数据分析
- 山东大学软件学院考试总结(2019-2022)
- chrome 禁止 http 自动跳转 https
- 视频教程-[曾贤志]-Excel数据透视表与SQL技术-Office/WPS
- python-如何解决执行pip install 命令的时候出现 File stdin,line 1 pip install XXX的问题
- 天正服务器不显示,天正软件客户端修改服务器地址
- 杀鸡用用牛刀 scrapy框架爬取豆瓣电影top250信息