z = pmdemod(y,Fc,Fs,phasedev)
z = pmdemod(y,Fc,Fs,phasedev,ini_phase)
z = pmdemod(y,Fc,Fs,phasedev) demodulates
the phase-modulated signal
y at the carrier frequency
the carrier signal have sampling rate
Fs must be at least
phasedev argument is the phase deviation of
the modulated signal, in radians.
z = pmdemod(y,Fc,Fs,phasedev,ini_phase) specifies
the initial phase of the modulated signal, in radians.
Set the sample rate. To plot the signals, create a time vector.
fs = 50; t = (0:2*fs+1)'/fs;
Create a sinusoidal input signal.
x = sin(2*pi*t) + sin(4*pi*t);
Set the carrier frequency and phase deviation.
fc = 10; phasedev = pi/2;
Modulate the input signal.
tx = pmmod(x,fc,fs,phasedev);
Pass the signal through an AWGN channel.
rx = awgn(tx,10,'measured');
Demodulate the noisy signal.
y = pmdemod(rx,fc,fs,phasedev);
Plot the original and recovered signals.
figure; plot(t,[x y]); legend('Original signal','Recovered signal'); xlabel('Time (s)') ylabel('Amplitude (V)')