Гость
Войти | Профиль | Очистить
Целевая тема:
Создать новую тему:
Автор:
Действия ...
Форумы / C++ [игнор отключен] [закрыт для гостей] / Метод вращения для решения систем уравнений. / 5 сообщений из 5, страница 1 из 1
07.04.2013, 02:32
    #38215804
kab18
kab18
Гость
Скрыть профиль Поместить в игнор-лист Сообщения автора в теме
Метод вращения для решения систем уравнений.
Пытаюсь самостоятельно написать по методичке, функцию для решения систем линейных уравнений методом вращений. Помогите найти ошибку в коде, потому что получаемые корни значительно отличаются от тех что я получил методом Гауса, в котором я уверен.
Код: plaintext
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
29.
30.
double __fastcall Rotations(double**A,double*X,unsigned int i)
{
//i-порядок системы уравнений
double C,S;
for(unsigned int n=1;n<i;n++)
{
C=A[0][0]/(sqrt((A[0][0]*A[0][0])+(A[n][0]*A[n][0])));
S=A[n][0]/(sqrt((A[0][0]*A[0][0])+(A[n][0]*A[n][0])));
for(unsigned int m=0;m<i+1;m++)
{
A[0][m]=C*A[0][m]+S*A[n][m];
A[n][m]=-S*A[0][m]+C*A[n][m];
}
}
X[i-1]=A[i-1][i]/A[i-1][i-1];
for(int n=i-2;n>=0;n--)
{
X[n]=A[n][i];
for(int m=i;m>=0;m--)
{
X[n]-=A[n][m]*X[m];
}
X[n]/=A[n][n];
}
for(unsigned int k=0;k<i;k++)
{
cout<<'\n'<<X[k];
}
return 0;
}
...
Рейтинг: 0 / 0
| Ответить | Цитировать | Написать  
08.04.2013, 10:41
    #38216685
pirovindos
pirovindos
Участник
Скрыть профиль Поместить в игнор-лист Сообщения автора в теме
Метод вращения для решения систем уравнений.
Расставьте, пожалуйста, по коду комментарии, в которых будет отражен смысл выполняемого шаг метода вращений.
Я этот метод в коде не нашел.
...
Рейтинг: 0 / 0
| Ответить | Цитировать | Написать  
08.04.2013, 12:55
    #38216997
mayton
mayton
Участник
Скрыть профиль Поместить в игнор-лист Сообщения автора в теме
Метод вращения для решения систем уравнений.
kab18
Код: plaintext
1.
2.
C=A[0][0]/(sqrt((A[0][0]*A[0][0])+(A[n][0]*A[n][0])));
S=A[n][0]/(sqrt((A[0][0]*A[0][0])+(A[n][0]*A[n][0])));


Тут избыточно как-то. Лишние расчёты. Нет это не страшно просто видно что формула
переписана в лоб. Бездумно как-то.
...
Рейтинг: 0 / 0
| Ответить | Цитировать | Написать  
11.04.2013, 20:42
    #38222683
kab18
kab18
Гость
Скрыть профиль Поместить в игнор-лист Сообщения автора в теме
Метод вращения для решения систем уравнений.
Вообщем переделал код, теперь находит корни коректно. Может у кого-то будут советы по его оптимизации.
Код: plaintext
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
29.
int __fastcall Rotations(double** A,double* X,unsigned int i)
{
double C,S,A1;
for(unsigned int k=0;k<i;k++)
{
for(unsigned int n=k+1;n<i;n++)
{
C=A[k][k]/(sqrt((A[k][k]*A[k][k])+(A[n][k]*A[n][k])));
S=A[n][k]/(sqrt((A[k][k]*A[k][k])+(A[n][k]*A[n][k])));
for(unsigned int m=0;m<i+1;m++)
{
A1=A[k][m];
A[k][m]=C*A1+S*A[n][m];
A[n][m]=-S*A1+C*A[n][m];
}
}
}
X[i-1]=A[i-1][i]/A[i-1][i-1];
for(int n=i-2;n>=0;n--)
{
X[n]=A[n][i];
for(unsigned int m=n+1;m<i;m++)
{
X[n]-=A[n][m]*X[m];
}
X[n]/=A[n][n];
}
return 0;
}
...
Рейтинг: 0 / 0
| Ответить | Цитировать | Написать  
23.04.2013, 22:39
    #38237411
mayton
mayton
Участник
Скрыть профиль Поместить в игнор-лист Сообщения автора в теме
Метод вращения для решения систем уравнений.
Как-то так вобщем. Проверяй сам. Мне лениво.
Код: plaintext
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
int __fastcall Rotations(double** A,double* X,unsigned int i)
{
  double C,S,A1;
  for(unsigned int k=0;k<i;k++){
   for(unsigned int n=k+1;n<i;n++){
     FuckenDenominator=(sqrt((A[k][k]*A[k][k])+(A[n][k]*A[n][k])));
     C=A[k][k]/FuckenDenominator;
     S=A[n][k]/FuckenDenominator;
     for(unsigned int m=0;m<i+1;m++){
       A1=A[k][m];
       FuckenMatrixElement=A[n][m];
       A[k][m]=C*A1+S*FuckenMatrixElement;
       A[n][m]=-S*A1+C*FuckenMatrixElement;
     }
   }
  }
  X[i-1]=A[i-1][i]/A[i-1][i-1];
  for(int n=i-2;n>=0;n--){
    //X[n]=A[n][i];
    FuckenWriteCachedElement=A[n][i];
    for(unsigned int m=n+1;m<i;m++){
       FuckenWriteCachedElement-=A[n][m]*X[m];
    }
    FuckenWriteCachedElement/=A[n][n];
    X[n]=FuckenWriteCachedElement;
  }
  return 0;
}
...
Рейтинг: 0 / 0
| Ответить | Цитировать | Написать  
Форумы / C++ [игнор отключен] [закрыт для гостей] / Метод вращения для решения систем уравнений. / 5 сообщений из 5, страница 1 из 1
Найденые пользователи ...
Разблокировать пользователей ...
Читали тему (0):
Читали форум (0):
Пользователи онлайн (0):
x
x
Закрыть


Просмотр
0 / 0
Close
Debug Console [Select Text]