Page 51 - 363_

52

також розміром вихідного вектора y). Наведемо приклад:
>> t =0 : 0.01: 0.2;
>> x = sin(2*pi*5*t) + cos(2*pi*12*t)
x =
Columns 1 through 7
1.0000 1.0380 0.6506 0.1716 -0.0411 0.1910 0.7637
Columns 8 through 14
1.3448 1.5564 1.1853 0.3090 -0.7348 -1.5176 -1.7388
Columns 15 through 21
-1.3768 -0.6910 -0.0747 0.1596 -0.0520 -0.4964 -0.8090
>> у = fft(x, 21)
y =
Columns 1 through 4
0.8378 2.3267-8.7884i 0.3399+6.4658i 1.2252-7.0142i
Columns 5 through 8
0.9466-2.3043i 0.8985-1.3424i 0.8802-0.8922i 0.8712-0.6097i
Columns 9 through 12
0.8664-0.4010i 0.8638-0.2289i 0.8626-0.0746i 0.8626+0.0746i
Columns 13 through 16
0.8638+0.2289i 0.8664+0.4010i 0.8712+0.6097i 0.8802+0.8922i
Columns 17 through 20
0.8985+1.3424i 0.9466+2.3043i 1.2252+7.0142i 0.3399-6.4658i
Column 21
2.3267+8.7884i
>> z = ifft(y, 21)
z =
Columns 1 through 4
1.0000+0.0000i 1.0380+0.0000i 0.6506+0.0000i 0.1716+0.0000i
Columns 5 through 8
-0.0411+0.0000i 0.1910+0.0000i 0.7637+0.0000i 1.3448-0.0000i
Columns 9 through 12
1.5564-0.0000i 1.1853-0.0000i 0.3090-0.0000i -0.7348-0.0000i
Columns 13 through 16
-1.5176 - 0.0000i -1.7388 - 0.0000i -1.3768+0.0000i -0.6910+0.0000i
Columns 17 through 20
-0.0747+0.0000i 0.1596+0.0000i -0.0520+0.0000i -0.4964+0.0000i
Column 21
-0.8090 + 0.0000i
Розглядаючи уважніше формулу дискретного перетворення Фур'є можна
помітити:

а) номер m відповідає моменту часу t т, в який виміряний вхідний

сигнал х(т);
б) номер k - це індекс значення частоти f k, якому відповідає

знайдений елемент y(k) дискретного перетворення Фур'є;

в) щоб перейти від індексів до часової та частотної областей
потрібно знати значення h дискрету (кроку) за часом, через який

вимірюється вхідний сигнал x(t) і повний діапазон T часу, протягом