Отправлено fontp 15 июля 2004 г. 12:44 В ответ на: в типичном случае отправлено fontp 15 июля 2004 г. 12:26

вот не поленился и нашёл

#ifndef M_PI
#define M_PI 3.14159265358979323846
#endif

static inline double arctg1 (double x)
//-------------------------------------------
// 2
// arctg(x) = A*x - B*x , if (x >= 0)
//
// 2
// arctg(x) = A*x + B*x , if (x < 0)
//
// error <= 0.004
//-------------------------------------------
{
if (x < 0) {
return x*(1.0597 + 0.27131*x);
}else {
return x*(1.0597 - 0.27131*x);
}
}/*arctan1*/

double arctg2 (double Rs, double Rc)
//-------------------------------------------
// 0 <= arctg2 <= 2*M_PI
//-------------------------------------------
{
double ret = 0.0;
if (Rc > +1.0E-10) {
if (Rs > +Rc) {
ret = (+M_PI/2.0 - arctg1(Rc/Rs));
}else if (Rs < -Rc) {
ret = (+3.0*M_PI/2.0 - arctg1(Rc/Rs));
}else { // |Rs| < |Rc|
if (Rs >= 0.0) {
ret = arctg1(Rs/Rc);
}else { // Rs < 0
ret = 2.0*M_PI + arctg1(Rs/Rc);
}
}
}else if (Rc < -1.0E-10) {
if (Rs < +Rc) { // Rs < 0, |Rs| > |Rc|
ret = (+3.0*M_PI/2.0 - arctg1(Rc/Rs));
}else if (Rs > -Rc) { // Rs > 0, |Rs| > |Rc|
ret = (+M_PI/2.0 - arctg1(Rc/Rs));
}else { // |Rs| < |Rc|
ret = (M_PI + arctg1(Rs/Rc));
}
}else {
if (Rs > 0.0) {
ret = +M_PI/2.0;
}else {
ret = +3.0*M_PI/2.0;
}
}
if (ret < 0.0) {
ret += 2*M_PI;
}else if (ret >= 2*M_PI) {
ret -= 2*M_PI;
}
return ret;
}/*arctg2*/

в коротком разложении понятно не ряд Тэйлора, а чебышевская аппроксимация

• Спасибо, арктангенс есть в dsplib 55xx, можно взять за основу. Видимо без честного деления и арктангенса не обойтись. Делаю демодуляцию DPSK с мягкими решниями, 00-strong 0 , 01- soft 0, 10 - soft 1, 11- strong 1. Мне нужно проверить разницу фаз c несколькими дискретными уровненями. — AntZ   (15.07.2004 13:25, пустое)
  • А зачем для демодулятора(DPSK) с мягким решением углы и вообще арктангенсы? — alostap   (15.07.2004 15:56, пустое)
    • Вообще для декодера достаточно метрик (мягких решений) = — alostap   (15.07.2004 16:11, 158 байт)
  • тогда лучше орктангенс не вычислять, а сделать типа векторного квантования — fontp   (15.07.2004 13:37, 225 байт)
    • Очень понравилось выражение — AntZ   (15.07.2004 15:29, 118 байт)
      • у вас есть плоскость (x, y) — fontp   (15.07.2004 15:43, 808 байт)
        • конфа не ест математической нотации x меньше a[i] умножить на y ;) — fontp   (15.07.2004 15:47, пустое)
          • Никаких проблем - меньше амперсандlt - <, больше амперсандgt - > — AntZ   (15.07.2004 16:47, пустое)
          • Съест если знак меньше написать по HTML-правилам — SM   (15.07.2004 15:58, пустое)
    • Согласен, еще использовал для расчета частоты опорника при устранении частотной ошибки. По-моему на каждый бит такое считать - слишком дорого. — andy_P   (15.07.2004 14:56, пустое)