|
|
#pragma once
|
|
|
|
|
|
#include "math.h"
|
|
|
#include "zxHugeNum.h"
|
|
|
|
|
|
// 贝塞尔函数
|
|
|
class I_ALGMATH_EXPORT zxBessel
|
|
|
{
|
|
|
public:
|
|
|
|
|
|
zxBessel();
|
|
|
~zxBessel();
|
|
|
|
|
|
public:
|
|
|
|
|
|
// TODO 注意，对于matlab的bessel函数，其实是有 besseli(0, q1, 1) 这种写法的
|
|
|
// 即：存在 第三个参数 1，即修正bessel方式，
|
|
|
// 如果不输入这个参数（或输入 0）时，besseli 直接计算函数值。
|
|
|
// 当自变量 q1 较大时，函数值会急剧增大，很容易超出双精度浮点数的上限（即“上溢”），
|
|
|
// 导致计算结果变为 Inf 而失效
|
|
|
// 如果输入这个参数 1 之后，MATLAB 会计算 缩放后的值，也就是说，它会乘上一个指数衰减因子
|
|
|
// exp(-abs(x)) 来抵消函数本身的高速增长。
|
|
|
// 即：
|
|
|
// besseli(0, q1) → 等价于 besseli(0, q1, 0) → 不乘衰减因子，直接算原函数值。
|
|
|
// besseli(0, q1, 1) → 才乘衰减因子 exp(-abs(q1))
|
|
|
|
|
|
// 贝塞尔函数（基础函数部分）
|
|
|
bool besselI0(double x, double &b);
|
|
|
double getbesselI0(double x);
|
|
|
|
|
|
bool besselI1(double x, double &b);
|
|
|
double getbesselI1(double x);
|
|
|
|
|
|
bool besselK0(double x, double &b);
|
|
|
double getbesselK0(double x);
|
|
|
|
|
|
bool besselK1(double x ,double &b);
|
|
|
double getbesselK1(double x);
|
|
|
|
|
|
// 贝塞尔函数（组合模式）
|
|
|
// K1/I1
|
|
|
bool besselK1I1(double x, double &b);
|
|
|
|
|
|
// K0/I0
|
|
|
bool besselK0I0(double x, double &b);
|
|
|
|
|
|
// I0/I1
|
|
|
bool besselI0I1(double x, double &b);
|
|
|
|
|
|
// I0(x)/I1(y)
|
|
|
bool besselI0I1(double x, double y, double &b);
|
|
|
|
|
|
// I0(x)/I0(y)
|
|
|
bool besselI0I0(double x, double y, double &b);
|
|
|
|
|
|
// I1(x)/I1(y)
|
|
|
bool besselI1I1(double x, double y, double &b);
|
|
|
|
|
|
// K1(x)/I1(y)
|
|
|
bool besselK1I1(double x, double y, double &b);
|
|
|
|
|
|
// K1(x)/I0(y)
|
|
|
bool besselK1I0(double x, double y, double &b);
|
|
|
|
|
|
// K0(x)/I0(y)
|
|
|
bool besselK0I0(double x, double y, double &b);
|
|
|
|
|
|
// (I0(x)/I1(x))/(I0(y)/I1(y))，暂未启用
|
|
|
bool besselI0I1I0I1(double x, double y, double &b);
|
|
|
|
|
|
// (K0(x)/I0(x))/(K1(y)/I1(y))，暂未启用
|
|
|
bool besselK0I0K1I1(double x, double y, double &b);
|
|
|
|
|
|
// 大数模式
|
|
|
zxHugeNum besselI(int which, double x);
|
|
|
zxHugeNum besselK(int which, double x);
|
|
|
|
|
|
public:
|
|
|
|
|
|
////////////////////////////////////////////////////////////////
|
|
|
// 由于单个的超过数值范围，故一起在内部进行处理
|
|
|
// 仅仅测试
|
|
|
//double calDemoA(double d1, double d2);
|
|
|
//double calDemoB(double d1, double d2);
|
|
|
|
|
|
zxHugeNum besselDivide(bool byK1, int which1, double x1, \
|
|
|
bool byK2, int which2, double x2, \
|
|
|
bool* pOk = nullptr);
|
|
|
zxHugeNum besselMultiple(bool byK1, int which1, double x1, \
|
|
|
bool byK2, int which2, double x2, \
|
|
|
bool* pOk = nullptr);
|
|
|
|
|
|
protected:
|
|
|
zxHugeNum besselDivMul(bool bDivide, \
|
|
|
bool byK1, int which1, double x1, \
|
|
|
bool byK2, int which2, double x2, \
|
|
|
bool* pOk = nullptr);
|
|
|
|
|
|
double _besselByY(double y, bool bK, int which, bool bLower);
|
|
|
|
|
|
public:
|
|
|
|
|
|
void setComMode(bool b) { m_bComMode = b; }
|
|
|
bool isComMode() { return m_bComMode; }
|
|
|
double getResultI() { return m_dIndex; }
|
|
|
double getResultM() { return m_dMultipler; }
|
|
|
|
|
|
private:
|
|
|
|
|
|
bool m_bComMode; //对于大数，采取复合模式
|
|
|
double m_dIndex; //exp的次数
|
|
|
double m_dMultipler; //乘数，合起来是 exp(m_dIndex) * m_dMultipler;
|
|
|
|
|
|
////////////////////////////////////////////////////////////////
|
|
|
|
|
|
};
|
|
|
|