Как найти спектр сигнала в mathcad

Корпорация "Центр"

Исследование входного сигнала , страница 7

С помощью программы MathCAD построим график амплитудного спектра выходного сигнала и представим его на рисунке 18:

Рис. 18 График амплитудного спектра выходного сигнала

Определим графически ширину спектра. Для этого проведём горизонтальный

уровень 0,1. Ширина спектра приблизительно равна:

С помощью программы MathCAD построим график фазового спектра выходного сигнала и представим его на рисунке 19:

Рис.19 График фазового спектра выходного сигнала.

5.Выходной сигнал.

Для построения выходного сигнала необходимо с помощью обратного преобразования Лапласа найти оригинал выражения:

Оригинал дроби представляет собой сумму трёх вычетов из пункта 3, тогда

Оригинал выражения можно найти используя функции включения:

Оригинал найдём используя свойство временного сдвига:

Подставим в формулу (16) наши значения, получим следующее выражения:

С помощью программы MathCAD построим график выходного сигнала и представим его на рисунке 20:

Рис. 20 График выходного сигнала.

6.Проверка с помощью программы ElectronicsWorkbench.

Амплитудно-частотная и фазо-частотная характеристики:

Список использованной литературы:

1. Баскаков С.И. Радиотехнические цепи и сигналы: Учебник для вузов по специальности «радиотехника».-2-е изд., перераб. И доп. — М.: Высш.шк., 1988-448 с.: ил.

2. Гоноровский И.С. Радиотехнические цепи и сигналы: Учебник для вузов.-4-е изд., перераб. и доп. — М.: Радио и связь,1986.-512 с.: ил.

Как получить верный спектр сигнала

В общем задался я целью посмотреть как выглядят спектры для разных спектров. Делал так.
Задал функцию реализующую исследуемый сигнал и стал применять к ней функции разложения в ряд Фурье(fft, FFT, cfft, CFFT). Вот только почему то графики получились разные. Вернее для функций fft и cfft они были идентичны, а остальные сильно отличались (файл во вложении).
В связи с этим меня сильно интересует два вопроса:
1) Почему результат у функций разный?
2) Как-же получить ВЕРНЫЙ спектр некоторой функции?

Научился картинки в пост вставлять
На первой картинке показана функция сигнала и её график, спектр которой я хочу получить. А также результат применения к этой функции функций fft и FFT:

Вложения

Сравнение алгоритмов Фурье.zip (62.7 Кб, 231 просмотров)

Как получить график Частотно Импульсной Модуляции (ЧИМ) импульсного сигнала и его спектр
Как сделать не знаю. Какие формулы вводить, значения какие брать, вообще не понятно. В инете ничего.

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

Провести спектральный анализ сигнала. Отобразить амплитудный и фазовый спектр сигнала
Так же разработать процедуру для вычисления дискретного преобразования Фурье заданного сигнала.

Вложения

Спектр гармонического колебания через ряды Фурье.rar (37.7 Кб, 66 просмотров)

Как построить спектр сигнала в MathCAD 15
Использую встроенную функцию cfft, но чем больше точек я учитываю, тем больше максимальное значение.

Как строить спектр чм сигнала в маткаде
Доброго времени суток. Подскажите, пожалуйста, как в маткаде можно построить спектр ЧМ сигнала вида

Спектр сигнала
Здравствуйте. Подскажите пожалуйста как построить спектр одиночного сигнала. С помощью БПФ.

Построить спектр ФМ сигнала
Удалось построить ФМС, но маткад не хочет строить спектр. Пробовал через интегрально.

Амплитудный спектр АМ-сигнала
Проблема такая: амплитудный спектр АМ сигнала, маткад строит только 1 гармонику,а две боковые.

Амплитудый спектр сигнала
Есть ли какой-то простой способ вывести амплитудный спектр сигнала в MathCAD? Исходные данные -.

Как найти спектр сигнала в mathcad

СПЕКТРАЛЬНЫЙ АНАЛИЗ
В ПАКЕТЕ ПРОГРАММ MATHCAD

Лабораторная работа №1
Спектральный анализ в пакете программ Mathcad

1.1. изучение функций Mathcad для анализа спектра сигналов;
1.2. изучение функций Mathcad для работы с файлами данных;
1.3. исследование спектра сигналов.

2. КРАТКИЕ ПОЯСНЕНИЯ К ЛАБОРАТОРНОЙ РАБОТЕ

Корпорация "Центр"

2.1. Спектры периодических сигналов
Периодические (полигармонические) сигналы
X(t)=X(t+kT), (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=KFд / 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 с.

Корпорация "Центр"

Добавить комментарий

Ваш адрес email не будет опубликован. Обязательные поля помечены *