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

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

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

Используйте блок 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) повторяются, если периодов сигнала несколько).

E.A.BUCHINSKIY, Симулинк считаю вещь канешн замечательна, но когда примерно представляешь, что внутри блока и как он работает...

Если нам это оперное творчество шоб получить амплитуду и фазу сигнала, то вот можно как-то так в скрипте:

пример:

t=0:0.001:0.019;%период - для ДПФ ширина окна нужна минимум в период, ну или кратно
f=50;%что выделяем
F=4*sin(2*pi*f+pi/3)+1*sin(2*2*pi*f-pi/6);%исх. сигнал - добавил 2ую гармонику для демонстр
Res=fft(F);%могем сразу воспользоваться встроенной функцией, на выходе получаем массив компл чисел для гармоник, кртаных основной (50 Гц). 1ый отсчет - нул гармоника, поэтому берем второй оттсчет...
A=2*abs(Res(2))*f/(t(2)-t(1));
Phase=angle(Res(2))+pi/2

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

for k=1:end-N
F=Y(k:(k+19))
.....
end

Добавлено: 2013-11-06 22:52:21

Для интереса можно добавить к сигналу шум или экспоненту и поисследовать что будет :)

вбейте
help fft
это просто дискр. преобразование Фурье. Разница с БПФ есть.

У меня выдает ампл=4, фаза=3.14 для
этого примера

Да, канешн, опечатка...

В строке расчета амплитуды поменяйте / на х

этот результат мне больше нравится =D

наоборот

спасибо с этим разобрался.

Может быть кто-то подскажет как в симулинк(именно в симулинк) реализовать вот это:

http://rzia.ru/extensions/hcs_image_uploader/uploads/50000/8000/58324/thumb/p188oubssrifu1dhshau1k461ns31.JPG

X - входной сигнал
g - оконная функция

пока голову ломаю:

- оконную функцию легко реализовать блоком window
- входной сигнал тоже не проблема получить sine wave
- преобразовать в цифровой вид тоже понятно как

с помощью каких блоков можно получить функцию комплексной экспоненты и выполнить данное интегрирование?

для понимания буду признателен если кто-то сможет выслать пример использования S-функции - простейшей.

http://matlab.exponenta.ru/simulink/book1/16.php
Я учил S-функции по этому хелпу. Для понимания основ лучше всего скачать пример с функцией timestwo (вход умножается на 2 и подается на выход). Этот пример можно использовать как "болванку" для всего остального.

Основной принцип работы с ними: берется эта болванка, и меняется содержимое последней функции на свое:

Если надо, можно будет потом поменять кое-какие опции в файле (там написано). Ничего сложного.

ВСпомнил - для интегрирования можно воспользоваться блоком 1/s