Цифровая обработка сигналов ЦОС / Матлаб / Лаба 4 / C4_106_Izumtsev_RT-31
.pdf
ЛР 4 — Изюмцев Елисей РТ-31 — Вариант 106
Задание 1
clc clear all close all
format short g
N = 106;
tu = N * 10^-6;
f_nes = (1 + 1/N) * 10^6; f_0v5 = 1/tu;
f_dis = 16 * 10^6;
f_vec = [0 (f_nes/2)/(f_dis/2) (f_nes-f_0v5/2)/(f_dis/2) f_nes/(f_dis/2) (f_nes+f_0v5/2)/(f_dis/2) (3*f_nes/2)/(f_dis/2) 1];
m_vec = [0.1 0.1 0.707 1 0.707 0.1 0.1];
a = 1; n = 16;
b = fir2(n,f_vec,m_vec); [h, t] = impz(b,a); figure
stem (t, h) figure freqz(b,a)
N |
Импульсная характеристика |
Частотная характеристика |
|
|
|
8
16
32
64
128
clc clear all close all
format short g
N = 106;
tu = N * 10^-6;
f_nes = (1 + 1/N) * 10^6; f_0v5 = 1/tu;
f_dis = 16 * 10^6;
f_vec = [0 (f_nes/2)/(f_dis/2) (f_nes-f_0v5/2)/(f_dis/2) f_nes/(f_dis/2) (f_nes+f_0v5/2)/(f_dis/2) (3*f_nes/2)/(f_dis/2) 1];
m_vec = [0.1 0.1 0.707 1 0.707 0.1 0.1];
a = 1; n = 16;
b = fir2(n,f_vec,m_vec); t=-tu:1/f_dis:2*tu; A=rectpuls(t-tu/2, tu); u=A.*sin(2*pi*f_nes*t); y=filter(b, a, u); figure
subplot(2,1,1)
plot(u)
title('INPUT SIGNAL u(n)') xlabel('Number sample, n') ylabel('Voltage x(n)') grid on
subplot(2,1,2)
plot(y)
title(' OUTPUT SIGNAL y(n)') xlabel('Number sample,, n') ylabel('Voltage y(n)') grid on
Реакция фильтра на входное воздействие в виде прямоугольного радиоимпульса
Совпадает с полосой сигнала |
Больше на 30% чем полоса |
Меньше на 30% чем полоса |
по первым нулям его спектра |
сигнала |
сигнал |
|
|
|
|
|
|
Задание 2
clc clear all close all
format short g
N = 106;
f_nes = (0.3 + 0.05 * N) * 10^6; f_0v5 = (64 + 2 * N) * 10^3; f_dis = 25 * 10^6;
f_vec = [0 (f_nes/2)/(f_dis/2) (f_nes-f_0v5/2)/(f_dis/2) f_nes/(f_dis/2) (f_nes+f_0v5/2)/(f_dis/2) (3*f_nes/2)/(f_dis/2) 1];
m_vec = [1 1 0.0001 0 0.0001 1 1];
a = 1; n = 256;
b = fir2(n,f_vec,m_vec); [h, t] = impz(b,a); figure
stem (t, h) figure freqz(b,a)