Технические разделы

Профилактика пожара

Приложение 3

Программный код и результаты моделирования для тестового сигнала

%Исходные данные= 96e3;= 1.5e3;

f2 = 20e3;= 25e3;= 1/Fs;= 1/f1;= 0: dt: T1*50;_shift = 10000;= (j*sin (2*pi*f1*t) + cos (2*pi*f1*t)) + (j*sin (2*pi*f2*t) + cos (2*pi*f2*t)) + (j*sin (2*pi*f3*t) + cos (2*pi*f3*t));

% shift to center frequnce for real filtering= S1. * exp (j*2*pi* (-df_shift) *t);

%Окно Хэмминга и построение спектра исходного сигнала

w = hamming (2048). ';= [];i=1: floor (length (S1) /2048)= [S2 S1 ( (i-1) *2048+1: i*2048). *w];=spectrum. periodogram;(1)= psd (Hs1,S2,'Fs',Fs); plot (h1)

%Кан. Фильтр order = 122

%fs = 96000 fpass = 10000 fstop = 12000 astop = 80= [zeros (1,length (k0) - 1) S1];= zeros (1, length (S1));i = length (k0): length (S1)n = 0: length (k0) - 1(i) = r0 (i) + S1 (i-n) *k0 (n+1);(1: length (k0) - 1) = [];on= psd (Hs1,r0,'Fs',Fs); plot (h1)

%фильтр перенос order = 33

%fs = 768000 fpass = 10000 fstop = 55000 astop = 60= upsample (r0,8);= [zeros (1,length (kn) - 1) B0];= zeros (1,length (B0));i = length (kn): length (B0)n = 0: length (kn) - 1(i) = r1 (i) +B0 (i-n) *kn (n+1);(1: length (kn) - 1) = [];(2)= psd (Hs1,r1,'Fs',768000); plot (h1)

%Первый FiR-фильтр order = 8

%fs = 1536000 fpass = 10000 fstop = 384000 astop = 80= upsample (r1,2);= [zeros (1,length (k1) - 1) B];= zeros (1,length (B));i = length (k1): length (B)n = 0: length (k1) - 1(i) = r2 (i) +B (i-n) *k1 (n+1);

%Второй FiR-фильтр order =8

%fs = 3072000 fpass = 10000 fstop = 768000 astop = 80= upsample (r2,2);(1: length (k2) - 1) = [];= zeros (1,length (B1));i = length (k2): length (r3)n = 0: length (k2) - 1(i) = r3 (i) +B1 (i-n) *k2 (n+1);(1: length (k2) - 1) = [];

%CIC= 64; % filter buffer depth= upsample (r3, 32);= zeros (1,D);= zeros (1,D);= zeros (1,D);= zeros (1,D);= 0;= 0;= 0;= 0;ii = 1: (length (B2) /10

%comb1= B2 (ii) - delayBuffer1 (end);(2: end) = delayBuffer1 (1: end-1);(1) = B2 (ii);

% integrator1= intOut1 + combOut1;(ii) = intOut1;

%comb2= r4 (ii) - delayBuffer2 (end);(2: end) = delayBuffer2 (1: end-1);(1) = r4 (ii);

% integrator2= intOut1 + combOut2;(ii) = intOut2;

%comb3= r5 (ii) - delayBuffer3 (end);(2: end) = delayBuffer3 (1: end-1);(1) = r5 (ii);

% integrator3= intOut2 + combOut3;(ii) = intOut3;

%comb4= r6 (ii) - delayBuffer4 (end);(2: end) = delayBuffer4 (1: end-1);(1) = r6 (ii);

% integrator3= intOut3 + combOut4;7 (ii) = intOut4;

end

%Спектр выходного сигнала

figure (3)= psd (Hs1,r7,'Fs',98304000); plot (h1)('Power Spectrum of a Sine Wave'); % Вывод названия графика.

xlabel ('Frequency (Hz) '); % Вывод величины и размерности по x.('Power'); % Вывод величины по y.

Спектр исходного сигнала (рис. 19).

Рис. 19

Спектр сигнала после всех фильтров (рис. 20)

Рис. 20

Перейти на страницу: 1 2 3 

Еще статьи по технике и технологиям

Расчет электрофизических характеристик структуры метал-диэлектрик-полупроводник
Основной задачей курсовой работы по дисциплине «Физика полупроводников и диэлектриков» является выяснение физической сущности явлений и процессов, лежащих в основе работы различных полупроводниковых структур и дискретных полупроводнико ...

Характеристика сетей и технологий Integrated Services Digital Network
В данной работе рассматривается тема "Характеристика сетей и технологий ISDN". ISDN - цифровая сеть комплексного обслуживания. В соответствии с определением МСЭ ISDN - это сеть, обеспечивающая полностью цифровые соединения м ...

© 2019 | www.techexpose.ru