clcFs = 16000; % Sampling rate of signalFc = 4000; % Carrier frequencyFt = 100;t = [0:2000]'/Fs; % Sampling timess = sin(2*pi*Ft*t);%調制信號a = modulate(s,Fc,Fs,'fm');%調制1dev = 500; % Frequency deviation in modulated signalb = fmmod(s,Fc,Fs,dev); %fmmod 調制2z = fmdemod(b,Fc,Fs,dev); %fmdemod 解調subplot(2,1,1);plot(z);subplot(2,1,2);plot(z);zz = fft(b, 1024);%fftf = 1000*(0:512)/1024;yy = zz.* conj(zz) / 1024;plot(f,yy(1:513))