精华内容
下载资源
问答
  • java 求解线性方程组 高斯列主元消去 迭代法 LU二分法 import java.util.Scanner; public class Gauss { /** * @列主元高斯消去法 */ static double a[][]; static double b[]; static double x[]; ...
  • java线性方程组 Java线性方程组 方法一:高斯消去法 import java.util.Scanner; public class Gauss { /*** @ 列主元高斯消去法*/ static double a[][]; static double b[]; static double x[]; static int n; ...

    41528d3028836879cd698677c3999917.gifjava解线性方程组

    Java 解线性方程组 方法一:高斯消去法 import java.util.Scanner; public class Gauss { /*** @ 列主元高斯消去法*/ static double a[][]; static double b[]; static double x[]; static int n; static int n2; //记录换行的次数public static void Elimination(){ //消元 for(int k=1;k=1;i--)x[i]=(b[i]-jisuan(i))/a[i][i];}public static double jisuan(int i){double he=0.0;for(int j=i+1;jmax){n1=i;max=Math.abs(a[i][k]);}}if(n1!=k){n2++;System.out.println(“当k=“+k+“ 时,要交换的行是:“+k+“ 和“+n1);for(int j=k;j=1;i--) x[i]=(x[i]-jisuan3(i))/a[i][i]; PrintX(1); } public static void HL(){//计算系数矩阵A 的行列式 double cj=1.0; for(int i=1;i<=n;i++) cj=cj*a[i][i]; System.out.println(“A的行列式:|A| = “+cj); } public static double jisuan(int r,int i){ double Sum=0.0; for(int k=1;k<=r-1;k++) Sum=Sum+a[r][k]*a[k][i]; return Sum; } public static double jisuan1(int r,int i){ double Sum=0.0; for(int k=1;k<=r-1;k++) Sum=Sum+a[i][k]*a[k][r]; return Sum; } public static double jisuan2(int i){ double Sum=0.0; for(int k=1;k<=i-1;k++) Sum=Sum+a[i][k]*x[k];

    展开全文
  • 求解线性方程组的解——java实现

    热门讨论 2009-11-30 19:32:08
    高斯列主消元法 LU分解法 迭代法求解线性方程组 高斯列主消元法 LU分解法 迭代法求解线性方程组 高斯列主消元法 LU分解法 迭代法求解线性方程组 高斯列主消元法 LU分解法 迭代法求解线性方程组 高斯列主消元法 LU...
  • 使用java实现直接法求解线性方程组,使用者可以根据自己需要进行修改
  • 问题补充:class Bird{{System.out.println("b1");}public Bird(){System.out.println("b2");}}class Raptor extends Bird{static{System.out.println("r1");}public Raptor(){System.out.println("r3");...

    问题补充:

    class Bird{

    {

    System.out.println("b1");

    }

    public Bird(){

    System.out.println("b2");

    }

    }

    class Raptor extends Bird{

    static{

    System.out.println("r1");

    }

    public Raptor(){

    System.out.println("r3");

    }

    {

    System.out.println("r2");

    }

    static{

    System.out.println("r4");

    }

    }

    class Hawk extends Raptor{

    public static void main(String[] args) {

    // TODO Auto-generated method stub

    System.out.println("pre");

    new Hawk();

    System.out.println("hawk");

    }

    }

    结果:

    r1

    r4

    pre

    b1

    b2

    r2

    r3

    hawk

    我想问一下r2 和 r3的输出顺序是怎么回事??为什么r3在前面反而是要先输出r2呢?

    展开全文
  • 用高斯消去法求解线性方程组JAVA程序实现
  • 数值分析实验之求解线性方程组实现 java实现的
  • 求解满足以下积分方程的f(kx,ky),k是2d中的向量方程:我将f(kx,ky)映射到2d N * N网格并求解了总N ^ 2个方程组 .from scipy.optimize import newton_krylovimport numpy as nL=10.0 # cut off of integralN=100...

    我求解满足以下积分方程的f(kx,ky),k是2d中的向量

    方程:

    d9c3c3ff-e6ba-4dde-be3e-1be644d501fb.png

    我将f(kx,ky)映射到2d N * N网格并求解了总N ^ 2个方程组 .

    from scipy.optimize import newton_krylov

    import numpy as n

    L=10.0 # cut off of integral

    N=100 # number of devision from 0 to L

    delta=L/(N-1)

    def f0(f): # define a larger matrix used in summation

    f1=n.ones((2*N,2*N))

    for i in range(0, N):

    for j in range(0, N):

    f1[i][j]=f[i][j]

    return f1

    def equation(f):

    f2=n.empty((N,N))

    for i in range(0, N): # i,j index of k

    for j in range(0, N):

    sum=0

    for a in range(-N+1, N): # a,b index of p

    for b in range(-N+1, N):

    if a!=0 or b!=0:

    for c in range(-N+1, N): # c,d index of q

    for d in range(-N+1, N):

    if c!=0 or d!=0:

    if a!=c or b!=d:

    if c!=-i or d!=-j:

    sum=sum-delta**2/(delta**2*(a**2+b**2)\

    +f0(f)[abs(a)][abs(b)])/(delta**2*((a-c)**2+(b-d)**2)\

    +f0(f)[abs(a-c)][abs(b-d)])*(1/(delta**2*((c+i)**2+(d+j)**2)\

    +f0(f)[abs(c+j)][abs(d+j)])-1/(delta**2*(c**2+d**2)\

    +f0(f)[abs(c)][abs(d)]))

    f2[i][j]=sum-f[i][j]

    return f2

    guess =n.ones((N, N), float)

    sol = newton_krylov(equation, guess, method='lgmres', verbose=1)

    print(sol)

    我首先尝试使用非常小的参数l = 1,N = 4 . 解决方案的发现非常缓慢 .

    0: |F(x)| = 23.6576; step 1; tol 0.487986

    .....

    519: |F(x)| = 21.4812; step 3.15471e-09; tol 0.9

    .....

    957: |F(x)| = 20.4225; step 2.03396e-06; tol 0.899997

    .....

    并且还在运行 .

    最初的猜测远非实际的解决方案,也许这就是原因 . 我还想知道在我的情况下使用newton_krylov是否是一个好方法?

    在实际解决方案中,注意f(0)也应该等于零 . 将奇点作为积分排除,求和是否合理?

    还有一个问题,如何在解决方案中强制f(0,0)始终等于零?

    展开全文
  • 《数值分析5-用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组》由会员分享,可在线阅读,更多相关《数值分析5-用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组(6页珍藏版)》请在人人文库网上搜索。1、作业六:...

    《数值分析5-用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组》由会员分享,可在线阅读,更多相关《数值分析5-用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组(6页珍藏版)》请在人人文库网上搜索。

    1、作业六:分别编写用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组Ax=B的标准程序,并求下列方程组的解。可取初始向量 X(0) =(0,0,0);迭代终止条件|x(k+1)-x(k)|efor i=1:nd=A(i,i);if abs(d)100warning(不收敛);endendx=x0;程序结果(1)(2)Gauss-Seidel迭代法:程序clear;clc;%A=8,-1,1;2,10,01;1,1,-5;%b=1;4;3;A=5,2,1;-1,4,2;2,-3,10;b=-12;20;3;m=size(A);if m(1)=m(2)error(矩阵A不是方阵);e。

    2、ndn=length(b);%初始化N=0;%迭代次数L=zeros(n);%分解A=D+L+U,D是对角阵,L是下三角阵,U是上三角阵U=zeros(n);D=zeros(n);G=zeros(n);%G=-inv(D+L)*Ud=zeros(n,1);%d=inv(D+L)*bx=zeros(n,1);for i=1:n%初始化L和Ufor j=1:nif ijU(i,j)=A(i,j);endendendfor i=1:n%初始化DD(i,i)=A(i,i);endG=-inv(D+L)*U;%初始化Gd=(D+L)b;%初始化d%迭代开始x1=x;x2=G*x+d;while norm(x2-x1,inf)10(-6)x1=x2;x2=G*x2+d;N=N+1;endx=x2;程序结果(1) (2。

    展开全文
  • JAVA线性方程组

    2013-12-05 23:23:38
    本文档利用Java编程语言求解线性方程组,不是原创,但很有用
  • 线性方程组问题可以利用矩阵变换求解。利用高斯消元法,将矩阵转换成一个行阶梯矩阵,最后得到一个简化行阶梯矩阵,就是方程的解。参考资料(高斯消元法) Java代码 public class FunctionResolver { public ...
  • 19 std::cout 除了colPivHouseholderQr、LLT、LDLT,还有以下的函数可以求解线性方程组,请注意精度和速度:解小矩阵(4*4)基本没有速度差别 1 //Solve Ax = b. Result stored in x. Matlab: x = A \ b. 2 x = A.ldlt...
  • 线性方程组的迭代求解java.docx 线性方程组的迭代求解线性方程组的迭代求解摘要迭代法是一种逐次逼近方法,在使用迭代法解方程组时,其系数矩阵在计算过程中始终不变。它利用计算机运算速度快、适合做重复性操作的...
  • /**求解线性方程组 * Ax = b ==> PAx = Pb ==>LUx = Pb ==> Ly = Pb ==> Ux = y */ public static float[] LUP_Solve(float[][] L,float[][] U,int[]P,float[] b){ int n = L[0].length,i = 0,j = 0; float...
  • GPS中用于定位的核心算法之一(牛顿迭代+最小二乘),最近在和同事交流时编写了相关的程序,其原理是通过泰勒展开的方式,将非线性方程线性化以方便求解。首先是简单的一元非线性方程:设f(x)为关于x的非线性函数,且...
  • 线性方程组的迭代求解 PAGE 线性方程组的迭代求解 摘要 迭代法是一种逐次逼近方法在使用迭代法解方程组时其系数矩阵在计算过程中始终不变它利用计算机运算速度快适合做重复性操作的特点让计算机对一组指令或一定步骤...
  • 数值分析实验之求解线性方程组java实现的
  • 该楼层疑似违规已被系统折叠隐藏此楼查看此楼请问牛顿迭代法求解三元非线性方程组,不同迭代初值貌似有很多不同结果,如何求解到满足要求的解,0FindRoot[{-6.565159793181527`*^-6 (0.008850684113090817` -0....
  • 二分法根据零点定理判断均分点的选择,依次分下去,能够求解到对应的值优点:计算简单缺点:收敛太慢不动点的迭代法及其收敛性1.不动点及其不动点迭代基本思想:将隐式的方程化成显式的计算公式几何上的理解:(额额...
  • 矩阵分解求解线性方程组:,其解可以表示为.为了提高运算速度,节约存储空间,通常会采用矩阵分解的方案,常见的矩阵分解有LU分解、QR分解、Cholesky分解、Schur分解、奇异分解等。这里简单介绍几种。LU分解:如果...
  • /*使用高斯赛德尔迭代法求解线性方程组*/ public class Gauss_Seidel_Iterate { /*求下三角*/ private static float[][] find_lower(float data[][],int k){ int length=data.length; float data2[][]=new ...
  • package linear_equation;.../*使用高斯赛德尔迭代法求解线性方程组*/public class Gauss_Seidel_Iterate {/*求下三角*/private static float[][] find_lower(float data[][],int k){int length=data.length;float ...
  • 方程组如下图,编写一个可以求解的应用程序,C,C++,java等都可以![图片说明](https://img-ask.csdn.net/upload/201601/10/1452397366_191324.jpg)
  • 单纯形法是解决线性规划问题的一个有效的算法。线性规划就是在一组线性约束条件下,求解目标函数最优解的问题。资源中包含了实现计算单纯形表的Java代码。
  • package linear_equation; import java.util..../*使用高斯赛德尔迭代法求解线性方程组*/ public class Gauss_Seidel_Iterate { /*求下三角*/ private static float[][] find_lower(float data[][],int k){ i
  • **编写程序,提示用户输入a,b,c,e和f,然后...应用Cramer规则求解线性方程组** public class linearEquation { public static void main(String[] args) { java.util.Scanner input = new java.util.Scanner(Syst...
  • 数值计算——一维非线性方程求解

    千次阅读 2017-04-21 23:16:26
    数值计算——一维非线性方程求解 1、二分法  把函数f(x)的零点所在的区间[a,b](满足f(a)●f(b) ...import java.text.... * 二分法求解线性方程组 * @author KeXin * */ public class Bisection { //
  • java写二元一次方程组求解

    千次阅读 2017-02-22 11:20:52
    1.运用克莱姆法则 ...先举一个实例,首先先运用行列式的性质进行适当的化简。...进而得出系数行列式的值不等于零,这就说明这个线性方程组有而且有唯一解。  同样的方法再处理第二行。
  • Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 2354 Accepted Submission(s): 862 Problem Description One day I was shoppi

空空如也

空空如也

1 2 3 4
收藏数 80
精华内容 32
关键字:

java求解线性方程组

java 订阅