调频FM仿真程序Word文档格式.docx

1、 %设定调制信号频率mt=am*cos（2*pi*fm*t）; %生成调制信号fc=50; %设定载波频率ct=cos（2*pi*fc*t）; %生成载波kf=10; %设定调频指数int_mt（1）=0;for i=1:length（t）-1 int_mt（i+1）=int_mt（i）+mt（i）*dt; %求信号m（t）的积分end %调制，产生已调信号sfm=am*cos（2*pi*fc*t+2*pi*kf*int_mt）; %已调信号disp（按任意键可以看到原调制信号、载波信号和已调信号的曲线）pause%*figure（1）*figure（1）subplot（3,1,1）;plot

2、（t,mt）; %绘制调制信号的时域图xlabel（时间t）;title（调制信号的时域图subplot（3,1,2）;plot（t,ct）; %绘制载波的时域图载波的时域图subplot（3,1,3）;plot（t,sfm）; %绘制已调信号的时域图已调信号的时域图%*二：%添加高斯白噪声sn1=10; %设定信躁比（小信噪比）sn2=30; %设定信躁比（大信噪比）sn=0; %设定信躁比（无信噪比）db=am2/（2*（10（sn/10）; %计算对应的高斯白躁声的方差n=sqrt（db）*randn（size（t）; %生成高斯白躁声nsfm=n+sfm; %生成含高斯白躁声的已调信号

3、（信号通%过信道传输）%FM解调length（t）-1 %接受信号通过微分器处理diff_nsfm（i）=（nsfm（i+1）-nsfm（i）./dt;enddiff_nsfmn = abs（hilbert（diff_nsfm）; %hilbert变换，求绝对值得到瞬时%幅度（包络检波）zero=（max（diff_nsfmn）-min（diff_nsfmn）/2;diff_nsfmn1=diff_nsfmn-zero;%时域到频域转换ts=0.001; %抽样间隔fs=1/ts; %抽样频率df=0.25; %所需的频率分辨率，用在求傅里叶变换%时，它表示FFT的最小频率间隔%对调制信号m（

4、t）求傅里叶变换m=am*cos（2*pi*fm*t）; %原调信号if nargin=2n1=0;elsen1=fs/df;n2=length（m）;n=2（max（nextpow2（n1）,nextpow2（n2）;M=fft（m,n）;m=m,zeros（1,n-n2）;df1=fs/n; %以上程序是对调制后的信号u求傅里变换M=M/fs; %缩放，便于在频谱图上整体观察f=0:df1:df1*（length（m）-1）-fs/2; %时间向量对应的频率向量%对已调信号u求傅里变换n2=length（sfm）;U=fft（sfm,n）;u=sfm,zeros（1,n-n2）; %以上是

5、对已调信号u求傅里变换U=U/fs; %缩放按任意键可以看到原调制信号和已调信号在频域内的图形%*figure（2）* figure（2）subplot（2,1,1）plot（f,abs（fftshift（M） %fftshift:将FFT中的DC分量移到频谱中心频率f原调制信号的频谱图subplot（2,1,2）plot（f,abs（fftshift（U）已调信号的频谱图三：按任意键可以看到原调制信号、无噪声条件下已调信号和解调信号的曲线%*figure（3）* figure（3） %绘制已调信号的时域图无噪声条件下已调信号的时域图nsfm=sfm; diff_nsfm（i）=（nsfm（i

6、+1）-nsfm（i）./dt; %hilbert变换，求绝对值得到%瞬时幅度（包络检波） %绘制无噪声条件下解调信号的时域图plot（1:length（diff_nsfmn1）./1000,diff_nsfmn1./400,r无噪声条件下解调信号的时域图四：按任意键可以看到原调制信号、小信噪比高斯白噪声条件下已调信号和解调信号已调信号的曲线%*figure（4）*figure（4）db1=am2/（2*（10（sn1/10）; %计算对应的小信噪比高斯白躁声的方差n1=sqrt（db1）*randn（size（t）;nsfm1=n1+sfm; diff_nsfm1（i）=（nsfm1（i+1

7、）-nsfm1（i）./dt;diff_nsfmn1 = abs（hilbert（diff_nsfm1）;diff_nsfmn1=diff_nsfmn1-zero;plot（1:length（diff_nsfm）,diff_nsfm）; %绘制含小信噪比高斯白噪声已%调信号的时域图含小信噪比高斯白噪声已调信号的时域图 %绘制含小信噪比高斯白噪声解调信号的时域图含小信噪比高斯白噪声解调信号的时域图%*五：按任意键可以看到原调制信号、大信噪比高斯白噪声条件下已调信号和解调信号已调信号的曲线%*figure（5）*figure（5）db1=am2/（2*（10（sn2/10）; %计算对应的大信噪比高斯白躁声的方差 %生成含高斯白躁声的已调信号（信号通过信道传输） %hilbert变换，求绝对值%得到瞬时幅度（包络检波）length（diff_nsfm1）,diff_nsfm1）; %绘制含大信噪比高斯%白噪声已调信号的时域图含大信噪比高斯白噪声已调信号的时域图 %绘制含大信噪比高斯白噪声解调信号%的时域图含大信噪比高斯白噪声解调信号的时域图