fmmode.m

fm调制和解调

clc
Fs = 16000; % Sampling rate of signal
Fc = 4000; % Carrier frequency
Ft = 100;
t = [0:2000]'/Fs; % Sampling times
s = sin(2*pi*Ft*t);%调制信号

a = modulate(s,Fc,Fs,'fm');%调制1

dev = 500; % Frequency deviation in modulated signal

b = fmmod(s,Fc,Fs,dev); %fmmod 调制2

z = fmdemod(b,Fc,Fs,dev); %fmdemod 解调

subplot(2,1,1);
plot(z);
subplot(2,1,2);
plot(z);
zz = fft(b, 1024);%fft
f = 1000*(0:512)/1024;
yy = zz.* conj(zz) / 1024;
plot(f,yy(1:513))