基于MATLAB的FM调制实现.pdf

《基于MATLAB的FM调制实现.pdf》由会员分享，可在线阅读，更多相关《基于MATLAB的FM调制实现.pdf（11页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、FM MATLAB 建模与仿真 一、FM 调制原理：FM 属于角度调制，角度调制与线性调制不同，已调信号频谱不再是原调制信号频谱的线性搬移，而是频谱的非线性变换，会产生与频谱搬移不同的新的频率成分，故又称为非线性调制。在本实验中使用正弦信号作为基带信号进行调制的分析.频率调制的一般表达式1为：FMs FM 调制是相位偏移随（）的积分呈线性变化。FM 调制模型的建立 图 1 FM 调制模型 其中，()m t为基带调制信号，设调制信号为()cos(2)mmtAf t 设正弦载波为()cos(2)cc tf t 信号传输信道为高斯白噪声信道，其功率为2。图 2 总体模型 二 调制过程的分析：在调制时

2、，调制信号的频率去控制载波的频率的变化，载波的瞬时频偏随调制信号()m t成正比例变化，即 发送m(t)信息源)(tSFM)(tSFM)(0tmFM 调制器 信道 FM 解调器 接收 FM MATLAB 建模与仿真 ()()fdtK mtdt 式中，fK为调频灵敏度（()rads V）。这时相位偏移为()()ftKmd 则可得到调频信号为()cos()FMcfstAt K m d FM 调制 1.对 FM 调制信号的频谱分析 clear all ts=0.00125;%信号抽样时间间隔 t=0:ts:10-ts;%时间向量 am=10;fs=1/ts;%抽样频率 df=fs/length(t)

3、;%fft的频率分辨率 msg=am*cos(2*pi*10*0:0.01:0.99);msg1=msg*ones(1,fs/10);%扩展成取样信号形式 msg2=reshape(msg1.,1,length(t);Pm=fft(msg2);%求消息信号的频谱 f=-fs/2:df:fs/2-df;subplot(3,1,1)plot(t,fft(abs(Pm)title(消息信号频谱)m=fft(msg,1024);%对 msg进行傅利叶变换 N=(0:length(m)-1)*fs/length(m)-fs/2;subplot(3,1,2)plot(N,abs(m);%调制信号频谱图 t

4、itle(调制信号频谱)int_msg(1)=0;%消息信号积分 for ii=1:length(t)-1 int_msg(ii+1)=int_msg(ii)+msg2(ii)*ts;end kf=50;fc=250;%载波频率 Sfm=am*cos(2*pi*fc*t+2*pi*kf*int_msg);%调频信号 FM MATLAB 建模与仿真 Pfm=fft(Sfm)/fs;%FM信号频谱 subplot(3,1,3);plot(f,fftshift(abs(Pfm)%画出已调信号频谱 title(FM信号频谱)Pc=sum(abs(Sfm).2)/length(Sfm)%已调信号功率 P

5、s=sum(abs(msg2).2)/length(msg2)%消息信号功率 fm=50;betaf=kf*max(msg)/fm%调制指数 W=2*(betaf+1)*fm%调制信号带宽 用 FFT函数进行傅利叶变换，进行傅立叶变化便分别得到调制信号与调制之后的 FM信号的频谱图如下：012345678910-202x 105消 息 信 号 频 谱-500-400-300-200-10001002003004005000200400600调 制 信 号 频 谱-500-400-300-200-10001002003004005000510FM信 号 频 谱 图 2-5 通过频谱图的对照比较我

6、们可以看出FM调制并不是使原正弦信号的频谱在原来位置上通过移动得到调制波形，调制后的波形与调制前的完全不同，这证明 FM调制并不是线性的，而是非线性的。由图像可知fft变换出来的图形是左右FM MATLAB 建模与仿真 对称的。2.改变采样频率对调制信号进行分析 采样频率分别为100HZ，200HZ，400HZ，800HZ 时的时域信号波形 00.10.20.30.40.50.60.70.80.91-1-0.500.51采 样 频 率 为 100HZ 时，时 域 信 号 的 波 形t(s)时域信号幅度 00.10.20.30.40.50.60.70.80.91-1-0.500.51采 样 频

7、率 为 200HZ时，时 域 信 号 的 波 形t(s)时域信号幅度 00.10.20.30.40.50.60.70.80.91-1-0.500.51采 样 频 率 为 400HZ 时，时 域 信 号 的 波 形t(s)时域信号幅度 FM MATLAB 建模与仿真 00.10.20.30.40.50.60.70.80.91-1-0.500.51采 样 频 率 为 800HZ 时，时 域 信 号 的 波 形t(s)时域信号幅度 采样频率分别为100HZ，200HZ，400HZ，800HZ 时的频域信号波形-50-40-30-20-100102030405001020304050100Hz 采 样

8、 频 率 FM信 号 的 频 域 频 谱频 率 分 量进过FM调制信号功率谱-100-80-60-40-2002040608010001020304050200Hz采 样 频 率 FM信 号 的 频 域 频 谱频 率 分 量进过FM调制信号功率谱 FM MATLAB 建模与仿真 -200-150-100-5005010015020001020304050400Hz采 样 频 率 FM信 号 的 频 域 频 谱频 率 分 量进过FM调制信号功率谱-400-300-200-100010020030040001020304050800Hz采 样 频 率 FM信 号 的 频 域 频 谱频 率 分 量进

9、过FM调制信号功率谱 载波频率不变，改变采样频率后，如果将采样频率取的太小，到小于两倍的载波频率，MATLAB 软件将不能输出波形，说明采样频率只有在大于两倍载波 在频域内，信号的能量主要集中在采样频率为|100|HZ 的范围内，可以看出采用等位置采样方法进行周期采样,频率随机波动的大小不会产生谱分析误差，当频率随机波动水平较大时,不仅会导致功率谱峰值的显著降低,。3加入噪声，绘制出加入噪声后的FM 信号的时频图，并对其进行分析噪声对信号的影响。设其载波频率fc=150;采样频率fs=800;调制信号x=sin(2*pi*30*t);FM 信号y=modulate(x,fc,fs,FM);利用

10、y1=y+awgn(y,1,0);加入噪声。得到其时域波形图如下：FM MATLAB 建模与仿真 FM 信号加入噪声的时域波形图 通过b=fft(y1,1024)；f=(0:length(b)-1)*fs/length(b)-fs/2;对 y 进行傅利叶变换，得到加入噪声后的FM 信号频谱为：FM 信号加入噪声的频谱图 通过对函数y1=y+awgn(y,10,0)改变信噪比后，信噪比由1 变为10，得到信噪比为10 的时域图为：信噪比改为10 的时域图 频谱图为：FM MATLAB 建模与仿真 信噪比改为10 的频域图 通过对函数y1=y+awgn(y,20,0)改变信噪比后，信噪比由10 变

11、为20，得到信噪比为20 的时域图为：信噪比改为20 的时域图 频域图为：信噪比改为20 的频域图 加入噪声后时域波形与原来的时域波形相比，波形明显失真，波形不仅不如FM MATLAB 建模与仿真 原本波形般规则，而且曲线之间还出现了为链接在一起的断裂，但随着信噪比的增大，与原有的波形的相似度也增大了，说信噪比越大，噪声对信号的影响也变小了。从加入噪声的图形与未加入噪声的对比中我们还可以看出噪声对时域图的变化明显比频域图的变化更为突出，白噪声在整个频谱内每个频点的能量为常数，且基本恒定，所以他对于时域的影响更大。通过对噪声的理解我们可以知道对于调频系统来说，增加传输带宽就可以改善抗噪声性能。调

12、频方式的这种以带宽取信噪比的特性是十分有益的 4 分析信噪声比大于多少的情况下，噪声会完全淹没在信号中，在接收端，解调出来的信号和没有加入噪声时传输到接收端的信号几乎完全一致 设置信噪比在10db，15db，25db，30db，40db，50db，80db-1000-800-600-400-20002004006008001000-20020原 始 消 息 信 号-1000-800-600-400-20002004006008001000-20020无 噪 声 时 的 解 调 信 号-1000-800-600-400-20002004006008001000-20020信 噪 比 为 10dB

13、时 的 解 调 信 号-1000-800-600-400-20002004006008001000-20020原 始 消 息 信 号-1000-800-600-400-20002004006008001000-20020无 噪 声 时 的 解 调 信 号-1000-800-600-400-20002004006008001000-20020信 噪 比 为 15dB时 的 解 调 信 号 FM MATLAB 建模与仿真 -1000-800-600-400-20002004006008001000-20020原 始 消 息 信 号-1000-800-600-400-2000200400600800

14、1000-20020无 噪 声 时 的 解 调 信 号-1000-800-600-400-20002004006008001000-20020信 噪 比 为 25dB时 的 解 调 信 号 012345678910-20020原 始 消 息 信 号012345678910-20020无 噪 声 时 的 解 调 信 号012345678910-20020信 噪 比 为 30dB时 的 解 调 信 号-1000-800-600-400-20002004006008001000-20020原 始 消 息 信 号-1000-800-600-400-20002004006008001000-20020无

15、 噪 声 时 的 解 调 信 号-1000-800-600-400-20002004006008001000-20020信 噪 比 为 40dB时 的 解 调 信 号 FM MATLAB 建模与仿真 -1000-800-600-400-20002004006008001000-20020原 始 消 息 信 号-1000-800-600-400-20002004006008001000-20020无 噪 声 时 的 解 调 信 号-1000-800-600-400-20002004006008001000-20020信 噪 比 为 50dB时 的 解 调 信 号-1000-800-600-400-20002004006008001000-20020原 始 消 息 信 号-1000-800-600-400-20002004006008001000-20020无 噪 声 时 的 解 调 信 号-1000-800-600-400-20002004006008001000-20020信 噪 比 为 80dB时 的 解 调 信 号 结论；随着信噪比的增加，在接收端解调出来的信号，失真会越来越小，并且在信噪比大于40db 以上，接收端接收到的信号几乎没有失真