MATLAB.Exponenta
MATLAB и Simulink на русском
Технологии разработки и отладки
		сложных технических систем

Обработка сигналов и изображений\Signal Processing Toolbox

Список функций Signal Processing Toolbox: Генерация сигналов

  В оглавление \ К следующему разделу \ К предыдущему разделу

CHIRP

Генерация сигнала с изменяющейся частотой

Синтаксис:

y = chirp(t,f0,t1,f1)
y = chirp(t,f0,t1,f1,'method')
y = chirp(t,f0,t1,f1,'method',phi)
y = chirp(t,f0,t1,f1,'quadratic',phi,'shape')

Описание:

y = chirp(t,f0,t1,f1)

Генерирует отсчеты сигнала с линейной частотной модуляцией для моментов времени, заданных входным параметром t. Входной параметр f0 задает мгновенную частоту сигнала в нулевой момент времени, а параметр f1 задает мгновенную частоту сигнала для момента времени t1. Оба параметра f0 и f1 задаются в герцах. По умолчанию используются следующие значения: f0 = 0, t1 = 1, f1 = 100.

y = chirp(T,F0,T1,F1,'method')

В данном варианте синтаксиса дополнительный входной строковый параметр 'method' позволяет задать закон изменения мгновенной частоты. Возможны следующие значения этого параметра:

  • 'linear' - линейный закон. Зависимость мгновенной частоты от времени описывается следующей формулой:



    где



    Такой выбор параметра обеспечивает в момент t1 значение мгновенной частоты, равное f1.

  • 'quadratic' - квадратичный закон. Зависимость мгновенной частоты от времени описывается следующей формулой:



    где



    Если f0 > f1 (уменьшение частоты), зависимость мгновенной частоты от времени по умолчанию является выпуклой (convex). Если f0 < f1 (увеличение частоты), зависимость мгновенной частоты от времени по умолчанию является вогнутой (concave).

  • 'logarithmic' - логарифмический закон. Несмотря на название метода, мгновенная частота в данном случае зависит от времени не по логарифмическому, а по экспоненциальному закону:



    где



    В данном случае должно выполняться неравенство f1 > f0.

Все три закона изменения частоты могут задаваться сокращенными строками 'li', 'q' и 'lo' соответственно.

y = chirp(t,f0,t1,f1,'method',phi)

Дополнительно позволяет задать начальную фазу колебания (в градусах) с помощью входного параметра phi. По умолчанию предполагается, что phi = 0. Значения по умолчанию подставляются вместо пустых ([]) и неуказанных параметров.

y = chirp(t,f0,t1,f1,'quadratic',phi,'shape')

Дополнительно позволяет задать форму квадратичной зависимости с помощью строкового параметра 'shape'. Возможны два значения этого параметра - concave (вогнутый) и convex (выпуклый). Они описывают форму параболы частотного закона в области положительных значений времени. Если параметр 'shape' опущен, по умолчанию форма параболы является выпуклой при уменьшении частоты (f0 > f1) и вогнутой при увеличении частоты (f0 < f1).

Примеры

Пример 1

Построим спектрограмму для сигнала с линейной частотной модуляцией:

t = 0:0.001:2;              % 2 секунды с частотой 1 кГц
y = chirp(t,0,1,150);       % частота меняется от 0 (при t=0) до 150 Гц (при t=1 с)
specgram(y,256,1e3,256,250) % спектрограмма 
title('Linear chirp')
Пример 2

Построим спектрограмму для сигнала с квадратичной частотной модуляцией:

t = -2:0.001:2;                     % ±2 секунды с частотой 1 кГц
y = chirp(t,100,1,200,'quadratic'); % при t=0 частота 100 Гц,
                                    % при t=1 частота 200 Гц
specgram(y,128,1e3,128,120)         % спектрограмма
title('Quadratic chirp')
Пример 3

Построим спектрограмму для сигнала с выпуклым законом квадратичной частотной модуляции:

t= -1:0.001:1;                    % ±1 секунда с частотой 1 кГц
fo=100; f1=400;                   % при t=0 частота 100 Гц, 
                                  % при t=1 частота 400 Гц
y=chirp(t,fo,1,f1,'q',[],'convex');
specgram(y,256,1000)              % спектрограмма
title('Convex chirp')
Пример 4

Построим спектрограмму для сигнала с вогнутым законом квадратичной частотной модуляции:

t= -1:0.001:1;                    % ±1 секунда с частотой 1 кГц
fo=100; f1=400;                   % при t=0 частота 100 Гц, 
                                  % при t=1 частота 400 Гц
y=chirp(t,fo,1,f1,'q',[],'concave');
specgram(y,256,1000)              % спектрограмма 
title('Concave chirp')

Сопутствующие функции: cos, diric, gauspuls, pulstran, rectpuls, sawtooth, sin, sinc, square, tripuls

  В оглавление \ К следующему разделу \ К предыдущему разделу


Поиск по сайту:

Система Orphus

Яндекс.Метрика