| Наблюдение эффектов размножения и наложения спектров.
 
 
 
 Исходные данные: Х1=2, Х2=0.8, Х3=0.8, ϕ1=0, ϕ2= π, ϕ3=0, F1=1МГц, Fд=8МГц, ∆t=1/512 мкс. Требуется: 1.Определить амплитудный спектр аналогового сигнала где  2.Определить амплитудный спектр дискретного сигнала в интервале частот от 0 до 2FД+3F1, где FД – частота дискретизации. 3.Выполнить моделирование аналогового сигнала, его дискретизации и восстановления аналогового сигнала из дискретного при двух значениях частоты дискретизации: первое значение FД= 8 МГц, второе значение в два раза меньше первого FД= 4 МГц. Моделирование выполняется по программе «Diskret». 4.Сравнить аналоговый сигнал, восстановленный из дискретного, с исходным аналоговым сигналом на входе дискретизатора. Сравнение выполнить при двух значениях частоты дискретизации.   Расчет частот спектральных составляющих дискретного сигнала по формуле: 
 k = 0 1. F=F1, f1 = F1 = 1 МГц, 2. F=2F1, f2 = 2F1 = 2 МГц, 3. F=3F1, f3 = 3F1 = 3 МГц.   k=1 4. F=3F1, f4 = FД – 3F1 = 8 – 3 = 5 МГц, 5. F=2F1, f5 = FД – 2F1 = 8 – 2 = 6 МГц, 6. F=F1, f6 = FД – F1 = 8 – 1 = 7 МГц, 7. F=F1, f7 = FД + F1 = 8 + 1 = 9 МГц, 8. F=2F1, f8 = FД + 2F1 = 8 + 2 = 10 МГц, 9. F=3F1, f9 = FД + 3F1 = 8 + 3 = 11 МГц.   k=2 10. F=3F1, f10 = 2FД – 3F1 = 16 – 3 = 13 МГц, 11. F=2F1, f11 = 2FД – 2F1 = 16 – 2 = 14 МГц, 12. F=F1, f12 = 2FД – F1 = 16 – 1 = 15 МГц, 13. F=F1, f13 = 2FД + F1 = 16 + 1 = 17 МГц, 14. F=2F1, f14 = 2FД + 2F1 = 16 + 2 = 18 МГц, 15. F=3F1, f15 = 2FД + 3F1 = 16 + 3 = 19 МГц.   Построение спектральных диаграмм аналогового и дискретного сигналов.   
 Рисунок 1 – Амплитудный спектр аналогового сигнала 
 Рисунок 2 – Амплитудный спектр дискретного сигнала Результат моделирования по программе «Diskret A» в виде временных и спектральных диаграмм аналогового и дискретного сигнала приведён на рисунке 3.   
 Рисунок 3. Временные и спектральные диаграммы при FД > 2F3 Вывод: При увеличении дискретизации возникает эффект размножения спектра и восстановленный аналоговый сигнал совпадает по форме с исходным, так как соблюдается теорема Котельникова.     Повторим эксперимент при уменьшенной в два раза частоты дискретизации. 
 Рисунок 4. Временные и спектральные диаграммы при FД < 2F3 Вывод: При уменьшении дискретизации в 2 раза, восстановленный и исходный сигнал не совпадает по форме, так как не соблюдается теорема Котельникова и происходит эффект наложения. Задание №1б.Моделирование процесса дискретизации аналогового сигнала, модулированного по амплитуде, и восстановления аналогового сигнала из дискретного.   Исходные данные: X0=2, f0=13, F=0.20, FД=4, Δt=1/256 Задание 1б выполняется по программе «Diskret_B» (Приложение Б). Эта программа отличается от программы «Diskret_A» только видом функции, описывающей исходный аналоговый сигнал. Результатом выполнения программы являются временные и спектральные диаграммы аналогового и дискретного сигналов (рисунок 5). 
 Рисунок 5 – Временные и спектральные диаграммы аналогового и дискретного сигналов при частоте несущей f0=13 МГц, частоте модуляции 0.20 МГц и частоте дискретизации FД=4МГц Вывод: восстановленный аналоговый сигнал идентичен сигналу на входе дискретизатора. Задание №2.Определение спектра восьмиточечной последовательности отсчетов сигнала с использованием алгоритмов БПФ с прореживанием во времени и с прореживанием по частоте и выполнение обратного преобразования.   2.1. Определение спектра последовательности xn, где n = 0,1,2..7, методом прямого дискретного преобразования Фурье 
 где  , N=8, k=0,1,2..7. 2.2. Определение последовательности отсчетов сигнала xn, где n = 0,1,2..7, по известной последовательности отсчетов спектра Sk, где k = 0,1,2..7, методом обратного дискретного преобразова- ния Фурье 
 где n = 0,1,2..7.   2.3. Определение спектра последовательности xn, где n = 0,1,2..7, методом БПФ с прореживанием во времени. 2.4. Определение последовательности отсчетов сигнала xn, где n=0,1,2..7, по известной последовательности отсчетов спектра Sk, где k=0,1,2..7, методом БПФ с прореживанием во времени. 2.5. Определение спектра последовательности xn, где n = 0,1,2..7, методом БПФ с прореживанием по частоте. 2.6. Определение последовательности отсчетов сигнала xn, где n=0,1,2..7, по известной последовательности отсчетов спектра Sk, где k=0,1,2..7, методом БПФ с прореживанием по частоте.                                           | Номер варианта | Последовательность x0 x1 x2 x3 x4 x5 x6 x7 |                       |  | 3.000 2.414 1.000 -0.414 -1.000 -0.414 1.000 2.414 |    Определим спектр заданной последовательности по программе «Прямое ДПФ»:  ,
 Спектр амплитуд: 
 Спектр фаз: 
 Спектры амплитуд и фаз приведены на рисунках 2.1 и 2.2 соответственно.   
 Рисунок 2.1 – Спектр амплитуд.     
 Рисунок 2.2 – Спектр фаз.   2. Определим последовательность отсчетов дискретного сигнала по её спектру.  .
 Воспользовавшись программой «Обратное ДПФ», получим: 
 3. Воспользуемся алгоритмом прямого БПФ с прореживанием во времени для исходной последовательности. 
 
 Рисунок 2.3   Выполним перестановку членов этой последовательности с учетом двоичной инверсии индексов (рисунок 2.3). В результате получим:   
 Воспользовавшись алгоритмом прямого БПФ с прореживанием во времени (рисунок 2.4), выполним операции «бабочка»: на первом уровне – 4 двухточечных «бабочки», на втором уровне – 2 четырёхточечных, на третьем уровне – одна восьмиточечная. 
   Рисунок 2.4 – Алгоритм прямого БПФ с прореживанием во времени.   4. Выполним перестановку членов полученной последовательности отсчетов спектра с учетом двоичной инверсии индексов. В результате получим: 
 Подадим эту последовательность на входы двухточечных бабочек рисунка 2.5 и выполним вычисления.     
 Примечание: результаты выполнения бабочки делятся на N=8   Рисунок 2.5 – Алгоритм обратного БПФ с прореживанием во времени.   5. Воспользуемся алгоритмом БПФ с прореживанием по частоте для определения прямого преобразования последовательности.   
   Согласно алгоритму прямого БПФ с прореживанием по частоте (рисунок 2.6), выполним операции «бабочка»: на первом уровне – одна восьмиточечная «бабочка», на втором уровне – 2 четырёхточечных, на третьем уровне – четыре двухточечных.   
 Рисунок 2.6 – Алгоритм прямого БПФ с прореживанием по частоте.   Выполним перестановку членов полученной последовательности с учетом двоичной инверсии индексов (рисунок 2.3). В результате получим: 
 6. По алгоритму рисунка 2.7 выполним обратное преобразование последовательности отсчетов спектра.  .
 
 Примечание: результаты на выходах двухточечных бабочек делятся на N=8 Рисунок 2.7 – Алгоритм обратного БПФ с прореживанием по частоте. Результатом этой операции является последовательность отсчетов дискретного сигнала .   
   Выполнив перестановку членов этой последовательности с учетом двоичной инверсии индексов, получим.   
   
 
 
 
 |