Добавил:
debilX2
t.me
Фулл всегда есть! А если нет, то Вы плохо его ищите! ИиКГ, СКДИКТ, ОКИТПЭС и тд https://t.me/whitedevil752rn Так же веду разработку КД (конструкторской документации) согласно ГОСТ. Имеется опыт работы при производстве на одном из ведущих в области радиэлектроники предприятии. Пишите)
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:pr1 задачка с доски / pr1
.m clear
clc
close all
% шум
N = 100;
m=8;
last = N;
load('Num.mat');
noize = dsp.ColoredNoise(1, N, 'Color','Brown');
x = noize();
% фильтр
y = filter(Num,1,x);
% задержка на 1 такт
x_delayed(1:last-1) = y(2:last);
% суммирование
x_trans = transpose(x_delayed);
x_dopoln = [x_trans; 0] ;
x_sum = x_dopoln - y;
% Нормирование
Max = max(x_sum);
x_norm=x_sum./Max;
d = 1/(2^(m-1));
part = [-1:d:1-d];
code = [-1:d:1];
[~,q_sig,distortion] = quantiz(x_norm,part,code);
plot(y)
hold on
c=cumsum(q_sig);
h = filter(Num,1,c);
h_domn = h.* (-1);
plot(h)
grid
clc
close all
% шум
N = 100;
m=8;
last = N;
load('Num.mat');
noize = dsp.ColoredNoise(1, N, 'Color','Brown');
x = noize();
% фильтр
y = filter(Num,1,x);
% задержка на 1 такт
x_delayed(1:last-1) = y(2:last);
% суммирование
x_trans = transpose(x_delayed);
x_dopoln = [x_trans; 0] ;
x_sum = x_dopoln - y;
% Нормирование
Max = max(x_sum);
x_norm=x_sum./Max;
d = 1/(2^(m-1));
part = [-1:d:1-d];
code = [-1:d:1];
[~,q_sig,distortion] = quantiz(x_norm,part,code);
plot(y)
hold on
c=cumsum(q_sig);
h = filter(Num,1,c);
h_domn = h.* (-1);
plot(h)
grid