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

Вычисление площади криволинейной трапеции - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 23, средняя оценка - 4.65
Peroksid123
0 / 0 / 0
Регистрация: 15.03.2010
Сообщений: 7
15.03.2010, 20:41     Вычисление площади криволинейной трапеции #1
Ребят,помогите с курсачём.Договорился с одним из участников форума,он сделал курсовую и пропал,а препод требует описание алгоритма.Помогите пожалуйста кому не лень.

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
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
/*вычислить с точностью Е площадь криволинейной трапеции ограниченной осью , Ох, прямыми х=а х=B , и кривой 
y=f(x)=(1-0,25sin^2tx)^1\2 . параметр "t"- корень уравнения t^3-0,39t^2-10,5t+11=0,
принад. отрезку [c, d] и отрезку с точ. Е
 
исх данные:
a=0 b=П\2
с=1
d=2
E=0,001
*/
#include <iostream>
#include <cmath>
#define _PI 3.1415926535897932384626433832795
using std::cin;
using std::cout;
using std::endl;
 
 
//////////////////////////////////////////////////////////////////////////////////////
//      уравнение трапеции(интегрируем)                                             //
//      вход:                                                                       //
//          X-координата  по оси абцисс                                             //
//          T-параметр является корнем уравнения f0                                 //
//      выход:                                                                      //
//          значение уравнения                                                      //
//////////////////////////////////////////////////////////////////////////////////////
float f1(float x, float t);
 
//////////////////////////////////////////////////////////////////////////////////////
//      уравнение (ищем корень)                                                     //
//      вход:                                                                       //
//          T-переменная                                                            //
//      выход:                                                                      //
//          значение уравнения в точке t                                            //
//////////////////////////////////////////////////////////////////////////////////////
float f0(float t);
 
//////////////////////////////////////////////////////////////////////////////////////
//      ПРЯМОУГОЛЬНИКИ (правые) - интегрирование                                    //
//      вход:                                                                       //
//          a-начало интервала интегрирования                                       //
//          б-конец интервала интегрирования                                        //
//          h-точность интегрирования(совпадает с шагом интегрирования              //
//          t-параметр из решения уравнения                                         //
//      выход:                                                                      //
//          интеграл функции f1  сточностью h на интервале [a,b]                    //
//////////////////////////////////////////////////////////////////////////////////////
float rectangle(float a,float b,float h,float t);
 
//////////////////////////////////////////////////////////////////////////////////////
//      ПОЛОВИННОЕ ДЕЛЕНИЕ - поиск корня                                            //
//      вход:                                                                       //
//          a-начало интервала поиска                                               //
//          б-конец интервала поиска                                                //
//          h-точность поиска                                                       //
//      выход:                                                                      //
//          интеграл функции f0  сточностью h на интервале [a,b]                    //
//          если корня нет то значение не попадающее в интервал поиска              //
//////////////////////////////////////////////////////////////////////////////////////
float separ(float a,float b,float e);
 
//////////////////////////////////////////////////////////////////////////////////////
//      главная функция                                                             //
//      вход:                                                                       //
//      выход:                                                                      //
//          0-нет ошибок                                                            //
//          1-ошибка нет корня уравнения fo                                         //
//////////////////////////////////////////////////////////////////////////////////////
int main(){
    float a,b,c,d,e,tmp;
    //a=0,b=_PI/2.,c=1,d=2,e=0.001;
    //cout<<rectangle(a,b,e,separ(c,d,e));
    cout<<"a: "; cin>>a;
    cout<<"b: "; cin>>b;
    cout<<"c: "; cin>>c;
    cout<<"d: "; cin>>d;
    cout<<"e: "; cin>>e;
    if((tmp=separ(c,d,e))<c)
        return 1;
    cout<<rectangle(a,b,e,tmp)<<endl;
    system("pause");
    return 0;
}
 
//////////////////////////////////////////////////////////////////////////////////////
//      уравнение трапеции(интегрируем)                                             //
//      вход:                                                                       //
//          X-координата  по оси абцисс                                             //
//          T-параметр является корнем уравнения f0                                 //
//      выход:                                                                      //
//          значение уравнения                                                      //
//////////////////////////////////////////////////////////////////////////////////////
float f1(float x, float t){ 
    return sqrt((1-0.25*sin(t*x)*sin(t*x)));
}
 
//////////////////////////////////////////////////////////////////////////////////////
//      уравнение (ищем корень)                                                     //
//      вход:                                                                       //
//          T-переменная                                                            //
//      выход:                                                                      //
//          значение уравнения в точке t                                            //
//////////////////////////////////////////////////////////////////////////////////////
float f0(float t){          
    return t*t*t-0.39*t*t-10.5*t+11;
}
 
//////////////////////////////////////////////////////////////////////////////////////
//      ПРЯМОУГОЛЬНИКИ (правые) - интегрирование                                    //
//      вход:                                                                       //
//          a-начало интервала интегрирования                                       //
//          б-конец интервала интегрирования                                        //
//          h-точность интегрирования(совпадает с шагом интегрирования              //
//          t-параметр из решения уравнения                                         //
//      выход:                                                                      //
//          интеграл функции f1  сточностью h на интервале [a,b]                    //
//////////////////////////////////////////////////////////////////////////////////////
float rectangle(float a,float b,float h,float t){
    float s,x1,x2;
    for(s=0,x1=a,x2=b;x1<x2;x1+=h)  //площадь=0,x1=начало интервала,х2=конец,пока х1 < х2
        s+=f1(x1,t);                //увелививаем х1 на величину шага и к площади прибавляем значение
                                    //функции в точке x1
    return s*h;                     //чтобыфиз суммызначений функции получить площадь умножим на величину шага
}
//////////////////////////////////////////////////////////////////////////////////////
//      ПОЛОВИННОЕ ДЕЛЕНИЕ - поиск корня                                            //
//      вход:                                                                       //
//          a-начало интервала поиска                                               //
//          б-конец интервала поиска                                                //
//          h-точность поиска                                                       //
//      выход:                                                                      //
//          интеграл функции f0  сточностью h на интервале [a,b]                    //
//          если корня нет то значение не попадающее в интервал поиска              //
//////////////////////////////////////////////////////////////////////////////////////
float separ(float a,float b,float e){
    float x1=a;
    float x2=b;
    float x;
    if(f0(a)*f0(b)>0)       //функция не пересекает ось X
        return a-1;         //если корня нет то возвращаем значение не попадающее в интервал поиска
    do{
        x=(x1+x2)/2;        //находим середину интервала поиска
        if(f0(x1)*f0(x)>0)  //если значения начала интервала и середины имеют разные знаки
            x1=x;           //теперь середина интервала будет являться началом интервала поиска
        else                //иначе корень на второй половине
            x2=x;           //теперь середина интервала будет являться концом интервала поиска
    }while(fabs(f0(x))>e);
    return x;
}
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
vet
 Аватар для vet
175 / 175 / 13
Регистрация: 08.04.2009
Сообщений: 1,309
16.03.2010, 05:31     Вычисление площади криволинейной трапеции #2
Цитата Сообщение от Peroksid123 Посмотреть сообщение
описание алгоритма.
Описание? Вы про блок - схему может говорите? Потому что комментов в проге даже больше чем нужноО_о
Peroksid123
0 / 0 / 0
Регистрация: 15.03.2010
Сообщений: 7
16.03.2010, 08:00  [ТС]     Вычисление площади криволинейной трапеции #3
Описание? Вы про блок - схему может говорите?
Возможно и про блок-схему,но препод называл это описание алгоритма
Zanoza
 Аватар для Zanoza
9 / 9 / 1
Регистрация: 16.03.2010
Сообщений: 44
17.03.2010, 12:49     Вычисление площади криволинейной трапеции #4
ну так это и есть блок схема,то есть описание работы программы
Peroksid123
0 / 0 / 0
Регистрация: 15.03.2010
Сообщений: 7
17.03.2010, 14:37  [ТС]     Вычисление площади криволинейной трапеции #5
ну так это и есть блок схема,то есть описание работы программы
нет,мне нужно чтобы было расписано по пунктам как эта программа вычисляет площадь т.е. например данной командой считается интеграл и.т.д.
Yandex
Объявления
17.03.2010, 14:37     Вычисление площади криволинейной трапеции
Ответ Создать тему
Опции темы

Текущее время: 19:12. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru