Моделирование в РЗА (Страница 11 из 15)

Посмотрите в симулинке "Phasor Solution Method". Если использовать этот способ, модель должна выдать синхронизированные векторы и не нужно будет городить АЦП, БПФ и пр.Ну а если Continuous или Discrete тогда нужна хотя бы простенькая модель измерительного тракта (непрерывный сигнал->дискретизатор->БПФ).На первой картинке subsystem, судя по входам-выходам, выполняет преобразование Фурье. Затем разложение комплексов на амплитуду, угол,преобразование угла в градусы и объединение всех полученных значений в один вектор (массив).Погрешность прибора, шум и синхронизацию в первом приближении, на мой взгляд, можно не делать, а внести в модель позже.Погрешность можно смоделировать ограничением точности расчетов - заданием типа сигнала, например int16 или fixdt(1,16).Шум - добавлением в полезный сигнал случайного шума.Несинхронность можно смоделировать блоками задержки (transport delay). Тогда надо будет еще моделировать схемы синхронизации или компенсации несинхронности.как-то так, если я правильно понял вопрос.

• E.A.BUCHINSKIY
• Пользователь
• Неактивен

• Откуда: Иркутск
• Зарегистрирован: 2011-05-25
• Сообщений: 358
• Репутация : [ 0 | 0 ]

Коллеги,в данный момент изучаю основы цифровой обработки сигналов, в т.ч. и средствами Simulink. Некоторые моменты прояснились, а кое-что дается тяжело.Буду благодарен если кто-то поделится моделью-примером в симулинк оконного преобразования Фурье синусоидального сигнала. Функция окна не имеет значения - любая. Пример необходим просто для понимания. Выполнил его сам блоком Short-Time FFT, но пока не могу понять детально что к чему. Заранее благодарен за помощь.

• Диоген
• Пользователь
• Неактивен

• Зарегистрирован: 2012-01-16
• Сообщений: 185
• Репутация : [ 0 | 0 ]

E.A.BUCHINSKIY, именно в Симулинке?

• E.A.BUCHINSKIY
• Пользователь
• Неактивен

• Откуда: Иркутск
• Зарегистрирован: 2011-05-25
• Сообщений: 358
• Репутация : [ 0 | 0 ]

можно и в самом Matlab - не так важно.

• retriever
• Пользователь
• Неактивен

• Зарегистрирован: 2012-11-26
• Сообщений: 2,622
• Репутация : [ 12 | 0 ]

Небольшой вопрос из любопытства: а для чего именно вы хотите использовать оконное преобразование Фурье? Его результатом будет спектр. а зачем вам спектр? Я так понял по вашим постам в этой теме, что вы занимаетесь синхронизированными векторными измерениями. Для этого достаточно выделить частоту 50 Гц и все. (Или вы хотите сделать особую, уличную магию?)Я это к тому, что выделить частоту 50 Гц - это гораздо проще, чем получить спектр.

• E.A.BUCHINSKIY
• Пользователь
• Неактивен

• Откуда: Иркутск
• Зарегистрирован: 2011-05-25
• Сообщений: 358
• Репутация : [ 0 | 0 ]

Небольшой вопрос из любопытства: а для чего именно вы хотите использовать оконное преобразование Фурье?

также из любопытства. изучаю программу, ее возможности. нашел прекрасную статью, где все проиллюстрировано и все понятно. Решил попробовать получить Рисунок 1 (нижний график). Пока не особо получилось. Поэтому интересно было бы посмотреть как это выполняется - особенно в симулинк. Как выбираются координаты для графика.

(Или вы хотите сделать особую, уличную магию?)

программу только начал осваивать поэтому получение магии не исключено

Я это к тому, что выделить частоту 50 Гц - это гораздо проще, чем получить спектр.

был бы признателен если бы Вы проиллюстрировали на синусоидальном сигнале.

• retriever
• Пользователь
• Неактивен

• Зарегистрирован: 2012-11-26
• Сообщений: 2,622
• Репутация : [ 12 | 0 ]

нашел прекрасную статью, где все проиллюстрировано и все понятно

ммм. какую именно?

был бы признателен если бы Вы проиллюстрировали на синусоидальном сигнале.

Используйте блок Fourier (на вход - Sine Wave, на выходах - два блока Display). Делает ДПФ, при этом для первого периода он замещает недостающие выборки какими-то числами (вроде нулями). Поэтому в течение первого периода -какая-то ерунда, далее нормально.

Можно также использовать Discrete Fourier, я так понимаю, что он вместо "ерунды" подставляет цифры, прописанные в поле Initial input. Поэтому у него такой ступенчатый график на выходе..

В матлабе для получения ДПФ можно воспользоваться формулами:N=10;w=2*pi*50;t=0:0.02/N:(0.02-0.02/N);%Это массив времени исходного сигналаu=sin(w*t);%Это исходный сигнал

T=t;%Это "окно". Выделил его здесь потому, что для получения изменяющегося во времени значения его нужно будет двигать. Число выборок в окне всегда равно N. Для первого периода можно либо не делать вообще, либо довтыкать нули вместо недостающих выборок. U=(2i/N)*sum(u.*exp(-1i*w*T))Итого будет комплекс. Вроде бы можно оптимизировать скрипт (ибо exp(-1i*w*t) повторяются, если периодов сигнала несколько).