电脑知识|欧美黑人一区二区三区|软件|欧美黑人一级爽快片淫片高清|系统|欧美黑人狂野猛交老妇|数据库|服务器|编程开发|网络运营|知识问答|技术教程文章 - 好吧啦网

您的位置:首頁技術文章
文章詳情頁

如何基于java實現Gauss消元法過程解析

瀏覽:6日期:2022-08-22 17:50:15

補充知識:

正定矩陣

奇異矩陣

嚴格對角占優

要理解Gauss消去法,首先來看一個例子:

如何基于java實現Gauss消元法過程解析

從上例子可以看出,高斯消去法實際上就是我們初中學的階二元一次方程組,只不過那里的未知數個數$n=2$

$n>2$時,Gauss消去法的思路實際上和解二元一次方程組是一樣的,方法如下:

將n方程組中的n−1個方程通過消元,形成一個與原方程組等價的一個新方程組,新方程組中的n−1個方程僅包含n−1個未知數。 故問題就轉化為了求解n−1元的方程組,這樣我們可以繼續消元,以次類推,直到最后一個方程組為一元一次方程組 從最后一個一元一次方程組求解出最后一個未知量,然后逐步回代入之前的方程組,從而得到所有的未知數。 我們可以看到Gauss實際上就分為兩步:消去和回代

下面通過一般化得到Gauss消元法的求解過程

如何基于java實現Gauss消元法過程解析

如何基于java實現Gauss消元法過程解析

如何基于java實現Gauss消元法過程解析

以上就是Gauss消去法的基本步驟,我們再回過頭看看有沒有什么問題?

我們在求比例$l_{ik}= frac{a_{ik}^{left (k-1 right )}}{a_{kk}^{left (k-1 right )}}$時,如果分母很小,即:

如何基于java實現Gauss消元法過程解析

$l_{ik}rightarrow infty$,那么

如何基于java實現Gauss消元法過程解析

總結一下,能否使用Gauss消元法的情況

如何基于java實現Gauss消元法過程解析

為了解決這個問題,我們可以使用列主元Gauss消元法。

如何基于java實現Gauss消元法過程解析

參考了一些網上的代碼,這里給出Gauss的Java實現

package peterxiazhe;import java.util.Scanner;public class Gauss { /** * 列主元高斯消去法 */ static double A[][]; static double b[]; static double x[]; static int n; //n表示未知數的個數 static int n_2; //記錄換行的次數 public static void main(String[] args) { System.out.println('--------------輸入方程組未知數的個數---------------'); Scanner sc = new Scanner(System.in); n = sc.nextInt();A = new double[n][n]; b = new double[n]; x = new double[n];System.out.println('--------------輸入方程組的系數矩陣A:---------------'); for(int i = 0; i < n; i++) { for(int j = 0; j < n; j++) {A[i][j] = sc.nextDouble(); } }System.out.println('--------------輸入方程組的常量向量b:---------------'); for(int i = 0; i < n; i++) {b[i] = sc.nextDouble(); }Elimination(); BackSubstitution(); PrintRoot(); } //消元法 public static void Elimination() { PrintA(); for(int k = 0; k < n; k++) { WrapRow(k); for(int i = k+1; i < n; i++) {double l = A[i][k] / A[k][k];A[i][k] = 0;for(int j = k+1; j < n; j++) { A[i][j] = A[i][j] - l * A[k][j];}b[i] = b[i] - l * b[k]; } //System.out.println('第' + k + '次消元后:'); //PrintA(); } } //回代法 public static void BackSubstitution() { x[n-1] = b[n-1] / A[n-1][n-1]; for(int i = n - 2; i >= 0; i--) { x[i] = (b[i] - solve(i)) / A[i][i]; } } public static double solve(int i) { double result = 0.0; for(int j = i; j < n; j++) result += A[i][j] * x[j]; return result; } //輸出方程組的根 public static void PrintRoot() { System.out.println('--------------方程組的根為---------------'); for(int i = 0; i < n; i++) { System.out.println('x' + (i+1) + ' = ' + x[i]); } } //交換Swap函數??? public static void Swap(double[] ar, int x, int y) { Double tmp = ar[x]; ar[x] = ar[y]; ar[y] = tmp; } public static void PrintA() { //輸出A的增廣矩陣 //System.out.println('--------------增廣矩陣---------------'); for(int i = 0; i < n; i++) { for(int j = 0; j < n; j++) {System.out.print(A[i][j] + ' '); } System.out.println(b[i]); } } //交換矩陣的行 public static void WrapRow(int k) { //k表示第k+1輪消元 double maxElement = Math.abs(A[k][k]);int WrapRowIndex = k; // 記住要交換的行 for(int i = k + 1; i < n; i++) { if (Math.abs(A[i][k]) > maxElement) {WrapRowIndex = i;maxElement = A[i][k]; } } if (WrapRowIndex != k) { //交換求得最大主元 n_2 += 1; System.out.println('k = ' + k + '時,' + '要交換的行為' + k + '和'+ WrapRowIndex); //先交換A for(int j = k; j < n; j++) {double[] arr = {A[k][j], A[WrapRowIndex][j]};Swap(arr, 0, 1);A[k][j] = arr[0]; A[WrapRowIndex][j] = arr[1];//double tmp = A[k][j];//A[k][j] = A[WrapRowIndex][j];//A[WrapRowIndex][j] = tmp; } //再交換b double[] arr = {b[k], b[WrapRowIndex]}; Swap(arr, 0, 1); b[k] = arr[0]; b[WrapRowIndex] = arr[1];// double tmp = b[k];// b[k] = b[WrapRowIndex];// b[WrapRowIndex] = tmp; System.out.println('--------------交換后---------------'); PrintA(); } }}

注意:由于Java不支持對基本數據類型的引用傳遞,這里使用了一個小技巧

java中交換兩個基本數據類型的變量函數swap(int[] source,int i,int j)

java中函數的參數傳遞機制是:基本數據類型采用值傳遞,對象采用傳引用。因此,如果要寫一個交換兩個int型變量數值的函數,還真是有點不方便,必須采用一個數組對象來作為輔助,具體實現如下:

//交換兩個整數 private static void swap(int[] source, int i, int j) { int temp = source[i]; source[i] = source[j]; source[j] = temp; }

以上就是本文的全部內容,希望對大家的學習有所幫助,也希望大家多多支持好吧啦網。

標簽: Java
相關文章:
主站蜘蛛池模板: 中医中药治疗血小板减少-石家庄血液病肿瘤门诊部 | 运动木地板_体育木地板_篮球馆木地板_舞台木地板-实木运动地板厂家 | 空气能暖气片,暖气片厂家,山东暖气片,临沂暖气片-临沂永超暖通设备有限公司 | 南京试剂|化学试剂|分析试剂|实验试剂|cas号查询-专业60年试剂销售企业 | 氟塑料磁力泵-不锈钢离心泵-耐腐蚀化工泵厂家「皖金泵阀」 | 工程管道/塑料管材/pvc排水管/ppr给水管/pe双壁波纹管等品牌管材批发厂家-河南洁尔康建材 | 胶水,胶粘剂,AB胶,环氧胶,UV胶水,高温胶,快干胶,密封胶,结构胶,电子胶,厌氧胶,高温胶水,电子胶水-东莞聚力-聚厉胶粘 | 板式换网器_柱式换网器_自动换网器-郑州海科熔体泵有限公司 | 卡诺亚轻高定官网_卧室系统_整家定制_定制家居_高端定制_全屋定制加盟_定制家具加盟_定制衣柜加盟 | 臭氧老化试验箱,高低温试验箱,恒温恒湿试验箱,防水试验设备-苏州亚诺天下仪器有限公司 | 深圳昂为官网-气体分析仪,沼气分析仪,动态配气仪,气体传感器厂家 | 高空重型升降平台_高空液压举升平台_高空作业平台_移动式升降机-河南华鹰机械设备有限公司 | 灰板纸、灰底白、硬纸板等纸品生产商-金泊纸业 | 面粉仓_储酒罐_不锈钢储酒罐厂家-泰安鑫佳机械制造有限公司 | 基业箱_环网柜_配电柜厂家_开关柜厂家_开关断路器-东莞基业电气设备有限公司 | Safety light curtain|Belt Sway Switches|Pull Rope Switch|ultrasonic flaw detector-Shandong Zhuoxin Machinery Co., Ltd | [品牌官网]贵州遵义双宁口腔连锁_贵州遵义牙科医院哪家好_种植牙_牙齿矫正_原华美口腔 | 提升海外网站流量,增加国外网站访客UV,定制海外IP-访客王 | 天津货架厂_穿梭车货架_重型仓储货架_阁楼货架定制-天津钢力仓储货架生产厂家_天津钢力智能仓储装备 | 缓蚀除垢剂_循环水阻垢剂_反渗透锅炉阻垢剂_有机硫化物-郑州威大水处理材料有限公司 | 不锈钢拉手厂家|浴室门拉手厂家|江门市蓬江区金志翔五金制品有限公司 | 钢格栅板_钢格板网_格栅板-做专业的热镀锌钢格栅板厂家-安平县迎瑞丝网制造有限公司 | 河南不锈钢水箱_地埋水箱_镀锌板水箱_消防水箱厂家-河南联固供水设备有限公司 | 切铝机-数控切割机-型材切割机-铝型材切割机-【昆山邓氏精密机械有限公司】 | 缝纫客| 浙江华锤电器有限公司_地磅称重设备_防作弊地磅_浙江地磅售后维修_无人值守扫码过磅系统_浙江源头地磅厂家_浙江工厂直营地磅 | 石牌坊价格石牌坊雕刻制作_石雕牌坊牌楼石栏杆厂家_山东嘉祥石雕有限公司 | 点胶机_点胶阀_自动点胶机_智能点胶机_喷胶机_点胶机厂家【欧力克斯】 | 河南不锈钢水箱_地埋水箱_镀锌板水箱_消防水箱厂家-河南联固供水设备有限公司 | 广东佛电电器有限公司|防雷开关|故障电弧断路器|智能量测断路器 广东西屋电气有限公司-广东西屋电气有限公司 | 垃圾处理设备_餐厨垃圾处理设备_厨余垃圾处理设备_果蔬垃圾处理设备-深圳市三盛环保科技有限公司 | 油漆辅料厂家_阴阳脚线_艺术漆厂家_内外墙涂料施工_乳胶漆专用防霉腻子粉_轻质粉刷石膏-魔法涂涂 | 桑茶-七彩贝壳桑叶茶 长寿茶| 电磁铁_推拉电磁铁_机械手电磁吸盘电磁铁厂家-广州思德隆电子公司 | 同学聚会纪念册制作_毕业相册制作-成都顺时针宣传画册设计公司 | 工业CT-无锡璟能智能仪器有限公司| 绿萝净除甲醛|深圳除甲醛公司|测甲醛怎么收费|培训机构|电影院|办公室|车内|室内除甲醛案例|原理|方法|价格立马咨询 | 高速龙门架厂家_监控杆_多功能灯杆_信号灯杆_锂电池太阳能路灯-鑫世源照明 | 西点培训学校_法式西点培训班_西点师培训_西点蛋糕培训-广州烘趣西点烘焙培训学院 | 西点培训学校_法式西点培训班_西点师培训_西点蛋糕培训-广州烘趣西点烘焙培训学院 | 无痕胶_可移胶_无痕双面胶带_可移无痕胶厂家-东莞凯峰 |