Сглаживание графика

@Robolightning · Регистрация: 02.07.2021

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

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

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
void smooth(std::vector<buf>& b, double g_smooth, int N)
{
    double dif1, dif2, acc = 1.0 + 10.0 / g_smooth;;
    bool end = true;
    int i;
    while (end)
    {
        end = false;
        dif1 = sum(b[N / 2].p, - b[N / 2 + 1].p);
        for (i = N / 2; i < N - 1; i++)
        {
            dif2 = sum(b[i].p, - b[i + 1].p);
            if (abs(dif1 / dif2) > acc)
            {
                dif1 = sum(b[i].p, - b[i - 1].p);
                b[i].p = sum(b[i].p, - dif1 / 3.0);
                b[i - 1].p = sum(b[i - 1].p, dif1 / 3.0);
                end = true;
            }
            else if (abs(dif2 / dif1) > acc)
            {
                dif1 = sum(b[i + 1].p, - b[i].p);
                b[i].p = sum(b[i].p, dif1 / 3.0);
                b[i + 1].p = sum(b[i + 1].p, -dif1 / 3.0);
                end = true;
            }
            dif1 = dif2;
        }
        if (!end) { break; }
        dif1 = sum(b[N / 2].p, - b[N / 2 - 1].p);
        for (i = N / 2; i > 0; i--)
        {
            dif2 = sum(b[i].p, - b[i - 1].p);
            if (abs(dif1 / dif2) > acc)
            {
                dif1 = sum(b[i].p, - b[i + 1].p);
                b[i].p = sum(b[i].p , -dif1 / 3.0);
                b[i + 1].p = sum(b[i + 1].p, dif1 / 3.0);
                end = true;
            }
            else if (abs(dif2 / dif1) > acc)
            {
                dif1 = sum(b[i - 1].p, - b[i].p);
                b[i].p = sum(b[i].p, dif1 / 3.0);
                b[i - 1].p = sum(b[i - 1].p, -dif1 / 3.0);
                end = true;
            }
            dif1 = dif2;
        }
    }
}

Алексей1153 · 28.07.2021, 11:36

первым делом исправить это

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

abs

->
std::abs

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

sum

- это откуда и что?

вторым делом - подать на вход заданные вручную тестовые данные из нескольких точек. После этого можно продолжить отладку

Добавлено через 52 секунды
неплохо бы компилируемый код в пост приложить. А то многое не определено в примере

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

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

dif2 = sum(b[i].p, - b[i + 1].p);
if (abs(dif1 / dif2) > acc)
{
dif1 = sum(b[i].p, - b[i - 1].p);
b[i].p = sum(b[i].p, - dif1 / 3.0);
b[i - 1].p = sum(b[i - 1].p, dif1 / 3.0);
end = true;
}
else if (abs(dif2 / dif1) > acc)
{
dif1 = sum(b[i + 1].p, - b[i].p);
b[i].p = sum(b[i].p, dif1 / 3.0);
b[i + 1].p = sum(b[i + 1].p, -dif1 / 3.0);
end = true;
}
dif1 = dif2;

dif1 не зависимо от ветвлений будет заполнен значением dif2

Добавлено через 1 минуту
в циклах повторяющийся код - его можно в лямбду оформить

@Robolightning · 28.07.2021, 14:13 **[ТС]**

sum это функция точного сложения, в новом варианте кода, приложенном ниже, я заменил её на обычное сложение, чтобы было понятнее. abs заменил, но это ничего не изменило. Весь код программы слишком большой и находится в нескольких файлах, не вижу смыла его кидать. На этот раз заполнил каждую двадцатую ячейку вектора значением 10, остальные -- нулями. Результат получился довольно интересный, однако это не меняет того факта, что с конечными значениями графика происходит нечто непонятное. Логично, что в теории должна получаться вертикальная прямая линия. С лямбдой не разобрался, использую два цикла по причине того, что один цикл сносил картину куда-то вбок, а покрывать каждым циклом весь диапазон почему-то не выходит: функция зацикливается. Потому и поставил до половины, так хоть работает. (Первый стриншот -- картина до сглаживания, второй -- после (коэффициент сглаживания g_smooth беру равным 10)).
P.S. То, что линии на первом скриншоте неровные, по всей видимости ошибки отображения, в векторе всё верно.

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
56
57
58
void smooth(std::vector<buf>& b, double g_smooth, int N)
{
    double dif1, dif2, acc = 1.0 + 10.0 / g_smooth;;
    bool end = true;
    int i;
 
    for (i = 0; i < N; i++)
    {
        if (i % 20 == 0) b[i].p = 10;
        else b[i].p = 0;
    }
 
    while (end)
    {
        end = false;
        dif1 = b[N / 2].p - b[N / 2 + 1].p;
        for (i = N / 2; i < N - 1; i++)
        {
            dif2 = b[i].p - b[i + 1].p;
            if (std::abs(dif1 / dif2) > acc)
            {
                dif1 = b[i].p - b[i - 1].p;
                b[i].p = b[i].p - dif1 / 3.0;
                b[i - 1].p = b[i - 1].p + dif1 / 3.0;
                end = true;
            }
            else if (std::abs(dif2 / dif1) > acc)
            {
                dif1 = b[i + 1].p - b[i].p;
                b[i].p = b[i].p + dif1 / 3.0;
                b[i + 1].p = b[i + 1].p - dif1 / 3.0;
                end = true;
            }
            dif1 = dif2;
        }
        if (!end) { break; }
        dif1 = b[N / 2].p - b[N / 2 - 1].p;
        for (i = N / 2; i > 0; i--)
        {
            dif2 = b[i].p - b[i - 1].p;
            if (std::abs(dif1 / dif2) > acc)
            {
                dif1 = b[i].p - b[i + 1].p;
                b[i].p = b[i].p - dif1 / 3.0;
                b[i + 1].p = b[i + 1].p + dif1 / 3.0;
                end = true;
            }
            else if (std::abs(dif2 / dif1) > acc)
            {
                dif1 = b[i - 1].p - b[i].p;
                b[i].p = b[i].p + dif1 / 3.0;
                b[i - 1].p = b[i - 1].p - dif1 / 3.0;
                end = true;
            }
            dif1 = dif2;
        }
    }
}

Новые блоги и статьи Все статьи Все блоги /
сукцессия микоризы: основная теория в виде двух уравнений. anaschu 11.01.2026 https:/ / rutube. ru/ video/ 7a537f578d808e67a3c6fd818a44a5c4/	WordPad для Windows 11 Jel 10.01.2026 WordPad для Windows 11 — это приложение, которое восстанавливает классический текстовый редактор WordPad в операционной системе Windows 11. После того как Microsoft исключила WordPad из. . .	Classic Notepad for Windows 11 Jel 10.01.2026 Old Classic Notepad for Windows 11 Приложение для Windows 11, позволяющее пользователям вернуть классическую версию текстового редактора «Блокнот» из Windows 10. Программа предоставляет более. . .	Почему дизайн решает? Neotwalker 09.01.2026 В современном мире, где конкуренция за внимание потребителя достигла пика, дизайн становится мощным инструментом для успеха бренда. Это не просто красивый внешний вид продукта или сайта — это. . .
Модель микоризы: классовый агентный подход 3 anaschu 06.01.2026 aa0a7f55b50dd51c5ec569d2d10c54f6/ O1rJuneU_ls https:/ / vkvideo. ru/ video-115721503_456239114	Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ФедосеевПавел 06.01.2026 Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .	Модель микоризы: классовый агентный подход 2 anaschu 06.01.2026 репозиторий https:/ / github. com/ shumilovas/ fungi ветка по-частям. коммит Create переделка под биомассу. txt вход sc, но sm считается внутри мицелия. кстати, обьем тоже должен там считаться. . . .	Расчёт токов в цепи постоянного тока igorrr37 05.01.2026 / * Дана цепь постоянного тока с сопротивлениями и источниками (напряжения, ЭДС и тока). Найти токи и напряжения во всех элементах. Программа составляет систему уравнений по 1 и 2 законам Кирхгофа и. . .