Қазақстан Республикасы ғылымжәне жоғары білім министрлігі



Дата26.10.2023
өлшемі377,46 Kb.
#188280
Байланысты:
Омарбек Айдана СРС
951 tarih-12-3, «Ôèçèêà» ï í³íåí 4-òî ñàí à àðíàë àí æèûíòû áà àëàóäû òàïñûðìà, Нургазы Айда фил срс., Sabey fizo 10

Қазақстан Республикасы ғылымжәне жоғары білім министрлігі
Әл - Фараби атындағы Қазақ Ұлттық Университеті



СӨЖ
Факультеті: « Механика-Математика »
Кафедра: « Есептеу ғылымдары және статистика »
Тақырыбы: Ритц және Галеркин әдісі

Орындаған: ЕҒжС-21-3-А


мамандығының 3-курс студенті
Омарбек Айдана Ғаниқызы

Тексерген: ЕҒжС кафедрсының аға оқытушысы


Касенов Сырым Еркинович
Алматы, 2023ж





Код:
#include


#include
#include
// кітапханалар
using namespace std;
// чтоб стд-ны каждый раз жазбас үшін
double phi1(double x) {
return sin(x);
}

double phi2(double x) {


return cos(x);
}
double f(double x) {
return x * x ;
} int main(){
double a = 0.5;
double b = 0.02;
double c = 0.001;
double interval_start = 0.0;
double interval_end = 1.0;
int num_basis_functions = 2;
int num_nodes = 100;
//Галеркин әдісімен шешім:
vector coefficients( num_nodes, 0.0);
//цикл
for (int i = 0; i < num_basis_functions ; ++i) {
for (int j = 0; j < num_basis_functions; ++j) {

double integral = 0.0;

for (int k = 0; k < num_nodes; ++k) {
double x_k = interval_start + k * (interval_end - interval_start) / (num_nodes - 1);
integral += (a * phi1(x_k) * phi1(x_k) + b * phi1(x_k) * phi2(x_k) + c * phi2(x_k) * phi2(x_k));
}

coefficients[i] += integral;


}
}
//негізгі функциялардың интегралдарын есептеу арқылы Галеркин әдісі үшін есептеулер жүргізеді
double solution = 0.0;

for (int i = 0; i < num_basis_functions; ++i) {


solution += coefficients[i] * phi1(interval_start + i * (interval_end - interval_start) / (num_nodes - 1));
}
//шығару
cout << "Solution: " << solution << endl;
return 0;
}

Есептеу нәтижесі экранға шығарылады.





Достарыңызбен бөлісу:




©engime.org 2024
әкімшілігінің қараңыз

    Басты бет