级数学习笔记

一、数学基础

1. 数项级数（Number Series）

数项级数是指形如：

∑(n=1 to ∞) aₙ = a₁ + a₂ + a₃ + ...

的无穷和。

1.1 收敛性判别法

比较判别法
比值判别法
根值判别法
积分判别法
莱布尼茨判别法（交错级数）

2. 幂级数（Power Series）

幂级数是指形如：

∑(n=0 to ∞) aₙ(x-x₀)ⁿ = a₀ + a₁(x-x₀) + a₂(x-x₀)² + ...

的级数。

2.1 收敛半径

对于幂级数，存在收敛半径R，使得：

当|x-x₀| < R时，级数绝对收敛
当|x-x₀| > R时，级数发散
当|x-x₀| = R时，需要单独判断

3. 泰勒级数（Taylor Series）

函数f(x)在点x₀处的泰勒级数展开：

f(x) = ∑(n=0 to ∞) [f⁽ⁿ⁾(x₀)/n!] (x-x₀)ⁿ

4. 傅里叶级数（Fourier Series）

周期函数f(x)的傅里叶级数展开：

f(x) = a₀/2 + ∑(n=1 to ∞) [aₙcos(nx) + bₙsin(nx)]

其中：

aₙ = (1/π)∫(-π to π) f(x)cos(nx)dx
bₙ = (1/π)∫(-π to π) f(x)sin(nx)dx

二、代码实现

1. 头文件 (series.h)

#ifndef SERIES_H
#define SERIES_H#include <functional>
#include <vector>
#include <cmath>
#include <stdexcept>/*** @class Series* @brief 级数计算类* * 用于计算各种级数的和、收敛性等。*/
class Series {
public:// 计算数项级数的部分和static double partialSum(const std::function<double(int)>& term,int n);// 判断数项级数的收敛性static bool isConvergent(const std::function<double(int)>& term,double epsilon = 1e-6,int max_terms = 1000);// 计算幂级数的收敛半径static double convergenceRadius(const std::function<double(int)>& coefficient,double epsilon = 1e-6,int max_terms = 100);// 计算幂级数的和static double powerSeriesSum(const std::function<double(int)>& coefficient,double x,double x0,int n);// 计算泰勒级数展开static std::vector<double> taylorSeries(const std::function<double(double)>& func,const std::function<double(double, int)>& derivative,double x0,int n);// 计算傅里叶级数系数static std::pair<std::vector<double>, std::vector<double>> fourierCoefficients(const std::function<double(double)>& func,int n,double period = 2 * M_PI);// 计算傅里叶级数的和static double fourierSeriesSum(const std::vector<double>& a,const std::vector<double>& b,double x,int n);private:// 辅助函数：计算阶乘static double factorial(int n);// 辅助函数：计算组合数static double combination(int n, int k);// 辅助函数：计算积分static double integrate(const std::function<double(double)>& func,double a,double b,int n = 1000);
};#endif // SERIES_H

2. 实现文件 (series.cpp)

#include "series.h"
#include <algorithm>
#include <numeric>double Series::partialSum(const std::function<double(int)>& term,int n
) {double sum = 0.0;for (int i = 1; i <= n; ++i) {sum += term(i);}return sum;
}bool Series::isConvergent(const std::function<double(int)>& term,double epsilon,int max_terms
) {double sum = 0.0;double prev_sum = 0.0;for (int i = 1; i <= max_terms; ++i) {sum += term(i);if (std::abs(sum - prev_sum) < epsilon) {return true;}prev_sum = sum;}return false;
}double Series::convergenceRadius(const std::function<double(int)>& coefficient,double epsilon,int max_terms
) {double r = 0.0;double prev_r = 0.0;for (int n = 1; n <= max_terms; ++n) {double a_n = coefficient(n);double a_n_plus_1 = coefficient(n + 1);if (std::abs(a_n_plus_1) < epsilon) {r = std::numeric_limits<double>::infinity();break;}r = std::abs(a_n / a_n_plus_1);if (std::abs(r - prev_r) < epsilon) {break;}prev_r = r;}return r;
}double Series::powerSeriesSum(const std::function<double(int)>& coefficient,double x,double x0,int n
) {double sum = 0.0;double power = 1.0;for (int i = 0; i <= n; ++i) {sum += coefficient(i) * power;power *= (x - x0);}return sum;
}std::vector<double> Series::taylorSeries(const std::function<double(double)>& func,const std::function<double(double, int)>& derivative,double x0,int n
) {std::vector<double> coefficients(n + 1);for (int i = 0; i <= n; ++i) {coefficients[i] = derivative(x0, i) / factorial(i);}return coefficients;
}std::pair<std::vector<double>, std::vector<double>> Series::fourierCoefficients(const std::function<double(double)>& func,int n,double period
) {std::vector<double> a(n + 1);std::vector<double> b(n + 1);// 计算a₀a[0] = integrate(func, -period/2, period/2) / period;// 计算aₙ和bₙfor (int i = 1; i <= n; ++i) {auto cos_term = [&](double x) {return func(x) * std::cos(2 * M_PI * i * x / period);};auto sin_term = [&](double x) {return func(x) * std::sin(2 * M_PI * i * x / period);};a[i] = 2 * integrate(cos_term, -period/2, period/2) / period;b[i] = 2 * integrate(sin_term, -period/2, period/2) / period;}return {a, b};
}double Series::fourierSeriesSum(const std::vector<double>& a,const std::vector<double>& b,double x,int n
) {double sum = a[0] / 2;for (int i = 1; i <= n; ++i) {sum += a[i] * std::cos(2 * M_PI * i * x) +b[i] * std::sin(2 * M_PI * i * x);}return sum;
}double Series::factorial(int n) {if (n < 0) {throw std::invalid_argument("Factorial is not defined for negative numbers");}double result = 1.0;for (int i = 2; i <= n; ++i) {result *= i;}return result;
}double Series::combination(int n, int k) {if (k < 0 || k > n) {return 0.0;}double result = 1.0;for (int i = 1; i <= k; ++i) {result *= (n - k + i) / i;}return result;
}double Series::integrate(const std::function<double(double)>& func,double a,double b,int n
) {double h = (b - a) / n;double sum = (func(a) + func(b)) / 2;for (int i = 1; i < n; ++i) {sum += func(a + i * h);}return h * sum;
}

3. 示例程序 (series_demo.cpp)

#include "series.h"
#include <iostream>
#include <iomanip>
#include <cmath>// 打印向量
void printVector(const std::vector<double>& v, const std::string& name) {std::cout << name << " = [";for (size_t i = 0; i < v.size(); ++i) {std::cout << std::fixed << std::setprecision(6) << v[i];if (i < v.size() - 1) std::cout << ", ";}std::cout << "]" << std::endl;
}int main() {try {// 1. 数项级数示例// 计算调和级数的部分和auto harmonic_term = [](int n) { return 1.0 / n; };double harmonic_sum = Series::partialSum(harmonic_term, 10);std::cout << "调和级数的前10项和 = " << harmonic_sum << std::endl;// 判断调和级数的收敛性bool is_harmonic_convergent = Series::isConvergent(harmonic_term);std::cout << "调和级数是否收敛: " << (is_harmonic_convergent ? "是" : "否") << std::endl;// 2. 幂级数示例// 计算几何级数的收敛半径auto geometric_coef = [](int n) { return 1.0; };double radius = Series::convergenceRadius(geometric_coef);std::cout << "几何级数的收敛半径 = " << radius << std::endl;// 计算几何级数的和double geometric_sum = Series::powerSeriesSum(geometric_coef, 0.5, 0.0, 10);std::cout << "几何级数在x=0.5处的和 = " << geometric_sum << std::endl;// 3. 泰勒级数示例// 计算e^x的泰勒级数展开auto exp_func = [](double x) { return std::exp(x); };auto exp_derivative = [](double x, int n) { return std::exp(x); };auto taylor_coef = Series::taylorSeries(exp_func, exp_derivative, 0.0, 5);printVector(taylor_coef, "e^x的泰勒级数系数");// 4. 傅里叶级数示例// 计算方波的傅里叶级数auto square_wave = [](double x) {return std::abs(std::fmod(x, 2 * M_PI)) < M_PI ? 1.0 : -1.0;};auto [a, b] = Series::fourierCoefficients(square_wave, 5);printVector(a, "方波的傅里叶系数a");printVector(b, "方波的傅里叶系数b");// 计算方波的傅里叶级数和double fourier_sum = Series::fourierSeriesSum(a, b, M_PI/4, 5);std::cout << "方波在x=π/4处的傅里叶级数和 = " << fourier_sum << std::endl;} catch (const std::exception& e) {std::cerr << "错误: " << e.what() << std::endl;return 1;}return 0;
}

4. Makefile

CXX = g++
CXXFLAGS = -std=c++17 -Wall -Wextra -O2all: series_demoseries_demo: series_demo.o series.o$(CXX) $(CXXFLAGS) -o $@ $^series_demo.o: series_demo.cpp series.h$(CXX) $(CXXFLAGS) -c $<series.o: series.cpp series.h$(CXX) $(CXXFLAGS) -c $<clean:rm -f *.o series_demo.PHONY: all clean

三、使用说明

1. 编译

使用提供的Makefile编译项目：

make

2. 运行示例

./series_demo

3. 基本用法

// 计算数项级数的部分和
double sum = Series::partialSum(term_func, n);// 判断级数的收敛性
bool is_convergent = Series::isConvergent(term_func);// 计算幂级数的收敛半径
double radius = Series::convergenceRadius(coef_func);// 计算幂级数的和
double sum = Series::powerSeriesSum(coef_func, x, x0, n);// 计算泰勒级数展开
auto coef = Series::taylorSeries(func, derivative, x0, n);// 计算傅里叶级数系数
auto [a, b] = Series::fourierCoefficients(func, n);// 计算傅里叶级数的和
double sum = Series::fourierSeriesSum(a, b, x, n);