软件无线电实验matlab程序
基于matlab的软件无线电滤波器设计文献综述
毕业设计（论文）文献综述软件无线电中的FIR滤波器设计1软件无线电的发展与应用软件无线电是最近几年在无线通信领域提出的一种新的通信系统体系结构，它的基本思想是以开放性、可扩展、结构最简的硬件为通用平台，把尽可能多的通信功能用可升级、可替换的软件来实现。
软件无线电是具有可重配置硬件平台的无线设备，可以跨多种通信标准。
因为具有更低的成本、更大的灵活性和更高的性能，软件无线电已迅速成为军事、公共安全和商用无线领域的事实标准。
SDR成为商用流行的主要原因之一是它能够对多种波形进行基带处理和数字中频(IF)处理。
IF处理将数字信号处理的领域从基带扩展到RF。
支持基带和中频处理的能力增加了系统灵活性，同时减小了制造成本。
无线标准不断地发展，通过先进的基带处理技术如自适应调制编码、空时编码(STC)、波束赋形和多入多出(MIMO)天线技术，支持更高的数据速率。
基带信号处理器件需要巨大的处理带宽，以支持这些技术中大计算量的算法。
例如，美国军事联合战术无线系统(JTRS)定义了军事无线中20多种。
而数字滤波器在软件无线电中的有着非常重要的地位。
调谐本振分波段滤波器放大放大DDCA/DDSP软件图1 宽带中频数字化软件无线电接收机框图2 数字滤波器国内外发展现状数字滤波在DSP（数字信号处理）中占有重要地位。
数字滤波器按实现的网络结构或者从单位脉冲响应，分为IIR(无限脉冲响应）和FIR（有限脉冲响应）滤波器。
如果IIR滤波器和FIR滤波器具有相同的性能，那么通常IIR滤波器可以用较低的阶数获得高的选择性，执行速度更快，所有的储存单元更少，所有既经济又高效。
数字滤波器精确度高、使用灵活、可靠性高，具有模拟设备所没有的许多优点，已广泛地应用于各个科学技术领域，例如数字电视、语音、通信、雷达、声纳、遥感、图像、生物医学以及许多工程应用领域。
随着信息时代数字时代的到来，数字滤波技术已经成为一门及其重要的科学和技术领域。
以往的滤波器大多采用模拟电路技术，但是，模拟电路技术存在很多难以解决的问题，例如，模拟电路元件对温度的敏感性，等等。
基于Matlab设计的多信道接收器——电子科技大学《软件无线电》课程设计
经多相滤波后，每个子信道的带宽为3.2MHz，它们的频谱图分别如下：
分析：图二为子信道y0-y3的频谱图，对应频率范围：-24.0MHz——-11.2MHz， 图二中没有输出信号分量；图三为子信道y4-y7的频谱图，对应输出频率范围为11.2MHz——1.6MHz，其中y4中的输出分量是-9.1MHz和-8.3MHz，y5的输出分量是 -7.7MHz和-6.9MHz，y6的输出分量是-4.75MHz（不明显）、-4.0MHz和-3.25MHz， y7的输出分量是-0.5MHz和0.5MHz；
%%%%%%%%%%%将输入x分为M个部分,逆序排列%%%%%%%%%%%%
x0=downsample(x,M,0); x1=downsample(x,M,15); x2=downsample(x,M,14); x3=downsample(x,M,13); x4=downsample(x,M,12); x5=downsample(x,M,11); x6=downsample(x,M,10); x7=downsample(x,M,9); x8=downsample(x,M,8); x9=downsample(x,M,7); x10=downsample(x,M,6); x11=downsample(x,M,5); x12=downsample(x,M,4); x13=downsample(x,M,3); x14=downsample(x,M,2); x15=downsample(x,M,1);
令：
则转移函数为：
其中，
称为
的多相分量；
若将
中的z换成������������������，
哈工大软件无线电实验报告
一．实验要求
本实验为演示实验，需观察实验现象，实验要求： 1. 掌握采样原理； 2. 对仿真信号绘制波形图，得到信号频谱。
二．实验原理——SDR 采样理论
1. 低通采样（Shannon Sampling Theory）
f s 2 f max
2. 带通采样（Bandpass Sampling）
实验二 信道化发射信号仿真实验
一．实验要求
本实验为仿真实验，实验要求： 1. 掌握信道化滤波器设计与使用； 2. 使用 Matlab 对信道化原理仿真； 3. 掌握多相滤波器结构； 4. 讨论多信道方法的频带使用率以及原理中快速算法的实现。
二．实验理
多频带的信道化发射以其优越的频谱利用率， 广泛应用于雷达和众多电子工程中， 多相 滤波器结构在实时滤波、并串转换等中也得到广泛的应用。 信道化发射信号原理：对于输入的多个频带信号，基本信道化的结构如图 1 所示。
图 1 基本信道化结构
图 2 变化后的信道化处理
基于多相滤波器结构，变化后的信道化处理如图 2 所示。
三．实验内容及结果分析
1. 产生 16 个复信号：1~16Hz
2. 设计滤波器
应用 Parks-McClellan optimal FIR filter order estimation 设计滤波器，其中参数为：内插 I=16，内插前采样率 50Hz，滤波器通带截止频率 16Hz, 阻带起始频率 25Hz；通带和阻带期 望的幅度分别为 1 和 0，起伏为 1dB 和 40dB。 利用 Parks-McClellan 方法得到的频率向量 fo，幅度向量 ao 和权值 w 设计最终使用的滤 波器系数（可以使用 remez 方法，得到指定阶数的滤波器系数） 。 显然，滤波器系数个数应为 16 的整数倍。 实验中所设计的滤波器特性如下：
adalm-pluto matlab 例程-概述说明以及解释
adalm-pluto matlab 例程-概述说明以及解释1.引言1.1 概述概述部分的内容可以写作如下：引言部分是文章的开端，用于介绍论文的主题和意义。
本文旨在探讨Adalm-Pluto与Matlab的连接以及Adalm-Pluto Matlab例程的使用方法。
Adalm-Pluto是一种可编程的软件无线电（SDR）平台，广泛应用于无线通信、无线电频谱监测、物联网和无线电教育等领域。
Matlab作为一种强大的数学软件工具，具有丰富的数据处理、算法设计和可视化分析能力。
将Adalm-Pluto与Matlab相结合，可以进一步发挥两者的优势，提高数据处理和系统设计的效率。
本文的主要结构包括引言、正文和结论三部分。
引言部分首先概述了文章的背景和目的。
接下来的正文部分将从Adalm-Pluto的介绍开始，介绍其硬件特点和功能。
然后，详细介绍了Matlab与Adalm-Pluto的连接方式，包括通过USB和无线连接两种方法。
最后，重点介绍了Adalm-Pluto Matlab例程的使用方法，包括数据采集、信号处理和系统设计等方面。
文章的结论部分总结了本文的主要内容，强调了Adalm-Pluto与Matlab相结合在无线通信和无线电教育等领域的应用前景。
同时，还对Adalm-Pluto Matlab例程的未来发展方向进行了展望。
最后，文章以结束语作为结尾，希望本文能够为读者提供一些有关Adalm-Pluto与Matlab的实践经验和参考。
1.2文章结构1.2 文章结构本文主要介绍了adalm-pluto与Matlab的连接以及Matlab上adalm-pluto的例程的使用方法。
具体结构如下：第2节将详细介绍adalm-pluto的概述，包括其基本特点、硬件配置和工作原理。
我们将了解adalm-pluto作为软件定义无线电（SDR）平台的重要性，以及其在射频领域的广泛应用。
第3节将重点讲解Matlab与adalm-pluto的连接方法。
软件无线电实验matlab程序
附录(1)、SDR低通采样理论%parametersfs = 5e4;%采样频率注意fs至少是f的2倍或2倍以上f = 5e3;%信号的频率N = 1024;n = 1:1024;t= n/fs;deta_t = 1/fs;%采样间隔%signalssn = sin(2*pi*f*t);%时域采样后的信号w = hanning(N);%加汉宁窗sn1 = sn.*w';%加窗减少频率泄露Sn_fft = fft(sn1);%频域信号Sn = abs(Sn_fft);%figuresfigure(1);plot(t,sn);title('s(n)的时域波形');grid on;xlabel('t');xlim([0 0.021]);ylabel('幅度');figure(2);xax = [-N/2:N/2-1]/N*fs/1000;%将f轴单位变为kHzplot(xax,20*log10(fftshift(Sn)));title('s(n)的频域波形');grid on;xlabel('f(kHz)');ylabel('幅度(dB)');ylim([-200 100]);(2)、SDR带通采样理论%parametersf = 1.3e3;%基频fo = 100e6;%载波频率N = 1024;n = 1:1024; fs = 4e3;%采样频率t= 0:0.001:2;st = cos(2*pi*(f+fo)*n/fs);%带通抽样信号St1 = fft(st);St = abs(St1);st0 = cos(2*pi*(f+fo)*t);St0 = fft(st0);St2 = abs(St0);figure(1);plot(n,st);title('带通抽样信号时域图像');xlabel('n');xlim([0 1026]);ylabel('幅度');grid on;figure(2);plot(t,st0);title('带通原始信号时域图像');xlabel('t');ylabel('幅度');grid on;figure(3);xax =[-N/2:N/2-1]/N*fs/1000;plot(xax,20*log10(fftshift(St)));title('带通信号抽样后频域图像');xlabel('f(kHz)');ylabel('幅度(dB)');grid on;(3)、频率调制信号%频率调制%parametersfs = 1e3;%抽样频率ts = 1/fs;%采样率N = 1024;n = 0:1:N-1;t = n/fs;kf = 100.1;A = 100;Aw = 10;fc = 10000;%载波频率fm = 10;%调制信号频率mf = kf*Aw/(2*pi*fm);%signalss1 = A*cos(2*pi*fc*t).*cos(mf*sin(2*pi*fm*t)); s2 = A*sin(2*pi*fc*t).*sin(mf*sin(2*pi*fm*t));s = s1-s2;%频率调制信号S0 = fft(s);S1 = fftshift(S0);S_0 = abs(S1);%figuresfigure(1);plot(t,s);title('频率调制信号时域图像');xlim([0,0.5]);xlabel('t(s)');ylabel('幅度');grid on;figure(2);xax = [-N/2:N/2-1]/N*fs/1000;plot(xax,20*log10(S_0));title('频率调制信号频域图像');xlabel('f(kHz)');ylabel('幅度(dB)');grid on;(4)、幅度调制信号、DSB信号、USB信号%parametersclear all;clc;close all;fs = 1e3;%抽样频率ts = 1/fs;%采样间隔N = 10000;n = 1:N;t = n/fs;ma = 0.5;fc = 200;%载波频率x = cos(2*pi*fc*t);%生成载波fm = 1;%调制信号频率m = cos(2*pi*fm*t);%基带信号%signalss = (1+ma*m).*x;%幅度调制信号s1= m.*x;%双边带信号DSBs2 = cos(2*pi*(fc+fm)*t);%SSB信号中的USBS1 = fft(s);S_1 = abs(S1);S2 = fft(s1);S_2 = abs(S2);S3 = fft(s2);S_3 = abs(S3);%figuresfigure(1);plot(t,s);title('幅度调制信号时域图像');xlim([0 1.5]);xlabel('t(s)');ylabel('幅度');grid on;figure(2);xax = [-N/2:N/2-1]/N*fs/1000;plot(xax,20*log10(fftshift(S_1)));title('幅度调制信号频域图像');xlim([0.19 0.21]);xlabel('f(kHz)');ylabel('幅度(dB)');grid on;figure(3);plot(t,s1);title('双边带调制信号时域图像');xlim([0 1.5]);xlabel('t(s)');ylabel('幅度');grid on;figure(4);xax1 = [-N/2:N/2-1]/N*fs/1000;plot(xax1,20*log10(fftshift(S_2)));title('双边带调制信号频域图像');xlim([0.19 0.21]);xlabel('f(kHz)');ylabel('幅度(dB)');grid on;figure(5);plot(t,s2);title('USB调制信号时域图像');xlim([0 1.0]);xlabel('t(s)');ylabel('幅度');grid on;figure(6);xax1 = [-N/2:N/2-1]/N*fs/1000;plot(xax1,20*log10(fftshift(S_3)));title('USB调制信号频域图像');xlim([0.19 0.21]);xlabel('f(kHz)');ylabel('幅度(dB)');grid on;(5)、2ASK信号%2ASK信号的时域和频域分析clear all;clc;close all;%清除一切。
软件无线电Labview仿真实验报告
《通信系统原理实验》课程研究性学习手册一、实验任务：1、在LabVIEW 平台上完成一个AM 演示程序，实现简单的AM 调制。
2、实现一个基于LabVIEW 和NI-USRP 平台的调频收音机，并正确接收空中的调频广播电台信号。
二、理论分析：1.幅度调制幅度调制（Amplitude Modulation ，AM ）是一种模拟线性调制方法。
频域上，已调信号频谱是基带调制信号频谱的线性位移；在时域上，已调波包络与调制信号波形呈线性关系。
AM 调制的载波信号通常是高频正弦波，作为载体来传递信源信号中的信息。
调制结果是一个双边带信号，中心是载波频率，带宽是原始信号的两倍。
调制信号的数学表达式为：()()()()()()000c o s c o sθωθω+++==t t f t A t c t m t s c c AM （1.1） 式中，)(t m 是调制信号，其直流分量为0A ，交流分量为；)(t c 是载波信号，其为角频率为c ω、初始相位为0θ的余弦信号。
从式1.1我们能够得出幅度调制的已调信号就是是)(t m 和)(t c 的乘积。
为了实现)(t m 能够对载波信号的幅度实现线性调制，)(t m 应该包含直流分量来保证0)(≥t m ，也就是 ()0m a x A t f ≤ （1.2）这样的话才能够保证()t s AM 的包络完全在时间轴上方，如图1所示。
根据式（1.2），为避免产生“过调幅”现象而导致包络检波的结果严重失真，兹定义一个重要参数：10≤=A A mAM β （1.3）式中，称AM β为调幅指数，或调幅深度；m A 代表信源信号()t f 的最大幅值。
一般AM β不超过0.8。
下面对AM 调制在频域上进行分析。
对于式1.1，我们能够直接通过傅里叶变换得到其频域表达式，如式1.4所示。
()()()[]()()[]22220000000θθωωωωδπωωωωδπωj j AM e F A e F A S -+-++++=-（1.4）频谱如图2所示：图2 调幅信号频谱由于软件无线电的核心思想是对天线感应的射频模拟信号尽可能地直接数字化，将其变为适合于数字信号处理器（DSP ）或计算机处理的数据流，然后由软件（算法）来完成各种各样的功能，使其具有更好的可扩展性和应用环境适应性，故而对信源信号的各种调制与解调过程都是在数字域实现的。
软件无线电基本试验
实验一软件无线电基础一、实验目的熟悉软件无线电实验平台，了解软件无线电平台的软硬件处理通信任务的过程，学会软件无线电的基本设计方法和开发工具软件使用方法。
二、实验内容用软件无线电实验平台和LabVIEW软件创建一个调频无线接收器；创建一个自定义LabVIEW用户界面，配置USRP，用LabVIEW设计无线通信系统原型。
三、实验仪器1 USRP实验平台一台2计算机一台四、实验原理1软件无线电平台原理无线通信测试创新论坛对软件无线电(SDR)的定义：〃无线电的一些或全部的物理层功能由软件定义。
〃软件无线电参考了这样一个技术：在通用硬件平台上运行软件模块，用于实现无线通信功能。
结合USRP通用软件无线电硬件和模块化软件的优势，提供了满足多功能需求且灵活性强的快速通信原型平台，适用于物理层设计、算法验证、多标准无线系统、无线信号录制与回放、通信情报等应用。
图1.软件无线电平台构架2软件无线电实现的数字通信系统2.1 典型的数字通信系统一个典型的数字通信系统包括：发射机、接收机和通信信道。
图3展示了一个数字 通信系统的通用组件。
放在第一行是发射机，包含信源编码、信道编码、调制、上变频 模块。
第二行是接收机由下变频、匹配滤波器、均衡器、解调、信道译码和信源译码模 块组成。
数字处理.I合成船1 __________________控制软件接收器控制器□ RF/IF □混含信号 口数亭部分 敬宇处理骐劭器1 放大可配置 H 接收器T N 口D/A可配置 发射端BEF E图2数字通信系统框图酬三P8 -KU-PCI3 Zu 」ncl6-SUUE J=Jj_n r BitsKEP£a白Bu」LTLT日均u.EH-np岂=O_£^CQWCSQQ2.2NI USRP无线通信实验系统图3 NI-USRP无线实验系统硬件、软件平台1） NI USRP硬件平台图4 NI-USRP硬件平台前面板射频信号输入到SMA连接器，USRP硬件平台通过直接变频接收机中的混频操作，产生同相正交（I/Q）基带信号，再经过一个2通道，速率为100 MS/s的14位模数转换器 依口0采样。
matlab操作实验报告
matlab操作实验报告一、实验目的1、学会使用matlab建立.m文件。
2、学会二机五节点的潮流计算计算原理。
3、学会使用matlab命令窗编写程序。
4、学会潮流计算的matlab的程序。
5、学会matlab中simulink模块库的模型用法。
二、实验原理潮流分布（1）描述电力系统运行状态的技术术语，它表明电力系统在某一确定运行方式和接线方式下，系统从电源经网络到负荷各处的电压、电流、功率的大小和方向的分布情况。
（2）电力系统潮流分布主要取决于负荷的分布、电力网参数、以及和供电电源之间的关系。
潮流计算的方法1）建立描述电力系统运行状态的数学模型；2）确定解算数学模型的方法；3）制定程序框图，编写计算机计算程序，并进行计算；4）对计算结果进行分析。
对图1所示电力系统拓扑图进行潮流计算该拓扑图为2机5节点的系统应用Matpower计算潮流技巧的核心在于输入好三个矩阵和部分参数，清晰的知道输入参数、矩阵中每一个元素的含义。
参数一、MATPOWER CaseFormat:Version2mpc.version='2';解释：目前普遍采用2形式的算法。
参数二、system MVA base mpc.baseMVA=100;解释：采用有铭值图1mpc.baseMVA=100；（Matpower只能计算有铭值得网络）矩阵一%%bus data %bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin解释：bus data母线参数也就是我们所说的节点参数，下面逐条注释：1bus number(positive integer)：第一列表示节点的编号（括号里面注释正整数）； 2bus type：第二列表示节点的类型，一般只用得到1、2、3三种节点类型，4类型的节点目前没有接触到。
PQ bus=1PV bu=2reference bus=3isolated bu=43Pd,(MW)：表示负荷所需要的有功功率（所有数据都是正数）（有铭值）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
附录(1) 、SDR 低通采样理论%parametersfs = 5e4;%采样频率注意或 2 倍以上f = 5e3;%信号的频率N = 1024;n = 1:1024;t= n/fs;deta_t = 1/fs;% 采样间隔%signalssn = sin(2*pi*f*t);% 时域采样后的信号w = hanning(N);% 加汉宁窗sn1 = sn.*w';% 加窗减少频率泄露Sn_fft = fft(sn1);% 频域信号Sn = abs(Sn_fft);%figuresfigure(1);plot(t,sn);title('s(n)的时域波形');grid on;xlabel('t');xlim([0 0.021]);ylabel('幅度');figure(2);xax = [-N/2:N/2-1]/N*fs/1000;% 将f 轴单位变为kHz plot(xax,20*log10(fftshift(Sn))); title('s(n)的频域波形');grid on;xlabel('f(kHz)'); ylabel('幅度(dB)');ylim([-200 100]);(2) 、SDR 带通采样理论%parametersf = 1.3e3;% 基频fo = 100e6;% 载波频率N = 1024;fs = 4e3;%采样频率t= 0:0.001:2; st = cos(2*pi*(f+fo)*n/fs);% 带通抽样信号St1 = fft(st);St = abs(St1);st0 = cos(2*pi*(f+fo)*t);St0 = fft(st0);St2 = abs(St0);figure(1);plot(n,st);title('带通抽样信号时域图像’)；xlabel('n');xlim([0 1026]);ylabel('幅度');grid on;figure(2);plot(t,st0);title('带通原始信号时域图像');xlabel('t');ylabel('幅度');grid on;figure(3);xax =[-N/2:N/2-1]/N*fs/1000;plot(xax,20*log10(fftshift(St)));title('带通信号抽样后频域图像’)；xlabel('f(kHz)');ylabel('幅度(dB)');grid on;(3) 、频率调制信号%频率调制%parametersfs = 1e3;%抽样频率ts = 1/fs;%采样率N = 1024;n = 0:1:N-1;t = n/fs;kf = 100.1;A = 100;Aw = 10;fc = 10000;% 载波频率fm = 10;% 调制信号频率mf = kf*Aw/(2*pi*fm);%signalss1A*cos(2*pi*fc*t).*cos(mf*sin(2*pi*fm*t)); s2fs 至少是f 的2 倍n = 1:1024;A*sin(2*pi*fc*t).*sin(mf*sin(2*pi*fm*t)); s = s1-s2;% 频率调制信号S0 = fft(s);S1 = fftshift(S0);S_0 = abs(S1);%figuresfigure(1);plot(t,s);title(' 频率调制信号时域图像');xlim([0,0.5]);xlabel('t(s)');ylabel('幅度');grid on;figure(2);xax = [-N/2:N/2-1]/N*fs/1000;plot(xax,20*log10(S_0));title(' 频率调制信号频域图像');xlabel('f(kHz)');ylabel('幅度(dB)');grid on;(4) 、幅度调制信号、DSB 信号、USB 信号%parametersclear all;clc;close all;fs = 1e3;%抽样频率ts = 1 /fs ;%采样间隔N = 10000;n = 1:N;t = n/fs;ma = 0.5;fc = 200;%载波频率x = cos(2*pi*fc*t);% 生成载波fm = 1 ;%调制信号频率m = cos(2*pi*fm*t);% 基带信号%signalss = (1+ma*m).*x;% 幅度调制信号s1= m.*x;% 双边带信号DSBs2 = cos(2*pi*(fc+fm)*t);%SSB 信号中的USB S1 = fft(s);S_1 = abs(S1);S2 = fft(s1); S_2 = abs(S2);S3 = fft(s2);S_3 = abs(S3);%figures figure(1); plot(t,s);title('幅度调制信号时域图像'); xlim([0 1.5]);xlabel('t(s)');ylabel('幅度');grid on;figure(2);xax = [-N/2:N/2-1]/N*fs/1000; plot(xax,20*log10(fftshift(S_1))); title('幅度调制信号频域图像'); xlim([0.19 0.21]);xlabel('f(kHz)');ylabel('幅度(dB)');grid on;figure(3);plot(t,s1);title('双边带调制信号时域图像'); xlim([0 1.5]);xlabel('t(s)');ylabel('幅度');grid on;figure(4);xax1 = [-N/2:N/2-1]/N*fs/1000; plot(xax1,20*log10(fftshift(S_2))); title(' 双边带调制信号频域图像'); xlim([0.19 0.21]);xlabel('f(kHz)');ylabel('幅度(dB)');grid on;figure(5);plot(t,s2);title('USB 调制信号时域图像'); xlim([0 1.0]);xlabel('t(s)'); ylabel('幅度');grid on;figure(6);xax1 = [-N/2:N/2-1]/N*fs/1000;plot(xax1,20*log10(fftshift(S_3))); title('USB 调制信号频域图像'); xlim([0.19 0.21]);xlabel('f(kHz)'); ylabel('幅度(dB)');grid on;(5) 、2ASK 信号%2ASK 信号的时域和频域分析clear all;clc;close all;% 清除一切。
m = [1 1 1 0 0 0 1 0 1 1 0 1];% 数字信号的序列(自己随便产生的)Lm = length(m);% 序列的长度F = 200 ;%数字信号的带宽('1'或'0'的带宽是200Hz)f = 800;% 正弦载波信号的频率A = 1;% 载波的幅度Q = f/F;%频率比，即一个码元宽度中的正弦周期个数%signalscarry1 = repmat(A*sin(2*pi*f*t),1,Q);% 一个码元宽度内的正弦载波信号Lcarry1 = length(carry1);% 一个码元宽度内的信号长度carry2 = kron(ones(size(m)),carry1);% 载波信号ask = kron(m,carry1);% 调制后的ASK 信号N = length(ask);%ASK 信号的长度tau = (0:N-1)/(M-1)/f;% 时间Tmin = min(tau);% 最小时间Tmax = max(tau);% 最大时间T = on es(size(carry1));% —个数字信号’1' dsig = kron(m,T);Ask = fft(ask); Ask0 = abs(Ask);%imagesfigure(1);plot(tau,dsig);grid on;axis([Tmin Tmax -0.2 1.2]); title('数字信号序列的图像');figure(2); plot(tau,carry2);xlim([0 0.06]);grid on;title('载波时域图像');figure(3);plot(tau,ask);grid on;axis([Tmin Tmax -1.2*A 1.2*A]); title('2ASK 信号的时域图像');M = 500;% 一个正弦周期内的采样点数t = (0:M-1)/M/f;% 一个正弦信号周期内的时间figure(4);xax = [-N/2:N/2-1]/N*Lm*F;plot(xax,20*log10(fftshift(Ask0)));xlim([-300 300]);title('2ASK 信号的频域图像');grid on;(6) 、16 路信道化发射机%parameters clc;close all;clear all;a = 1;% 幅度fs = 50;% 内插前采样频率k = 1:16;%16 路信号I = 16;% 内插系数N = 1024;fk = 1:1:N;sk = zeros(16,N);% 创建16 行N 列的零矩阵用来存放16 路时域信号Sk = zeros(16,N);% 创建16 行N 列的零矩阵用来存放16 路频域信号for k = 1:16sk(k,:) = a*(1+0.5*exp(j*2*pi*k*fk/fs));% 输入16 路时域信号endfor k = 1:16Sk(k,:) = fft(sk(k,:));%16 路信号频域endfigure(1);xax = [-N/2:N/2-1]/N*fs;plot(xax,20*log10(abs(fftshift(Sk(:,1:N))))); title('16 原始复信号时域叠加图');xlabel('f(Hz)');xlim([-5 25]);ylabel('幅度');figure(2);h0 = fir1(265,0.5);HW = fft(h0,N);xax1 = [-N/2:N/2-1]/N*fs;plot(xax1,10*log10(abs(fftshift(HW)))); title(' 滤波器响应');xlabel('f(Hz)'); ylabel('幅度');%I 倍内插，即16 倍内插，此时fs 变大N_1 = N*16;st = zeros(16,N_1);St = zeros(16,N_1);for k = 1:16for g = 0:N-1st(k,g*I+k) = sk(k,g+1);endendsum = zeros(1,N_1);for k = 1:16sum = sum+st(k,1:N_1);endfor k = 1:16St(k,:) = fft(sum(1,1:N_1)); endfor k = 1:16St(k,1:N) = St(k,1:N).*HW; endfigure(3);xax3 = [-N_1/2:N_1/2-1]/N_1*fs*16;plot(xax3,sum(1,:));title('内插后信号时域图像’)；xlim([-60 60]);xlabel('t')； ylabel('幅度');figure(4)；xax2 = [-N/2:N/2-1]/N*fs*16;plot(xax2,db(abs(fftshift(St(:,1:N)))));title('内插后信号频域图像’)；xlim([-10 200]);xlabel('f(Hz)'); ylabel('幅度');。