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

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

Восстановить пароль Регистрация
Другие темы раздела
C++ Разработать подпрограмму обработки данного файла http://www.cyberforum.ru/cpp-beginners/thread83734.html
2) Создать файл начисления студентам стипендии. Студент, получивший все оценки “5” и активно участвующий в общественной работе, зачисляется на повышенную стипендию (доплата 50%), не активно участвующие – доплата 25%, студенты, получившие “4” и “5” зачисляются на обычную стипендию. Студент, получивший одну оценку “3”, но активно занимающийся общественной работой, также зачисляется на стипендию. В...
C++ копирование файлов Зраствуйте можете подсказать алгоритм копирования файлов любого типу(разширения)?? http://www.cyberforum.ru/cpp-beginners/thread83733.html
Число n C++
Дано натуральное число n. Получить последовательность b1,...,bn, где при i=1,2,...,n значение bi равно: p.s зааранее all thanks
C++ Что такое дескиптор?
По понятнее пожалуйста на википедиа смотрел ни черта не понял.
C++ рисование в с++ http://www.cyberforum.ru/cpp-beginners/thread83703.html
можно ли рисовать в borland c++ (если можно то как)reen"]
C++ НАСЛЕДОВАНИЕ И ВИРТУАЛЬНЫЕ ФУНКЦИИ Определить иерархию классов (игрушка, продукт, товар, молочный продукт). Определить в классе статическую компоненту - указатель на начало связанного списка объектов и статическую функцию для просмотра списка. Реализовать классы. Написать демонстрационную программу, в которой создаются объекты различных классов и помещаются в список, после чего список просматривается. Сделать соответствующие... подробнее

Показать сообщение отдельно
AraxniD
0 / 0 / 0
Регистрация: 08.01.2010
Сообщений: 7
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;
}
 
Текущее время: 17:24. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru