Апроксимация функции с фурье, компилятор ошибок не находит, но программа работает не верно.

@AraxniD · Регистрация: 08.01.2010

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

Надо апроксимировать функцию..грубо говоря, я раскладываю функцию в ряд фурье и считаю значение функции в каждой точке х. После этого, по этим значениям строю график, но это уже касается другой программы.
Проблема в том, что вместо циферок на экран выводится белеберда. Где проблема не вижу.
Писал на вижуал с++ 2008.
Буду очень благодарен за помощь.

2*l(L)=1.
T=6
самое уравнение "раскладывание"
f(x)=Ao/2+(ряд от к=1 до м) (Ak*cos(k*w*x)+Bk*sin(k*w*x))
Ao=1/l (интеграл от -l до l) f(x)dx
Ak=1/l (интеграл от -l до l) f(x)cos(k*w*x)dx
Bk аналогчно Ak, вместо косинуса синус.
W=pi/l

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
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
//
 
#include "stdafx.h"
#include <iostream>
#include <math.h>
using namespace std;
double a_nol(double,int,int,double );
double a_k(double,int,int,int,double );
double b_k(double,int,int,int,double );
double fun(double);
double integral(int,int,double (*f)(double),int,int,double);
void main(void)
{
    int o;
    double Ao,Ak,Bk;
    double l=0.5;
    int n=10;
    double eps=0.01;
    double epsx=0.05;
    int k=1;
    int m=10;
    double x=-1; // Дальше я считаю значение функции в каждой точки х. Пока это я выкинул
    double S1,S2=0;
    while (k<=m) // Считается ряд. От к=1 до м, где м грубо говоря точность подсчета
        {
            Ao=a_nol(x,l,n,eps);
            Ak=a_k(x,k,l,n,eps);
            Bk=b_k(x,k,l,n,eps);
            S1=Ak*cos(k*3.14*x/0.5)+Bk*sin(k*3.14*x/0.5);
            S2+=S1;
            k++;
        }
    S2=Ao+S2;
    cout<<S2<<endl;
    cout<<Ao<<Ak<<Bk<<endl;
    cin>>o;
}
 
double a_nol(double x,int l,int n,double eps)
{
    int k=0;
    double I1,I2;
    int z=0;
    I2=integral(l,n,fun,z,k,x);
    do
    {
        I1=I2;
        n=n*2;
        I2=integral(l,n,fun,z,k,x);
    }
    while (fabs(I1-I2)>eps);
    return 2*I2;
}
double a_k(double x,int k,int l,int n,double eps)
{
    double I1,I2;
    int z=1;
    I2=integral(l,n,fun,z,k,x);
    do
    {
        I1=I2;
        n=n*2;
        I2=integral(l,n,fun,z,k,x);
    }
    while (fabs(I1-I2)>eps);
    return 2*I2;
}
double b_k(double x,int k,int l,int n,double eps)
{
    double I1,I2;
    int z=2;
    I2=integral(l,n,fun,z,k,x);
    do
    {
        I1=I2;
        n=n*2;
        I2=integral(l,n,fun,z,k,x);
    }
    while (fabs(I1-I2)>eps);
    return 2*I2;
}
double fun(double x)
{
    if (x<=0)
        return 1.5;
    else
        return -1.5;
 
}
double integral(int l,int n,double (*f)(double),int z,int k,double x) //Метод прямоугольников
{
    double h1,S=0;                                                       //Возможн ошибка здесь, но не вижу.
    h1=l/n;
    for (int i=1;i<=n;i++)
    {
        switch(z)
        {
        case 0:
            S+=f(0+i*h1);
            break;
        case 1:
            S+=f(0+i*h1)*cos(k*3.14*x/l);
            break;
        case 2:
            S+=f(0+i*h1)*sin(k*3.14*x/1);
            break;
        }
    }
    return h1*S;
}

@AraxniD · 11.01.2010, 22:19 **[ТС]**

Значит все исправил, прога работает, но вместо нормального результата, координаты сильно смещены вверх\вниз. В чем ошибка найти не могу.
Хелп!

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
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
#include "stdafx.h"
#include <iostream>
#include <math.h>
using namespace std;
double anol(double,int,double (*f)(double,double,int,int),int);
double bkoef(double,int,double (*f)(double,double,int,int),int,int);
double akoef(double,int,double (*f)(double,double,int,int),int,int);
double fun(double,double,int,int);
double integral(double,int,double (*f)(double,double,int,int),int,int,int);
void main(void)
{
    double Ak,Ao,Bk;
    double p;
    double l=0.5;
    int n=5;
    int k=1,tochnost=10;
    double S=0;
    double shag=0.1;
    for (double x=-0.5;x<=-0.4;x+=shag)
    {
        Ao=anol(l,n,fun,x);
        while (k<=tochnost)
        {
          Ak=akoef(l,n,fun,x,k);
          Ak=Ak*cos(k*3.14*x/0.5);
          Bk=bkoef(l,n,fun,x,k);
          Bk=Bk*sin(k*3.14*x/0.5);
          S=Ak+Bk;
          cout<<S<<" ";
          k++;
        }
        S=Ao+S;
        k=1;
        cout<<S<<" "<<x<<endl;
        S=0;
    }
    cin>>p;
}
double akoef(double l,int n,double (*f)(double,double,int,int),int x,int k)
{
    double eps=0.001;
    double I1,I2;
    int c=1;
    I2=integral(l,n,fun,x,k,c);
    do
    {
        I1=I2;
        n=n*2;
        I2=integral(l,n,fun,x,k,c);
    }
    while (fabs(I1-I2)>eps);
    return 2*I2/0.5;
}
double fun(double z,double x,int k,int c)
{
    switch(c)
    {
        case 0:
            if (z>=0)
                return 1.5;
            else
                return -1.5;
            break;
        case 1:
            if (z>=0)
                return 1.5*cos(k*3.14*x/0.5);
            else
                return -1.5*cos(k*3.14*x/0.5);
            break;
        case 2:
            if (z>=0)
                return 1.5*sin(k*3.14*x/0.5);
            else
                return -1.5*sin(k*3.14*x/0.5);
            break;
    }
}
double integral(double l,int n,double (*f)(double,double,int,int),int x,int k,int c)
{
    double h,s=0;
    h=l/n;
    for (int i=1;i<=n;i++)
        s+=f(0+i*h,x,k,c);
    return h*s;
}
double bkoef(double l,int n,double (*f)(double,double,int,int),int x,int k)
{
    double eps=0.001;
    double I1,I2;
    int c=2;
    I2=integral(l,n,fun,x,k,c);
    do
    {
        I1=I2;
        n=n*2;
        I2=integral(l,n,fun,x,k,c);
    }
    while (fabs(I1-I2)>eps);
    return 2*I2/0.5;
}
double anol(double l,int n,double (*f)(double,double,int,int),int x)
{
    double eps=0.001;
    double I1,I2;
    int k=0;
    int c=0;
    I2=integral(l,n,fun,x,k,c);
    do
    {
        I1=I2;
        n=n*2;
        I2=integral(l,n,fun,x,k,c);
    }
    while (fabs(I1-I2)>eps);
    return 2*I2/0.5;
}

Новые блоги и статьи Все статьи Все блоги /
Новый ноутбук 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	Расскажи мне о Мире, бродяга kumehtar 12.11.2025 — Расскажи мне о Мире, бродяга, Ты же видел моря и метели. Как сменялись короны и стяги, Как эпохи стрелою летели. - Этот мир — это крылья и горы, Снег и пламя, любовь и тревоги, И бескрайние. . .