Сигнальные процессоры. Аксенов В.П. - 126 стр.

z22(m) = Fs*acos(z1)/(2*pi); % расчет частоты выходного сигнала
z3=sqrt(abs(s1(m+1)^2-s1(m)*s1(m+2)));
 if abs(z22(m)) > 0.02 z34(m)=z3*Fs/(2*pi*z22(m));
 else z34(m) = z34(m-1) ;% расчет амплитуды выходного сигнала
end
end
m=1:maxt;
subplot(1,2,2)
plot(m,z34)              % график мгновенной амплитуды А_ВЫХ (m)
grid on                      % время t = m/FS
title ('амплитуда выходного сигнала от m');


                                                                      ПРИЛОЖЕНИЕ 9

              Синтез КИХ-фильтра НЧ при фазовой манипуляции


%         КИХ-фильтр нижних частот, фазовая манипуляция
%         окно Кайзера
clear all; % уничтожаем все переменные
A = 5; Fs =32000; V = 600; % амплитуда, частота дискрет - ии, скорость
Rp = 1.5; Rs = 25; % пульсации и минимальное ослабление в децибелах
Kp =(10^(Rp/10)-1)^0.5; Km = 10^(-Rs/20); % пересчет Rp, Rs в абсолют-ные значения
f1 = 4.0*V; f2 = f1; % f1 - частота входного сигнала
f3 = 1.20*f2; f4=1.65*f2;         % точки АЧХ фильтра НЧ
[n,Wn,beta,ftype] = kaiserord ([f3,f4], [1 0], [Kp,Km], Fs);
% определен порядок n фильтра с окном Кайзера
b = fir1(n, Wn, ftype, kaiser (n+1, beta), 'noscale');   % расчет коэфф-ов
[h,f]=freqz(b);                % расчет АЧХ фильтра
subplot(1,2,1) % 1-ый график АЧХ с автоматическим выбором масштаба
plot(f/pi,20*log10(abs(h))) % график АЧХ в логарифмическом масштабе
grid on           % включение сетки на графике
title ('АЧХ в децибелах');
subplot(1,2,2) %                 2-ой график АЧХ для определения пульсаций
plot(f/pi,20*log10(abs(h))) % график АЧХ в логарифмическом масштабе
ylim([-Rp Rp])      % диапазон Rp в децибелах


                                                126