fm2.asv

这是关于模拟调制的MATLAB源代码

% fm2.m
% 频率调制的Matlab源程序。
% 被调信号m(t)=sinc(100t)
t0=0.2;
ts=0.001; 
fc=250; 
kf=100;
fs=1/ts;  
df=0.3;   
t=[-t0/2:ts:t0/2];
%% 信息信号
m=sinc(100*t);
int_m(1)=0;
for i=1:length(t)-1  % m的积分
    int_m(i+1)=int_m(i)+m(i)*ts;
end
[M,m,df1]=fftseq(m,ts,df);
M=M/fs;
f=[0:df1:df1*(length(m)-1)]-fs/2; 
u=cos(2*pi*fc*t+2*pi*kf*int_m); % 调制信号
[U,u,df1]=fftseq(u,ts,df);
U=U/fs;
[v,phase]=env_phas(u,ts,250);
phi=unwrap(phase); %恢复原相位
dem=(1/(2*pi*kf))*(diff(phi)/ts);
%% 画图
subplot(3,2,1);
plot(t,m(1:length(t)))
axis([0 0.15 -2.1 2.1])
title('被调信号m(t)')
subplot(3,2,2);
plot(t,u(1:length(t)))
axis([0 0.15 -2.1 2.1])
title('调制信号u(t)')
subplot(3,2,3);
plot(f,abs(fftshift(M)))
title('被调信号M(f)')
subplot(3,2,4);
plot(f,abs(fftshift(U)))
title('调制信号U(f)')
subplot(3,2,5);
plot(t,m(1:length(t)))
axis([0 0.15 -2.1 2.1])
title('被调信号m(t)')
subplot(3,2,6);
plot(t,dem(1:length(t)))
axis([0 0.15 -2.1 2.1])
title('调制信号dem(t)')