华北电力大学通信系统仿真综合实验报告.doc

综合实验报告 ( 2015--2016 年度 第 一 学期) 名 称： 通信系统仿真 题 目： FM的调制与解调 院 系： 信息工程系 班 级： 学 号： 学生姓名： 指导教师： 孙景芳、王雅宁 设计周数： 1 成 绩： 日期： 2015 年 11 月 10 日实验名称 实验一：1-7 实验环境 Matlab仿真平台 实 验 目 的 1、能够熟练掌握和综合运用通信领域中的基本理论和专业知识； 2、能够完成通信基本理论的仿真FM信号的频谱有如下特点： （1）以载频为中心，由无穷多对以调制信号频率F为间隔的边频分量组成，各分量幅值取决于Bessel函数，且以对称分布； （2）载波分量并不总是最大，有时为零； （3）FM信号的功率大部分集中在载频附近； （4）频谱结构与F密切相关； 调频波解调又称鉴频，其中一种方法为将输入调频信号进行特定波形变换，使变换后波形包含反映瞬时频率变化的平均分量。然后通过低通滤波器就能输出所需的解调电压。 加入高斯白噪声后的已调波会随输入噪声信噪比的不同呈现不同的情况。 当信噪比比较大时，频率调制系统的调制制度增益很高，即抗噪性能好。 当信噪比较小时，频率调制系统的调制制度增益较低，当（S/N）低于一定数值时，解调器的输出信噪比就会急剧恶化，这种现象称为解调信号解调的门限效应。 实 验 步 骤 function[]=jietiao() t0=0.2; tz=0.0001; %设定时间步长 fz=1/tz; %设定抽样频率 t=-t0:tz:t0; %产生时间向量 kf=100; %设定调频指数 fc=1000; %设定载波频率 kd=0.5; %设定鉴频增益/鉴频器灵敏度 df=0.05; %设定分辨率 m_fun=cos(400*pi*t); int_m=zeros(1,length(t)-1); int_m(1)=0; %对m_fun 积分 for i=1:length(t)-1 int_m(i+1)=int_m(i)+m_fun(i)*tz; end x=sin(2*pi*fc*t+2*pi*kf*int_m); %调制信号 noisycode=awgn(x,2);%加入噪声 u9=m_fun+noisycode; y=u9.*kf; %解调信号 [M,m_fun,df1]=fftseq(m_fun,tz,df); %对原始信号快速傅里叶变换 M=M/fz; f=[0:df1:df1*(length(m_fun)-1)]-fz/2; [X,x,df1]=fftseq(x,tz,df); %对已调信号快速傅里叶变换 X=X/fz; [Y,y,df1]=fftseq(y,tz,df); %对解调信号快速傅里叶变换 Y=Y/fz; figure(1); %生成原始信号的时域图形 plot(t,m_fun(1:length(t)),linewidth,3); axis([-0.01 0.01 -1.5 1.5]); title(原始信号的时域图形); xlabel(时间); legend(m(t)) figure(2); %生成原始信号的频域图形 plot(f,abs(fftshift(M)),linewidth,3); axis([-400 400 -0.01 0.1]); title(原始信号的频域图形); xlabel(频率); legend(M(f)); figure(3); %生成已调信号的时域图形 plot(t,x(1:length(t)),linewidth,3); axis([-0.015 0.015 -1.5 1.5]); title(已调信号的时域图形); xlabel(时间); legend(x(t)); 实验步骤 figure(4); %生成已调信号的频域图形 plot(f,abs(fftshift(X)),linewidth,3); axis([-1500 1500 0 0.1]); title(已调信号的频域图形); xlabel(频率); legend(X(f)); figure(5); %生成解调信号的时域图形 plot(t,y(1:length(t)),linewidth,3); axis([-0.01 0.01 -100 100]); title