解线性方程组(一)——克拉默法则求解(C++)

克拉默法则

解线性方程组最基础的方法就是使用克拉默法则,需要注意的是,该方程组必须是线性方程组。
假设有方程组如下:
{ a 11 x 1 + a 12 x 2 + ⋯ + a 1 n x n = b 1 a 21 x 1 + a 22 x 2 + ⋯ + a 2 n x n = b 2 ⋯ ⋯ ⋯ a n 1 x 1 + a n 2 x 2 + ⋯ + a n n x n = b n \begin{cases} a_{11}x_1+a_{12}x_2+\cdots+a_{1n}x_n=b_1\\ a_{21}x_1+a_{22}x_2+\cdots+a_{2n}x_n=b_2\\ \cdots \qquad \qquad\cdots \qquad \qquad \cdots \\ a_{n1}x_1+a_{n2}x_2+\cdots+a_{nn}x_n=b_n\\ \end{cases} a11x1+a12x2++a1nxn=b1a21x1+a22x2++a2nxn=b2an1x1+an2x2++annxn=bn
将其转换为矩阵形式
A x ⃗ = b ⃗ A\vec{x}=\vec{b} Ax =b
[ a 11 a 12 ⋯ a 1 n a 21 a 22 ⋯ a 2 n ⋮ ⋮ ⋱ ⋮ a m 1 a m 2 ⋯ a m n ] [ x 1 x 2 ⋮ x n ] = [ b 1 b 2 ⋮ b n ] \begin{bmatrix} {a_{11}}&{a_{12}}&{\cdots}&{a_{1n}}\\ {a_{21}}&{a_{22}}&{\cdots}&{a_{2n}}\\ {\vdots}&{\vdots}&{\ddots}&{\vdots}\\ {a_{m1}}&{a_{m2}}&{\cdots}&{a_{mn}}\\ \end{bmatrix} \begin{bmatrix} {x_{1}}\\ {x_{2}}\\ {\vdots}\\ {x_{n}}\\ \end{bmatrix}= \begin{bmatrix} {b_{1}}\\ {b_{2}}\\ {\vdots}\\ {b_n} \end{bmatrix} a11a21am1a12a22am2a1na2namn x1x2xn = b1b2bn
根据克拉默法则有:
x i = ∣ A i ∣ ∣ A ∣ x_i=\frac{|A_i|}{|A|} xi=AAi
A i = [ a 11 ⋯ a 1 i − 1 b 1 a 1 i + 1 ⋯ a 1 n a 21 ⋯ a 2 i − 1 b 2 a 2 i + 1 ⋯ a 2 n ⋮ ⋮ ⋱ ⋮ a n 1 ⋯ a n i − 1 b n a n i + 1 ⋯ a n n ] A_i= \begin{bmatrix} {a_{11}}&{\cdots}&{a_{1i-1}}&{b1}&{a_{1i+1}}&{\cdots}&{a_{1n}}\\ {a_{21}}&{\cdots}&{a_{2i-1}}&{b2}&{a_{2i+1}}&{\cdots}&{a_{2n}}\\ {\vdots}&{\vdots}&{\ddots}&{\vdots}\\ {a_{n1}}&{\cdots}&{a_{ni-1}}&{bn}&{a_{ni+1}}&{\cdots}&{a_{nn}}\\ \end{bmatrix} Ai= a11a21an1a1i1a2i1ani1b1b2bna1i+1a2i+1ani+1a1na2nann

代码

实现方法可以参考行列式求值(C++)、插值(一)——多项式插值(C++)

//克拉默法则求解线性方程组
/*
5x1+2x2-2x3=1
x1+4x2+x3=2
x1-2x2+4x3=-1
*/
#include<iostream>
#include<cmath>
//使用代数余子式进行求解
double determinant_value(double **D,int n)
{if(n==1){return D[0][0];}else if(n==0){throw "error: determinant is empty";}char flag = 1;//符号位double D_value = 0.0;double **D_temp = new double*[n];for (int i = 0; i < n;i++){D_temp [i]= new double[n];for(int j=0;j<n;j++){D_temp[i][j]=D[i][j];}}// 转为上三角for (int i = 0; i < n; i++){if (D_temp[i][i] == 0){int j = i + 1;// 将主对角线上的值转为非0值for (; j < n; j++){if (D_temp[j][i] != 0){double temp;for (int k = 0; k < n; k++){temp = D_temp[i][k];D_temp[i][k] = D_temp[j][k];D_temp[j][k] = temp; // 交换两行}flag = -flag;break;}}if (j == n){return D_value;}}// 将主对角线下面的值转为0for (int j = i + 1; j < n; j++){double temp = D_temp[j][i] / D_temp[i][i];for (int k = i; k < n; k++){D_temp[j][k] -= temp * D_temp[i][k];}}}// 计算行列式的值D_value = 1.0;for (int i = 0; i < n; i++){D_value *= D_temp[i][i];}D_value*=flag;for (int i = 0; i < n;i++){delete[] D_temp[i];}delete [] D_temp;return D_value;
}
//克拉默法则求解线性方程组
void Kramer(double **A,double *b,double *x,int n)
{double **A_i=new double*[n];double A_value=determinant_value(A,n);if(A_value==0){std::cout<<"该方程组不是线性方程组"<<std::endl;exit(0);}for (int i = 0; i < n; i++){A_i[i] = new double[n];for (int j = 1; j < n; j++){A_i[i][j] = A[i][j];}}//下面是为了求D_i,每次只需要修改两列数据for (int i = 0; i < n;i++){if(i==0){for(int j = 0;j < n;j++){A_i[j][0]=b[j];}}else{for(int i2 = 0;i2 < n;i2++){A_i[i2][i-1]=A[i2][i-1];A_i[i2][i]=b[i2];}}// for(int i2 = 0;i2 < n;i2++)// {//     for (int j2 = 0;j2 < n;j2++)//     {//         std::cout<<A_i[i2][j2]<<" ";//     }//     std::cout<<std::endl;// }//求多项式系数x[i] = determinant_value(A_i, n) / A_value;}for (int i = 0; i < n; i++){delete[] A_i[i];}delete[] A_i;
}
int main()
{int n;//矩阵维度std::cout<<"请输入矩阵维度:";std::cin>>n;double *x=new double[n];double **A = new double *[n];std::cout<<"请输入矩阵A的元素:"<<std::endl;for(int i = 0;i < n;i++){A[i] = new double [n];for(int j = 0;j < n;j++){std::cin>>A[i][j];}}double *b = new double [n];std::cout<<"请输入矩阵b的元素:"<<std::endl;for(int i = 0;i < n;i++){std::cin>>b[i];}Kramer(A, b, x, n);for (int i = 0; i < n; i++){std::cout<<"x"<<i+1<<"="<<x[i]<<" ";}for (int i = 0; i < n; i++){delete[] A[i];}delete[] A;delete [] b;delete[] x;return 0;
}

结果分析

运行结果如下:
在这里插入图片描述

matlab结果如下:
在这里插入图片描述
此方法虽然可以求解,但是存在一个问题,就是当矩阵维度较大时,程序运行时间会很长,所以在一些精度要求不高的地方,不需要使用此方法。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/685488.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

[前端开发] CSS基础知识 [下]

上篇:CSS 基础知识 [上] CSS基础知识 [下] CSS 新特性媒体查询雪碧图字体图标 CSS 新特性 圆角 (border-radius) 通过 border-radius 属性为元素添加圆角。类型: border-radius: a b c d: 四个值分别为左上|右上|右下|左下角border-radius: a b c : 三个值分别为左上|右上和左下…

re:从0开始的CSS之旅 18. z-index

1. z-index z-index 属性用于设置元素的叠放层次&#xff0c;属性值可以为一个整数&#xff0c;整数值越大越优先显示 注意&#xff1a;z-index只对开启了定位的元素有效 示例如下&#xff1a; <!DOCTYPE html> <html lang"en"><head><meta…

洛谷C++简单题小练习day13—文字处理软件

day13--文字处理软件--2.16 习题概述 题目描述 你需要开发一款文字处理软件。最开始时输入一个字符串作为初始文档。可以认为文档开头是第 0 个字符。需要支持以下操作&#xff1a; 1 str&#xff1a;后接插入&#xff0c;在文档后面插入字符串 str&#xff0c;并输出文档的…

ubuntu20修改xorg.conf实现双屏幕输出

如果显卡工作正常是不需要自己手动编写xorg.conf的&#xff08;这个文件一般不存在或者是空的&#xff09;&#xff0c;系统会根据xorg.conf的缺省自动设置屏幕。 但有时候有的屏幕输出不对&#xff0c;想手动固定一下配置。比如我的ROG想设置内屏用intel驱动&#xff08;集显…

HTTP请求的构造方式

前言&#xff1a; 在详解完HTTP协议的请求和响应格式以后&#xff0c;对HTTP协议就会有更深层次的了解。接下来就要了解HTTP协议的具体的一些用法基础。 在HTTP协议中有两个核心的作用&#xff1a;&#xff08;1&#xff09;如何让客户端构造一个HTTP请求&#xff1b;&#xff…

macOS 安装 conda

macOS 安装 conda 安装 conda参考 Conda是一个开源的软件包管理系统和环境管理系统&#xff0c;用于安装和管理软件包和其依赖项。 安装 conda mkdir miniconda3 cd miniconda3 bash Miniconda3-latest-MacOSX-x86_64.sh$ conda list参考 macOS 安装 conda开始使用conda

OpenAI Sora是世界模型?

初见Sora&#xff0c;我被OpenAI的野心震撼了。 他们不仅想教会AI理解视频&#xff0c;还要让它模拟整个物理世界&#xff01;这简直是通用人工智能的一大飞跃。 但当我深入了解后&#xff0c;我发现Sora比我想象的更复杂、更强大。 Sora不是简单的创意工具&#xff0c;而是…

8086指令小结

所有指令 &#xff08; 1 &#xff09;立即数不能作为目的操作数。 &#xff08; 2 &#xff09;不能在 2 个存储单元之间直接进行操作&#xff08;串操作除外&#xff09; 。 &#xff08; 3 &#xff09; MOV 指令和堆栈指令是惟一能对段寄存器进行操作的指令。 &…

【图像分割 2023】BRAU-Net++

【图像分割 2023】BRAU-Net 论文题目&#xff1a;BRAU-Net: U-Shaped Hybrid CNN-Transformer Network for Medical Image Segmentation 中文题目&#xff1a; 论文链接&#xff1a;[2401.00722] BRAU-Net: U-Shaped Hybrid CNN-Transformer Network for Medical Image Segment…

人工智能专题:通过AI转变保险(英译中)

今天分享的是人工智能系列深度研究报告&#xff1a;《人工智能专题&#xff1a;通过AI转变保险&#xff08;英译中&#xff09;》。 &#xff08;报告出品方&#xff1a;VIEWPOINT&#xff09; 在新时代释放数据的力量 在数据和人工智能 &#xff08; AI &#xff09; 融合的…

Leetcode 42. 接雨水

题意理解&#xff1a; 给定 n 个非负整数表示每个宽度为 1 的柱子的高度图&#xff0c;计算按此排列的柱子&#xff0c;下雨之后能接多少雨水。 左边的柱子和右边的柱子形成围栏&#xff0c;可以使中间能够积水 求最大的积水面积。h*w 解题思路&#xff1a; 1.横向求解 这里的单…

vscode写MATLAB配置

vscode写MATLAB python下载 官网说明Versions of Python Compatible with MATLAB Products by Release - MATLAB & Simulink 不确定这三列都表示什么意思&#xff0c;尽量安装这三列都有的python版本吧&#xff0c;我安装的 MATLAB R2023b,python选择的是3.11.5 …

政安晨:【示例演绎】【Python】【Numpy数据处理】快速入门(三)—— 数组的操作

准备 这是Numpy数据处理的示例演绎系列文章的第三篇&#xff0c;我的前两篇文章为&#xff1a; 政安晨&#xff1a;【示例演绎】【Python】【Numpy数据处理】快速入门&#xff08;一&#xff09;https://blog.csdn.net/snowdenkeke/article/details/136125773 政安晨&#x…

【前端工程化面试题目】webpack 的热更新原理

可以在顺便学习一下 vite 的热更新原理&#xff0c;请参考这篇文章。 首先有几个知识点需要明确 热更新是针对开发过程中的开发服务器的&#xff0c;也就是 webpack-dev-serverwebpack 的热更新不需要额外的插件&#xff0c;但是需要在配置文件中 devServer属性中配置&#x…

计算机设计大赛 深度学习YOLO安检管制物品识别与检测 - python opencv

文章目录 0 前言1 课题背景2 实现效果3 卷积神经网络4 Yolov55 模型训练6 实现效果7 最后 0 前言 &#x1f525; 优质竞赛项目系列&#xff0c;今天要分享的是 &#x1f6a9; **基于深度学习YOLO安检管制误判识别与检测 ** 该项目较为新颖&#xff0c;适合作为竞赛课题方向&…

17.3.1.5 二值化(黑白)

版权声明&#xff1a;本文为博主原创文章&#xff0c;转载请在显著位置标明本文出处以及作者网名&#xff0c;未经作者允许不得用于商业目的。 二值化的算法主要有以下两种&#xff1a; 1、彩色图像灰度化后&#xff0c;根据灰度和阈值来确定颜色是黑色还是白色。通常情况下使…

STM32

Unix时间戳 BKP 读写备份寄存器 main.c #include "stm32f10x.h" // Device header #include "Delay.h" #include "OLED.h" #include "Key.h"uint8_t KeyNum;uint16_t ArrayWrite[] {0x1234,0x2134}; uint16_t Arr…

Nvidia 携手 RTX 推出的本地运行 AI 聊天机器人

每周跟踪AI热点新闻动向和震撼发展 想要探索生成式人工智能的前沿进展吗&#xff1f;订阅我们的简报&#xff0c;深入解析最新的技术突破、实际应用案例和未来的趋势。与全球数同行一同&#xff0c;从行业内部的深度分析和实用指南中受益。不要错过这个机会&#xff0c;成为AI领…

2024.2.16日总结(小程序开发8)

数据监听器 监听对象属性的变化 数据监听器支持监听对象中单个或多个属性的变化 纯数据字段 什么是纯数据字段 纯数据字段指的是哪些不用于页面渲染的data字段 应用场景:例如有些情况下&#xff0c;某些 data 中的字段既不会展示在界面上&#xff0c;也不会传递给其他组件…

html从零开始9:javaScript简介,语句、标识符,变量,JavaScript引入到文件【搬代码】

javaScript简介 javaScript语句、标识符 变量 var num 10; var就是固定声明,num就是变量名&#xff0c;10就是变量&#xff1b;<!DOCTYPE html> <html lang"en"> <head><meta charset"UTF-8"><meta http-equiv"X-UA-Comp…