ReSQL.ru
     
powered by simpleCommunicator - 2.0.59     © 2025 Programmizd 02
Мобильная версия    Контакт    Правила    FAQ    Помощь
Целевая тема:
Создать новую тему:
Автор:
Закрыть
Цитировать
Доб. в избранное | Игнор. тему | Прикреп. тему | Пометить прочит. / непрочит. | Фильтр
Форумы / C++ [игнор отключен] [закрыт для гостей] / Метод вращения для решения систем уравнений.
5 сообщений из 5, страница 1 из 1
Метод вращения для решения систем уравнений.
    #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
07.04.2013, 02:32
| Ответить | Цитировать | Написать  
Метод вращения для решения систем уравнений.
    #38216685
pirovindos
Скрыть профиль Поместить в игнор-лист Сообщения автора в теме
pirovindos
Участник
Расставьте, пожалуйста, по коду комментарии, в которых будет отражен смысл выполняемого шаг метода вращений.
Я этот метод в коде не нашел.
...
Рейтинг: 0 / 0
08.04.2013, 10:41
| Ответить | Цитировать | Написать  
Метод вращения для решения систем уравнений.
    #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
08.04.2013, 12:55
| Ответить | Цитировать | Написать  
Метод вращения для решения систем уравнений.
    #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
11.04.2013, 20:42
| Ответить | Цитировать | Написать  
Метод вращения для решения систем уравнений.
    #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
23.04.2013, 22:39
| Ответить | Цитировать | Написать  
5 сообщений из 5, страница 1 из 1
Форумы / C++ [игнор отключен] [закрыт для гостей] / Метод вращения для решения систем уравнений.
Найденые пользователи ...
Разблокировать пользователей ...
Читали тему (0):
Читали форум (0):
Пользователи онлайн (0):
x
x
Закрыть


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