Амплитудный ограничитель

>> f0=1e3;
>> fs=100e3;
>> fd=2e6;
>> t=0:1/fd:0.01;
>> x=cos(2*pi*f0*t); % сигнал
>> dw=fs/fd*2*pi/max(max(x));%задаём интенсивность флуктуаций мгновенной чатоты(при девиации частоты 100кГц
>> dm=0.03*dw;%Задаём частоту девиации, равную 3кГц
>> [y,t]=modulate(x,fs,fd,'fm', dm);%модулируем сигнал
>> n=wgn(1,length(t),10,50,'dBm','real');%создаём белый шум
>> kn=0.1*n;% делаем отношение с/ш=100
>> sign =y+kn; %   Добавляем к модулированному сигналу шумы
yf=filter(PF30e3,y);%пропускаем сигнал через фильтр входной цепи
fg=555e3; %частота гетеродина
>> sg=cos(2*pi*fg*t); %сигнал гетеродина
sm=sg.*yf;%смеситель
>> fsm=filter(PF45510e3,sm);%пропускаем сигнал через ПФ с центральной частотой 455кГц
d=demod(fsm,455e3,fd,'fm',dm);%демодулируем сигнал
>> lpf=filter(LPF5e3, d);%пропускаем через ФНЧ с частотой среза 5кГц
subplot(411);plot(t, y);subplot(412);plot(t, yf);subplot(413);plot(t, fsm); subplot(414);plot(t, lpf)