Этот пример генерирует последовательность импульсов, используя прямоугольный импульс по умолчанию модуля ширины. Частота повторения составляет 0,5 Гц, длина сигнала составляет 60 с, и частота дискретизации составляет 1 кГц. Коэффициент усиления является синусоидой частоты 0,05 Гц.

t = 0:1/1e3:60;
d = [0:2:60;sin(2*pi*0.05*(0:2:60))]';
x = @rectpuls;
y = pulstran(t,d,x);

plot(t,y)
hold off
xlabel('Time (s)')
ylabel('Waveform')

Этот пример генерирует асимметричную пилообразную форму волны с частотой повторения 3 Гц. Пилюля имеет ширину 0,2 с и коэффициент наклона -1. Длина сигнала составляет 1 с, а частота дискретизации - 1 кГц. Постройте график последовательности импульсов.

fs = 1e3;
t = 0:1/1e3:1;
d = 0:1/3:1;   
x = tripuls(t,0.2,-1);
y = pulstran(t,d,x,fs);

plot(t,y)
hold off
xlabel('Time (s)')
ylabel('Waveform')

Постройте график 10 кГц по Гауссову РЧ импульсу с 50% полосой пропускания, дискретизированной со скоростью 10 МГц. Обрезайте импульс, где огибающая падает на 40 дБ ниже пика.

fs = 1e7;
tc = gauspuls('cutoff',10e3,0.5,[],-40); 
t = -tc:1/fs:tc; 
x = gauspuls(t,10e3,0.5); 

plot(t,x)
xlabel('Time (s)')
ylabel('Waveform')

Частота повторения импульсов составляет 1 кГц, частота дискретизации - 50 кГц, а длина последовательности импульсов - 25 мс. Коэффициент усиления является синусоидой частоты 0,1 Гц.

ts = 0:1/50e3:0.025;
d = [0:1/1e3:0.025;sin(2*pi*0.1*(0:25))]';
y = pulstran(ts,d,x,fs);

Постройте график периодической последовательности импульсов.

plot(ts,y)
xlim([0 0.01])
xlabel('Time (s)')
ylabel('Waveform')

Напишите функцию, которая генерирует пользовательские импульсы, состоящие из синусоиды, демпфированной экспоненциальной. Импульс является нечетной функцией времени. Функция генерации имеет второй входной параметр, который задает одно значение для синусоидальной частоты и коэффициента затухания. Отобразите сгенерированный импульс, дискретизированный на частоте 1 кГц в течение 1 секунды, со значением частоты и демпфирования, оба равными 30.

fnx = @(x,fn) sin(2*pi*fn*x).*exp(-fn*abs(x));

ffs = 1000;
tp = 0:1/ffs:1;

pp = fnx(tp,30);

plot(tp,pp)
xlabel('Time (s)')
ylabel('Waveform')

Используйте pulstran функция для генерации train пользовательских импульсов. Дискретизацию train производят при частоте дискретизации 2 кГц в течение 1,2 с. Импульсы происходят каждую треть секунды и имеют экспоненциально уменьшающиеся амплитуды.

Первоначально задайте сгенерированный импульс в качестве прототипа. Включите частоту выборки прототипа в вызов функции. В этом случае pulstran повторяет импульсы в заданных местоположениях.

fs = 2e3;
t = 0:1/fs:1.2;

d = 0:1/3:1;
dd = [d;4.^-d]';

z = pulstran(t,dd,pp,ffs);

plot(t,z)
xlabel('Time (s)')
ylabel('Waveform')

Сгенерируйте последовательность импульсов снова, но теперь используйте сгенерирующую функцию как входной параметр. Включите частоту и параметр демпфирования в вызов функции. В этом случае pulstran генерирует импульс так, чтобы он был центрирован около нуля.

y = pulstran(t,dd,fnx,30);

plot(t,y)
xlabel('Time (s)')
ylabel('Waveform')

Написание функции, которая генерирует пользовательскую экспоненциально затухающую пилообразную форму волны с частотой 0,25 Гц. Функция генерации имеет второй входной параметр, который задает одно значение для частоты пилообразного дерева и коэффициента затухания. Отобразите сгенерированный импульс, дискретизированный с частотой 0,1 кГц в течение 1 секунды, со значением частоты и демпфирования, равными 50.

fnx = @(x,fn) sawtooth(2*pi*fn*0.25*x).*exp(-2*fn*x.^2);

fs = 100;
t = 0:1/fs:1;

pp = fnx(t,50);

plot(t,pp)

Используйте pulstran функция для генерации train пользовательских импульсов. Дискретизация train составляет 0,1 кГц в течение 125 секунд. Импульсы происходят каждые 25 секунд и имеют экспоненциально уменьшающиеся амплитуды.

Задайте сгенерированный импульс в качестве прототипа. Сгенерируйте три последовательности импульсов с помощью метода линейной интерполяции по умолчанию, интерполяции по ближайшему соседу и кусочно-кубической интерполяции. Сравните последовательности импульсов на одном графике.

d = [0:25:125; exp(-0.015*(0:25:125))]';
ffs = 100;
tp = 0:1/ffs:125;

r = pulstran(tp,d,pp);
y = pulstran(tp,d,pp,'nearest');
q = pulstran(tp,d,pp,'pchip');


plot(tp,r)
hold on
plot(tp,y)
plot(tp,q)
xlim([0 125])
legend('Linear interpolation','Nearest neighbor interpolation','Piecewise cubic interpolation')
hold off