Привет.
Вопрос - как сгенерировать сигнал 50 Гц с гармониками со 2-ой по 100-ую (т.е. до 5кГц) в Симулинке без применения сотни генераторов синусов. С различными амплитудами и фазами. Частота дискретизации(sample time) 20мкс.

Интересует нетрадиционное решение. Можно из области математики. Я помню в DSP генератор синуса можно было реализовать очень прикольненько и простенько на Си.

Сгенерировать Фурье-спектр с псевдо-дельтафункциями на нужных частотах и сделать обратное БПФ. Я так студентам "синтетические спектры" для обработки делал. Потом еще awgn добавить, и будет нормальный псевдосигнал.

Дело в том, что "сделать обратное БПФ" я только генераторами синусов и умею.

пропустить через нелинейность?

Многоканальный генератор на матлабе - это три строки


Вырианты, если три строки - это слишком мало и/или не достаточно наукоемко:

> DDS. 1 таблица с периодом синусоиды с нужной частотой дискретизации. Обновляете 100 фазовых аккумуляторов с кратным инкрементом. Частота дискретизации периода в таблице определяет SFDR.

> генератор - нелинейность - фильтр. Будет сложновато построить фильтр, обеспечивающий заданный набор амплитуд и фаз на выходе.

> ОБПФ, как вам писали, но там видится одна проблема - сшивка окон, если возникнет требование изменять частоту основной гармоники. Придется выкручиваться с интерполяцией на выходе обпф. Или напрямую менять частоту отсчетов цап.

Можно подробнее, как это работает?

Мама дорогая! Пишем:

Симулинк с векторами работает также как и с одномерными сигналами, ваш один простенький DSP генератор также может генерировать 100 синусов, т. е. не требуется рисовать 100 генераторов, они автоматически размножаются под размерность сигнала.

Тогда это должна быть нелинейность, которая описывается полиномом 100-го порядка)))

LPF{ ( 1 + cos(wt) ) * sign(cos( wt )) }

естественно Embedded Block и немного кода


cos(x)^2 = 1/2*(1+cos(2x));

Какова же будет точность представления амплитуды 100-ой гармоники, которая пропорциональна ~(100!)^-1?

Различные фазы так не получить..

Ровно такая, на какую у вас хватит фантазии
https://gmplib.org/
http://www.mpfr.org/

Понятно, за границы double пока не выходил, небыло необходимости