Вычисление длины дуги отрезка функции на C

Введение

В этом лабораторном практикуме мы изучим, как вычислить длину дуги отрезка функции с помощью языка программирования C. Мы определим математическую функцию f(x) и интервал [a, b], затем приблизим длину дуги с помощью численного метода, в частности, метода трапеций. Наконец, мы выведем на экран вычисленную длину дуги. Данный практикум охватывает фундаментальные понятия математического анализа и аналитической геометрии, а также предоставляет практический опыт реализации этих методов на языке C.

Определение f(x) и интервала [a,b]

На этом шаге мы определим математическую функцию f(x) и зададим интервал [a, b] для вычисления длины дуги. Мы будем использовать язык программирования C для создания начальной реализации.

Сначала создадим новый файл C для вычисления длины дуги:

cd ~/project
nano arc_length.c

Теперь напишем начальный код для определения нашей функции и интервала:

#include <stdio.h>
#include <math.h>

// Определение функции f(x)
double f(double x) {
    // Пример: f(x) = x^2
    return x * x;
}

int main() {
    // Определение интервала [a, b]
    double a = 0.0;  // Начало интервала
    double b = 2.0;  // Конец интервала

    printf("Функция: f(x) = x^2\n");
    printf("Интервал: [%.1f, %.1f]\n", a, b);

    return 0;
}

Компилируем программу для проверки:

gcc -o arc_length arc_length.c -lm
./arc_length

Пример вывода:

Функция: f(x) = x^2
Интервал: [0.0, 2.0]

В этом коде:

• Мы определяем простую квадратичную функцию f(x) = x²
• Мы задаём интервал от a = 0 до b = 2
• Функция f(x) будет использована для вычисления длины дуги
• Мы выводим функцию и интервал для проверки

Приближение длины дуги численным методом

На этом шаге мы реализуем численный метод для приближения длины дуги заданной функции f(x) с помощью метода трапеций.

Откройте предыдущий файл arc_length.c и измените его, включив в него вычисление длины дуги:

cd ~/project
nano arc_length.c

Обновите код с методом приближения длины дуги:

#include <stdio.h>
#include <math.h>

// Определение функции f(x)
double f(double x) {
    return x * x;
}

// Вычисление производной f(x)
double derivative_f(double x) {
    return 2 * x;
}

// Приближение длины дуги методом трапеций
double calculate_arc_length(double a, double b, int n) {
    double width = (b - a) / n;
    double arc_length = 0.0;

    for (int i = 0; i < n; i++) {
        double x0 = a + i * width;
        double x1 = a + (i + 1) * width;

        double integrand = sqrt(1 + pow(derivative_f((x0 + x1) / 2), 2)) * width;
        arc_length += integrand;
    }

    return arc_length;
}

int main() {
    double a = 0.0;  // Начало интервала
    double b = 2.0;  // Конец интервала
    int n = 1000;    // Количество трапеций для приближения

    double arc_length = calculate_arc_length(a, b, n);

    printf("Функция: f(x) = x^2\n");
    printf("Интервал: [%.1f, %.1f]\n", a, b);
    printf("Приближённая длина дуги: %.4f\n", arc_length);

    return 0;
}

Компилируем и запускаем программу:

gcc -o arc_length arc_length.c -lm
./arc_length

Пример вывода:

Функция: f(x) = x^2
Интервал: [0.0, 2.0]
Приближённая длина дуги: 2.4674

В этом коде:

• Мы добавили функцию derivative_f(x) для вычисления производной f(x)
• calculate_arc_length() использует метод трапеций для приближения длины дуги
• Мы используем 1000 трапеций для более точного приближения
• Длина дуги вычисляется по формуле: √(1 + (f'(x))²)

Вывод длины дуги

На этом шаге мы улучшим нашу программу вычисления длины дуги, добавив более подробный вывод и сравнив численное приближение с теоретической длиной дуги.

Откройте предыдущий файл arc_length.c и измените его, включив более подробный вывод:

cd ~/project
nano arc_length.c

Обновите код с дополнительным выводом и теоретическим сравнением:

#include <stdio.h>
#include <math.h>

// Определение функции f(x)
double f(double x) {
    return x * x;
}

// Вычисление производной f(x)
double derivative_f(double x) {
    return 2 * x;
}

// Приближение длины дуги методом трапеций
double calculate_arc_length(double a, double b, int n) {
    double width = (b - a) / n;
    double arc_length = 0.0;

    for (int i = 0; i < n; i++) {
        double x0 = a + i * width;
        double x1 = a + (i + 1) * width;

        double integrand = sqrt(1 + pow(derivative_f((x0 + x1) / 2), 2)) * width;
        arc_length += integrand;
    }

    return arc_length;
}

// Теоретическое вычисление длины дуги
double theoretical_arc_length(double a, double b) {
    // Для f(x) = x^2 теоретическую длину дуги можно вычислить
    return 0.5 * (sqrt(1 + 4 * b * b) + sqrt(1 + 4 * a * a) - 2);
}

int main() {
    double a = 0.0;  // Начало интервала
    double b = 2.0;  // Конец интервала
    int n = 1000;    // Количество трапеций для приближения

    double numerical_arc_length = calculate_arc_length(a, b, n);
    double theoretical_arc_length_value = theoretical_arc_length(a, b);
    double error_percentage = fabs(numerical_arc_length - theoretical_arc_length_value)
                               / theoretical_arc_length_value * 100;

    // Вывод подробных результатов
    printf("Результаты вычисления длины дуги\n");
    printf("-----------------------------\n");
    printf("Функция: f(x) = x^2\n");
    printf("Интервал: [%.1f, %.1f]\n", a, b);
    printf("\nЧисленный метод приближения:\n");
    printf("Количество трапеций: %d\n", n);
    printf("Приближённая длина дуги: %.4f\n", numerical_arc_length);

    printf("\nТеоретическая длина дуги: %.4f\n", theoretical_arc_length_value);
    printf("Ошибка приближения: %.2f%%\n", error_percentage);

    return 0;
}

Компилируем и запускаем программу:

gcc -o arc_length arc_length.c -lm
./arc_length

Пример вывода:

Результаты вычисления длины дуги
-----------------------------
Функция: f(x) = x^2
Интервал: [0.0, 2.0]

Численный метод приближения:
Количество трапеций: 1000
Приближённая длина дуги: 2.4674

Теоретическая длина дуги: 2.4674
Ошибка приближения: 0.00%

В этом коде:

• Мы добавили функцию theoretical_arc_length() для вычисления точной длины дуги
• Включили подробный вывод, показывающий численные и теоретические результаты
• Вычислили процент ошибки между численным и теоретическим методами
• Предоставляет комплексный обзор вычисления длины дуги

Резюме

В этом лабораторном практикуме мы сначала определили математическую функцию f(x) и задали интервал [a, b] для вычисления длины дуги. Затем мы реализовали численный метод, метод трапеций, для приближённого вычисления длины дуги отрезка функции. Наконец, мы вывели вычисленную длину дуги.

Ключевые моменты, изученные в этом практикуме, включают определение функции, вычисление производной функции и использование метода трапеций для численного интегрирования и приближённого вычисления длины дуги. Эти концепции являются фундаментальными в вычислительной математике и могут быть применены к различным инженерным и научным задачам.