Форум программистов, компьютерный форум CyberForum.ru
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Polynomyj
0 / 0 / 0
Регистрация: 19.12.2011
Сообщений: 5
#1

Метод релаксации неправильно работает - C++

15.04.2012, 11:52. Просмотров 959. Ответов 0
Метки нет (Все метки)

Мне необходимо решить уравнение Пуассона методом релаксации для области во вложении. (Она целиком лежит в квадрате 3х3). Я написал код для метода релаксации, однако он даёт неправильный результат. Помогите найти ошибку.

N - количество частей, Omega - параметр релаксации, eps - точность, phi - матрица начальных приближений, f - функция распределения, U - нужная матрица.

C++
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
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
void Relax(int N, double Omega, double eps, double phi[LEN][LEN], PF f, double U[LEN][LEN])
{
 double h = 3/(double)N;
 double ro[LEN][LEN];
 double c[LEN][LEN];
 for (int i = 0; i < N; i++)
  {
   double x = i * h;
   for (int j = 0; j < N; j++)
    {
      double y = j * h;
      if (x >= 1 && x <= 2)
       {
        if (y >= 1 && y <= 2)
         ro[i][j] = f(x,y);
        else
         ro[i][j] = 0;
       }
      else ro[i][j] = f(x,y);
    }
  }
 
  for (int i = 1; i < N-1; i++)
    for (int j = 1; j < N-1; j++)
     {
      if (i * h >= 1 && i * h <= 2)
       if (j*h >= 1 && j*h <= 2)
        U[i][j] = (1-Omega) * phi[i][j] + (Omega * (phi[i+1][j] + phi[i-1][j] + phi[i][j+1] + phi[i][j-1] + h*h*ro[i][j]))/4;
       else U[i][j] = 0;
      else  U[i][j] = (1-Omega) * phi[i][j] + (Omega/4) * (phi[i+1][j] + phi[i-1][j] + phi[i][j+1] + phi[i][j-1] + h*h*ro[i][j]);
     }
  Diff(U,phi,c,N);
 
  int k = 0;
 
  while (Norm1(c,N) >= eps)
  {
   for (int i = 0; i < N; i++)
    for (int j = 0; j < N; j++)
     phi[i][j] = U[i][j];
 
  for (int i = 1; i < N-1; i++)
    for (int j = 1; j < N-1; j++)
     {
      if (i * h >= 1 && i * h <= 2)
       if (j*h >= 1 && j*h <= 2)
        U[i][j] = (1-Omega) * phi[i][j] + (Omega/4) * (phi[i+1][j] + phi[i-1][j] + phi[i][j+1] + phi[i][j-1] + h*h*ro[i][j]);
       else U[i][j] = 0;
      else  U[i][j] = (1-Omega) * phi[i][j] + (Omega/4) * (phi[i+1][j] + phi[i-1][j] + phi[i][j+1] + phi[i][j-1] + h*h*ro[i][j]);
     }
 
   Diff(U,phi,c,N);
   Form1->ProgressBar1->Position = ++k;
  }
}
0
Изображения
 
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
15.04.2012, 11:52
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Метод релаксации неправильно работает (C++):

неправильно работает метод класа - C++
подскажите пожалуйста написал часть кода, он ище далек от окончания, но что у меня щас уже не нравильно работает, мой метод класса...

Метод Зейделя, работает неправильно - C++
сообственно проблема в том что считает немного неточно. и какую функцию можно использовать для закрытия программы?(exit не работает) ...

Проверьте задачку по циклам, неправильно работает. [думаю что неправильно] - C++
Спасибо что решили зайти. Задание выгладит так: http://*******/PW95p А результат выплнения: http://*******/KwhuS #include...

Метод золотого сечения считается неправильно - C++
#include &quot;stdafx.h&quot; #include &quot;stdio.h&quot; #include &quot;conio.h&quot; #include &quot;math.h&quot; #include &lt;iostream&gt; float Afunc(float x1) { ...

Метод простой итерации (неправильно выдает корни уравнения) - C++
Привет всем)) Помогите, пожалуйста) В системе уравнений неправильно выдает корни, уже и не знаю, что делать((( Почему-то, вместо корней...

Условие работает неправильно - C++
string getNumberEmployeeEtype(string x, int temp){ /*первым аргументом является строка должности сотрудника, вторым - номер сотрудника*/ ...

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
15.04.2012, 11:52
Привет! Вот еще темы с ответами:

Неправильно работает программа - C++
файл complex.h #include &lt;iostream&gt; #include &lt;conio.h&gt; #include &lt;math.h&gt; using namespace std; class Complex // класс...

Цикло do while неправильно работает - C++
В общем, проблема в том, что с цикла начинает работать непонятно как, и я не понимаю почему. При чём, если в task_1() ответить...

Неправильно работает функция - C++
Проблема заключается в том, что в случае, если пользователь вводит числа от 1 до 9 все работает правильно, а если от 10 и далее, то...

программа работает неправильно - C++
Дана символьная матрица размера 13 на 18. Найти: номер первой по порядку строки, содержащий наибольшее число цифр; #include...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru