added plot

This commit is contained in:
ada-dmitry
2023-12-07 12:51:59 +03:00
parent 084d0ae413
commit 9c3ba250d4
8 changed files with 96 additions and 126 deletions
+1 -56
View File
@@ -94,7 +94,7 @@ void choose_path() // Функция выбора действия после в
printf("\nПродолжить поиск корня?\n"); printf("\nПродолжить поиск корня?\n");
printf("\nC - продолжить, A - закончить работу программы, \ printf("\nC - продолжить, A - закончить работу программы, \
R - начать поиск на другом отрезке, M - выбрать другой способ решения: \n"); R - начать поиск на другом отрезке: \n");
do do
{ {
@@ -146,61 +146,6 @@ R - начать поиск на другом отрезке, M - выбрать
} }
} }
return; return;
case 'M':
case 'm':
__fpurge(stdin);
// Включение символа ECHO
term.c_lflag |= ECHO;
// Применение новых настроек терминала
tcsetattr(0, TCSANOW, &term);
printf("Выберите метод вычисления корня:\
0 - метод вилки, 1 - метод хорд, 2 - метод касательных\n(По умолчанию: 0)\n");
*mth = getchar();
printf("Введите новые значения границ a и b: ");
scanf("%lf%lf", pa, pb);
if (f(*pa) * f(*pb) > 0)
{
printf("В указанных границах нуля функции не существует.\n");
printf("Нажмите Ctrl+C и поменяйте границы\n");
}
if (*pa == *pb)
{
if (f(*pa) == 0)
{
printf("\nКорень уравнения: %lf\n", *pa);
printf("Работа программы завершена.\n");
exit(0);
}
else
{
printf("Нажмите Ctrl+C и поменяйте границы\n");
}
}
if (*mth == 0)
{
printf("Fork\n");
fork_method();
}
else if (*mth == 1)
{
printf("Chord\n");
chord_method();
}
else if (*mth == 2)
{
printf("Newton\n");
newton_method();
}
else
{
printf("Fork\n");
fork_method();
}
return;
default: default:
printf("Ошибка ввода, повторная попытка...\n"); printf("Ошибка ввода, повторная попытка...\n");
Binary file not shown.
+3 -2
View File
@@ -6,19 +6,20 @@ CFLAGS= -c -Wall
default: build default: build
general: main.o functions.o general: main.o functions.o
killall -9 function_zeros_calculator
$(CC) ./out/main.o ./out/functions.o declare.h -o ./out/function_zeros_calculator -lm
$(CC) ./out/main.o ./out/functions.o -o ./out/function_zeros_calculator -lm
main.o: main.c main.o: main.c
$(CC) $(CFLAGS) main.c -o ./out/main.o $(CC) $(CFLAGS) main.c -o ./out/main.o
functions.o: functions.c functions.o: functions.c
$(CC) $(CFLAGS) functions.c -o ./out/functions.o $(CC) $(CFLAGS) functions.c -o ./out/functions.o
clean: clean:
rm -rf ./out/main.o ./out/functions.o rm -rf ./out/main.o ./out/functions.o
build: general clean build: general clean
run: run:
View File
@@ -13,3 +13,4 @@ void chord_method();
void newton_method(); void newton_method();
double numerical_derivative(double x); double numerical_derivative(double x);
void exit_777(); void exit_777();
void show_grapgh();
+84 -62
View File
@@ -9,7 +9,36 @@
double f(double x) // Функция, заданная пользователем, для которой происходит поиск корней. double f(double x) // Функция, заданная пользователем, для которой происходит поиск корней.
{ {
return x * x * x - 125; return x - 125;
}
void show_grapgh()
{
struct termios term;
tcgetattr(STDIN_FILENO, &term);
term.c_lflag &= ~(ECHO | ICANON);
tcsetattr(STDIN_FILENO, TCSANOW, &term);
char func[100] = "x - 125";
// Подключение GNU Plot
FILE *gp = popen("gnuplot -persist", "w");
if (gp == NULL)
{
printf("GNUerror\n");
exit(3);
}
fprintf(gp, "set xrange [-500:500]\n");
fprintf(gp, "set yrange [0:500]\n");
fprintf(gp, "plot %s\n", func);
fflush(gp); // очищение буфера вывода
printf("Введите любую букву, чтобы продолжить...\n");
getchar();
fprintf(gp, "exit\n");
term.c_lflag |= (ECHO | ICANON);
tcsetattr(STDIN_FILENO, TCSANOW, &term);
return;
} }
// Приближенное вычисление производной по методу конечных разностей // Приближенное вычисление производной по методу конечных разностей
@@ -20,38 +49,34 @@ double numerical_derivative(double x)
void newton_method() void newton_method()
{ {
sig_flag = 0; *pb = *pa;
double left = *pa, right = *pb;
right = left;
do do
{ {
printf("%lf %lf\n", left, right); printf("%lf %lf\n", *pa, *pb);
left = right; *pa = *pb;
right = left - f(left) / numerical_derivative(left); *pb = *pa - f(*pa) / numerical_derivative(*pa);
// Используем численно вычисленную производную // Используем численно вычисленную производную
*px = right; *px = *pb;
sleep(1); sleep(1);
} while (fabs(right - left) > eps && f(x) != 0); } while (fabs(*pb - *pa) > eps && f(x) != 0);
return; return;
} }
void chord_method() void chord_method()
{ {
sig_flag = 0; *px = *pa;
double left = *pa, right = *pb; while ((*pb - *pa) > eps && f(x) != 0)
*px = left;
while ((right - left) > eps && f(x) != 0)
{ {
// Находим точку пересечения хорды с осью x // Находим точку пересечения хорды с осью x
*px = (left * f(right) - right * f(left)) / (f(right) - f(left)); *px = (*pa * f(*pb) - *pb * f(*pa)) / (f(*pb) - f(*pa));
printf("%lf %lf\n", left, right); printf("%lf %lf\n", *pa, *pb);
// Проверяем знаки функции в точках a и x // Проверяем знаки функции в точках a и x
if (f(*px) == 0) if (f(*px) == 0)
break; break;
else if (f(*px) * f(left) < 0) else if (f(*px) * f(*pa) < 0)
right = *px; *pb = *px;
else else
left = *px; *pa = *px;
sleep(1); sleep(1);
} }
return; return;
@@ -59,25 +84,24 @@ void chord_method()
void fork_method() // Функция поиска корней, использующая метод вилки. void fork_method() // Функция поиска корней, использующая метод вилки.
{ {
sig_flag = 0;
double left = *pa, right = *pb;
do do
{ {
*px = (left + right) / 2; // метод деления отрезка пополам *px = (*pa + *pb) / 2; // метод деления отрезка пополам
// Печать границ для отладки и для отслеживания хода вычислений // Печать границ для отладки и для отслеживания хода вычислений
printf("%lf %lf\n", left, right); printf("%lf %lf\n", *pa, *pb);
// Проверка, какой из отрезков [a, x] или [x, b] подходит под условие существования корня, // Проверка, какой из отрезков [a, x] или [x, b] подходит под условие существования корня,
// т.е. различие знаков на границе. // т.е. различие знаков на границе.
if (f(left) * f(x) <= 0)
right = x; if (f(*pa) * f(x) <= 0)
else if (f(right) * f(x) <= 0) *pb = x;
left = x; else if (f(*pb) * f(x) <= 0)
*pa = x;
sleep(1); // Принудительное ожидание программой, для возможности прерывания пользователем. sleep(1); // Принудительное ожидание программой, для возможности прерывания пользователем.
} while (fabs(right - left) > eps && f(x) != 0); } while (fabs(*pb - *pa) > eps && f(x) != 0);
return; return;
} }
@@ -85,16 +109,19 @@ void choose_path() // Функция выбора действия после в
{ {
sig_flag = 1; sig_flag = 1;
char choice, *ch = &choice; char choice;
struct termios original_term, modified_term; // Задание структуры, определяющей терминал Linux. struct termios default_term, modif_term; // Задание структуры, определяющей терминал Linux.
tcgetattr(STDIN_FILENO, &original_term); // Получение текущих настроек терминала
modified_term = original_term; tcgetattr(0, &default_term);
modified_term.c_lflag &= ~(ICANON | ECHO); // отключить канонического режим и отображение ввода modif_term = default_term;
tcsetattr(STDIN_FILENO, TCSANOW, &modified_term); // Отключение символа ECHO
modif_term.c_lflag &= ~ECHO;
// Применение новых настроек терминала
tcsetattr(0, TCSANOW, &modif_term);
printf("\nПродолжить поиск корня?\n"); printf("\nПродолжить поиск корня? Напечатайте букву и нажмите Enter\n");
printf("\nC - продолжить, A - закончить работу программы, \ printf("\nC - продолжить, A - закончить работу программы, \
R - начать поиск на другом отрезке: \n"); R - начать поиск на другом отрезке: \n");
@@ -102,45 +129,26 @@ R - начать поиск на другом отрезке: \n");
{ {
// Очистка стандартного потока ввода от "лишних" // Очистка стандартного потока ввода от "лишних"
__fpurge(stdin); __fpurge(stdin);
*ch = getchar(); // Запись буквы, введенной пользователем в переменную choice. scanf("%c", &choice); // Запись буквы, введенной пользователем в переменную choice.
switch (*ch) switch (choice)
{ {
case 'C': // Выбор вида "Продолжить вычисление". case 'C': // Выбор вида "Продолжить вычисление".
case 'c': case 'c':
tcsetattr(STDIN_FILENO, TCSANOW, &original_term); sig_flag = 0;
if (*mth == 0)
{
printf("Fork\n");
fork_method();
}
else if (*mth == 1)
{
printf("Chord\n");
chord_method();
}
else if (*mth == 2)
{
printf("Newton\n");
newton_method();
}
else
{
printf("Fork\n");
fork_method();
}
return; // Выход из функции без каких-либо изменений. return; // Выход из функции без каких-либо изменений.
case 'A': // Выбор вида "Завершить вычисление". case 'A': // Выбор вида "Завершить вычисление".
case 'a': case 'a':
printf("\nКорень уравнения: %lf\n", *px); printf("\nКорень уравнения: %lf\n", *px);
printf("Работа программы завершена.\n"); printf("Работа программы завершена.\n");
exit(0); // Принудительное завершение программы с выводом текущего приближения корня. exit_777(); // Принудительное завершение программы с выводом текущего приближения корня.
case 'R': // Выбор вида "Продолжить вычисления с новым отрезком". case 'R': // Выбор вида "Продолжить вычисления с новым отрезком".
case 'r': case 'r':
tcsetattr(STDIN_FILENO, TCSANOW, &original_term); // Возврат к обычному выводу
tcsetattr(0, TCSANOW, &default_term);
printf("Введите новые значения границ a и b: "); printf("Введите новые значения границ a и b: ");
@@ -157,22 +165,26 @@ R - начать поиск на другом отрезке: \n");
{ {
printf("\nКорень уравнения: %lf\n", *pa); printf("\nКорень уравнения: %lf\n", *pa);
printf("Работа программы завершена.\n"); printf("Работа программы завершена.\n");
exit(0); exit_777();
} }
else else
{ {
printf("Нажмите Ctrl+C и поменяйте границы\n"); printf("Нажмите Ctrl+C и поменяйте границы\n");
} }
} }
sig_flag = 0;
__fpurge(stdin);
return; return;
default: default:
printf("Ошибка ввода, повторная попытка...\n"); printf("Ошибка ввода, повторная попытка...\n");
__fpurge(stdin); __fpurge(stdin);
} }
} while (*ch != 'A' || *ch != 'R' || *ch != 'C' || *ch != 'a' || *ch != 'r' || *ch != 'c'); } while (choice != 'A' || choice != 'R' || choice != 'C' || choice != 'a' || choice != 'r' || choice != 'c' || choice != 'M' || choice != 'm');
tcsetattr(STDIN_FILENO, TCSANOW, &original_term); // Применение новых настроек терминала
tcsetattr(0, TCSANOW, &default_term);
sig_flag = 0; sig_flag = 0;
} }
@@ -184,4 +196,14 @@ void ctrlc_handler(int signum) // Функция, заменяющая стан
printf("\nТекущее приближение: %lf\n", *px); printf("\nТекущее приближение: %lf\n", *px);
choose_path(); choose_path();
} }
return;
}
void exit_777()
{
struct termios term;
tcgetattr(0, &term);
term.c_lflag |= ~ECHO;
tcsetattr(0, TCSANOW, &term);
exit(1);
} }
+4 -3
View File
@@ -11,7 +11,6 @@ short method, *mth, sig_flag;
int main() int main()
{ {
double tmp; double tmp;
sig_flag = 0; sig_flag = 0;
method = 0; method = 0;
@@ -32,7 +31,9 @@ int main()
return 1; return 1;
} }
printf("Введите интервал [a, b] и точность eps: "); show_grapgh();
printf("Введите интервал [a, b] и точность eps и нажмите Enter: ");
scanf("%lf %lf %lf", pa, pb, &eps); scanf("%lf %lf %lf", pa, pb, &eps);
__fpurge(stdin); __fpurge(stdin);
@@ -65,7 +66,7 @@ int main()
{ {
printf("\nКорень уравнения: %lf\n", *pa); printf("\nКорень уравнения: %lf\n", *pa);
printf("Работа программы завершена.\n"); printf("Работа программы завершена.\n");
exit(0); exit_777();
} }
else else
{ {
BIN
View File
Binary file not shown.