Правильно ли реализовано вычисление членов последовательности

@Qazan · Регистрация: 30.04.2013

Студворк — интернет-сервис помощи студентам

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
void SUR(double **a,int n,double w,double *f,double *x_app)
{
 
    //-------------------------------------------
    double * fut = new double[n];
    //-------------------------------------------
    for(int i = 0;i < n;i++)
    {
        fut[i] = 0;//x_app[i];
    }
    //-------------------------------------------
//  for(int k = 0;k < 6;k++)
//  {
        for(int i = 0;i < n;i++)
        {
            fut[i] = 0;//x_app[i];
        }
        for(int i = 0;i <n;i++)
        {
            double t = 0;
            for(int j = i + 1;j < n;j++)
            {
                t += a[i][j]*x_app[j];
            }
 
            double temp = 0;
            for(int j = 0;j < i;j++)
            {
                temp += a[i][j]*fut[j];
            }
            double qs = 1./a[i][i]; 
            fut[i] = (1. - w)*x_app[i] + w*qs*(f[i] - temp - t);
        }
 
        for(int i = 0;i <n;i++)
        {
            x_app[i] = fut[i];
        }
//  }
}

$https://www.cyberforum.ru/cgi-bin/latex.cgi? x^{(k+1)}_i = (1-\omega)x^{(k)}_i + \frac{\omega}{a_{ii}} \left(b_i - \sum_{j<i} a_{ij}x^{(k+1)}_j - \sum_{j>i} a_{ij}x^{(k)}_j \right),\quad i=1,2,\ldots,n.$

Добавлено через 1 минуту
Правильно ли реализовал, математическую запись ?

Мне кажеться ,что все норм ,однако считает ужасно

@outoftime · 20.05.2014, 18:30

Qazan, Что такое k, k+1? Что такое i, j?

@Qazan · 20.05.2014, 19:23 **[ТС]**

outoftime, k, k + 1 это соответственно текущий и предыдущий вектор

Это итерационный метод решения,в начале задаем приближенное значение ,а потом приближаем дальше

i,j - индексы вектора

формируется вектор приближенных значений постепенно

Добавлено через 1 минуту
$https://www.cyberforum.ru/cgi-bin/latex.cgi? a_{ij}$ Элементы матрицы A.

@outoftime · 20.05.2014, 20:09

Сообщение от Qazan

k, k + 1 это соответственно текущий и предыдущий вектор

А теперь распиши что из себя представляют переменные.

@Qazan · 20.05.2014, 21:30 **[ТС]**

outoftime, переменный это эле менты вектора тобишь просто числа

Добавлено через 32 секунды
Successive_over-relaxation

Добавлено через 2 минуты
Successive_over-relaxation

@outoftime · 21.05.2014, 00:35

Qazan, на самом деле, на вики ниже есть алгоритм, и даже там вместо двух циклов используется один с условием (кода на 1 операцию больше, зато читабельней).

Я только не понял как выбирается начальное значение и что значит "повторять пока не сойдется", сам алгоритм выдает результат схождением (по аналогии с разложением функции в ряд быстрого схождения). У меня только одна мысль по этому поводу - делать по аналогии с рядами измеряя эпсилон как разницу между значениями на смежных итерациях.

Если после этого у вас возникают вопросы не касающиеся кода - можете начинать писать код. Если у вас есть явные ошибки в коде - форумчане помогут. Но копаться в этом методе решения линейных уравнений я не хочу (остальные, скорей всего, меня поддержат).

@Qazan · 21.05.2014, 05:37 **[ТС]**

outoftime, про реализацию это да, там есть псевдокод

но мне он не очень ясен ;

Хотелось бы расписать все так ,как описано математический;

Добавлено через 1 минуту
Про остановку пока не сойдется это понятно,
там нужно какое то условие остановки

например, как вы и сказали , разность последующей и предыдущей итераций ,только норма от этого дела

@outoftime · 21.05.2014, 11:33

Сообщение от Qazan

но мне он не очень ясен

Это плохо, псевдокод на вики, обычно, лучше всего описывает реализацию (оптимизации только language specific можно проводить).

@Qazan · 21.05.2014, 15:20 **[ТС]**

outoftime, Обычно мне было понятно , но в данном случае мне не ясно

$https://www.cyberforum.ru/cgi-bin/latex.cgi? \sigma \leftarrow \sigma + a_{ij} \phi_j$

Вот ,что означает такая запись ?

@outoftime · 21.05.2014, 16:46

Qazan, эквивалент на плюсах

C++
1
sigma = sigma + a[i][j] * phi[j];

Где sigma, a, phi - переменные соответствующих типов и для тих определены данные операции.

@Qazan · 21.05.2014, 18:02 **[ТС]**

outoftime, а дак стрелка это присваивание .Ясно

@outoftime · 21.05.2014, 20:23

Не по теме:

Qazan, стандартная нотация, кажись, в Кормене такая же. Советую почитать.

@Qazan · 21.05.2014, 21:28 **[ТС]**

outoftime, Пасибо ,почитаю

Новые блоги и статьи Все статьи Все блоги /
Access VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .	Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .
Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .	Создание Single Page Application на фреймах krapotkin 16.11.2025 Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .	Фото: Daniel Greenwood kumehtar 13.11.2025

@outoftime ║XLR8║ 1212 / 909 / 270 Регистрация: 25.07.2009 Сообщений: 4,360 Записей в блоге: 5
	20.05.2014, 18:30
	Qazan, Что такое k, k+1? Что такое i, j? 0

@Qazan 215 / 63 / 25 Регистрация: 30.04.2013 Сообщений: 865 Записей в блоге: 10
	20.05.2014, 19:23 [ТС]
	outoftime, k, k + 1 это соответственно текущий и предыдущий вектор Это итерационный метод решения,в начале задаем приближенное значение ,а потом приближаем дальше i,j - индексы вектора формируется вектор приближенных значений постепенно Добавлено через 1 минуту $https://www.cyberforum.ru/cgi-bin/latex.cgi? a_{ij}$ Элементы матрицы A. 0

@Qazan 215 / 63 / 25 Регистрация: 30.04.2013 Сообщений: 865 Записей в блоге: 10
	20.05.2014, 21:30 [ТС]
	outoftime, переменный это эле менты вектора тобишь просто числа Добавлено через 32 секунды Successive_over-relaxation Добавлено через 2 минуты Successive_over-relaxation 0

@outoftime ║XLR8║ 1212 / 909 / 270 Регистрация: 25.07.2009 Сообщений: 4,360 Записей в блоге: 5
	21.05.2014, 00:35
	Qazan, на самом деле, на вики ниже есть алгоритм, и даже там вместо двух циклов используется один с условием (кода на 1 операцию больше, зато читабельней). Я только не понял как выбирается начальное значение и что значит "повторять пока не сойдется", сам алгоритм выдает результат схождением (по аналогии с разложением функции в ряд быстрого схождения). У меня только одна мысль по этому поводу - делать по аналогии с рядами измеряя эпсилон как разницу между значениями на смежных итерациях. Если после этого у вас возникают вопросы не касающиеся кода - можете начинать писать код. Если у вас есть явные ошибки в коде - форумчане помогут. Но копаться в этом методе решения линейных уравнений я не хочу (остальные, скорей всего, меня поддержат). 0

@Qazan 215 / 63 / 25 Регистрация: 30.04.2013 Сообщений: 865 Записей в блоге: 10
	21.05.2014, 05:37 [ТС]
	outoftime, про реализацию это да, там есть псевдокод но мне он не очень ясен ; Хотелось бы расписать все так ,как описано математический; Добавлено через 1 минуту Про остановку пока не сойдется это понятно, там нужно какое то условие остановки например, как вы и сказали , разность последующей и предыдущей итераций ,только норма от этого дела 0

@Qazan 215 / 63 / 25 Регистрация: 30.04.2013 Сообщений: 865 Записей в блоге: 10
	21.05.2014, 15:20 [ТС]
	outoftime, Обычно мне было понятно , но в данном случае мне не ясно $https://www.cyberforum.ru/cgi-bin/latex.cgi? \sigma \leftarrow \sigma + a_{ij} \phi_j$ Вот ,что означает такая запись ? 0

@Qazan 215 / 63 / 25 Регистрация: 30.04.2013 Сообщений: 865 Записей в блоге: 10
	21.05.2014, 18:02 [ТС]
	outoftime, а дак стрелка это присваивание .Ясно 0

@outoftime
	21.05.2014, 20:23
	Не по теме: Qazan, стандартная нотация, кажись, в Кормене такая же. Советую почитать. 0

@Qazan 215 / 63 / 25 Регистрация: 30.04.2013 Сообщений: 865 Записей в блоге: 10
	21.05.2014, 21:28 [ТС]
	outoftime, Пасибо ,почитаю 0