当前位置：首页 > news >正文

＜C++学习＞C++ Boost 数学与科学计算教程

news 来源：原创 2025/5/7 0:01:38

C++ Boost 数学与科学计算教程

Boost 提供了强大的数学与科学计算功能模块，包括常用数学函数、特殊函数、矩阵运算、随机数生成器、统计分析工具等。这些工具覆盖了从基本数学操作到复杂科学计算的广泛场景。

1. Boost.Math 简介

Boost.Math 提供了丰富的数学功能，主要包括：

特殊函数：如伽玛函数、贝塞尔函数等。
数值工具：如区间算术、根求解等。
分布和统计：如正态分布、泊松分布等。
精度管理：支持高精度浮点计算。

头文件：

#include <boost/math/constants/constants.hpp>
#include <boost/math/special_functions.hpp>
#include <boost/math/distributions.hpp>

2. 特殊函数

Boost 提供了常见的特殊数学函数，包括伽玛函数、贝塞尔函数、误差函数等。

2.1 伽玛函数

示例：计算伽玛函数

#include <boost/math/special_functions/gamma.hpp>
#include <iostream>int main() {double x = 5.0;// 计算 Gamma(x) 和 ln(Gamma(x))double gamma_value = boost::math::tgamma(x);double log_gamma_value = boost::math::lgamma(x);std::cout << "Gamma(" << x << ") = " << gamma_value << std::endl;std::cout << "ln(Gamma(" << x << ")) = " << log_gamma_value << std::endl;return 0;
}

输出

Gamma(5) = 24
ln(Gamma(5)) = 3.17805

2.2 贝塞尔函数

示例：计算贝塞尔函数

#include <boost/math/special_functions/bessel.hpp>
#include <iostream>int main() {double nu = 0.0; // 阶数double x = 1.0;  // 自变量double bessel_value = boost::math::cyl_bessel_j(nu, x); // 一类贝塞尔函数std::cout << "J_" << nu << "(" << x << ") = " << bessel_value << std::endl;return 0;
}

输出

J_0(1) = 0.765198

2.3 误差函数

示例：计算误差函数

#include <boost/math/special_functions/erf.hpp>
#include <iostream>int main() {double x = 1.0;double erf_value = boost::math::erf(x);double erfc_value = boost::math::erfc(x); // 互补误差函数std::cout << "erf(" << x << ") = " << erf_value << std::endl;std::cout << "erfc(" << x << ") = " << erfc_value << std::endl;return 0;
}

输出

erf(1) = 0.842701
erfc(1) = 0.157299

3. 数学常量

Boost 提供了一些常用的数学常量，如 π（pi）、e（自然对数的底）。

示例：使用数学常量

#include <boost/math/constants/constants.hpp>
#include <iostream>int main() {constexpr double pi = boost::math::constants::pi<double>();constexpr double e = boost::math::constants::e<double>();std::cout << "Pi = " << pi << std::endl;std::cout << "e = " << e << std::endl;return 0;
}

输出

Pi = 3.14159
e = 2.71828

4. 分布与统计

Boost 提供了对概率分布和统计分析的支持，包括正态分布、泊松分布等。

4.1 正态分布

示例：正态分布概率密度函数

#include <boost/math/distributions/normal.hpp>
#include <iostream>int main() {boost::math::normal_distribution<> normal_dist(0.0, 1.0); // 均值 0，标准差 1double x = 1.0;double pdf = boost::math::pdf(normal_dist, x); // 概率密度函数值double cdf = boost::math::cdf(normal_dist, x); // 累积分布函数值std::cout << "PDF at x = " << x << " : " << pdf << std::endl;std::cout << "CDF at x = " << x << " : " << cdf << std::endl;return 0;
}

输出

PDF at x = 1 : 0.241971
CDF at x = 1 : 0.841345

4.2 泊松分布

示例：泊松分布的概率计算

#include <boost/math/distributions/poisson.hpp>
#include <iostream>int main() {boost::math::poisson_distribution<> poisson_dist(3.0); // 平均值 λ = 3double k = 2.0;double pmf = boost::math::pdf(poisson_dist, k); // 概率质量函数double cdf = boost::math::cdf(poisson_dist, k); // 累积分布函数std::cout << "PMF at k = " << k << " : " << pmf << std::endl;std::cout << "CDF at k = " << k << " : " << cdf << std::endl;return 0;
}

输出

PMF at k = 2 : 0.224042
CDF at k = 2 : 0.42319

5. 精度管理

Boost 支持高精度浮点运算，可以指定计算结果的精度。

示例：高精度计算 π 的平方

#include <boost/multiprecision/cpp_dec_float.hpp>
#include <boost/math/constants/constants.hpp>
#include <iostream>int main() {using namespace boost::multiprecision;cpp_dec_float_50 pi = boost::math::constants::pi<cpp_dec_float_50>();cpp_dec_float_50 pi_squared = pi * pi;std::cout << "Pi squared (50 digits precision): " << pi_squared << std::endl;return 0;
}

输出（部分）

Pi squared (50 digits precision): 9.869604401089358618834490999876151135313699407240790626413

6. 区间算术

Boost.Math 提供区间算术支持，可用于误差传播和范围计算。

示例：区间相加

#include <boost/numeric/interval.hpp>
#include <iostream>int main() {using namespace boost::numeric;interval<double> a(1.0, 2.0);interval<double> b(3.0, 4.0);interval<double> c = a + b;std::cout << "Interval result: [" << c.lower() << ", " << c.upper() << "]" << std::endl;return 0;
}

输出

Interval result: [4, 6]

7. 综合示例：正态分布与误差函数结合

以下示例计算正态分布在某个区间内的概率，结合误差函数进行精确计算。

#include <boost/math/distributions/normal.hpp>
#include <boost/math/special_functions/erf.hpp>
#include <iostream>int main() {double mean = 0.0;double stddev = 1.0;boost::math::normal_distribution<> dist(mean, stddev);double lower = -1.0;double upper = 1.0;double probability = boost::math::cdf(dist, upper) - boost::math::cdf(dist, lower);std::cout << "Probability in range [" << lower << ", " << upper << "] = " << probability << std::endl;return 0;
}

输出

Probability in range [-1, 1] = 0.682689

8. 学习建议

熟悉特殊函数：
- 学习常用函数如伽玛函数、贝塞尔函数。
- 掌握误差函数在概率和统计中的应用。
了解分布与统计分析：
- 学会使用正态分布、泊松分布等常见分布的计算。
尝试高精度计算：
- 在科学计算中，探索 boost::multiprecision 的高精度运算。
参考文档：
- Boost.Math 官方文档

通过系统学习和实践，Boost 的数学与科学计算工具可以大大提升你在复杂场景下的开发效率！

＜C++学习＞C++ Boost 数学与科学计算教程

C Boost 数学与科学计算教程 Boost 提供了强大的数学与科学计算功能模块，包括常用数学函数、特殊函数、矩阵运算、随机数生成器、统计分析工具等。这些工具覆盖了从基本数学操作到复杂科学计算的广泛场景。 1. Boost.Math 简介 Boost.Math 提供了丰富的数学功能&a…...

编程日记 2025/5/7 0:01:38

Spring Framework 5.3.x源码构建（jdk-1.8, gradle 7.5.1, idea2024.3）

1、下载jdk安装并配置环境变量（自行百度） https://www.oracle.com/java/technologies/downloads/#java8 2、下载spring-framework源码，切换分支到5.3.x https://github.com/spring-projects/spring-framework.git 备用地址 https://gitco…...

编程日记 2025/5/3 21:13:23

有心力场的两体问题

有心力场中的两体问题是经典力学中的重要研究对象，中心力场问题通常涉及两个相互作用的物体（例如行星与恒星、电子与原子核等）。为了简化分析，问题往往可以转化为一个等效的单体问题。这种方法大大提高了问题的可解性，是解决两体和多体问题的基础步骤之一。将两体问题简化…...

编程日记 2025/5/5 18:03:28

【Pandas】pandas Series rdiv

Pandas2.2 Series Binary operator functions 方法描述Series.add()用于对两个 Series 进行逐元素加法运算Series.sub()用于对两个 Series 进行逐元素减法运算Series.mul()用于对两个 Series 进行逐元素乘法运算Series.div()用于对两个 Series 进行逐元素除法运算Series.true…...

编程日记 2025/5/3 21:30:03

Openstack持久存储-Swift,Cinder,Manila三者之间的区别

总结不易，给个三连吧！！！ 补充： 文件共享存储服务Manila 在OpenStack生态系统中，Cinder和Manila分别提供了两种不同类型的存储服务，类似于传统的SAN（存储区域网络）和NAS&…...

编程日记 2025/4/29 15:44:52

单片机（MCU）-简单认识

简介： 内部集成了CPU，RAM，ROM，定时器，中断系统，通讯接口等一系列电脑的常用硬件功能。单片机的任务是信息采集（依靠传感器），处理（依靠CPU）&…...

编程日记 2025/4/27 4:16:22

Redis常见知识点

1、什么是缓存穿透？ 缓存穿透是指查询一定某个key是否存在，每次不存在都查询DB会把DB压垮。解决方案的话，我们通常都会用布隆过滤器来解决。布隆过滤器：在查找一个数是否在一个集合中使用的，通过对数组长度取模&a…...

编程日记 2025/5/3 11:59:15

SQL进阶实战技巧：统计相同时刻多地登陆的用户？

目录 0 问题描述 1 数据准备 2 代码实现 3 问题拓展 3.1 查询每个用户登录日期的最大空档期...

编程日记 2025/5/3 23:31:16

【数据库初阶】表的增删改语句

🎉博主首页： 有趣的中国人 🎉专栏首页： 数据库初阶 🎉其它专栏： C初阶 | C进阶 | 初阶数据结构亲爱的小伙伴们，大家好！在这篇文章中，我们将深入浅出地为大家讲解 MySQL…...

编程日记 2025/5/6 20:28:36

水水水水水水

为了拿推广卷，但不想把我原本完整的文章拆成零散的多篇，只能出此下策随便发一篇，认真写的都笔记专栏里网络技术：数字时代的基础设施在当今社会，网络技术无疑是推动现代生活和经济发展的核心动力之一。从简单的信息传…...

编程日记 2025/4/30 4:39:09

基于STM32的智能电表可视化设计：ESP8266、AT指令集、python后端Flask（代码示例）

一、项目概述随着智能家居的普及，智能电表作为家庭用电管理的重要工具，能够实时监测电流、电压及功率，并将数据传输至后台进行分析和可视化。本项目以STM32C8T6为核心，结合交流电压电流监测模块、ESP8266 Wi-Fi模块、OLED显示屏…...

编程日记 2025/5/6 3:14:34

SpringBoot的@Scheduled和@Schedules有什么区别

Scheduled 的详细解析参数详解 cron: 使用Cron表达式来指定复杂的调度模式。Cron表达式的格式如下： 秒（0-59）分钟（0-59）小时（0-23）日（1-31）月（1-12 或 JAN-…...

编程日记 2025/5/6 11:33:48

Qiskit快速编程探索（进阶篇）

五、量子电路模拟：探索量子世界的虚拟实验室 5.1 Aer模拟器：强大的模拟引擎在量子计算的探索旅程中，Aer模拟器作为Qiskit的核心组件之一，宛如一座功能强大的虚拟实验室，为开发者提供了在经典计算机上模拟量子电路运行的卓越能力。它打破了硬件条件的限制，使得研究者无…...

编程日记 2025/5/6 3:59:51

【漫话机器学习系列】043.提前停止训练（Early Stopping）

提前停止训练（Early Stopping） 提前停止（Early Stopping） 是一种在训练机器学习模型（尤其是深度学习模型）时常用的正则化技术，用于防止过拟合并提升模型的泛化能力。它通过监控验证集的性能&am…...

编程日记 2025/5/6 15:07:39

Linux(上）：基本知识篇

提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档文章目录前言一、Linux初识1 Linux简介2 Linux学习环境配置（1）安装Linux（2）FinalShell远程连接Linux服务器二、Linux基础命令1 Linux目录结构，根目录 /2 Linux命令基础（1）什么是命令、命令行？（2）…...

编程日记 2025/5/1 1:31:03

Python爬虫与1688图片搜索API接口：深度解析与显著收益

在电子商务的浩瀚海洋中，数据是驱动业务决策的核心引擎。阿里巴巴旗下的1688平台，作为全球领先的B2B在线市场，不仅汇聚了海量的商品信息，还提供了丰富的API接口，为开发者提供了强大的数据获取工具。本文将深入探讨1688…...

编程日记 2025/5/4 8:37:04

生物医学信号处理--随机信号通过线性时不变系统

本章主要讨论 (1)输出、输入间自相关函数和功率谱的关系以及两者间的互相关函数和互谱，把它们和系统的冲激响应与频率特性联系起来。这些是用于随机问题的基本关系，将分别讨论连续时间和离散时间两种情况。 (2)初步介绍一些线性系统在处理随机信号时的应用。但是，应该指出，…...

编程日记 2025/5/6 4:29:59

《HeadFirst设计模式》笔记（上）

设计模式的目录： 1 设计模式介绍要不断去学习如何利用其它开发人员的智慧与经验。学习前人的正统思想。我们认为《Head First》的读者是一位学习者。一些Head First的学习原则： 使其可视化将文字放在相关图形内部或附近，而不是放在底部…...

编程日记 2025/5/5 3:46:57

Ubuntu更改内核

需求背景： 由于软件需要在较低版本或者指定版本才可以运行版本： 配置文件： vi /etc/default/grub 启动界面： 可运行版本： 解决方案： 方案1、更改启动顺序 sudo vi /etc/default/grub 方案2、调整启动顺…...

编程日记 2025/5/2 14:09:06

广告公司咋找客户？怎么获取目标企业的采购部联系方式

在广告行业，获取目标企业采购部的联系方式是开展业务、拓展客户的关键一步。分享一些实用的方法，希望能帮到正在为获取联系方式而发愁的广告行业的朋友们。一、利用官方网站和社交媒体平台 1. 官网大多数企业的官方网站都会提供一些联系方式&#xf…...

编程日记 2025/5/6 23:17:18

个人在技术领导力方面的自我反思与提升

大家好！我是 [数擎 AI]，一位热爱探索新技术的前端开发者，在这里分享前端和 Web3D、AI 技术的干货与实战经验。如果你对技术有热情，欢迎关注我的文章，我们一起成长、进步！ 开发领域：前端开发 | A…...

编程日记 2025/5/2 14:10:23

鸿蒙面试 2025-01-11

ArkTs 和TS的关系？ ArkTS（方舟开发语言）与 TypeScript（TS）存在紧密联系，同时也有显著区别： 联系语法基础：ArkTS 在语法层面大量借鉴了 TypeScript ，TypeScript 里诸如…...

编程日记 2025/5/4 22:17:31

Vim的使用方法

Vim的使用方法来自Linux 日常操作与基础知识 | archlinux 简明指南终端编辑器 vim 的使用我们需要掌握一个能在终端中进行文本编辑的软件，这里介绍 vim。创建并编辑名为 hello.txt 的文件： vim hello.txt此时可以看到进入了一个空的界面&#xf…...

编程日记 2025/5/4 8:08:05

什么是卷积网络中的平移不变性？平移shft在数据增强中的意义

今天来介绍一下数据增强中的平移shft操作和卷积网络中的平移不变性。 1、什么是平移 Shift 平移是指在数据增强（data augmentation）过程中，通过对输入图像或目标进行位置偏移（平移），让目标在图像中呈现出…...

编程日记 2025/5/2 13:37:04

Java基础：equals()方法与==的区别

1、超类Object的equals()底层原理： 在Object超类中已经提供了equals()方法，源码如下： public boolean equals(Object obj) { return (this obj); } 所有的对象都拥有标识（内存地址）和状态（数据&a…...

编程日记 2025/5/5 18:03:45

备战蓝桥杯队列和queue详解

目录队列的概念队列的静态实现总代码 stl的queue 队列算法题 1.队列模板题 2.机器翻译 3.海港双端队列队列的概念和栈一样，队列也是一种访问受限的线性表，它只能在表头位置删除，在表尾位置插入，队列是先进先出&…...

编程日记 2025/5/6 6:05:50

《分布式光纤测温：解锁楼宇安全的 “高精度密码”》

在楼宇建筑中，因其内部空间庞大，各类电器设施众多，如何以一种既高效又稳定，兼具低成本与高覆盖特性的方式，为那些关键线路节点开展温度监测，是目前在安全监测领域一项重点研究项目，而无锡布里渊…...

编程日记 2025/5/2 14:07:16

每日一题（一）：识别字符串中的字符是否唯一

目录一、前言二、题目三、题目分析 （一）明确需求 （二）分析思路四、将思路转变为一个程序 C代码 C代码注释： 五、总结一、前言作为博主在本系列文章中的第一篇文章，博主先对本系列做一个简单的介绍…...

编程日记 2025/5/4 5:52:07

【翻译】2025年华数杯国际赛数学建模题目+翻译pdf自取

保存至本地网盘链接：https://pan.quark.cn/s/f82a1fa7ed87 提取码：6UUw 2025年“华数杯”国际大学生数学建模竞赛比赛时间于2025年1月11日（周六）06:00开始，至1月15日（周三）09:00结束&#xff…...

编程日记 2025/5/4 15:43:52

如何稳定使用 O1 / O1 Pro，让“降智”现象不再困扰？

近期，不少朋友在使用 O1 或 O1 Pro 模型时，都会碰到“降智”或“忽高忽低”的智力波动，比如无法识图、无法生成图片、甚至回答准确度也不稳定。面对这些问题，你是不是也感到头疼呢？ 为了找到更可靠的解决办法&#xf…...

编程日记 2025/5/4 23:35:57

【Linux】4.Linux常见指令以及权限理解（2）

文章目录 3. Linux指令3.1 ls指令和rm指令补充3.2 man指令（重要）3.3cp指令（重要）输出重定向3.3.1ubuntu20.04如何安装tree 3.4 mv指令（重要）mv指令更改文件名mv指令更改目录名如何看待指令指令的重命名3.5…...

编程日记 2025/5/2 15:35:07

神州数码交换机和路由器命令总结

神州数码交换机和路由器命令总结一、神州数码交换机命令总结 1. 交换机恢复出厂设置及其基本配置. 1) //进入特权模式 2) del startup.cfg 2. Telnet方式管理交换机. 1) //进入全局配置模式 2) enable password 0 [密码] 3) Line 0 4 4) Password 0 [密码] 5) Login 3. 交换机…...

编程日记 2025/5/2 14:39:30

echarts: 双柱柱状图同时两侧双y轴问题

问： 现在又有一个问题：echarts双柱柱状图，同时左侧的y轴有刻度，右侧的y轴也有刻度，双y轴，y轴每一个刻度都对应虚线，这时候会有一个问题，左侧的y轴的虚线和右侧y轴的虚线没对应&…...

编程日记 2025/5/5 11:40:21

计算机网络八股文学习笔记

文章目录计算机网络基础网络分层模型OSI七层模型TCP/IP四层模型 HTTP从输入URL到页面展示到底发生了什么?(非常重要)HTTP状态码HTTP Header中常见的字段有哪些?HTTP和HTTPS有什么区别?(重要)HTTP/1.0和HTTP/1.1有什么区别?HTTP/1.1和HTTP/2.0有什么区别?HTTP/2.0和HTTP/3…...

编程日记 2025/5/4 1:11:16

java中json字符串键值获取

<dependency><groupId>com.alibaba</groupId><artifactId>fastjson</artifactId><version>1.2.83</version> </dependency>使用fastjson依赖 JSONObject jsonObject JSON.parseObject(s); 这个jsonObject本质就是一个map&…...

编程日记 2025/5/2 14:13:24

MySQL的安装

MySQL典型的关系型数据库（RDBMS）：oracle、MySQL、SqlServer MySQL的版本 5.5~5.7、8.0 MySQL的安装和配置下载地址： https://downloads.mysql.com/archives/community/ 安装包 (x86, 64-bit), MSI Installer 执行下一步即…...

编程日记 2025/5/2 13:57:54

Open FPV VTX开源之第一次出图

Open FPV VTX开源之第一次出图 1. 源由2. 连线2.1 飞控2.2 调试 3. serial3.1 启动log - uboot3.2 登录版本 - linux3.3 获取有线IP 4. ssh - linux5. PixelPilot出图6. 总结7. 参考资料8. 补充 - 8812AU网卡 1. 源由在《Open FPV VTX开源之硬件规格及组成》章节中&#xff0…...

编程日记 2025/5/2 14:06:11

R语言在森林生态研究中的魔法：结构、功能与稳定性分析——发现数据背后的生态故事！

森林生态系统结构、功能与稳定性分析与可视化研究具有多方面的重要意义，具体如下： 一、理论意义 ●深化生态学理论通过研究森林生态系统的结构、功能与稳定性，可以深化对生态系统基本理论的理解。例如，生物多样性与生态系统稳定性…...

编程日记 2025/5/3 17:35:15

qt QLabel QPushButton 控件重写paintEvent后控件消失

qt 继承自PushButton控件的类重写paintEvent后控件消失解决办法，在paintevent结尾加上这条语句：QPushButton::paintEvent(event); void MyButton::paintEvent(QPaintEvent *event) {QPushButton::paintEvent(event); } 这里QPushButton不能写成Q…...

编程日记 2025/5/4 13:38:24

安全基础-互联网技术基础

互联网技术基础概述：计算机网络、网络协议、HTTP协议、前端与后端技术、Web服务器、数据库以及浏览器等目录互联网技术基础前言一、计算机网络定义二、网络协议和协议分层1.OSI七层模型2.TCP/IP四层模型三、HTTP协议1、HTTP协议的特点2、HTTP请求3、HTTP响应4、…...

编程日记 2025/5/4 7:06:35

【Rust自学】11.9. 单元测试

喜欢的话别忘了点赞、收藏加关注哦（加关注即可阅读全文），对接下来的教程有兴趣的可以关注专栏。谢谢喵！(･ω･) 11.9.1. 测试的分类 Rust把测试分为两类，一个是单元测试，一个是集成…...

编程日记 2025/5/2 1:17:14

了解Node.js

Node.js是一个基于V8引擎的JavaScript运行时环境，它允许JavaScript代码在服务器端运行，从而实现后端开发。Node.js的出现，使得前端开发人员可以利用他们已经掌握的JavaScript技能，扩展技能树并成为全栈开发人员。本文将深入浅出地…...

编程日记 2025/5/4 9:56:12

二十三种设计模式-原型模式

原型模式（Prototype Pattern）是一种创建型设计模式，它通过拷贝现有的实例来创建新的实例，而不是通过新建实例。这种方式可以避免复杂的构造过程，同时还能保持对象的创建和使用分离，提高系统的灵活性和扩展性…...

编程日记 2025/5/2 14:03:02

蓝桥与力扣刷题（66 加一）

题目： 给定一个由整数组成的非空数组所表示的非负整数，在该数的基础上加一。最高位数字存放在数组的首位， 数组中每个元素只存储单个数字。你可以假设除了整数 0 之外，这个整数不会以零开头。示例 1： 输入…...

编程日记 2025/5/2 14:07:03

Markdown 分级标题

Markdown 分级标题 1. 前言2. 语法详解2.1 SeText风格2.2 Atx风格 3. 使用场景及实例4. 小结5. 其他风格5. 其他文章快来试试吧🎈 Markdown 分级标题 👈点击这里也可查看 1. 前言关于 Markdown 的分级标题，官方给出了如下定义： …...

编程日记 2025/5/5 23:37:30

回归预测 | MATLAB基于RF-Adaboost多输入单输出回归预测

回归预测 | MATLAB基于RF-Adaboost多输入单输出回归预测目录回归预测 | MATLAB基于RF-Adaboost多输入单输出回归预测预测效果基本介绍程序设计参考资料预测效果基本介绍回归预测 | MATLAB基于RF-Adaboost多输入单输出回归预测。 1.Matlab实现RF-Adaboost随机森林集成学习…...

编程日记 2025/5/2 14:34:07