C++实现第一类贝塞尔函数

最新推荐文章于 2022-12-22 14:37:21 发布
最新推荐文章于 2022-12-22 14:37:21 发布
阅读量3k 收藏 6
点赞数 1
分类专栏： 日常学习 文章标签： 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44332298/article/details/108112937
版权
本文介绍了在实习中遇到的复数域下半整阶第一类贝塞尔函数的实现问题。由于递推公式求解存在精度不足的问题，作者选择了直接实现贝塞尔函数的定义公式来确保计算精确。通过求导递推公式，确保了贝塞尔函数的精确计算，特别是对于特殊值J(1/2)和J(-1/2)。函数能够处理任意大小的n值，包括负数和小数，且对于z为复数的情况（将虚部设为0）也适用。在实际应用中，发现循环次数设为15即可达到收敛，优于原先的100次循环设定。
摘要由CSDN通过智能技术生成
第一类贝塞尔函数

这几天实习碰到了半整阶的第一类的贝塞尔函数，而且是复数域的。最开始在网上查资料是用递推公式求解的，因为在计算特殊值J(1/2)和J(-1/2)的过程中会省略掉很多项并不精确，且在递推过程中也会省略很多项，导致最终递推得到的结果并不精确，故而又改为实现贝塞尔函数的定义函数。
定义公式：
在这里插入图片描述
递推公式：
在这里插入图片描述

求导递推公式：
对贝塞尔函数求导是用的如下递推公式，需要保证贝塞尔函数一定要算的精确，求导结果才能精确。在这里插入图片描述
该函数的n可以取任意值，包括负数、小数。因为项目需要，z为复数，把虚部设为0，同样适用于实数。

std::complex<double
最低0.47元/天 解锁文章
IT圈的01羊
关注
专栏目录
第一类贝塞尔函数
09-23
第一类贝塞尔函数，使计算第一类贝塞尔函数更快更方便。
C/C++ 内使用贝塞尔函数
baidu_41901258的博客
03-17 7190
今天在捣鼓MMSE，倒在了贝塞尔函数上，翻找好久，原来C++的cmath里有。 官方文档在这里： https://docs.microsoft.com/en-us/cpp/standard-library/cmath?view=vs-2019#regular-modified-cylindrical-bessel-functions https://docs.microsoft.com/...
贝塞尔函数
08-15
贝塞尔函数贝塞尔函数的有关公式 贝塞尔方程 的持解Bp(z)为(柱)贝塞尔函数。有 第一类柱贝塞尔函数Jp(z) p为整数n时，Jn=(1) nJn； p不为整数时，Jp与Jp线性无关。 第二类柱贝塞尔函数N p(z)(柱诺依曼函数) n为整数时Nn=(1) nNn。 第三类柱贝塞尔函数Hp(z) (柱汉开尔函数)： 第一类柱汉开尔函数 Hp(1)(z)= Jp(z)+j N p(z) 第二类柱汉开尔函数 Hp(2)(z)= Jp(z)j N p(z)
贝塞尔函数求值 C++实现
A_Pointer的博客
08-16 1656
本实现完全基于一般意义下的贝塞尔函数定义，验证了部分阶数为正数的贝塞尔函数值。 此外，本实现需要利用我博客下的那个伽马函数实现。 #include <cmath> #include <complex> //任意阶三的第一类贝塞尔函数 //n 为阶数 //z 为复数（也满足实数求值，仅需设置虚部为0） std::complex<double> JFunction(const double n, const std::complex<double>&amp
贝塞尔函数 c语言程序,贝塞尔函数Bessel（原理与实现）
weixin_39777464的博客
05-20 1855
程序实现：// crt_bessel1.c#include #include int main( void ){double x = 2.387;int n = 3, c;printf( "Bessel functions for x = %f:\n", x );printf( " Kind Order Function Result\n\n" );printf( " First ...
c语言整数阶贝塞尔函数,第一类整数阶贝塞尔函数(贝塞尔函数)
weixin_31887861的博客
05-22 1210
#include #include double Bessel(int n,double x){int i,m;double t,y,z,p,q,s,b0,b1;static double a[6]={ 57568490574.0,-13362590354.0,651619640.7,-11214424.18,77392.33017,-184.9052456};static double b[6]...
第二类修正贝塞尔函数
haodawei123的博客
09-17 9651
1 递推算式 摘自《MATLAB语音常用程序集》 第二类修正贝塞尔函数需要用第一类修正贝塞尔函数来计算 2 vc编程实现 #include <stdio.h> #include <math.h> double first_xzbsl(int n,double x)//第一类修正的贝塞尔函数 { double a[7]={1.0,3.5156229,3.0899424,1.2067492,0.2659732,0.0360768,0.0045813}; double c[7]
Complex Order and Argument的贝塞尔函数：复阶和参数的贝塞尔函数-matlab开发
06-01
贝塞尔函数是一类与圆周率和阶数n有关的无穷级数，分为整数阶（通常称为第一类贝塞尔函数J_n(x)和第二类贝塞尔函数Y_n(x)）和半整数阶（通常称为奈曼函数H_n^(±1)(x)）。在MATLAB中，可以使用`besselj`, `bessely`,...
贝塞尔曲线
12-24
vs结合easyX实现贝塞尔曲线的绘制，并可以通过鼠标拖拽控制点实现曲线的变化
一类含贝塞尔函数积分的数值算法
03-04
一类含贝塞尔函数积分的数值算法，里面是具体数学算法，不是编好的程序，有难度，没信心看懂的绕道。
第一类整数阶贝塞尔函数(贝塞尔函数)
01-05 8373
#include #include double Bessel(int n,double x) { int i,m; double t,y,z,p,q,s,b0,b1; static double a[6]={ 57568490574.0,-13362590354.0, 651619640.7,-11214424.18,77392.33017, -184.9052456}; static double b[6]={ 57568490411.0,102
C++： 计算非整数阶的贝塞尔 J 函数(附完整源码)
希望我的博客，能帮上你解决学习中工作中所遇到的问题
12-22 310
C++： 计算非整数阶的贝塞尔 J 函数(附完整源码)
第一类修正贝塞尔函数的C语言实现
业精于勤，荒于嬉
09-20 3070
第一类修正贝塞尔函数公式 Iν(z)=(z2)ν∑(k=0)∞(z24)kk!Γ(ν+k+1)I_{\nu}(z)=\left(\frac{z}{2}\right)^{\nu} \sum_{(k=0)}^{\infty} \frac{\left(\frac{z^{2}}{4}\right)^{k}}{k ! \Gamma(\nu+k+1)}Iν​(z)=(2z​)ν(k=0)∑∞​k!Γ(ν+k+1)(4z2​)k​ 其中，Γ(a)\Gamma(a)Γ(a)是gamma函数。 常量和结构体定义 常量和结构体
第一类修正Bessel 函数
热门推荐
qq_36465370的博客
10-30 1万+
摘自图书：
C++简单实现Bézier曲线
weixin_52435933的博客
03-08 495
贝塞尔曲线有着很多特殊的性质在图形设计和路径规划中应用都非常广泛。贝塞尔曲线完全由其控制点决定其形状，n个控制点对应着n-1阶的贝塞尔曲线，并且可以通过递归来定义。本篇文章的重点在于对Bézier Curve的理解以及用OpenGL绘制Bézier Curve。 〄 Bézier Curve·1阶给定两个已知坐标的控制点，，那么1阶贝塞尔曲线可以用一个关于t的参数方程来描...
c语言贝塞尔函数J1,贝塞尔函数：_j0、_j1、_jn、_y0、_y1、_yn
weixin_36340890的博客
05-25 2835
贝塞尔函数：_j0、_j1、_jn、_y0、_y1、_yn4/2/2020本文内容计算第一种或第二种贝塞尔函数，顺序为 0、1 或 n。 贝赛耳函数通常用于电磁波理论的数学学科中。语法double _j0(double x);double _j1(double x);double _jn(int n,double x);double _y0(double x);double _y1(double ...
C++求两条贝塞尔曲线的交点
最新发布
05-20
求两条贝塞尔曲线的交点，可以通过求解方程组来实现。假设两条贝塞尔曲线分别为 $B_1(t)$ 和 $B_2(t)$，则它们的交点满足以下方程组： $$ \begin{cases} B_1(t_1) - B_2(t_2) = 0\\ B_1'(t_1) \cdot B_2''(t_2) - B_2'(t_2) \cdot B_1''(t_1) = 0 \end{cases} $$ 其中 $B_i(t)$ 表示第 $i$ 条贝塞尔曲线在参数 $t$ 处的点，$B_i'(t)$ 表示第 $i$ 条贝塞尔曲线在参数 $t$ 处的切向量，$B_i''(t)$ 表示第 $i$ 条贝塞尔曲线在参数 $t$ 处的二阶导向量。 具体实现可以采用牛顿迭代法来求解方程组的根。 以下是一个简单的 C++ 实现： ```cpp #include <iostream> #include <cmath> using namespace std; const double EPSILON = 1e-6; // 停止迭代的精度要求 struct Point { double x; double y; }; struct BezierCurve { Point p0; Point p1; Point p2; Point p3; }; Point bezier(BezierCurve curve, double t) { double u = 1 - t; double uu = u * u; double uuu = uu * u; double tt = t * t; double ttt = tt * t; Point p; p.x = uuu * curve.p0.x + 3 * uu * t * curve.p1.x + 3 * u * tt * curve.p2.x + ttt * curve.p3.x; p.y = uuu * curve.p0.y + 3 * uu * t * curve.p1.y + 3 * u * tt * curve.p2.y + ttt * curve.p3.y; return p; } Point bezierDerivative(BezierCurve curve, double t) { double u = 1 - t; double uu = u * u; double tt = t * t; Point p; p.x = 3 * uu * (curve.p1.x - curve.p0.x) + 6 * u * t * (curve.p2.x - curve.p1.x) + 3 * tt * (curve.p3.x - curve.p2.x); p.y = 3 * uu * (curve.p1.y - curve.p0.y) + 6 * u * t * (curve.p2.y - curve.p1.y) + 3 * tt * (curve.p3.y - curve.p2.y); return p; } Point bezierSecondDerivative(BezierCurve curve, double t) { double u = 1 - t; Point p; p.x = 6 * u * (curve.p2.x - 2 * curve.p1.x + curve.p0.x) + 6 * t * (curve.p3.x - 2 * curve.p2.x + curve.p1.x); p.y = 6 * u * (curve.p2.y - 2 * curve.p1.y + curve.p0.y) + 6 * t * (curve.p3.y - 2 * curve.p2.y + curve.p1.y); return p; } double distance(Point a, Point b) { double dx = a.x - b.x; double dy = a.y - b.y; return sqrt(dx * dx + dy * dy); } Point findIntersection(BezierCurve curve1, BezierCurve curve2) { double t1 = 0.5; // 初始化迭代参数 double t2 = 0.5; double dt = 1.0; while (dt > EPSILON) { // 迭代求解方程组 Point p1 = bezier(curve1, t1); Point p2 = bezier(curve2, t2); double d = distance(p1, p2); if (d < EPSILON) { // 当两条曲线已经相交时，直接返回交点 return p1; } Point v1 = bezierDerivative(curve1, t1); Point v2 = bezierDerivative(curve2, t2); Point a1 = bezierSecondDerivative(curve1, t1); Point a2 = bezierSecondDerivative(curve2, t2); double f1 = p1.x - p2.x; double f2 = p1.y - p2.y; double f3 = v1.x * a2.y - v1.y * a2.x - v2.x * a1.y + v2.y * a1.x; double f4 = a1.x * a2.y - a1.y * a2.x; double d1 = v1.x * a2.y - v1.y * a2.x; double d2 = v2.x * a1.y - v2.y * a1.x; double d3 = a1.x * a2.y - a1.y * a2.x; double h1 = f1 - d1 / d3 * f3; double h2 = f2 - d2 / d3 * f3; double h3 = f4; double det = h1 * h2 - h3 * h3; double t1_new = t1 - (h2 * f3 - h3 * f2) / det; double t2_new = t2 - (h1 * f3 - h3 * f1) / det; if (t1_new < 0 || t1_new > 1 || t2_new < 0 || t2_new > 1) { // 当迭代参数越界时，重新初始化 t1_new = 0.5; t2_new = 0.5; } dt = max(abs(t1_new - t1), abs(t2_new - t2)); t1 = t1_new; t2 = t2_new; } return bezier(curve1, t1); } int main() { BezierCurve curve1 = {{0, 0}, {1, 2}, {3, -1}, {4, 0}}; BezierCurve curve2 = {{2, 0}, {3, 1}, {1, 3}, {0, 2}}; Point p = findIntersection(curve1, curve2); cout << "Intersection point: (" << p.x << ", " << p.y << ")" << endl; return 0; } ``` 在上述代码中，`BezierCurve` 结构体表示一条贝塞尔曲线，包含四个控制点。`bezier` 函数用来计算曲线在参数 $t$ 处的点，`bezierDerivative` 和 `bezierSecondDerivative` 函数分别用来计算曲线在参数 $t$ 处的一阶导向量和二阶导向量。`distance` 函数用来计算两个点之间的距离。`findIntersection` 函数则用来求解两条贝塞尔曲线的交点，采用牛顿迭代法求解方程组的根。
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值