Вычислить Модуль Вектора |x|

Тема в разделе "C/C++/C#", создана пользователем Ruslan737, 24 май 2012.

  1. Ruslan737

    Ruslan737 New Member

    Регистрация:
    4 апр 2012
    Сообщения:
    4
    Симпатии:
    0
    Найти вектор X как решение системы уравнений A(n,n)X(n)=B(n). Вычислить мо-
    дуль вектора |X|.
    Задание: Необходимо написать и отладить программу с обязательным применением динамических матриц.
    Есть какие-то наброски,помогите исправить;)
    Код (C++):
    float* res(float* *a1, float* b1, int n, float tol)
    {
    float* *a=new float*[n];
    float* b=new float[n];
    for(int i=0; i<n; i++)
    {
    b[i]=b1[i];
    a[i]=new float[n];
    for(int j=0; j<n; j++)
    a[i][j]=a1[i][j];
    }
    for(int i=0; i<n-1; i++)//Прямой проход; a1 -- матрица, tol -- точность, b1 -- свободные члены
    {
    if(abs(a[i][i])>tol)
    {
    for(int j=i+1; j<n; j++)
    {
    float k1=a[j][i]/a[i][i];
    for(int k=0; k<n; k++)
    a[j][k]-=k1*a[i][k];
    b[j]-=k1*b[i];
    }
    }
    else
    {
    bool flag=false;
    for(int j=i+1; j<n; j++)
    if(abs(a[j][i])>tol)
    {
    flag=true;
    float k1;//Буффер
    for(int k=0; k<n; k++)
    {
    k1=a[i][k];
    a[i][k]=a[j][k];
    a[j][k]=k1;
    }
    k1=b[i];
    b[i]=b[j];
    b[j]=k1;
    break;
    }
    if(flag==false)
    {
    return NULL;//Решений нет
    }
    i--;
    }
    }
    if(abs(a[n-1][n-1])<tol)
    return NULL;
    for(int i=n-1; i>0; i--)//Обратный проход
    {
    for(int j=i-1; j>=0; j--)
    {
    float k1=a[j][i]/a[i][i];
    for(int k=0; k<n; k++)
    a[j][k]-=k1*a[i][k];
    b[j]-=k1*b[i];
    }
    }
    for(int i=0; i<n; i++)//Привели матрицу к единичной
    b[i]/=a[i][i];
    delete a;
    return b;
    }
     
Загрузка...

Поделиться этой страницей