lab5_variant2

k = 0:1:(10^5-1);
A = 1.02;
w = 0.72;
phi = 4.99;
x = A*cos(w*k+phi);
%% Шаг квантования 256
x_q = round(x*256)/256; %квантование
delta_x = x - x_q;
figure(1);
plot(k,delta_x);
%Шум квантования в 200 отсчетов
figure(2)
plot(delta_x(1:200))
%Гистограмма
hist(delta_x, 100)
%xcorr
[Rx, dk] = xcorr(delta_x, 100, 'unbiased');
figure(3)
plot(dk,Rx);
%Для расчета оценки и построения графика СПМ шума квантования в дан-
%ной работе используется функция pwelch
Lb = 256;
figure(4)
pwelch(delta_x, Lb);
%% Шаг квантования 16
x_q = round(x*16)/16; %квантование
delta_x = x - x_q;
figure(1);
plot(k,delta_x);
%Шум квантования в 200 отсчетов
figure(2)
plot(delta_x(1:200))
%Гистограмма
hist(delta_x, 100)
%xcorr
[Rx, dk] = xcorr(delta_x, 100, 'unbiased');
figure(3)
plot(dk,Rx);
%Для расчета оценки и построения графика СПМ шума квантования в дан-
%ной работе используется функция pwelch
Lb = 256;
figure(4)
pwelch(delta_x, Lb);
%% Шаг квантования 1
x_q = round(x*1)/1; %квантование
delta_x = x - x_q;
figure(1);
plot(k,delta_x);
%Шум квантования в 200 отсчетов
figure(2)
plot(delta_x(1:200))
%Гистограмма
hist(delta_x, 100)
%xcorr
[Rx, dk] = xcorr(delta_x, 100, 'unbiased');
figure(3)
plot(dk,Rx);
%Для расчета оценки и построения графика СПМ шума квантования в дан-
%ной работе используется функция pwelch
Lb = 256;
figure(4)
pwelch(delta_x, Lb);
%% дискретный белый гауссов шум (БГШ)
sn = randn(1, 10^5);
sn2 = sn/max(abs(sn));
x_q = round(sn*256)/256; %квантование
delta_x = sn2 - x_q;
figure(1);
%plot(k,delta_x);
plot(k,sn2);
%Шум квантования в 200 отсчетов
figure(2)
plot(delta_x(1:200))
%Гистограмма
hist(delta_x, 100)
%xcorr
[Rx, dk] = xcorr(delta_x, 100, 'unbiased');
figure(3)
plot(dk,Rx);
%Для расчета оценки и построения графика СПМ шума квантования в дан-
%ной работе используется функция pwelch
Lb = 256;
figure(4)
pwelch(delta_x, Lb);
%% Речевой сигнал
load mtlb;
x = mtlb/max(abs(mtlb));
x_q = round(x*256)/256; %квантование
delta_x = x - x_q;
figure(1);
%plot(k,delta_x);
plot(x);
%Шум квантования в 200 отсчетов
figure(2)
plot(delta_x(1:200))
%Гистограмма
hist(delta_x, 100)
%xcorr
[Rx, dk] = xcorr(delta_x, 100, 'unbiased');
figure(3)
plot(dk,Rx);
%Для расчета оценки и построения графика СПМ шума квантования в дан-
%ной работе используется функция pwelch
Lb = 256;
figure(4)
pwelch(delta_x, Lb);