通信原理实验模拟调制系统(AM,FM)实现方法

%对已调信号u求傅里叶变换 fs=1/ts;

if nargin==2 n1=0; else

n1=fs/df; end

n2=length(sfm);

n=2^(max(nextpow2(n1),nextpow2(n2))); U=fft(sfm,n);

u=[sfm,zeros(1,n-n2)];

df1=fs/n;

U=U/fs; %缩放

figure(2)

subplot(2,1,1)

plot(f,abs(fftshift(M))) ?tshift:将FFT中的DC分量移到频谱中心 xlabel('频率f')

title('原调制信号的频谱图') subplot(2,1,2)

plot(f,abs(fftshift(U))) xlabel('频率f')

title('已调信号的频谱图')

figure(3)

subplot(3,1,1);plot(t,mt); %绘制调制信号的时域图 title('调制信号的时域波形'); subplot(3,1,2);plot(t,sfm); %绘制已调信号的时域图 title('已调信号的时域波形');

nsfm=sfm;

for i=1:length(t)-1 %接受信号通过微分器处理 diff_sfm(i)=(sfm(i+1)-sfm(i))./dt; end

diff_sfmn = abs(hilbert(diff_sfm)); %hilbert变换，求绝对值得到瞬时幅度（包络检波） zero=(max(diff_sfm)-min(diff_sfm))/2; diff_sfm1=diff_sfm-zero;

subplot(3,1,3); %绘制解调信号的时域图 plot((1:length(diff_sfm1))./1000,diff_sfm1./400,'r'); title('解调信号的时域波形');

生成图像如下：

－ 13 －

放大后可以看到，调制信号频谱在6HZ处有冲击，符合要求；已调信号频谱宽度较宽，由BFM=2*（mf+1）*fm=2*(10+1)*6=132HZ,由图中可见，理论与实验结果一致。

－ 14 －

从图中可以看出，解调信号的包络与原调制信号是一致的。 五、实验心得

这次实验主要做了关于标准调幅波信号AM，抑制载波的双边带信号DSB，单边带信号SSB以及调频FM的仿真实验。通过实验，自己进一步掌握了AM,DSB,SSB及FM的联系与区别，理解了课本上的理论知识。通过实验的形式将理论知识加以应用，得出的实验结果又进一步论证了理论的正确性。

－ 15 －