Простой анализ спектра простого сигнала Опции

[Punk]

Здравствуйте мегаматематики. Попросили(добровольно-принудительно) меня в технаре написать програмку на камп чтобы строила 2 графика гармонических функций а ниже нарисовала спектр сигнала, являющегося суммой этих двух гармоник.
Синус с косинусом нарисовать и результирующий сигнал найти непроблема, а вот спектр мозгов не хватает, ПОМОЖИТЕ ЛЮДИ ДОБРЫЕ !!! Препод понимает C, Pascal, qBasic объектноориентированное непонимает.Заранее спасибо за любую инфу

[Punk]

я всеже хочу найти простой метод без библиотек в котором смогу самостоятельно разобраться.

методом перерывания инфы вот че нашел из "ИСКУССТВО СХЕМОТЕХНИКИ" в главе посвященной цифровой фильтрации,..код на фортране

A=EXP(-TS/TC)
B=1.-A
DO 10I=2,N
10X(I)=A*X(I-1)+B*X(I)

где Х-массив данных TS-интервал времени между отсчетами TC-выбранная постоянная времени фильтра

я сравнил эту формулу с формулами из литературы посвященной БПФ


и пришел к выводу что спектрограмму можно построить при помощи неё, я прав ?

[SM]

я всеже хочу найти простой метод без библиотек в котором смогу самостоятельно разобраться.

методом перерывания инфы вот че нашел из "ИСКУССТВО СХЕМОТЕХНИКИ" в главе посвященной цифровой фильтрации,..код на фортране

A=EXP(-TS/TC)
B=1.-A
DO 10I=2,N
10X(I)=A*X(I-1)+B*X(I)

где Х-массив данных TS-интервал времени между отсчетами TC-выбранная постоянная времени фильтра

я сравнил эту формулу с формулами из литературы посвященной БПФ


и пришел к выводу что спектрограмму можно построить при помощи неё, я прав ?

Нет, это цифровой фильтр НЧ первого порядка. К анализу спектра отношения не имеет. А ДПФ для точки номер "k" рассчитывается так:

Код

a=b=0.0;
for (i=0; i<N; i++) {
  a += fcos(2*pi*i*k/N) * signal[i];
  b += fsin(2*pi*i*k/N) * signal[i];
}

[Punk]

Нет, это цифровой фильтр НЧ первого порядка. К анализу спектра отношения не имеет. А ДПФ для точки номер "k" рассчитывается так:

Код

a=b=0.0;
for (i=0; i<N; i++) {
  a += fcos(2*pi*i*k/N) * signal[i];
  b += fsin(2*pi*i*k/N) * signal[i];
}

Из этого цикла выходят два коэффициента a и b че сопсна с имя делать, и как определить число k, равно ли оно N ? если че простите за тупизну и настойчивость.

[SM]

Из этого цикла выходят два коэффициента a и b че сопсна с имя делать, и как определить число k, равно ли оно N ? если че простите за тупизну и настойчивость.

Тяжелый случай... k это номер выборки спектральной ф-ции в частотной области. Вычисляется для всех k от 0 до N-1. Коэфф. a и b дают на выходе X(k) = a+j*b где j сами попробуйте догадаться что.

[Punk]

Тяжелый случай... k это номер выборки спектральной ф-ции в частотной области. Вычисляется для всех k от 0 до N-1. Коэфф. a и b дают на выходе X(k) = a+j*b где j сами попробуйте догадаться что.

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

В общем я разобрался всем спасибо.

P.S. а j это наверно мнимая еденица ,...может быть