Отчет по семинару №2.

Выполнил: Краснов Д.С., МП-47.

1. Номер варианта и условие (исходные данные) задания (задач).

21 вариант (1 по модулю)

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

С помощью графических средств MATLAB, например, функции stem или plot, отобразите дискретную последовательность . Примените ДПФ к и также с помощью функции stem или plot отобразите амплитудный и фазовый спектр (амплитудные и фазовые спектральные отсчеты).

Определите следующие характеристики исследуемых сигналов – ширину спектра по уровню 0,707 и по ближайшим нулям, частоту на которой сигнал достигает максимума, значение трех ближайших по уровню боковых компонент спектра и заполните таблицу 1.

1. Одиночный видеоимпульс прямоугольной формы длительностью  = 2 мкс.

2. Последовательность 4-х видеоимпульсов треугольной формы длительностью  = 6 мкс с передним фронтом бесконечной крутизны и длительностью паузы между ними 12 мкс.

3. Одиночный Гауссов радиоимпульс на несущей частоте =3 МГц, относительной шириной спектра 10% от несущей на уровне -10 дБ.

Задание 1.

Одиночный видеоимпульс прямоугольной формы длительностью  = 2 мкс.

Код программы:

clear; % очищаем память

t0=0; % начало импульса 0 с.

Fs=48e6; % частота дискретизации 48 МГц

taui=2e-6; % длительность импульса 2 мкс

t=t0-taui/4:1/Fs:t0+2*taui; % вектор дискретных значений времени t(n)

u=rectpuls(t-taui/2,taui)

axis tight; % масштаб такой, чтобы вошел весь диапазон значений по обеим осям

plot(t,u); % график сигнала

title('Signal u(n)'); % заголовок графика

xlabel('time, t=nT=n/Fs (n max=N), [s]'); % надписи по горизонтальной оси

ylabel('u(n), [V]'); % надписи по вертикальной оси

ylim([-0.1 1.1]) %масштаб по оси ординат (минимальное и максимальное значение)

Nfft=2^nextpow2(length(u)); % ближайшая к числу отсчетов степень двойки

scompl=fft(u,Nfft); % комплексный спектр сигнала

s=abs(scompl); % амплитудный спектр сигнала

phi=180/pi*angle(scompl); % фазовый спектр сигнала

figure; % создание нового окна

f =0:Nfft-1; % вектор отсчетов нормированной частоты для построения ДПФ

f=(f/(Nfft-1))*Fs; % вектор отсчетов частоты для ДПФ

stem(f,s); % отображение амплитудного спектра

title('Spectre abs(scompl(f))');

xlabel('Number, n');

ylabel('amplitude s(f), [V]');

figure;

phi=unwrap(phi); % сглаживаем график фазового спектра

plot(phi); % построение фазового спектра

title('Spectre angle(scompl(f))');

xlabel('Number, n');

ylabel('angle s(f), [degree]');