This repository was archived by the owner on Jan 3, 2023. It is now read-only.
-
Notifications
You must be signed in to change notification settings - Fork 40
/
Copy pathSource.cpp
67 lines (46 loc) · 1.41 KB
/
Source.cpp
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
#include<iostream>
#include<Windows.h>
#include<fstream>
#include<cmath>
#include<vector>
using namespace std;
/*Íà÷àëüíàÿ òåìïåðàòóðà*/
const double INITIAL_TEMPERATURE = 18;
int main()
{
SetConsoleCP(1251);
SetConsoleOutputCP(1251);
vector <double> linearTemperature(1, INITIAL_TEMPERATURE);
vector <double> nonLinearTemperature(1, INITIAL_TEMPERATURE);
int iterNum;
cout << "Êîëè÷åñòâî èòåðàöèé: "; cin >> iterNum;
/*Êîýôôèöèåíòû äëÿ ëèíåéíîãî óðàâíåíèÿ*/
double a = 1, b_lin = 0.05;
/*Êîëè÷åñòâî òåïëîòû*/
double warm;
cout << "Êîëè÷åñòâî òåïëîòû: "; cin >> warm;
/*Ôàéëû, äëÿ çàïèñè çíà÷åíèé ôóíêöèé*/
ofstream linOutput("linear.txt");
ofstream nonLinOutput("nonlinear.txt");
for (int i = 0; i < iterNum; i++) {
linearTemperature.push_back(a * linearTemperature[i] + b_lin * warm);
}
cout << "Ëèíåéíàÿ ôóíêöèÿ:\n";
for (int i = 0; i < iterNum; i++)
{
cout << ". " << " " << linearTemperature[i] << endl;
linOutput << i << "\t" << linearTemperature[i] << endl;
}
/*Êîýôôèöèåíòû äëÿ íåëèíåéíîãî óðàâíåíèÿ*/
double b_nonLin = 0.0001;
for (int i = 0; i < iterNum; i++)
{
nonLinearTemperature.push_back(linearTemperature[i] - b_nonLin * pow(1.005, 2) + 1.69 * warm + 1.25 * sin(0.95));
}
cout << "Íåëèíåéíàÿ ôóíêöèÿ:\n";
for (int i = 0; i < iterNum; i++)
{
cout << i << "s\t" << nonLinearTemperature[i] << endl;
nonLinOutput << i << "\t" << nonLinearTemperature[i] << endl;
}
}