ограничение амплитуды, для квадр. с-ла Опции

[shf_05]

необходимо ограничить по амплитуде синус(затем и произвольный сигнал) (например по уровню X), далее он фильтруется и проч., в итоге получаем его огибающую.
- это предельно прозрачно, если сигнал один, скажем sin(1000Гц) или случ. сигнал, как быть для квадратурного с-ла?
можно для sin ограничить отдельно Re и Im, вроде в этом случае все будет корректно, будет ли все корректно для произвольного сигнала?
можно конечно повыводить формулы, помоделить, надеюсь кто-нить скажет "с ходу", поделитесь опытом.

[petrov]

необходимо ограничить по амплитуде синус(затем и произвольный сигнал) (например по уровню X), далее он фильтруется и проч., в итоге получаем его огибающую.
- это предельно прозрачно, если сигнал один, скажем sin(1000Гц) или случ. сигнал, как быть для квадратурного с-ла?
можно для sin ограничить отдельно Re и Im, вроде в этом случае все будет корректно, будет ли все корректно для произвольного сигнала?
можно конечно повыводить формулы, помоделить, надеюсь кто-нить скажет "с ходу", поделитесь опытом.

Нет не корректно ограничивать квадратуры отдельно, надо переводить в полярные координаты, заменять амплитуду на константу, затем обратно в декартовы.

[shf_05]

т.е. сделать sat(Ampl)*(cos(phi)+i*sin(phi))?
можн как нибудь попроще- ведь на дсп вычислять tan, atan, sin, cos оччень нехочется.

[petrov]

т.е. сделать sat(Ampl)*(cos(phi)+i*sin(phi))?
можн как нибудь попроще- ведь на дсп вычислять tan, atan, sin, cos оччень нехочется.

Не вижу никаких сложностей в вычислении аргумента комплексного числа, кстати для каких вам целей, может и нет нужды обратно в квадратуры переходить?

[shf_05]

соложность в том, что вычислять аргумент надо в режиме РВ, на медленном проце.
далее квадратуры фильтруются, и вычисляется их огибающая.

[alex_os]

соложность в том, что вычислять аргумент надо в режиме РВ, на медленном проце.
далее квадратуры фильтруются, и вычисляется их огибающая.

Если фазу считать не хочется можно так :

Код

s = x+1j*y; // входной сигнал
tmp = x^2 + y^2;
if( tmp>A^2 ) // A - порог ограничения
{
    y = s * A/sqrt(tmp)  // Ограниченный выходной сигнал
}
else
{
   y = s;
}

[shf_05]

+1 )))
это эквивалентно ограничению амплитуды и умножению на фазу

зы: раньше я и делить не хотел теперь все-таки видимо придется