数字系统设计实习报告(滤波器 FM调制算法的实现)
FM调制与解调系统课程设计报告
FM调制与解调系统的设计摘要：调频和调相是广泛采用的两种调角的基本调制方式。
其中调频（FM）是载波信号的频率按调制信号的规律变化；调相（PM）是载波信号的相位按调制信号的规律变化。
两种调制方式都表现为信号的瞬时相位受到调变。
调频波的解调称为鉴频；调相波的解调称为鉴相。
在掌握模拟系统FM和PM调制与解调原理和设计方法的基础上，可以通过MATLAB进行编程仿真实现对系统的时域、频域特性分析，可以通过Simulink动态建模和Labview虚拟仪器对系统进行仿真，检测所设计系统的功能，还可以通过GUI设计实现针对该系统的图形用户界面。
关键词：调制，解调，系统，仿真一、课题的目的本课程设计课题主要研究FM 调制与解调模拟系统的理论设计和软件仿真方法。
通过完成本课题的设计，拟主要达到以下几个目的：1．掌握模拟系统FM 调制与解调的原理。
2．掌握模拟系统FM 调制与解调的设计方法；3．掌握应用MATLAB分析系统时域、频域特性的方法，进一步锻炼应用Matlab进行编程仿真的能力；4．熟悉基于Simulink的动态建模和仿真的步骤和过程；5．了解基于LabVIEW虚拟仪器的特点和使用方法，熟悉采用LabVIEW进行仿真的方法。
二、课题任务设计FM调制与解调模拟系统，仿真实现相关功能。
包括: 可实现单音调制的FM调制及解调、PM调制及解调的系统设计及仿真，要求给出系统的设计框图、源程序代码及仿真结果，并要求给出程序的具体解释说明，记录系统的各个输出点的波形和频谱图。
具体内容为：(1）设计FM调制与解调、PM调制与解调的模拟系统，给出系统的原理框图,对系统的主要参数进行设计说明。
(2）采用Matlab语言设计相关程序,实现系统的功能，要求采用两种方式进行仿真，即直接采用Matlab语言编程的静态仿真方式、采用Simulink进行动态建模和仿真的方式。
要求采用两种以上调制信号源进行仿真，并记录系统的各个输出点的波形和频谱图。
大学数字实验报告
实验名称：数字信号处理实验实验日期：2023年3月15日实验地点：大学计算机实验室实验目的：1. 理解数字信号处理的基本概念和原理。
2. 掌握数字滤波器的设计和实现方法。
3. 学会使用数字信号处理软件进行实验和分析。
4. 培养实验操作能力和数据分析能力。
实验原理：数字信号处理（Digital Signal Processing，DSP）是利用计算机或专用处理硬件对数字信号进行操作的一门技术。
它包括信号的采样、量化、滤波、变换、压缩、解压缩等处理过程。
本实验主要涉及数字滤波器的设计和实现。
实验仪器：1. 实验计算机2. 数字信号处理软件（如MATLAB）3. 示波器4. 音频播放器实验内容：一、实验一：数字滤波器的基本概念1. 实验目的：理解数字滤波器的基本概念，掌握滤波器的类型和特点。
2. 实验步骤：a. 在MATLAB中创建一个简单的数字滤波器模型，例如低通滤波器。
b. 使用MATLAB的内置函数进行滤波器的设计，并观察滤波器的频率响应。
c. 分析滤波器的性能，包括通带、阻带、过渡带等。
3. 实验结果：a. 设计了一个低通滤波器，其截止频率为1000Hz。
b. 频率响应显示，在截止频率以下，滤波器对信号有较好的抑制效果；在截止频率以上，滤波器对信号有较好的通过效果。
c. 分析结果表明，该滤波器满足实验要求。
二、实验二：FIR滤波器的设计1. 实验目的：掌握FIR滤波器的设计方法，学会使用MATLAB进行FIR滤波器的设计和实现。
2. 实验步骤：a. 设计一个具有线性相位特性的FIR滤波器，例如汉明窗设计。
b. 使用MATLAB的内置函数进行滤波器的设计，并观察滤波器的频率响应。
c. 对滤波器进行时域和频域分析，评估滤波器的性能。
3. 实验结果：a. 设计了一个具有线性相位特性的低通FIR滤波器，其截止频率为1000Hz。
b. 频率响应显示，该滤波器具有较好的线性相位特性，且在截止频率以下对信号有较好的抑制效果。
设计滤波器实验报告
设计滤波器实验报告设计滤波器实验报告引言：滤波器是信号处理中常用的工具，它可以通过选择性地传递或抑制特定频率的信号，对信号进行滤波。
本实验旨在设计并实现一个滤波器，通过对不同类型的信号进行滤波，验证滤波器的性能和效果。
一、实验目的本实验的主要目的是：1. 了解滤波器的基本原理和分类；2. 掌握滤波器的设计方法和实现技巧；3. 验证滤波器的性能和效果。
二、实验原理滤波器根据其频率响应特性可分为低通、高通、带通和带阻滤波器。
低通滤波器能够通过低频信号，抑制高频信号。
高通滤波器则相反，能够通过高频信号，抑制低频信号。
带通滤波器则能够通过一定范围内的频率信号，抑制其他频率信号。
带阻滤波器则相反，能够抑制一定范围内的频率信号，通过其他频率信号。
三、实验步骤1. 确定滤波器类型和频率响应特性；2. 根据所选滤波器类型和频率响应特性，设计滤波器的传递函数；3. 根据传递函数，计算滤波器的电路参数；4. 根据计算结果，搭建滤波器电路；5. 连接信号源和示波器，输入信号；6. 调节信号源的频率，并观察示波器上的输出信号；7. 对比输入信号和输出信号的频谱特性，验证滤波器的性能和效果。
四、实验结果与分析在实验中，我们设计了一个低通滤波器，频率响应特性为通过0-1 kHz的低频信号，抑制1 kHz以上的高频信号。
通过计算和搭建电路，我们成功实现了滤波器的设计。
在实验中，我们输入了不同频率的信号，并观察了输出信号的频谱特性。
结果显示，当输入信号的频率低于1 kHz时，输出信号基本保持不变；当输入信号的频率高于1 kHz时，输出信号的幅度逐渐减小，直至完全抑制。
通过对比输入信号和输出信号的频谱特性，我们可以清楚地看到滤波器对高频信号的抑制效果。
这表明我们设计的滤波器能够有效地滤除高频噪声，保留低频信号。
五、实验总结本实验通过设计滤波器并验证其性能，使我们更加深入地了解了滤波器的原理和应用。
通过实际操作，我们掌握了滤波器的设计方法和实现技巧。
调试滤波实习报告
一、实习目的本次实习旨在通过实际操作，掌握滤波电路的基本原理，了解滤波器在电子系统中的应用，并学会调试滤波电路，提高自己的实践操作能力。
二、实习时间2023年3月1日至2023年3月10日三、实习地点某电子科技公司实验室四、实习内容1. 学习滤波电路的基本原理，包括低通、高通、带通和带阻滤波器等。
2. 分析滤波电路的组成，了解各种滤波元件的特性。
3. 学习滤波电路的调试方法，包括电路参数的调整、滤波效果的测试等。
4. 实际操作，调试滤波电路，解决实际问题。
五、实习过程1. 学习滤波电路的基本原理在实习初期，我通过查阅资料、课堂学习等方式，了解了滤波电路的基本原理。
滤波电路主要利用电容、电感等元件对信号频率进行选择性过滤，从而实现信号的平滑处理。
根据信号频率的不同，滤波电路可分为低通、高通、带通和带阻滤波器。
2. 分析滤波电路的组成在了解了滤波电路的基本原理后，我开始分析滤波电路的组成。
滤波电路主要由滤波元件（电容、电感）、放大器、反馈网络等组成。
通过调整滤波元件的参数，可以实现不同频率的信号过滤。
3. 学习滤波电路的调试方法在实习过程中，我学习了滤波电路的调试方法。
主要包括以下几个方面：（1）调整滤波元件参数，实现不同频率的信号过滤；（2）通过测试滤波效果，判断滤波电路的性能；（3）根据实际需求，对滤波电路进行优化。
4. 实际操作，调试滤波电路在掌握了滤波电路的调试方法后，我开始实际操作，调试滤波电路。
首先，根据设计要求，选择合适的滤波元件；然后，按照电路图组装滤波电路；最后，通过测试滤波效果，对滤波电路进行调整，以达到预期的滤波效果。
六、实习收获通过本次实习，我收获颇丰：1. 掌握了滤波电路的基本原理和组成；2. 学会了滤波电路的调试方法，提高了自己的实践操作能力；3. 深入了解了滤波器在电子系统中的应用，为今后从事相关工作打下了基础。
七、实习建议1. 在实习过程中，要注重理论与实践相结合，多动手操作，提高自己的实践能力；2. 加强与指导老师的沟通交流，及时解决实习过程中遇到的问题；3. 积极参与实习项目，积累实际工作经验，为今后的发展奠定基础。
fm调制与解调实验心得
fm调制与解调实验心得
进行FM调制与解调的实验是一个有趣且富有挑战性的过程。
在这个实验中，我学到了关于调制和解调的基本原理，并且亲自动手进行了实践和验证。
通过实验，我深刻理解了FM调制的原理。
调制过程中，通过改变信号的频率，将音频信号转换为一个频率可变的载波信号，从而实现信号的传输。
在实验中，我使用了一个音频信号发生器来产生音频信号，将其输入到调制器中。
通过调整调制器的调制指数，我成功地调制出了频率可变的载波信号。
解调过程同样具有重要的意义。
通过解调，我们可以将调制过的信号恢复回原始的音频信号。
在实验中，我使用了一个解调器和一个带通滤波器来进行解调。
通过将调制后的信号输入到解调器中，并结合合适的解调器参数和带通滤波器的设置，我成功地将信号解调并得到了原始的音频信号。
这个实验使我更加了解了FM调制与解调的原理和技术，同时也对信号传输和调制解调技术有了更深入的理解。
同时，通过实践，我也学会了如何正确地操作和调整相关的设备，以实现预期的调制和解调效果。
在实验过程中，我也遇到了一些挑战。
调制参数的选择和调整是一个关键的部分，需要进行反复尝试和调整。
另外，信号的噪声和干扰也会对调制和解调的效果产生影响，因此需要采取适当的措施来减少干扰，并保证信号的质量。
南理工电类综合实验报告FM调制解调的数字实现
《电类综合实验》仿真报告实验课题：FM调制解调的数字实现指导教师：刘光祖学生姓名：院系：电光学院专业：通信与信息系统实验时间：2016.05.09至2016.05.13一、实验背景1.FPGA简介FPGA（Field－Programmable Gate Array，现场可编程逻辑门阵列）是在PAL、GAL、CPLD等可编程逻辑器件的基础上进一步发展的产物，是专用集成电路领域一种半定制的集成数字芯片，其最大特点是现场可编程，既解决了全定制电路的不足，又克服了原有可编程逻辑器件门电路数有限的缺点。
FPGA的内部结构由CLB、RAM、DCM、IOB、Interconnect 等构成。
如下给出了FPGA的最典型的结构：FPGA开发的一般流程：1(a) 1(b)2.DE2-115开发板简介本实验中所用开发板为Altera公司的DE2-115。
如下为开发板的配置列表：•Altera Cyclone® IV 4CE115 FPGA 器件•Altera 串行配置芯片：EPCS64•USB Blaster 在线编程；也支持JTAG和AS可编程方式•2MB SRAM•两片64MB SDRAM•8MB Flash memory•SD卡插槽•4个按钮•18个滑动开关•18个红色LED•9个绿色LED•50M时钟源•24位音频编解码器，麦克风插孔•电视解码•RJ45 2G以太网接口•VGA连接器•含有USB_A和USB_B连接器的主从控制器•RS232收发器和9针连接器•PS/2鼠标和键盘连接器•红外接收器ControlPanel是开发板自带的一个工具软件，可以通过该软件提供的图形界面直接对FPGA上的各个外设进行操作。
通过该操作可以确认PC机与开发板的连接是否正确，开发板的硬件工作是否正常。
ControlPanel的安装过程如下：确保QUARTUS II 10.0 或以上版本能被成功安装；将开关RUN/PROG切换到RUN位置；将USB接线连接至USB驱动端口，供12V电源并打开开关；打开主机上的可执行文件DE2_115_ControlPanel.exe，controlPanel的用户界面如下：DE2_115_ControlPanel.exe一旦被启动，DE2_115_ControlPanel.sof程序流文件将会被自动加载；如果未连接，点击CONNECT，点sof文件将会重新加载到板子上；注意，控制面板将会占用一直到你关闭那个端口，除非你关闭USB端口，否则你不能使用QUARTUS II 来下载文件；控制面板现在可以使用了，通过设置一些LED灯ON/OFF的状态来观察DE2-115上的状态。
数字滤波器设计实验报告
数字滤波器设计实验报告刘古城65100609一、实验目的研究数字滤波器的设计思想，理解数字频域，模拟频域的关系，掌握数字系统处理模拟信号的方法。
FIR数字滤波器设计：掌握窗函数设计FIR数字滤波器的方法，理解FIR的意义：线性相位。
二、实验原理1、FIR的特点（1）系统的单位冲击响应在有限个n值处不为零。
（2）对于稳定系统，系统函数在| z |>0处收敛，极点全部在z=0处。
（3）结构上主要是非递归结构，没有输出到输入的反馈，但在个别结构中（如频率抽样结构）也包含反馈的递归部分‘2、FIR滤波器的优点（1）即具有严格的线性相位，又具有任意的幅度’（2）FIR滤波器的抽样响应是有限长的，因而滤波器的性能稳定。
（3）只要经过一定的延时，任何非因果的有限长序列都能变成有限长的因果的序列，因而能用因果系统来实现。
（4）FIR滤波器单位冲击响应是有限长的，因而可以进行快速傅立叶变换，提高运算效率。
3、用窗函数设计FIR数字滤波器对函数加窗处理，实际是用一个有限长函数来逼近原函数。
常用的窗函数有矩形窗、三角窗，汉宁窗、海明窗、布莱克曼窗、凯撒窗等。
三、实验要求1、设计FIR数字低通滤波器，要求在不同窗口长度（N=15，33）下，分别求出h(n)，画出相应的幅频特性和相频特性曲线，观察3dB带宽和20dB带宽，总结窗口长度N对滤波特性的影响。
2、对三个拟合三角函数进行滤波处理。
3、对含噪心电信号函数进行滤波处理。
四、实验内容1、不同窗函数长度对于滤波特性的影响fs=100,N=32;n=0:N-1;t=n/fs;f0=n*fs/N;y=exp(-2*t);z=fft(y);m=abs(z);w1=blackman(N);z1=w1'.*y;x1=fft(z1),mo1=abs(x1);subplot(1,2,1);plot(f0,m/fs);subplot(1,2,2);plot(f0,mo1/fs)运行结果改变N值，令N=14，得到结果2、对三个拟合三角函数进行滤波clear;fs=2000;t=(1:1000)/fs;x=10*cos(2*pi*30*t)+cos(2*pi*150*t)+5*cos(2*pi*600*t); L=length(x);N=2^(nextpow2(L));Hw=fft(x,N);figure(1);subplot(2,1,1);plot(t,x);grid on;title('滤波前信号x');xlabel('时间/s');% 原始信号subplot(2,1,2);plot((0:N-1)*fs/L,abs(Hw));% 查看信号频谱grid on;title('滤波前信号频谱图');xlabel('频率/Hz');ylabel('振幅|H(e^jw)|');%% x_1=10*cos(2*pi*30*t)Ap=1;As=60;% 定义通带及阻带衰减dev=[(10^(Ap/20)-1)/(10^(Ap/20)+1),10^(-As/20)];% 计算偏移量mags=[1,0];% 低通fcuts=[60,100];% 边界频率[N,Wn,beta,ftype]=kaiserord(fcuts,mags,dev,fs);% 估算FIR滤波器阶数hh1=fir1(N,Wn,ftype,kaiser(N+1,beta));% FIR滤波器设计x_1=filter(hh1,1,x);% 滤波x_1(1:ceil(N/2))=[];% 群延时N/2，删除无用信号部分L=length(x_1);N=2^(nextpow2(L));Hw_1=fft(x_1,N);figure(2);subplot(2,1,1);plot(t(1:L),x_1);grid on;title('x_1=10*cos(2*pi*30*t)');xlabel('时间/s');subplot(2,1,2);plot((0:N-1)*fs/L,abs(Hw_1));% 查看信号频谱grid on;title('滤波后信号x_1频谱图');xlabel('频率/Hz');ylabel('振幅|H(e^jw)|');%% x_2=cos(2*pi*150*t)Ap=1;As=60;% 定义通带及阻带衰减dev=[10^(-As/20),(10^(Ap/20)-1)/(10^(Ap/20)+1),10^(-As/20)];% 计算偏移量mags=[0,1,0];% 带通fcuts=[80,120,180,220];% 边界频率[N,Wn,beta,ftype]=kaiserord(fcuts,mags,dev,fs);% 估算FIR滤波器阶数hh2=fir1(N,Wn,ftype,kaiser(N+1,beta));% FIR滤波器设计x_2=filter(hh2,1,x);% 滤波x_2(1:ceil(N/2))=[];% 群延时N/2，删除无用信号部分L=length(x_2);N=2^(nextpow2(L));Hw_2=fft(x_2,N);figure(3);subplot(2,1,1);plot(t(1:L),x_2);grid on;title('x_2=cos(2*pi*150*t)');xlabel('时间/s');subplot(2,1,2);plot((0:N-1)*fs/L,abs(Hw_2));% 查看信号频谱grid on;title('滤波后信号x_2频谱图');xlabel('频率/Hz');ylabel('振幅|H(e^jw)|');%% x_3=5*cos(2*pi*600*t)Ap=1;As=60;% 定义通带及阻带衰减dev=[10^(-As/20),(10^(Ap/20)-1)/(10^(Ap/20)+1)];% 计算偏移量mags=[0,1];% 高通fcuts=[500,550];% 边界频率[N,Wn,beta,ftype]=kaiserord(fcuts,mags,dev,fs);% 估算FIR滤波器阶数hh2=fir1(N,Wn,ftype,kaiser(N+1,beta));% FIR滤波器设计x_3=filter(hh2,1,x);% 滤波x_3(1:ceil(N/2))=[];% 群延时N/2，删除无用信号部分L=length(x_3);N=2^(nextpow2(L));Hw_3=fft(x_3,N);figure(4);subplot(2,1,1);plot(t(1:L),x_3);grid on;title('x_3=5*cos(2*pi*600*t)');xlabel('时间/s');subplot(2,1,2);plot((0:N-1)*fs/L,abs(Hw_3));% 查看信号频谱grid on;title('滤波后信号x_3频谱图');xlabel('频率/Hz');ylabel('振幅|H(e^jw)|');运行结果3、对含噪心电信号函数进行滤波处理。
数字滤波器设计实验报告
数字滤波器设计实验报告实验目的：1.掌握数字滤波器的基本理论知识。
2.学习数字滤波器设计方法。
3.实现数字滤波器的设计与模拟。
实验原理：FIR滤波器的特点是稳定性好、相位响应线性和易于设计。
FIR滤波器的设计方法主要有窗函数法、频率采样法和最小最大化法等。
IIR滤波器的特点是具有较窄的通频带宽率、相位响应非线性和较高的处理效率。
IIR滤波器的设计方法主要有双线性变换法、脉冲响应不变法和双正交变换法等。
实验步骤：1.根据实验要求和给定的参数，选择适合的滤波器类型（FIR或IIR）。
2.根据滤波器的设计方法，计算滤波器的系数。
3.使用MATLAB或其他工具进行滤波器的设计和仿真。
4.分析仿真结果，评估滤波器的性能。
5.根据实际需求，进行滤波器参数的优化和调整。
6.进行实验数据的滤波处理，并比较滤波前后的信号质量。
7.总结实验结果，写出实验报告。
实验结果：根据实验要求，我们选择了FIR滤波器进行设计。
通过使用窗函数法和最小最大化法，计算得到了滤波器的系数。
将滤波器的设计结果导入MATLAB进行仿真，得到了滤波器的频率响应和时域波形。
通过分析仿真结果，发现滤波器的设计基本满足了要求，但仍存在一些性能方面的改进空间。
根据实验需求和实际情况，我们对滤波器的参数进行了优化和调整。
经过多次迭代和调试，最终得到了满意的结果。
将优化后的滤波器应用于实验数据的滤波处理，可以看到滤波效果明显，信号质量得到了显著提升。
实验结论：通过本次实验，我们学习并掌握了数字滤波器的基本理论知识和设计方法。
通过实际操作和实验仿真，对数字滤波器的设计和应用有了更深入的了解。
实验结果表明，数字滤波器可以有效地对信号进行滤波处理，提高信号质量和准确度。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
数字系统设计实习报告指导老师：李林、邹修国、徐友
专业：电子信息科学与技术
组员：3238212李丽林
3238215郝文星
3238207王雷鸣
3238214金煜
基础部分—滤波器
一、设计目的
1.掌握用Matlab设计FIR数字高通滤波器的方法；
2.掌握高通滤波器截止频率变换的原理和方法；
3.利用设计的FIR数字滤波器，检验、观察滤波效果；
二、设计原理
用sine wave 、highpass filter和scope模块搭建FIR数字低通率波器。
改变高通滤波器的滤波系数，通过波特图观察滤波效果。
三、设计结果
FIR滤波器的滤波系数以及Matlab仿真的波特图
f=1KHz
(1)截止频率
c
f=2KHz
(2)截止频率
c
f=3KHz (3)截止频率
c
f=4KHz (4)截止频率
c
f=5KHz (5)截止频率
c
f=6KHz (6)截止频率
c
f=7KHz (7)截止频率
c
f=8KHz (8)截止频率
c
f=9KHz (9)截止频率
c
f=10KHz (10)截止频率
c
四、心得体会
实习过程中我了解到要快而准备的完成设计，必须对低通滤波器有个深刻的了解，另外模块的搭建成功与否在于参数的设计，这要求我们需对各个模块有个大致的了解。
专业部分—FM 调制算法的实现
一、 设计目的
1. 设计FM 调制通信系统，并得出仿真结果。
2. 熟悉MATLAB 文件中M 文件的使用方法，包括函数、原理和方法的应用。
3. 增强在通信原理仿真方面的动手能力与自学能力。
二、 设计内容和实验要求
系统建模、确定仿真算法、建立仿真模型、设计仿真程序、运行仿真程序、输出仿真结果
三、 实验原理
FM 调制原理：
频率调制的一般表达式[1]
为：
（2-1）
FM 和PM 非常相似，如果预先不知道调制信号的具体形式，则无法判断已调信号是调频信号还是调相信号。
图 2-3
()()()FM m t dt PM S t →∙→→⎰
图 2-4
图（2-3）所示的产生调频信号的方法称为直接调频法，图（2-4）所示的产生调频信号
的方法称为间接调频法[4]。
由于实际相位调制器的调节范围不可能超出，因而间接调频的方法仅适用于相位偏移和频率偏移不大的窄带调制情形，而直接调频则适用于宽带调制情形。
根据调制后载波瞬时相位偏移的大小，可将频率调制分为宽带调频（WBFM ）与窄带调频（NBFM ）。
宽带与窄带调制的区分并无严格的界限，但通常认为由调频所引起的最大瞬时相位偏移远小于30°时，称为窄带调频。
否则，称为宽带调频。
（2-2）
为方便起见，无妨假设正弦载波的振幅A ＝1，则由式（2-1）调频信号的一般表达式，得
()cos[()]
t
FM c F S t t K m t d -∞
=ω+τ⎰=
cos cos[()]sin sin[()]t
t
c F c F t K m
d K m d -∞
-∞
ωττ-ωττ⎰⎰ （2-3）
通过化解，利用傅立叶变化公式可得NBFM 信号的频域表达式：
（2-4）在NBFM中，由于下边频为负，因而合成矢量不与载波同相，而是存在相位偏移，当最大相位偏移满足式（2-2）时，合成矢量的幅度基本不变，这样就形成了FM信号
图2-5 NBFM信号频谱
四、实验结果
MATLAB源程序如下：
dt=0.001; %设定时间步长
t=0:dt:1.5; %产生时间向量
am=5; %设定调制信号幅度
fm=5; %设定调制信号频率
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); %调制信号
ts=0.001; %抽样间隔
fs=1/ts; %抽样频率
df=0.25; %所需的频率分辨率，用在求傅里叶变换时，它表示FFT的最小频率间隔
m=am*cos(2*pi*fm*t); %原调信号
fs=1/ts;
if nargin==2
n1=0;
else
n1=fs/df;
end
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; %时间向量对应的频率向量
fs=1/ts;
if nargin==2
n1=0;
else
n1=fs/df;
end
n2=length(sfm);
n=2^(max(nextpow2(n1),nextpow2(n2)));
U=fft(sfm,n);
u=[sfm,zeros(1,n-n2)];
df1=fs/n; %以上是对已调信号u求傅里变换
WORD格式整理
U=U/fs; %缩放
subplot(3,1,1);plot(t,mt); %绘制调制信号的时域图
xlabel('时间t');
title('调制信号的时域图');
subplot(3,1,2);plot(t,ct); %绘制载波的时域图
xlabel('时间t');
title('载波的时域图');
subplot(3,1,3);
plot(t,sfm); %绘制已调信号的时域图
xlabel('时间t');
title('已调信号的时域图');
subplot(2,1,1)
plot(f,abs(fftshift(M))) %fftshift:将FFT中的DC分量移到频谱中心
xlabel('频率f')
title('原调制信号的频谱图')
subplot(2,1,2)
plot(f,abs(fftshift(U)))
xlabel('频率f')
title('已调信号的频谱图')
五、实验小结与体会
通过这一次课程设计，我了解很多关于专业的知识，以前每次学这些知识时，总是不知道这些东西具体拿来有什么用，现在才知道，几个短短输入信号，在有了一个简单的电路流程后，就能仿真成我们生活中很多常见的东西。
总的来说，这次课程设计过程还是比较愉快轻松的，虽然中间有过一些困难，但是在老师与同学的指点下我还是渡过了，在这里我要谢谢帮助我的老师和同学。
专业资料值得拥有。