t=0:0.001:18;
s=sin (2*pi*0.25*t);
thetaFM= -Df*(1/(2*pi*0.25))*cos(2*pi*0.25*t);
thetaPM = Dp*s;
fc = 1;
A = 1;
Dp = pi;
Df = 2*pi;

port = A*cos(2*pi*fc*t);
FM = A*cos(2*pi*fc*t + thetaFM);
PM = A*cos(2*pi*fc*t + thetaPM);

subplot(4,1,1);
plot(t, s); title('Información'); xlabel('Tiempo'); ylabel('Amplitud');

subplot(4,1,2);
plot(t, port); title('Portadora'); xlabel('Tiempo'); ylabel('Amplitud');

subplot(4,1,3);
plot(t, FM); title('FM'); xlabel('Tiempo'); ylabel('Amplitud');

subplot(4,1,4);
plot(t, PM); title('PM'); xlabel('Tiempo'); ylabel('Amplitud');