一、复数的表示
复数的两种表示:
gsl复数结构的声明和部分宏在gsl_complex.h中,方法的声明和另一部分宏在gsl_complex_math.h。复数的表示(结构)有三种,即float型、double型和long double型。dat[0]表示实部,dat[1]表示虚部。应注意的是,在gsl_complex_math.h中,所有的方法都是关于gsl_complex和double的。
1 typedef struct 2 { 3 long double dat[2]; 4 } 5 gsl_complex_long_double; 6 7 typedef struct 8 { 9 double dat[2]; 10 } 11 gsl_complex; 12 13 typedef struct 14 { 15 float dat[2]; 16 } 17 gsl_complex_float;
二、复数的创建
1 gsl_complex gsl_complex_polar (double r, double theta); /* r= r e^(i theta) */ 2 gsl_complex gsl_complex_rect (double x, double y); /* r= real+i*imag */
两个函数对应两种不同的表示方法,但返回的结构是一样的。在gsl_complex.h中定义了用于获取实部、获取虚部的宏,在gsl_complex_math.h声明了获取模和弧度的函数。
用于判断复数是否相等、设置复数、设置实部、设置虚部的宏定义在gsl_complex.h中。
#define GSL_REAL(z) ((z).dat[0]) #define GSL_IMAG(z) ((z).dat[1]) double gsl_complex_arg (gsl_complex z); /* return arg(z), -pi< arg(z) <=+pi */ double gsl_complex_abs (gsl_complex z); /* return |z| */
三、复数的运算
- 1个求弧函数
- 3个取模运算(|z|,|z|^2,log |z|)
- 12个四则运算(复数与复数,复数与实数,复数与虚数)
- 共轭
- 复数倒数
- -z(“复数相反数”)
- 2个开方运算(被开方数为复数或实数,结果都是复数)
- 2个乘方运算(指数为复数或实数)
- exp()
- 3个对数运算(底数为e、10、z)
- 三角函数反三角函数
- 双曲函数反双函数
2018-01-02