FM调制与解调系统课程设计报告

1、 43 / 43信号与线性系统课程设计报告FM调制与解调系统的设计摘要：调频和调相是广泛采用的两种调角的基本调制方式。其中调频（FM）是载波信号的频率按调制信号的规律变化；调相（PM）是载波信号的相位按调制信号的规律变化。两种调制方式都表现为信号的瞬时相位受到调变。调频波的解调称为鉴频；调相波的解调称为鉴相。在掌握模拟系统FM和PM调制与解调原理和设计方法的基础上， 可以通过MATLAB进行编程仿真实现对系统的时域、频域特性分析，可以通过Simulink动态建模和Labview虚拟仪器对系统进行仿真，检测所设计系统的功能，还可以通过GUI设计实现针对该系统的图形用户界面。关键词：调制，解调，系

2、统，仿真一、 课题的目的本课程设计课题主要研究FM 调制与解调模拟系统的理论设计和软件仿真方法。通过完成本课题的设计，拟主要达到以下几个目的：1掌握模拟系统FM 调制与解调的原理。2掌握模拟系统FM 调制与解调的设计方法； 3掌握应用MATLAB分析系统时域、频域特性的方法，进一步锻炼应用Matlab进行编程仿真的能力；4熟悉基于Simulink的动态建模和仿真的步骤和过程；5了解基于LabVIEW虚拟仪器的特点和使用方法，熟悉采用LabVIEW进行仿真的方法。二、 课题任务设计FM调制与解调模拟系统，仿真实现相关功能。 包括: 可实现单音调制的FM调制及解调、PM调制及解调的系统设计及仿真，

3、要求给出系统的设计框图、源程序代码及仿真结果，并要求给出程序的具体解释说明，记录系统的各个输出点的波形和频谱图。具体内容为：(1）设计FM调制与解调、PM调制与解调的模拟系统，给出系统的原理框图,对系统的主要参数进行设计说明。(2）采用Matlab语言设计相关程序,实现系统的功能，要求采用两种方式进行仿真，即直接采用Matlab语言编程的静态仿真方式、采用Simulink进行动态建模和仿真的方式。要求采用两种以上调制信号源进行仿真，并记录系统的各个输出点的波形和频谱图。(3)采用LabVIEW进行仿真设计,实现系统的功能,要求给出系统的前面板和框图，采用两种以上调制信号源进行仿真，并记录仿真结

4、果。(4)要求对系统的时域、频域特性进行分析，并与理论设计结果进行比较分析。(5)对系统功能进行综合测试，整理数据，撰写设计报告。三、 设计方案论证调频（FM）：载波信号的频率按调制信号的规律变化；调相（PM）：载波信号的相位按调制信号的规律变化；两种调制方式均表现为载波信号的瞬时相位收到调变。在设计调频（FM）和调相（PM）系统时，由于调频（FM）可以看成是相位按调制信号的时间积分值规律变化的的调相信号，故可以都借助调相系统，只是调频时，需先将调制信号积分以后再通过调相系统。设计调相系统时，在确定载波幅度和频率、调制信号幅度和频率、调频或调相比例常数后，将调制信号与比例常数的乘积（调相时）或

5、 调制信号积分以后与比例常数的乘积（调频时），即可完成调相功能。由于等幅调频波通过延时网络以后，在限制0/12(为调制信号的角频率，0为延时量)的条件下，输出调频波与输入时相比会产生附加相移，该附加相移反映了输入调频波的瞬时频率变化，即输出为调频-调相波，再将其通过鉴相器即可完成鉴频。此即相位鉴频器的理论设计模型。假设调相波和载波之间的相位差为，使调相波产生900固定相移后，和载波一起通过相乘器，再通过低通滤波器后，输出波在限制|/12的条件下，幅度反映了输入调相波的瞬时相位变化，即完成了鉴相功能。此即乘积型鉴相器的理论设计模型。故在设计鉴频（FM）和鉴相（PM）系统时，只需在确定系统参数的情

6、况下，按上述原理完成相位鉴频器和乘积型鉴相器的功能即可完成相应的解调功能。在MATLAB静态编程仿真部分，将上述原理转换成相应的数学关系，即可完成编程；在Simulink动态建模仿真部分，根据上述原理选择合适的模型，并做合适的参数设置，即可完成系统方框图的设计；在Labview仿真部分，根据上述原理选择合适的函数和控件，并做合适的参数设置，即可完成前面板和程序框图的设计；在GUI图形用户界面设计部分，在了解设计方法的基础上，根据上述原理给添加的按钮或文本框添加对应的代码，即可实现图形用户界面。四、 设计内容、步骤和要求4.1必选部分(1) 设计实现FM、PM调制与解调的模拟系统，给出系统的原理

7、框图,对系统的主要参数进行设计说明,具体参数包括：载波频率、调制信号频率、载波大小、调制信号大小、调制系数等参数。并对所设计的系统进行理论分析计算。(2) 根据所设计的FM、PM调制与解调的模拟系统，进行基于Matlab语言的静态仿真设计。分别实现单音调制的FM调制及解调、PM调制及解调的系统仿真设计，要求给出系统的Matlab编程仿真程序及结果，并要求写出程序的具体解释说明，记录系统的各个输出点的波形和频谱图。要求调制信号分别采用不同类型的信号进行仿真，至少给出两种以上调制信号源，具体参数自定。载波信号频率根据设计情况设定。(3) 根据所设计的FM、PM调制与解调的模拟系统，采用Simuli

8、nk进行动态建模仿真设计。分别实现单音调制的FM调制及解调、PM调制及解调的系统动态仿真设计，要求包括调制和解调的部分，并给出采用Simulink进行动态建模仿真的系统方框图，同时记录系统的各个输出点的波形和频谱图。要求采用两种以上调制信号源进行仿真，具体参数自定。载波信号频率根据设计情况设定。(4) 根据仿真结果,对系统的时域、频域特性进行分析，并与理论设计结果进行比较分析。4.2 拓展部分(1) 根据所设计的FM、PM调制与解调的模拟系统，说明具体的参数，进行基于LabVIEW环境的仿真，分别实现单音调制的FM调制及解调、PM调制及解调的系统仿真设计，要求包括调制和解调的部分，给出系统的前

9、面板和框图，并记录仿真结果。(2) 要求调制信号采用不同类型的信号源，进行进一步的仿真，给出系统的前面板和框图，并记录仿真结果，观察分析频谱的变化情况。(3) 比较分析采用以上两种软件环境：Matlab与LabVIEW，进行仿真的各自的特点，分析说明不同语言环境的各自优势。4.3 选作部分(1) 根据所设计的FM、PM调制与解调的模拟系统，采用Matlab语言设计相关程序，并且利用GUI设计图形用户界面，完成4.1必选部分所要求的功能。(2) 分析比较采用Matlab语言的GUI设计图形用户界面与采用LabVIEW进行设计的各自的特点。五、 实验具体过程、结果与分析实现FM、PM调制与解调的模

10、拟系统的原理框图：图3 相位鉴频器原理框图图1 间接调频原理框图图4 乘积型鉴相器原理框图图2 矢量合成法调相原理框图MATLAB静态编程仿真FM和PM的调制程序流程图：确定系统的主要参数构建调制信号和载波信号调频OR调相对调制信号积分与相应的比例常数相乘加到载波的相位上得到调频或调相信号调频调相编程仿真程序代码参数设置为：调制信号频率：10Hz；调制信号幅度：1V；载波信号频率：100Hz；载波信号幅度：1V；调频比例常数kf：2000；调相比例常数kp：100；采样频率fs: 500Hz。正弦波：T=0.002;%采样间隔fs=1/T;%采样频率N=256;%采样点数n=0:N-1;t=n

11、*T;vx=cos(2*pi*10*t);%调制信号int_vx(1)=0;for i=1:(length(t)-1) int_vx(i+1)=int_vx(i)+vx(i)*T;endvc=cos(2*pi*100*t);%载波信号kf=2000;%调频比例常数vf=cos(2*pi*100*t+2*pi*kf*int_vx);%调频信号VX=fft(vx,N);VC=fft(vc,N);VF=fft(vf,N);f=n*fs/N;figure(1)subplot(3,2,1)plot(vx);title('时域调制信号')grid onsubplot(3,2,2)plot(

12、f(1:N/2),abs(VX(1:N/2);title('频域调制信号')xlabel('Hz')grid onsubplot(3,2,3)plot(vc);title('时域载波信号')grid onsubplot(3,2,4)plot(f(1:N/2),abs(VC(1:N/2);title('频域载波信号')xlabel('Hz')grid onsubplot(3,2,5)plot(vf);title('时域调频信号')grid onsubplot(3,2,6)plot(f(1:N/2),a

13、bs(VF(1:N/2);title('频域调频信号')xlabel('Hz')grid onkp=100;%调相比例常数vp=cos(2*pi*100*t+kp*vx);%调相信号VP=fft(vp,N);figure(2)subplot(3,2,1)plot(vx);title('时域调制信号')grid onsubplot(3,2,2)plot(f(1:N/2),abs(VX(1:N/2);title('频域调制信号')xlabel('Hz')grid onsubplot(3,2,3)plot(vc);tit

14、le('时域载波信号')grid onsubplot(3,2,4)plot(f(1:N/2),abs(VC(1:N/2);title('频域载波信号')xlabel('Hz')grid onsubplot(3,2,5)plot(vp);title('时域调相信号')grid onsubplot(3,2,6)plot(f(1:N/2),abs(VP(1:N/2);title('频域调相信号')xlabel('Hz')grid on方波：T=0.002;%采样间隔fs=1/T;%采样频率N=256;%采

15、样点数n=0:N-1;t=n*T;vx=square(2*pi*10*t,50);%调制信号int_vx(1)=0;for i=1:(length(t)-1) int_vx(i+1)=int_vx(i)+vx(i)*T;endvc=cos(2*pi*100*t);%载波信号kf=2000;%调频比例常数vf=cos(2*pi*100*t+kf*int_vx);%调频信号VX=fft(vx,N);VC=fft(vc,N);VF=fft(vf,N);f=n*fs/N;figure(1)subplot(3,2,1)plot(vx);title('时域调制信号')axis(0,300,

16、-1.2,1.2)grid onsubplot(3,2,2)plot(f(1:N/2),abs(VX(1:N/2);title('频域调制信号')xlabel('Hz')grid onsubplot(3,2,3)plot(vc);title('时域载波信号')grid onsubplot(3,2,4)plot(f(1:N/2),abs(VC(1:N/2);title('频域载波信号')xlabel('Hz')grid onsubplot(3,2,5)plot(vf);title('时域调频信号')g

17、rid onsubplot(3,2,6)plot(f(1:N/2),abs(VF(1:N/2);title('频域调频信号')xlabel('Hz')grid onkp=100;%调相比例常数vp=cos(2*pi*100*t+kp*vx);%调相信号VP=fft(vp,N);figure(2)subplot(3,2,1)plot(vx);title('时域调制信号')axis(0,300,-1.2,1.2)grid onsubplot(3,2,2)plot(f(1:N/2),abs(VX(1:N/2);title('频域调制信号'

18、;)xlabel('Hz')grid onsubplot(3,2,3)plot(vc);title('时域载波信号')grid onsubplot(3,2,4)plot(f(1:N/2),abs(VC(1:N/2);title('频域载波信号')xlabel('Hz')grid onsubplot(3,2,5)plot(vp);title('时域调相信号')grid onsubplot(3,2,6)plot(f(1:N/2),abs(VP(1:N/2);title('频域调相信号')xlabel(&

19、#39;Hz')grid on运行结果：正弦波:方波：FM和PM的解调程序流程图：确定系统参数构建调制信号、载波信号和调角信号鉴频OR鉴相将调频信号延时，得调频-调相信号鉴频鉴相给调角信号加上固定相移载波信号通过相乘器得到解调信号编程仿真程序代码参数设置如下：调制信号频率：10Hz；调制信号幅度：1V；载波信号频率：100Hz；载波信号幅度：1V；调频比例常数kf：pi/2；调相比例常数kp：pi/20；采样频率fs：500Hz；低通滤波器参数：正弦波：fir1(80,0.1)，方波：fir1(80,0.35)。正弦波：T=0.002;%采样间隔fs=1/T;%采样频率N=256;%采

20、样点数n=0:N-1;t=n*T;kf=pi/2;%调频比例常数kp=pi/20;%调相比例常数vx=cos(2*pi*10*t);%调制信号int_vx =sin(20*pi*t)/(20*pi);%对调制信号积分vc=cos(2*pi*100*t);%载波信号vf=cos(2*pi*100*(t-0.002)+2*pi*kf*(sin(20*pi*(t-0.002)/(20*pi);%对调频FM信号延时,相位变化中含0.4*pi的恒定相移vf1=cos(2*pi*100*(t-0.002)+2*pi*kf*(sin(20*pi*(t-0.002)/(20*pi)-pi/10);%对调频FM

21、信号延时后，再加上pi/10的固定相移vp1=cos(2*pi*100*t+kp*vx-pi/2);%给调相信号加上pi/2的固定相移vip=vc.*vp1;%通过相乘器vif=vc.*vf1;%通过相乘器b=fir1(80,0.1);vof=filter(b,1,vif);vop=filter(b,1,vip);VOF=fft(vof,N);VOP=fft(vop,N);f=n*fs/N;figure(1)subplot(2,2,1)plot(vx);title('时域调制信号')grid onVX=fft(vx,N);subplot(2,2,2)plot(f(1:N/2),

22、abs(VX(1:N/2);title('频域调制信号')xlabel('Hz')grid onsubplot(2,2,3)plot(vof);title('时域解调FM信号')grid onsubplot(2,2,4)plot(f(1:N/2),abs(VOF(1:N/2);title('频域调制FM信号')xlabel('Hz')grid onfigure(2)subplot(2,2,1)plot(vx);title('时域调制信号')grid onVX=fft(vx,N);subplot(2,

23、2,2)plot(f(1:N/2),abs(VX(1:N/2);title('频域调制信号')xlabel('Hz')grid onsubplot(2,2,3)plot(vop);title('时域解调PM信号')grid onsubplot(2,2,4)plot(f(1:N/2),abs(VOP(1:N/2);title('频域调制PM信号')xlabel('Hz')grid on方波：T=0.002;%采样间隔fs=1/T;%采样频率N=256;%采样点数n=0:N-1;t=n*T;kf=pi/2;%调频比例常

24、数kp=pi/20;%调相比例常数vx=square(2*pi*10*t,50);%调制信号int_vx=25*sawtooth(2*pi*t*10,0.5);%调制信号vc=cos(2*pi*100*t);%载波信号vf=cos(2*pi*100*(t-0.002)+2*pi*kf*25*sawtooth(2*pi*(t-0.002)*10,0.5);%对调频FM信号延时,相位变化中含0.4*pi的恒定相移vf1=cos(2*pi*100*(t-0.002)+2*pi*kf*25*sawtooth(2*pi*(t-0.002)*10,0.5)-pi/10);%对调频FM信号延时后，再加上pi

25、/10的固定相移vp1=cos(2*pi*100*t+kp*vx-pi/2);%给调相信号加上pi/2的固定相移vip=vc.*vp1;%通过相乘器vif=vc.*vf1;%通过相乘器b=fir1(80,0.35);b1=fir1(80,0.2);vof=filter(b1,1,vip);vop=filter(b,1,vip);VOF=fft(vof,N);VOP=fft(vop,N);f=n*fs/N;figure(1)subplot(2,2,1)plot(vx);title('时域调制信号')axis(0,300,-1.2,1.2)grid onVX=fft(vx,N);s

26、ubplot(2,2,2)plot(f(1:N/2),abs(VX(1:N/2);title('频域调制信号')xlabel('Hz')grid onsubplot(2,2,3)plot(vof);title('时域解调FM信号')grid onsubplot(2,2,4)plot(f(1:N/2),abs(VOF(1:N/2);title('频域调制FM信号')grid onfigure(2)subplot(2,2,1)plot(vx);title('时域调制信号')axis(0,300,-1.2,1.2)gri

27、d onVX=fft(vx,N);subplot(2,2,2)plot(f(1:N/2),abs(VX(1:N/2);title('频域调制信号')xlabel('Hz')grid onsubplot(2,2,3)plot(vop);title('时域解调PM信号')axis(0,300,-0.1,0.1)grid onsubplot(2,2,4)plot(f(1:N/2),abs(VOP(1:N/2);title('频域调制PM信号')grid on运行结果：正弦波：方波：结果分析：Simulink动态建模仿真FM调制和解调系统

28、方框图：参数设置：正弦波：方波:仅改变signal generator的参数设置如下，其余参数设置与正弦波时相同。运行结果：正弦波：示波器（从上到下一次是：调制信号、调频信号、解调信号）：调制信号频谱:调频信号频谱：解调信号频谱：方波:示波器（从上到下一次是：调制信号、调频信号、解调信号）：调频信号频谱:解调信号频谱：调制信号频谱：PM的调制和解调系统方框图：参数设置：正弦波：方波：仅改变signal generator的参数设置如下，其余参数设置与正弦波时相同：运行结果：正弦波：示波器（从上到下一次是：调制信号、调相信号、解调信号）：调相信号频谱：调制信号频谱： 解调信号频谱：方波：示波器（

29、从上到下一次是：调制信号、调相信号、解调信号）：调相信号频谱：调制信号频谱：调制信号频谱：解调信号频谱：结果分析：Labview仿真PM调制前面板：程序框图：参数设置及运行结果：FM调制前面板：程序框图：参数设置及运行结果：结果分析GUI设计图形用户界面程序代码：function varargout = jiaodu_tiaozhiANDjietiao(varargin)% JIAODU_TIAOZHIANDJIETIAO M-file for jiaodu_tiaozhiANDjietiao.fig% JIAODU_TIAOZHIANDJIETIAO, by itself, creates

30、a new JIAODU_TIAOZHIANDJIETIAO or raises the existing% singleton*.% H = JIAODU_TIAOZHIANDJIETIAO returns the handle to a new JIAODU_TIAOZHIANDJIETIAO or the handle to% the existing singleton*.% JIAODU_TIAOZHIANDJIETIAO('CALLBACK',hObject,eventData,handles,.) calls the local% function named C

31、ALLBACK in JIAODU_TIAOZHIANDJIETIAO.M with the given input arguments.% JIAODU_TIAOZHIANDJIETIAO('Property','Value',.) creates a new JIAODU_TIAOZHIANDJIETIAO or raises the% existing singleton*. Starting from the left, property value pairs are% applied to the GUI before jiaodu_tiaozhiA

32、NDjietiao_OpeningFcn gets called. An% unrecognized property name or invalid value makes property application% stop. All inputs are passed to jiaodu_tiaozhiANDjietiao_OpeningFcn via varargin.% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one% instance to run (singleton)&qu

33、ot;.% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help jiaodu_tiaozhiANDjietiao% Last Modified by GUIDE v2.5 26-Dec-2011 15:32:02% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct('gui_Name', mfilename, . 'gui_Singleton&#

34、39;, gui_Singleton, . 'gui_OpeningFcn', jiaodu_tiaozhiANDjietiao_OpeningFcn, . 'gui_OutputFcn', jiaodu_tiaozhiANDjietiao_OutputFcn, . 'gui_LayoutFcn', , . 'gui_Callback', );if nargin && ischar(varargin1) gui_State.gui_Callback = str2func(varargin1);endif nargo

35、ut varargout1:nargout = gui_mainfcn(gui_State, varargin:);else gui_mainfcn(gui_State, varargin:);end% End initialization code - DO NOT EDIT% - Executes just before jiaodu_tiaozhiANDjietiao is made visible.function jiaodu_tiaozhiANDjietiao_OpeningFcn(hObject, eventdata, handles, varargin)% This funct

36、ion has no output args, see OutputFcn.% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% varargin command line arguments to jiaodu_tiaozhiANDjietiao (see VARARGIN)% Choose default command line outp

37、ut for jiaodu_tiaozhiANDjietiaohandles.output = hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes jiaodu_tiaozhiANDjietiao wait for user response (see UIRESUME)% uiwait(handles.figure1);% - Outputs from this function are returned to the command line.function varargout = jiao

38、du_tiaozhiANDjietiao_OutputFcn(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT);% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Get default command l

39、ine output from handles structurevarargout1 = handles.output;% - Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with hand

40、les and user data (see GUIDATA)T=0.002;%采样间隔fs=1/T;%采样频率N=256;%采样点数n=0:N-1;t=n*T;f=n*fs/N;kf=2000;%调频比例常数vx=cos(2*pi*10*t);%调制信号int_vx(1)=0;for i=1:(length(t)-1) int_vx(i+1)=int_vx(i)+vx(i)*T;endvc=cos(2*pi*100*t);%载波信号vf=cos(2*pi*100*t+2*pi*kf*int_vx);%调频信号VX=fft(vx,N);VC=fft(vc,N);VF=fft(vf,N);sub

41、plot(3,2,1)plot(vx);title('时域调制信号')grid onsubplot(3,2,2)plot(f(1:N/2),abs(VX(1:N/2);title('频域调制信号')xlabel('Hz')grid onsubplot(3,2,3)plot(vc);title('时域载波信号')grid onsubplot(3,2,4)plot(f(1:N/2),abs(VC(1:N/2);title('频域载波信号')xlabel('Hz')grid onsubplot(3,2,5

42、)plot(vf);title('时域调频信号')grid onsubplot(3,2,6)plot(f(1:N/2),abs(VF(1:N/2);title('频域调频信号')xlabel('Hz')grid on% - Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)% hObject handle to pushbutton2 (see GCBO)% eventdata reserved - t

43、o be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)T=0.002;%采样间隔fs=1/T;%采样频率N=256;%采样点数n=0:N-1;t=n*T;kf=pi/2;%调频比例常数vx=cos(2*pi*10*t);%调制信号int_vx =sin(20*pi*t)/(20*pi);%对调制信号积分vc=cos(2*pi*100*t);%载波信号vf=cos(2*pi*100*(t-0.002)+2*pi*kf*(sin(20*pi*(t-0

44、.002)/(20*pi);%对调频FM信号延时,相位变化中含0.4*pi的恒定相移vf1=cos(2*pi*100*(t-0.002)+2*pi*kf*(sin(20*pi*(t-0.002)/(20*pi)-pi/10);%对调频FM信号延时后，再加上pi/10的固定相移vif=vc.*vf1;%通过相乘器b=fir1(80,0.1);vof=filter(b,1,vif);VOF=fft(vof,N);f=n*fs/N;subplot(2,2,1)plot(vx);title('时域调制信号')grid onVX=fft(vx,N);subplot(2,2,2)plot(

45、f(1:N/2),abs(VX(1:N/2);title('频域调制信号')xlabel('Hz')grid onsubplot(2,2,3)plot(vof);title('时域解调FM信号')grid onsubplot(2,2,4)plot(f(1:N/2),abs(VOF(1:N/2);title('频域调制FM信号')xlabel('Hz')grid on% - Executes on button press in pushbutton3.function pushbutton3_Callback(hO

46、bject, eventdata, handles)% hObject handle to pushbutton3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)T=0.002;%采样间隔fs=1/T;%采样频率N=256;%采样点数n=0:N-1;t=n*T;kp=100;%调相比例常数f=n*fs/N;vx=cos(2*pi*10*t);%调制信号vc=cos(2*p

47、i*100*t);%载波信号VX=fft(vx,N);VC=fft(vc,N);vp=cos(2*pi*100*t+kp*vx);%调相信号VP=fft(vp,N);subplot(3,2,1)plot(vx);title('时域调制信号')grid onsubplot(3,2,2)plot(f(1:N/2),abs(VX(1:N/2);title('频域调制信号')xlabel('Hz')grid onsubplot(3,2,3)plot(vc);title('时域载波信号')grid onsubplot(3,2,4)plot(

48、f(1:N/2),abs(VC(1:N/2);title('频域载波信号')xlabel('Hz')grid onsubplot(3,2,5)plot(vp);title('时域调相信号')grid onsubplot(3,2,6)plot(f(1:N/2),abs(VP(1:N/2);title('频域调相信号')xlabel('Hz')grid on% - Executes on button press in pushbutton4.function pushbutton4_Callback(hObject,

49、 eventdata, handles)% hObject handle to pushbutton4 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)T=0.002;%采样间隔fs=1/T;%采样频率N=256;%采样点数n=0:N-1;t=n*T;kp=pi/20;%调相比例常数vx=cos(2*pi*10*t);%调制信号vc=cos(2*pi*100*t);%载波信

50、号vp1=cos(2*pi*100*t+kp*vx-pi/2);%给调相信号加上pi/2的固定相移vip=vc.*vp1;%通过相乘器b=fir1(80,0.1);vop=filter(b,1,vip);VOP=fft(vop,N);f=n*fs/N;subplot(2,2,1)plot(vx);title('时域调制信号')grid onVX=fft(vx,N);subplot(2,2,2)plot(f(1:N/2),abs(VX(1:N/2);title('频域调制信号')xlabel('Hz')grid onsubplot(2,2,3)pl

51、ot(vop);title('时域解调PM信号')grid onsubplot(2,2,4)plot(f(1:N/2),abs(VOP(1:N/2);title('频域调制PM信号')xlabel('Hz')grid on% - Executes on button press in pushbutton5.function pushbutton5_Callback(hObject, eventdata, handles)% hObject handle to pushbutton5 (see GCBO)% eventdata reserved

52、- to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)T=0.002;%采样间隔fs=1/T;%采样频率N=256;%采样点数n=0:N-1;t=n*T;vx=square(2*pi*10*t,50);%调制信号int_vx(1)=0;for i=1:(length(t)-1) int_vx(i+1)=int_vx(i)+vx(i)*T;endvc=cos(2*pi*100*t);%载波信号kf=2000;%调频比例常数vf=cos(2*p

53、i*100*t+kf*int_vx);%调频信号VX=fft(vx,N);VC=fft(vc,N);VF=fft(vf,N);f=n*fs/N;subplot(3,2,1)plot(vx);title('时域调制信号')axis(0,300,-1.2,1.2)grid onsubplot(3,2,2)plot(f(1:N/2),abs(VX(1:N/2);title('频域调制信号')xlabel('Hz')grid onsubplot(3,2,3)plot(vc);title('时域载波信号')grid onsubplot(3,2,4)plot(f(1:N/2),abs(VC(1:N/2);title('频域载波信号')xlabel('Hz')grid onsubplot(3,2,5)plot(vf);title('时域调频信号')grid onsubplot(3,2,6)plot(f(1:N/2)