matlab 频率调制

题目

代码

clear all;clc;
%FM调制与解调
t0=0.15;
ts=0.001;
Fc=100;
FM=50;
Fs=1/ts;
df=0.3;
t=[0:ts:t0];
m=[ones(1,t0/(3*ts)),2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)];%定义信号序列
%%%%%%%%%%%方法一%%%%%%%
% int_m(1)=0; %积分
% for i=1:length(t)-1
%    int_m(i+1)=int_m(i)+m(i)*ts;
% end
% int_m= cumsum(m)/Fs; %积分 
% u=cos(2*pi*Fc*t+2*pi*FM*int_m); %调制信号
% y=m.*FM; %解调信号
%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%方法二%%%%%%%
u=fmmod(m,Fc,Fs,FM);
y=fmdemod(u,Fc,Fs,FM);
%%%%%%%%%%%%%%%%%%%%%%%%%
[M,m,df1]=fft_seq(m,ts,df); %对原信号傅里叶变换
M=M/Fs;
[U,u,df1]=fft_seq(u,ts,df); %对已调信号傅里叶变换
U=U/Fs;
[Y,y,df1]=fft_seq(y,ts,df); %对解调信号傅里叶变换
Y=Y/Fs;
f=[0:df1:df1*(length(m)-1)] -Fs/2;

subplot(3,2,1);
plot(t,m(1:length(t)));
title('FM未调制信号');

subplot(3,2,2);
plot(t,u(1:length(t)));
title('FM已调制信号');

subplot(3,2,3);
plot(t,y(1:length(t)));
title('FM解调信号');

subplot(3,2,4);
plot(f,abs(fftshift(Y)));
title('FM解调信号频谱');

subplot(3,2,5);
plot(f,abs(fftshift(M)));
title('FM未调制信号频谱');

subplot(3,2,6);
plot(f,abs(fftshift(U))); 
title('FM已调制信号频谱');

效果