基于MATLAB的模拟调制实验报告.doc

基于matlab的模拟调制实验报告1、 实验目的 1.进一步学习调制的知识，掌握调频与调角两种模拟调制技术。 2.进一步学习matlab的编程，熟练使用matlab进行作图。2、 实验原理1.调制的概念 调制（modulation）就是对信号源的信息进行处理加到载波上，使其变为适 合 于信道传输的形式的过程，是使载波随信号而改变的技术。一般，用来传送消息的信号叫作载波或受调信号，代表所欲传送消息的信号叫作调制信号，调制后的信号叫作已调信号。用调制信号控制载波的某些参数，使之随而变化，就可实现调制。2.调制的目的 频谱变换 当所要传送的信号的频率或者太低，或者频带很宽，对直接采用电磁波的形 式进行发送很不利，需要的天线尺寸很大，而且发射和接受短的天线与谐振回路的参数变化范围很大。为了信息有效与可靠传输，往往需要将低频信号的基带频谱搬移到适当的或指定的频段。这样可以提高传输性能，以较小的发送功率与较短的天线来辐射电磁波。 实现信道复用 为了使多个用户的信号共同利用同一个有较大带宽的信道，可以采用各种复用技术。如模拟电话长途传输是通过利用不同频率的载波进行调制。将各用户话音每隔4 khz搬移到高频段进行传输。 提高抗干扰能力不同的调制方式，在提高传输的有效性和可靠性方面各有优势。如调频广播系统，它采用的频率调制技术，付出多倍带宽的代价，由于抗干扰性能强，其音质比只占10 khz带宽的调幅广播要好得多。扩频通信就是以大大扩展信号传输带宽，以达到有效抗拒外部干扰和短波信道多径衰落的特殊调制方式。3.调制的种类 根据和的不同类型和完成调制功能的调制器传递函数不同，调制分为以下多种方式：（1）按调制信号的类型分为：l 模拟调制：调制信号是连续变化的模拟量，如话音与图像信号。l 数字调制：调制信号是数字化编码符号或脉冲编码波形。（2）按载波信号的类型分：l 连续波调制：载波信号为连续波形，通常以正弦波作为载波。l 脉冲调制：载波信号是脉冲波形序列。（3）按调制器的不同功能分：l 幅度调制：以调制信号去控制载波的幅度变化，如模拟调幅，脉冲幅度调制（pam），幅移键控（ask）。l 频率调制：以调制信号去控制载波信号的频率变化，如模拟调频（fm），频移键控（fsk），脉宽调制（pdm）。l 相位调制：以调制信号去控制载波信号的相位变化，如模拟调相（pm），相移键控（psk），脉位调制（ppm）。（4）按调制器的传输函数分：l 线性调制：已调信号的频谱与调制信号频谱是线性的频谱位移关系。如各种幅度调制，幅移键控（ask）。l 非线性调制：已调信号的频谱与调制信号频谱没有线性关系，即调制后派生出大量不同于调制信号的新的频率成份。如调频（fm），调相（pm），频移键控（fsk）。3、 实验过程 这次实验主要利用matlab实现幅度调制与角度调制。 1.幅度调制原理 am调制 am是用调制信号去控制高频正弦载波的幅度，使其按调制信号的规律变化的过程。如图1所示。 图1 调制信号叠加直流后再与载波相乘，则输出的信号就是常规双边带（am）调幅波。其时域表达式为： am信号波形的包络与输入基带信号成正比，故用包络检波的方法很容易恢复原始调制信号。但为了保证包络检波时不发生失真，必须满足，否则将出现过调幅现象而带来失真。am信号的频谱是由载频分量和上、下两个边带组成。上边带的频谱与原调制信号的频谱结构相同，下边带是上边带的镜像。am信号是带有载波的双边带信号，它的带宽信号带宽的两倍。下面利用matlab实现am的调制。为方便观察，令输入信号为单一正弦波信号。 -am调制输出信号和频谱- dt=0.0001; %时间采样频率 fh=10; %调制信号最高频率 fc=1000; %载波中心频率 t=0.2; %信号时长 n=t/dt; t=0:n-1*dt; mt=sqrt(2)*cos(2*pi*fh*t); %调制信号 a=2; s_am=(a+mt).*cos(2*pi*fc*t); f,xf=fft_shift(t,s_am); %已调信号频谱 subplot(211); plot(t,s_am);hold on; %画出am信号波形 plot(t,a+mt,r-); %表示am上包络 plot(t,-(a+mt),r-); %表示am下包络 title(am调制信号及其包络); xlabel(时间t); subplot(212); %画出am波频谱 plot(f,abs(xf); axis(700 1300 0 max(xf); title(am信号功率谱); xlabel(频率f); 图2结果分析：从图2中频谱图可以看出，调幅过程实际就是频谱搬移，是一个线性过程。经调制后，调制信号的频谱被搬移到载频附近，成为下边频（带）和下边频（带）。信号功率正比于频谱的平方，所以可以看出载波信号的功率占整个调幅波功率的绝大部分。同样由公式计算可以得到 其中，为载波功率。所以在发送中，我们可以只发送上、下边带，此时为dsb调制；只发送上边带或者下边带，此时为ssb。 dsb调制 在幅度调制的一般模型中，若假设滤波器为全通网络（1），调制信号m(t)中无直流分量，则输出的已调信号就是无载波分量的双边带调制信号，或称抑制载波双边带（dsb-sc）调制信号，简称双边带（dsb）信号。dsb调制器模型如下图3所示。 图3下面利用matlab实现dsb的调制。为方便观察，令输入信号为单一正弦波信号。 -dsb调制信号输出和频谱- dt=0.0001; %时间采样频率 fh=10; %调制信号最高频率 fc=1000; %载波中心频率 t=0.1; %信号时长 n=t/dt; t=0:n-1*dt; mt=sqrt(2)*cos(2*pi*fh*t); %调制信号 s_dsb=mt.*cos(2*pi*fc*t); f,sf=fft_shift(t,s_dsb); %已调信号频谱 subplot(211) plot(t,s_dsb);hold on; %画出dsb信号波形 plot(t,mt,r-); %标示mt波形 plot(t,-mt,r-); title(dsb调制信号及其包络); xlabel(时间t); subplot(212) plot(f,psd); axis(750, 1250, 0 ,0.015); title(dsb信号频谱); xlabel(频率f); 图4结果分析：从图4可以看出，dsb信号时域图与频谱中不含有载波分量，仅有上下边带，发射功率大大降低。 ssb由于dsb信号的上、下两个边带是完全对称的，皆携带了调制信号的全部信息，因此，从信息传输的角度来考虑，仅传输其中一个边带就够了。用滤波法实现单边带调制的原理图如下图5所示 图5 下面利用matlab实现ssb的调制。为方便观察，令输入信号为单一正弦波信号。 -ssb调制信号输出和频谱- dt=0.0001; %时间采样频率 fh=10; %调制信号最高频率 fc=1000; %载波中心频率 t=0.1; n=t/dt; t=0:n-1*dt; mt=sqrt(2)*cos(2*pi*fh*t); %调制信号 s_ssb=real(hilbert(mt).*exp(j*2*pi*fc*t); f,sf=fft_shift(t,s_ssb); %单边带信号频谱 subplot(211) plot(t,s_ssb);hold on; %画出ssb信号波形 plot(t,mt,r-); %标示mt 的包络 plot(t,-mt,r-); %标示mt 的包络 title(ssb调制信号); xlabel(时间t); subplot(212) plot(f,abs(sf); axis(750 1250 0 max(abs(sf); title(ssb信号频谱); xlabel(频率f); 图6结果分析： 从图6的调制信号频谱可以看出，输出信号只有上（或者下）边带，大大减小了发射功率。而且ssb信号频带可节约一半，这对于日益拥挤的短波波段（330mhz）来说有重大意义，因为这样就能在同一波段中，使所容纳的频道数目增加一倍，大大提高了短波波段的利用率。总结： 在以上模拟线性调制中，am波发送信号的上下边带和载波，带宽为调信号的最大频率的两倍，发射功率较大；dsb波只发送信号的上下边带，发射功率大大减小，但带宽仍然为调信号的最大频率的两倍；ssb波只发送信号上边带或者下边带，带宽是am与dsb波的一半。但是，在解调过程中，am波采用非相干解调（包络检波），电路十分简单，采用二极管检波电路就可以恢复出调制信号。dsb、ssb需采用相干解调，电路较复杂。所以在广播通信中广泛采用am调制。2.角度调制 在调制时，若载波的频率随调制信号变化，称为频率调制，简称调频；若载波的相位随调制信号变化，称为相位调制，简称调相。在这两种调制中，载波的幅度保持恒定不变，而频率和相位的变化都表现为载波的瞬时相位变化，故把调频和调相统称为角度调制或调角。调角的一般表达式为 式中：a为载波的恒定振幅，为信号的瞬时相位，极为；为相对于载波的相位的瞬时相位偏移；是信号的瞬时角频率，记为；而成为相对于载频的瞬时频偏。 fm调频调频就是使瞬时频率偏移随调制信号成比例变化，即 为调制灵敏度。 所以频率调制的一般表达式为： 根据调制后载波瞬时相位偏移的大小，可将频率调制分为宽带调频（wbfm）与窄带调频（nbfm）。宽带与窄带调制的区分并无严格的界限，但通常认为由调频所引起的最大瞬时相位偏移远小于30时，即 称为窄带调频。否则，称为宽带调频。宽带调频一般用直接调频法，框图如图7。 图 7 直接调频法窄带调制一般使用间接间接调频法，框图如图8。 图 8 间接调频法 为方便起见，无妨假设正弦载波的振幅a1，则由调频信号的一般表达式，通过化解，利用傅立叶变化公式可得nbfm信号的频域表达式：在nbfm中，由于下边频为负，因而合成矢量不与载波同相，而是存在相位偏移，当最大相位偏移满足式时，合成矢量的幅度基本不变，这样就形成了fm信号。 图9 nbfm信号频谱根据卡森公式所以，在宽带调频中， 在窄带调皮中， 下面利用matlab实现fm的调制。 -fm调制信号输出和频谱- t0=0.2; %信号的持续时间,用来定义时间向量ts=0.001; %抽样间隔fs=1/ts; %抽样频率fc=300; %载波频率,fc可以任意改变t=-t0/2:ts:t0/2; %时间向量kf=100; %偏差常数df=0.25; %所需的频率分辨率，用在求傅里叶变换时，它表示 %fft的最小频率间隔m=sin(100*t); %调制信号，m(t)可以任意更改int_m(1)=0; %求信号m(t)的积分for i=1:length(t)-1 int_m(i+1)=int_m(i)+m(i)*ts;endm,m,df1=fftseq(m,ts,df); %对调制信号m(t)求傅里叶变换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=u/fs; %缩放figure(1)subplot(2,1,1) %子图形式显示结果plot(t,m(1:length(t) %现在的m信号是重新构建的信号 %因为在对m求傅里叶变换时axis(-0.1 0.1 -1 1) %定义两轴的刻度xlabel(时间t) title(原调制信号的时域图)subplot(2,1,2)plot(t,u(1:length(t)axis(-0.1 0.1 -1 1)xlabel(时间t)title(已调信号的时域图)figure(2)subplot(2,1,1)plot(f,abs(fftshift(m)%fftshift：将fft中的dc分量移到频谱中心axis(-600 600 0 0.04)xlabel(频率f)title(原调制信号的频谱图)subplot(2,1,2)plot(f,abs(fftshift(u)axis(-600 600 0 0.04)xlabel(频率f)title(已调信号的频谱图)function m,m,df=fftseq(m,ts,df)fs=1/ts;if nargin=2 n1=0; %nargin为输入参量的个数else n1=fs/df;endn2=length(m);n=2(max(nextpow2(n1),nextpow2(n2); %nextpow2(n)取n最接近的较大2次幂m=fft(m,n); %m为信号m的傅里叶变换，n为快速傅 里叶变换的点数，及基n-fft变换m=m,zeros(1,n-n2); %构建新的m信号df=fs/n; 图10 图11结果分析： 从图10已调信号的时域图中可以看出载波的频率随着调制信号的改变而改变，幅度不发生变化。从图11频谱图中可以得出，已调信号的带宽明显大于调制信号最高频率的两倍，并且产生了很多新的频率分量，是一种非线性调制。从频谱的幅值大小还可以看出，载波频谱的幅值明显大于已调波频谱的幅值，实际上，调频实现了载波的能量的搬移。 pm调制在模拟调制中，当幅度和频率保持不变时，改变载波的相位使之随未调信号的大小而改变，这就是调相。设高频载波为u=ucost，调制信号为u(t)，则调相信号的瞬时相位 (t)=+ku(t)瞬时角频率 (t)=+k调相信号 u=ucost+ku(t) 将信号的信息加在载波的相位上则形成调相信号，调相的表达式为： s(t)=acost+kf(t)+ 这里k称为相移指数，这种调制方式，载波的幅度和角频率不变，而瞬时相位偏移是调制信号f(t)的线性函数，称为相位调制。调相与调频有着相当密切的关系，我们知道相位与频率有如下关系式： =+kf(t) (t)=t+k所以在调相时可以先将调制信号进行微分后在进行频率调制，这样等效于调相，此方法称为间接调相，与此相对应，上述方法称为直接调相。调相信号的产生如图12所示： 图12 pm调相信号的产生实现相位调制的基本原理是使角频率为的高频载波u(t)通过一个可控相移网络, 此网络产生的相移受调制电压u(t)控制, 满足=ku(t)的关系。下面利用matlab实现pm的调制。-pm调制信号输出和频- t0=0.2; %信号的持续时间,用来定义时间向量ts=0.0001; %抽样间隔fs=1/ts; %抽样频率fc=300; %载波频率,fc可以任意改变t=-t0/2:ts:t0/2; %时间向量kf=100; %偏差常数df=0.25; %所需的频率分辨率，用在求傅里叶变换时，它表示fft的最小频率间隔m=sin(100*t); %调制信号，m(t)可以任意更改int_m(1)=0; %求信号m(t)的积分for i=1:length(t)-1 int_m(i+1)=int_m(i)+m(i)*ts;endm,m,df1=fftseq(m,ts,df); %对调制信号m(t)求傅里叶变换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=u/fs; %缩放 %通过调用子程序env_phas和loweq来实现解调功能v,phase=env_phas(u,ts,fc); %解调，求出u的相位phi=unwrap(phase); %校正相位角，使相位在整体上连续，便于%后面对该相位角求导dem=(1/(2*pi*kf)*(diff(phi)*fs); %对校正后的相位求导 %再经一些线性变换来恢复原调制信号 %乘以fs是为了恢复原信号，因为前面使用了缩放 figure(1)subplot(2,1,1) %子图形式显示结果plot(t,m(1:length(t) %现在的m信号是重新构建的信号 axis(-0.1 0.1 -1 1) %定义两轴的刻度xlabel(时间t) title(原调制信号的时域图)subplot(2,1,2)plot(t,u(1:length(t)axis(-0.1 0.1 -1.5 1.5)xlabel(时间t)title(已调信号的时域图)figure(2)subplot(2,1,1)plot(f,abs(fftshift(m) %fftshift：将fft中的dc分量移到频谱中心axis(-600 600 0 0.1)xlabel(频率f)title(原调制信号的频谱图)subplot(2,1,2)plot(f,abs(fftshift(u)axis(0 600 0 0.05)xlabel(频率f)title(已调信号的频谱图) 图13 图14结果分析： 从图15已调信号的时域图中可以看出载波的频率随着调制信号的改变而改变，幅度不发生变化。从图16频谱图中可以得出，已调信号的带宽明显大于调制信号最高频率的两倍，并且产生了很多新的频率分量，是一种非线性调制。总结：v 调频与调相的比较 fm、pm很相似，都是一种非线性调制，产生了新的频率分量，带宽明显增大。调频与调相还是有区别的，从公式可以看出，调频波的额最大频移与调制频率无关，最大相移与调制频率成反比，调相波的最大频移与调制频率成正比，最大相移则与调制频率无关。这是两种调制的根本区别。正是由于这一点，调频波的频谱宽度对