| СПЕКТРАЛЬНЫЙ АНАЛИЗ В ПАКЕТЕ ПРОГРАММ MATHCAD | |
| Автор: student | Категория: Технические науки / Автоматизация | Просмотров: 6056 | Комментирии: 0 | 29-12-2013 00:24 |
СПЕКТРАЛЬНЫЙ АНАЛИЗ
В ПАКЕТЕ ПРОГРАММ MATHCAD
Лабораторная работа №1
Спектральный анализ в пакете программ Mathcad
1. ЦЕЛЬ РАБОТЫ
1.1. изучение функций Mathcad для анализа спектра сигналов;
1.2. изучение функций Mathcad для работы с файлами данных;
1.3. исследование спектра сигналов.
2. КРАТКИЕ ПОЯСНЕНИЯ К ЛАБОРАТОРНОЙ РАБОТЕ
2.1. Спектры периодических сигналов
Периодические (полигармонические) сигналы
X(t)=X(t+kT), (1)
где k – целое число, Т – период повторения сигнала.
Для полигармонических сигналов вводится понятие спектра, т.е. совокупности простых гармонических составляющих, на которые можно разложить сигнал. Спектр сиг-нала выражает его частотный (спектральный) состав, т.е. распределение по частоте ам-плитуды и фазы гармоник.
Полигармонический сигнал может быть представлен дискретным рядом Фурье (дискретным спектром):
X(t)= , (2)
где , , , .
X(t)= , (3)
где – среднее значение, , .
Информативные параметры сигналов.
Xmax – максимальное значение;
– среднее значение (постоянная составляющая);
– средневыпрямленное значение;
– действующее значение (СКЗ);
= .
, – коэффициенты амплитуды и формы;
Прямоугольный импульс.
, bk=0, – среднее, – СКЗ, .
а) б)
Рис. 1. Прямоугольный импульс (а) и его дискретный спектр (б)
Треугольный импульс (симметричный).
, bk=0, – среднее, – СКЗ.
а) б)
Рис. 2. Треугольный импульс (симметричный) (а) и его спектр (б)
Треугольный импульс (пилообразный).
, аk=0, – среднее, – СКЗ.
а) б)
Рис. 3. Треугольный импульс (пилообразный) (а) и его спектр (б)
Полусинусоида.
, bk=0, – среднее, – СКЗ.
а) б)
Рис. 4. Полусинусоида (а) и ее спектр (б)
2.2. Моделирование спектра периодических сигналов в программе Mathcad
2.2.1. При моделировании спектра периодического сигнала в пакете программ Mathcad достаточно знать выражения для коэффициентов ряда Фурье и задать требуемое количество анализируемых гармоник.
Например, для прямоугольного сигнала требуется проанализировать первые 10 гармоник.
а) б)
Рис. 5. Прямоугольный периодический сигнал (а) и его спектр (б)
2.2.2. Более сложный способ вычисления спектра, но вместе с тем более универ-сальный, использует встроенные функции интегрирования Mathcad и позволяет вычислять коэффициенты ряда Фурье при произвольной форме периодических сигналов.
Например, для прямоугольного сигнала требуется проанализировать первые 10 гармоник.
Сигнал может быть задан симметрично (рис. 5.а) на интервале времени от –(T/2) до +(T/2). В этом случае изменятся пределы интегрирования при вычислении коэффициентов ряда, выражения 2–3.
Результат вычисления амплитудного спектра будет такой же, как на рис. 5.б.
2.3. Спектры непериодических сигналов
Отличие переходных (непериодических) сигналов от периодических в том, что их нельзя представить суммой гармоник, то есть в виде дискретного спектра. Однако такие сигналы можно представить в виде непрерывного спектра, получаемого преобразовани-ем Фурье
. (4)
Учитывая равенства
и
спектр сигнала вычисляется также в виде
. (5)
Экспоненциальный импульс.
(6)
а) б)
Рис. 6. Экспоненциальный импульс (а) и его спектр (б)
Затухающие колебания.
(7)
а) б)
Рис. 7. Затухающие колебания (а) и их спектр (б)
Прямоугольный импульс.
(8)
а) б)
Рис. 8. Прямоугольный импульс (а) и его спектр (б)
2.4. Моделирование спектра непериодических сигналов в программе Mathcad
Вычисление спектра непериодических сигналов рекомендуется производить в со-ответствии с выражениями (5). В этом случае удается повысить скорость вычислений по сравнению с прямым вычислением интегрального преобразования Фурье (4).
Далее приведен пример вычисления спектра прямоугольного импульса.
а) б)
Рис. 9. Прямоугольный импульс (а) и его спектр (б)
2.5. Анализ спектра с помощью Быстрого преобразования Фурье (БПФ)
Для анализа спектра сигналов, записанных в виде файлов данных, в пакете про-грамм Mathcad2001 имеются следующие функции:
READPRN("file”) – считывание массива данных из текстового файла;
WRITEPRN("file”) – запись данных из матрицы в текстовый файл;
READWAV(file) – считывание массива амплитуд из WAV-файла (колонки соот-ветствуют каналам данных, строки – моментам времени).
WRITEWAV(file,s,b) – запись матрицы в формате WAV-файла;
CFFT(A) – вычисление БПФ матрицы A (следует не путать с функцией cfft(A));
ICFFT(B) – вычисление обратного БПФ матрицы B.
Функции вычисления БПФ требуют, чтобы количество обрабатываемых отсчетов было кратно 2V, где V – целое число.
2.5.1. Например, требуется записать текстовый файл данных сигнала, состоящего из основной гармоники и некоторого равномерно распределенного шума.
Для записи не менее двух периодов основной гармоники fx при условии, что коли-чество отсчетов Nmax кратно 2V, Nmax должно определяться следующим образом:
В данном примере используются кратные основной частоте сигнала гармоники шума.
2.5.2. Далее производится запись полученных отсчетов сигнала в файл «DataX.prn»
2.5.3. Если требуется произвести считывание данных, например, записанных с по-мощью выше указанных процедур, выполняется обратная операция
а) б)
Рис. 10. Сигнал, записанный в файл данных «DataX.prn» (а) и его спектр (б)
2.5.4. Поскольку количество отсчетов в файле кратно 2V, то нет необходимости вы-полнять проверку и усечение количества обрабатываемых данных.
Далее выполняется БПФ и строится график амплитудного спектра сигнала.
2.5.5. По номеру можно определить соответствующую частоту гармоники в спек-тре.
В данном случае гармонике с номером 3 соответствует частота 96 Гц.
3. ПРОГРАММНОЕ ОБЕСПЕЧЕНИЕ
При выполнении лабораторной работы используется пакет программ Mathcad вер-сии 2000 и выше.
4. ПРОГРАММА ЛАБОРАТОРНОЙ РАБОТЫ
4.1. Изучить методы анализа спектра сигналов в пакете программ Mathcad.
4.2. Выполнить запись и считывание файла, содержащего отсчеты сигнала.
4.3. Рассчитать и проанализировать спектры сигналов с помощью быстрого преоб-разования Фурье (БПФ).
4.4. Рассчитать и проанализировать спектры непериодических сигналов с помощью преобразования Фурье.
4.5. Рассчитать и проанализировать спектры периодических сигналов с помощью дискретного ряда Фурье.
5. КОНТРОЛЬНЫЕ ВОПРОСЫ
5.1. Определение спектра сигнала.
5.2. Методы расчета спектра периодического сигнала.
5.3. Методы расчета спектра непериодического сигнала.
5.4. Определение среднего, среднеквадратического, средневыпрямленного значе-ний, коэффициентов амплитуды и формы сигнала.
5.5. Основные функции Mathcad для записи/считывания файлов данных и анализа спектра.
5.6. Как определить ближайшее меньшее число кратное 2V, где V – целое число.
5.7. Как задать в Mathcad прямоугольный импульс (треугольный симметричный).
5.8. Какие функции используются в Mathcad для записи/считывания текстовых файлов данных.
5.9. Как определить количество отсчетов в текстовом файле данных.
5.10. Какие функции Mathcad используются для вычисления БПФ, обратного БПФ.
6. ПОРЯДОК ВЫПОЛНЕНИЯ ЛАБОРАТОРНОГО ЗАДАНИЯ
6.1. Спектральный анализ с использованием БПФ
6.1.1. Подготовить и ввести в файл программы «lab2_1.mcd» исходные данные (п.2.5.1), необходимые для синтеза сигнала X: частота сигнала Fс= 100 Гц; ам-плитуда – 10; коэффициент гармоник KГ=20 % (L=10, первые 10 гармоник с рав-номерным распределением), частота дискретизации Fд=10 кГц. Количество от-счетов сигнала должно быть 2B , где B – целое число (рекомендуется записать не менее 2-х периодов сигнала).
6.1.2. Записать отсчеты полученного сигнала в файл «dataX.prn». Построить график сигнала (п. 2.5.2).
6.1.3. Создать второй программный файл «lab2_2.mcd», в котором произвести считы-вание данных (п. 2.5.3) из файла «dataX.prn» в матрицу Y.
6.1.4. Построить график сигнала Y (п. 2.5.3). Определить количество отсчетов сигнала в файле (считая, что количество заранее неизвестно).
6.1.5. Вычислить спектр сигнала методом БПФ (функция C=СFFT(Y)). Построить график амплитудного спектра сигнала (п. 2.5.4).
6.1.6. Проверить соответствие между номером гармоники в полученном спектре и действительной частотой гармоники в сигнале (п. 2.5.5). Для этого необходимо найти номер K максимальной гармоники в спектре – данная гармоника соответ-ствует основной частоте сигнала Fc=100 Гц. По формуле Fk=KFд / N вычисляет-ся частота гармоники (N – количество отсчетов анализируемого сигнала). В об-щем случае имеется отклонение между действительной частотой Fc и частотой Fk.
6.1.7. Повторить пп. 6.1.1 – 6.1.6 для гармонического сигнала с частотой Fc=27=128 Гц. (при частоте дискретизации Fд=10 кГц). Построить графики исходного сиг-нала X и его спектра C, а также проверить соответствие между номером гармо-ники в полученном спектре и действительной частотой гармоники в сигнале.
6.1.8. Повторить пп. 6.1.1 – 6.1.6 для гармонического сигнала с частотой Fc=100 Гц. (частота дискретизации при этом равна 2V, где V – целое число, например V=13). Построить графики исходного сигнала X и его спектра C, а также проверить соответствие между номером гармоники в полученном спектре и действительной частотой гармоники в сигнале.
6.1.9. Повторить пп. 6.1.1 – 6.1.6 для гармонического сигнала с частотой Fc=27=128 Гц. (частота дискретизации при этом равна 2V, где V – целое число). Построить графики исходного сигнала X и его спектра C, а также проверить соответствие между номером гармоники в полученном спектре и действительной частотой гармоники в сигнале.
6.2. Анализ спектра с использованием преобразования Фурье и ряда Фурье
6.2.1. Создать программный файл «lab2_3.mcd». Задать непериодический сигнал в ви-де функции времени (например, прямоугольный сигнал с амплитудой Ux=10 и длительностью =0.1 с), п. 2.4.
6.2.2. Вычислить спектр сигнала, используя преобразование Фурье (мнимую и дейст-вительную части рекомендуется вычислять раздельно). Построить график полу-ченного спектра сигнала (п. 2.4).
6.2.3. Определить коэффициенты ряда Фурье для выбранного сигнала в случае, когда сигнал является периодическим (период задать в два раза больше длительности сигнала). Построить график полученного спектра сигнала (представить его в ви-де дискретных гармоник), п. 2.2.2.
6.2.4. Сравнить графики полученных спектров.
7. СПИСОК ЛИТЕРАТУРЫ
1. Орнатский П.П. Теоретические основы информационно-измерительной техни-ки. – Киев: Вища школа, 1983. – 455 с.
2. Каганов В.И. Радиотехника+компьютер+Mathcad. – М.: Горячая линия - Теле-ком, 2001. – 416 с.
3. Иванов В.А. и др. Математические основы теории автоматического регулиро-вания. Учеб. пособие для вузов. / Под ред. Чемоданова Б.К. – М.: Высшая шко-ла, 1971. – 808 с.
4. Сергиенко А.Б. Цифровая обработка сигналов. – СПб.: Питер, 2003. – 604 с.