1

Пишу программу для решения системы уравнений с комплексными числами методом Гаусса. Сам алгоритм я реализовал. Встал вопрос с описанием комплексных чисел. Код программы прилагается.

public class Gauss{
static int i,j,n,k,h,c;
double m,s;
static boolean f;
static double[][] buf;
static double[][] matrix;
static double[] b,x;
static void Prisvoi(int n){
    buf = new double[n][n];
    matrix = new double[n][n];
    b = new double[n];  
    x = new double[n];
    k=0; c=0; f=false;
}
void MetodGaussa(){
    //Прямой ход
    for(h=1;h<n;h++){
        if(matrix[k][k]==0){
            for(i=1+k;i<n;i++){if(matrix[i][k]!=0){c=i;}
            if(c!=0){
                for(j=0;j<n;j++){
                    buf[0][j]=matrix[k][j];
                    matrix[k][j]=matrix[c][j];
                    matrix[c][j]=buf[0][j];
                    if(j==n-1){
                        buf[0][j+1]=b[k];
                        b[k]=b[c];
                        b[k]=buf[0][j+1];
                    }
                }
                //buf=nill;
                for(i=1+k;i<n;i++){
                    for(j=0+k;j<n;j++){
                        if(j==0+k){m=matrix[i][k]/matrix[k][k];}
                        matrix[i][j]=matrix[i][j]-matrix[k][j]*m;
                        if(j==n-1){b[i]=b[i]-b[k]*m;}
                    }
                    f=true;
                }
            }
            else{
                System.out.println("Нет максимума в системе уравнений");
                f=false;
            }
            }
        }
        else{
            for(i=1+k;i<n;i++){
                for(j=0+k;j<n;j++){
                    if(j==0+k){m=matrix[i][k]/matrix[k][k];}
                    matrix[i][j]=matrix[i][j]-matrix[k][j]*m;
                    if(j==n-1){b[i]=b[i]-b[k]*m;}
                }
            } 
            f=true;
        }
    }
    //Расчёт корней
    if(f==true){
        if(matrix[n-1][n-1]==0){
            x[n-1]=b[n-1]/matrix[n-1][n-1];
            for(i=n-2;i>-1;i--){
                s= 0;
                for(j=i+1;j<n;j++){s=s-matrix[i][j]*x[j];}
                x[i]=(b[i]+s)/matrix[i][i];
            }
        }
        else{
            if(b[n-1]==0){System.out.println("Система уравнений не имеет решения");}
            else System.out.println("Система уравений имеет бесконечное множество решений");
        }
    }
}
2
  • А в чём проблема? Напишите свой класс Complex, и дело с концом. Правда, Java не позволяет перегружать операторы, так что вместо + у вас будет .add() и т. д. Ну или если лень, возьмите любую из свободных реализаций. Список тут.
    – VladD
    31 мая 2014 в 10:09
  • Спасибо уже разобрался с вопросом
    – seg_pro
    2 июн 2014 в 6:37

1 ответ 1

1

Прекрасная реализация комплексных чисел есть в библиотеке Apache Math (javadoc, учебник, конфиг для maven и других систем, скачать исходники).

Кстати, там есть и реализации алгоритмов решения СЛАУ, причём несколько. Рекомендую ознакомиться хотя бы потому, что код, на мой взгляд, очень качественный.

Приводить код здесь не буду, т.к. кода много. Думаю, в обозримом будущем Apache никуда не денется.

Ваш ответ

By clicking “Отправить ответ”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Всё ещё ищете ответ? Посмотрите другие вопросы с метками или задайте свой вопрос.