clear; clc; clf; close, hold off;



Am=1;
Ac=1;
fm=1500;
fc=40000;
prob=16;
f1=1000;
k=10;

fp=prob*fc;
N=fp/500;
n=0:N-1;
Tp=1/fp;
t=n*Tp;
f=n/Tp/N;


ss=0.2*sin(2*pi*f1*t)+0.5*sin(2*pi*fm*t);
m=ss/Am;
s=Ac*cos(2*pi*fc*t).*(1+k.*m);

znak=(sign(s)+1)/2;
o=s.*znak;

r=2*fm/fp;
[B,A]=butter(6,r);

y=filter(B,A,o);

widmo=fft(o);
wa=abs(widmo); ef=angle(widmo);

subplot(512);plot(f,s);title('Sygnal zmodulowany');
subplot(513);plot(f,o);title('Wyprostowany sygnal zmodulowany');
subplot(514 );plot(n,y);title('Sygnal po demodulacji');
subplot(511);plot(t,ss);title('Sygnal przed modulacja');
subplot(515);stem(f,wa);title('Widmo wyprostowanego sygnalu zmodulowanego');