Регулирование мощности шума по принципам адаптивной фильтрации, отвечающей высоким требованиям к точности настройки
Дипломная работа - Безопасность жизнедеятельности
Другие дипломы по предмету Безопасность жизнедеятельности
вая обработка сигналов. - 2005. - №2. - с.11-14.
.Лайонс Р. Цифровая обработка сигналов: Второе издание. Пер. с англ. - М.: ООО Бином-Пресс, 2006 г. - 656 с.
.Электронный справочник по функциям Matlab (
.Стешенко В. Б. ПЛИС фирмы ALTERA: Элементная база, система проектирования и языки описания аппаратуры. - М.: Издательский дом Додэка - XXI, 2002.
.Электронный справочник по микросхемам фирмы Altera (www.altera.com).
.Грибов В.Д., Грузинов В.П. Экономика предприятия. - М.: Финансы и статистика, 2007. - 336 с.
.Сергеев И.В. Экономика организации предприятия. - М.: Финансы и статистика, 2007. - 574с.
.Григорьев Л.Н. Радиолокационные системы. Учебное пособие. МИРЭА (ТУ) - М,,2003.-144с.
.Г.М. Батишева, М.В. Забродина, Р.Р. Курганова, Н.Д. Николаева, И.А. Нуль, И.Р. Никорук, ".И. Солнцев, В.И. Солнцев. Экономика предприятия и основы предпринимательства - М.: МИРЭА, 2001. - 20с.
.ГОСТ 12.2.032-78 СС БТ. Рабочее место при выполнении работы сидя. Общие эргономические требования.
.ГОСТ 12.1.005-88 СС БТ. Воздух рабочей зоны. Общие санитарно-гигиенические требования
.Безопасность жизнедеятельности/Под ред. Русака О.Н.- С.-Пб.: ЛТА, 1996.
.СП 52.13330.2011 Естественное и искусственное освещение.
.Безопасность жизнедеятельности. Арустамов Э.А. - 2006
Приложение А
Программный код модели погрешности расчета коэффициента передачи
shary_Errors;all
% noise.txt содержит дискретные отсчеты шума в виде таблицы,
% в которой первый столбец - номер отсчета, второй - значение %шума.
IN = load(noise.txt);= IN(:,2);= cat(1,NOISE,NOISE,NOISE,NOISE,NOISE);= cat(1,NOISE,NOISE);
%параметры= 32;
M = 0.01;
razm = 16;= 70;
% расчет коэффициентов передачи для всех алгоритмов
g = 0;G = cat(2, 0.3:0.1:1, 1.3, 1.5, 1.7, 2, 2.5, 3 )= g+1;(g) = G;= NOISE.*G;
%%%%%% ско теор
[H_OUT_E1,Wgg_E1,CKO_RAS_E1] = shary_teor_do_kv(X,CKOt,0,K,M);
%%%%%% ско кв теор
[H_OUT_E2,Wgg_E2,CKO_RAS_E2] = shary_teor_posle_kv
(X,CKOt,0,K,M);
%%%%%% дисп теор
[H_OUT_ED,Wgg_ED,CKO_RAS_ED] = shary_teor_disp_do_kv_E
(X,CKOt,0,K,M);
%%%%%% ско
[H_OUT_1,Wgg_1,CKO_RAS_1] = shary_do_kv(X,CKOt,razm,K,M);
% %%%%% кв ско
[H_OUT_2,Wgg_2,CKO_RAS_2] = shary_posle_kv(X,CKOt,razm,K,M);
%%%%%% дисп
[H_OUT_D1,Wgg_D1,CKO_RAS_D1] = shary_disp_do_kv
(X,CKOt,razm,K,M);
%%%%%% кв дисп
[H_OUT_D2,Wgg_D2,CKO_RAS_D2] = shary_disp_posle_kv
(X,CKOt,razm,K,M);
%%%%% УСТОЯВШИЙСЯ КОЭФФИЦИЕНТ= Wgg_E1;
WE1_dB(g) = 20*log10(WE1);= Wgg_E2;_dB(g) = 20*log10(WE2);= Wgg_ED;_dB(g) = 20*log10(WED);= Wgg_D1;_dB(g) = 20*log10(WD1);= Wgg_D2;_dB(g) = 20*log10(WD2);= Wgg_1;_dB(g) = 20*log10(W1);= Wgg_2;_dB(g) = 20*log10(W2);
%%%% Вид шума на входе=NOISE(1:8000);
figure(1)(NOISE)on(2)
stem(NOISE(1000:1100))on
%%%% изменение коэффициента передачи
figure(3)(Gn,WE1_dB,Gn, WE2_dB, Gn, WED_dB,Gn, WD1_dB,Gn, WD2_dB,Gn, W1_dB,Gn, W2_dB)(W ско теор, W ско кв теор, W дисп теор,W дисп, W кв дисп, W ско, W кв ско)
grid on
%%%% отклонение от теоретического алгоритма
figure(4)(Gn,WE1_dB-WE1_dB,Gn, WE2_dB-WE1_dB, Gn, WED_dB-WE1_dB,Gn, WD1_dB-WE1_dB,Gn, WD2_dB-WE1_dB,Gn, W1_dB-WE1_dB,Gn, W2_dB-WE1_dB)(W ско теор, W ско кв теор, W дисп теор,W дисп, W кв дисп, W ско, W кв ско)on
%%%%%% Вспомогательные функцииY = hilb_tran(X)= length(X);mod(N,2)==1 N=N-1; end= 1:2:N-1;(:,1) = X(k);= 2:2:N;(:,2) = X(n);= 1:2:N/2-1;= 2:2:N/2;(n,1) = Y(n,1).*-1;(k,2) = Y(k,2).*-1;[H_OUT,W_out,CKO_OUT] = shary_teor_do_kv
(X,CKOt,razr,K,M);= length(X);_ocenok = floor(N/K);_razr = 2^razr;= N_razr;k=1:N_ocenok(((k-1)*K+1):((k-1)*K+K))=(X( ((k-1)*K+1):((k-1)*K+K)).*Wg./N_razr);= Z(((k-1)*K+1):((k-1)*K+K));= X(((k-1)*K+1):((k-1)*K+K));= sqrt(sum(Zn.^2)/K);= sqrt(sum(Xn.^2)/K);= 2*M*N_razr*(CKOt-CKOz);= Qm / CKOx;= Wg + Qd;_OUT(k) = CKOz;_OUT = hilb_tran(Z);_out = Wg/N_razr;[H_OUT,W_out, CKO_OUT] = shary_teor_posle_kv
(IN,CKOt,razr,K,M);= length(IN);= K*2;_ocenok =floor(N/Kh);_razr = 2^razr;= N_razr;k=1:1:Nh_ocenok_OUT(((k-1)*K+1):((k-1)*K+K),:) = hilb_tran
(floor(IN(((k-1)*Kh+1):(k-1)*Kh+Kh).*Wg./N_razr));= H_OUT(((k-1)*K+1):((k-1)*K+K),1);= H_OUT(((k-1)*K+1):((k-1)*K+K),2);_out = sqrt(sum(X.^2+Y.^2)./Kh);= IN(((k-1)*Kh+1):((k-1)*Kh+Kh));_in = sqrt(sum( Xin.^2)/Kh);= 2*M*N_razr*(CKOt-CKO_out);= Qm / CKO_in;= Wg + Qd;_OUT(k) = CKO_out;_out = Wg / N_razr;[H_OUT,W_out,CKO_OUT] = shary_teor_disp_do_kv_E
(X,CKOt,razr,K,M);= length(X);_ocenok = floor(N/K);_razr = 2^razr;= N_razr;k=1:N_ocenok(((k-1)*K+1):((k-1)*K+K)) =(X( ((k-1)*K+1):((k-1)*K+K)).*Wg./N_razr);= Z(((k-1)*K+1):((k-1)*K+K));= X(((k-1)*K+1):((k-1)*K+K));= sum(Zn.^2) / K;= sum(Xn.^2) / K;= 2*M*N_razr*(CKOt.^2-CKOz);=Qm / CKOx;= Wg + Qd;_OUT(k) = floor(sqrt(floor(sum(Zn.^2)/K)));_OUT = hilb_tran(Z);_out = Wg / N_razr;[H_OUT,W_out,CKO_OUT]=shary_do_kv(X,CKOt,razr,K,M);= length(X);_ocenok = floor(N/K);_razr = 2^razr;= N_razr;k=1:N_ocenok(((k-1)*K+1):((k-1)*K+K)) =(X( ((k-1)*K+1):((k-1)*K+K)).*Wg./N_razr);= Z(((k-1)*K+1):((k-1)*K+K));= X(((k-1)*K+1):((k-1)*K+K));= floor(sum(abs(Zn))/K);= floor(sum(abs(Xn))/K);= floor(2*M*N_razr)*(CKOt-CKOz);= floor(Qm/CKOx);= Wg + Qd;_OUT(k) = CKOz;_OUT = hilb_tran(Z);_out = Wg / N_razr;[H_OUT,W_out,CKO_OUT]=shary_posle_kv(IN,CKOt,razr,Kh,M);= length(IN);= Kh / 2;_ocenok = floor(N/Kh);_razr = 2^razr;= N_razr;k=1:Nh_ocenok_OUT(((k-1)*K+1):((k-1)*K+K),:)=hilb_tran
(floor(IN(((k-1)*Kh+1):(k-1)*Kh+Kh).*Wg./N_razr));= abs(H_OUT(((k-1)*K+1):((k-1)*K+K),1));= abs(H_OUT(((k-1)*K+1):((k-1)*K+K),2));_out=floor(sum(7/8*max(X,Y)+9/16*min(X,Y))./(K*sqrt(2)));= IN(((k-1)*Kh+1):((k-1)*Kh+Kh));_in = floor(sum( abs(Xn))./Kh) ;= floor(2*M*N_razr)*(CKOt-CKO_out);= floor(Qm/CKO_in);= Wg + Qd;_OUT(k) = CKO_out;_out=Wg / N_razr;[H_OUT,W_out,CKO_OUT]=shary_disp_do_kv(X,CKOt,razr,K,M);= length(X);_ocenok = floor(N/K);_razr = 2^razr;= N_razr;k=1:N_ocenok(((k-1)*K+1):((k-1)*K+K)) =(X( ((k-1)*K+1):((k-1)*K+K)).*Wg./N_razr);= Z(((k-1)*K+1):((k-1)*K+K));= X(((k-1)*K+1):((k-1)*K+K));= floor(sum(Zn.^2)/K);= floor(sum(Xn.^2)/K);= floor(2*M*N_razr)*(CKOt.^2-CKOz);= floor(Qm/CKOx);= Wg + Qd;_OUT(k) = floor(sqrt(floor(sum(Zn.^2)/K)));_OUT = hilb_tran(Z);_out = Wg / N_razr;[H_OUT,W_out,CKO_OUT] = shary_disp_posle_kv
(IN,CKOt,razr,Kh,M);= length(IN);= Kh/2;_ocenok = floor(N/Kh);_razr = 2^razr;= N_razr;k=1:Nh_ocenok_OUT(((k-1)*K+1):((k-1)*K+K),:) = hilb_tran
(floor(IN(((k-1)*Kh+1):(k-1)*Kh+Kh).*Wg./N_razr));= abs(H_OUT(((k-1)*K+1):((k-1)*K+K),1));= abs(H_OUT(((k-1)*K+1):((k-1)*K+K),2));_out = floor(sum(X.^2+Y.^2)./Kh);= IN(((k-1)*Kh+1):((k-1)*Kh+Kh));_in = floor(sum(Xn.^2)./Kh) ;= floor(2*M*N_razr)*(CKOt.^2-CKO_out);= floor(Qm/CKO_in);= Wg + Qd;_OUT(k) = sqrt(sum(H_OUT(((k-1)*K+1):((k-1)*K+K),1).^2
+ H_OUT(((k-1)*K+1):((k-1)*K+K),2).^2)./Kh);_out=Wg / N_razr;
Приложение Б
Программный код модели, демонстрирующий работу системы ШАРУ
function shary_disp_demoall
% noise.txt содержит дискретные отсчеты шума в виде таблицы,
% в которой первый столбец - номер отсчета, второй - значение %шума.
IN = load(noise.txt);= IN(:,2);
%параметры= 15;_razm = 2^razm; = 70;= 32;= 0.01;
% расчет СКО в процессе подстройки для 3х случаев
g=0;G = [0.5 1 2]= NOISE.*G;% при 1 cko примерно 70= g+1;
[Y(:,g), W_out, Wgp(:,g),CKO_OUT(:,g)] =_disp_do_kv(X,CKOt,razm,K,M);_ocenok = floor(length(NOISE)/K);k=1:K_ocenok_S = ((k-1)*K+1):((k-1)*K+K);_SHUM_IN(k,g) = sqrt(sum(X(D_S).^2)/K);
end
%изменение шума в процессе подстройки для G=2
OUT(1:1999) = 2*NOISE(1:1999);(2000:2000+length (Y(:,3))-1) = Y(:,3);= OUT(1:8000);(1)( OUT)
grid on
% графики изменения СКО в процессе подстройки для 3х случаев