北邮通信原理软件实验报告
通信原理软件实验报告
学 院:信息与通信工程学院 班 级: 班内序号: 学生姓名: 学 号:
实验二 时域仿真精度分析 一【实验目的】
1. 了解时域取样对仿真精度的影响 2. 学会提高仿真精度的方法
二 【实验原理】
一般来说,任意信号 s(t)是定义在时间区间(-无穷,+无穷)上的连续函数,但所有计算机的CPU 都只能按指令周期离散运行,同时计算机也不能处理这样一个时间段。为此 将把 s(t)按区间[-T/2 ,+T/2 ]截短为按时间间隔dert T均匀取样,得到的取样点数为N=T/t
仿真时用这个样值集合来表示信号 s(t)。t反映了仿真系统对信号波形的分辨率,t 越小则仿真的精确度越高。据通信原理所学,信号被取样以后,对应的频谱是频率的周
期函数,其重复周期是1/t。如果信号的最高频率为fH 那么必须有
fH
1 2t
才能保证不发生频域混叠失真,这是奈奎斯特抽样定理。设
则称BS为仿真系统的系统带宽。如果在仿真程序中设定的采样间隔是t,那么不能用 此仿真程序来研究带宽大于这的信号或系统。换句话说,就是当系统带宽 一定的情况下,信号的采样频率最小不得小于 2*Bs,如此便可以保证信号的不失真,在此基础上时域采样频率越高,其时域波形对原信号的还原度也越高,信号波形越平滑。也就是说,要保证信号的通信成功,必须要满足奈奎斯特抽样定理,如果需要观察时域波形的某些特性,那么采样点数越多,可得到越真实的时域信号。
三 【实验步骤】
1.将正弦波发生器模块、示波器模块、时钟模块按图 5.9 方式连接:
2.设置相关参数进行仿真。
四【实验结果与分析】
1.时钟周期设置为0.01S 信号周期为
1S
分析此时每个周期采样10个点显示在示波器上,能很好地显示原信号的波形,基本无失真 2.修改时钟周期为
0.3S
分析:此时每个信号周期只采样10/3个样点,显示在示波器上,波形出现了明显失真,但是由于符合奈奎斯特抽样定理,所以依旧能显示出大概波形 3.修改时钟周期为
0.5S
分析:此时取样周期为信号周期的一半,对于正弦信号而言,也就是每个周期只取样两个点,从零开始取样的话,刚好是所以取样点都为零,所以此时显示波形为一条直线。
实验三 频域仿真精度分析
一【实验目的】
(1)理解 DFT 的数学定义及物理含义。 (2)学会应用 FFT 模块进行频谱分析。
(3)进一步加深对计算机频域仿真基本原理以及方法的学习掌握。
二【实验原理】
在通信系统仿真中,经常要用有限长序列来模拟实际的连续信号,用有限长序列的 DFT 来近似实际信号的频谱。DFT 只适用于有限长序列,在进行信号的频谱分析时,它的处理 结果会含有一定的偏差。分析可知由于DFT是针对有限长序列所定义的变换方法,因此,
利用这一变换分析信号的频谱前,首先构造了一系列有限长序列xn,又该序列的DFT结果X(k)来表示原信号xat的频谱xaj,也就是用Xk逼近xaj。 上述原理可以用图表示:
基于上述方法所得到的结果只是对原信号频谱xaj的一种近似,也就是说,Xk同
xaj之间存在着幅度偏差,造成这一偏差的原因主要体现在两个方面:1.时域混叠与频
域混叠2.频谱泄露
三【实验步骤】
1、将正弦波发生器 (sinusoid generator)、触发时钟(CLOCK_c)和频谱示波器模块按图 5.13 所示连接。
2、设置相关参数进行仿真。
四【实验结果与分析】
利用这一变换分析信号的频谱前,首先构造了一系列有限长序列xn,又该序列的DFT结果X(k)来表示原信号xat的频谱xaj,也就是用Xk逼近xaj。
上述原理可以用图表示:
基于上述方法所得到的结果只是对原信号频谱xaj的一种近似,也就是说,Xk同xaj之间存在着幅度偏差,造成这一偏差的原因主要体现在两个方面:1.时域混叠与频域混叠2.频谱泄露
三【实验步骤】
1、将正弦波发生器 (sinusoid generator)、触发时钟(CLOCK_c)和频谱示波器模块按图 5.13 所示连接。
2、设置相关参数进行仿真。
四【实验结果与分析】
2.FFT size 40960 窗口类型:1 图二
分析:上述两图的主瓣宽度不同是由于FFT SIZE不同所致,FFT size 越大,得到的频谱的谐波分量越多,频谱主瓣变得很尖锐
图三相对于图二而言是由于窗函数不同,导致有一定的频谱泄露
4.FFT size 40960 窗口类型:2 图四
5.FFT size 40960 窗口类型:4 图五
分析:由图三图四图五可知,而FFT size一致的时候,窗口类型对频谱的影响不太大,主瓣宽度基本一致,幅度基本一样,谐波分量也基本一样。但是相对于图二而言这些都有不同程度的频谱泄露,只是加窗不同,对泄露的处理结果也就不同。
五【思考题】
1、对于同一正弦信号,观察图 5.14、图 5.15 中所示频谱图的不同,分析其原因。
答:这个主要是因为FFT size的不同引起的,窗口宽度加宽的时候,就不会有更多的谐波分量被滤掉,导致频谱高频谐波分量的增加。
2. 观察图 5.15、图 5.16 所示频谱图的不同,解释其原因。
答:频谱的主瓣宽度增加,高频谐波分量减少。原因就是,采用了不同的窗函数,不同的窗函数对信号的滤波特性是不一致的。
3.将 FFT 模块中的参数 Type of window 改成 2 和 4,观察仿真结果的变化,解释其原因。 答:换成窗口2之后频谱变得越平滑,主要是因为滤去了更多的谐波分量,窗口4相对于窗口3也变得相对光滑,但是没有窗口2明显。
实验五 取样和重建
一【实验目的】
(1)了解取样定理的原理,取样后的信号如何恢复原信号
(2)了解取样时钟的选取。
二【实验原理】
数字信号是通过对模拟信号进行采样、量化和编码得到的,模拟信号是时间和幅度都连续的信号,记作 x(t)。采样的结果是产生幅度连续而时间离散的信号,这样的信号常被称为采样数据信号。 原理如下:
低通采样定理:如果采样频率fS大于2fh(fh是带限信号的截止频率),那么带限信号就可以无差错地通过其采样信号恢复。
模型:
具体原理见讲义。
在满足采样定理条件的情况下,初始输入信号可以从这些抽样值中恢复出来
三【实验步骤】
1. 脉冲信号产生器(Pulse generator, 来自 Scicom_sources 元件库)、正弦波发生器
(sinusoid generator)、模拟低通滤波器(analog low pass filter)、直流发生器 DC、触发时钟 (CLOCK_c)、乘法器、示波器模块(MScope)、频谱示波器(FFT)模块按下图所示连接。
设置相关参数进行仿真。
四【实验结果与分析】
时域仿真波形:
FFT(2)重建信号的频谱:
FFT(1)取样信号频谱:
取样信号放大频谱图:
分析:由上述各图图可以看出,当取样脉冲宽度为0.00001,也就是非常小的时候,重建信号基本能够反映原信号,周期和原信号相同,只是幅值稍有减小,几乎没有衰减。
第二次验证:
实验参数的设置,脉冲发生器高电平时间0.1,常数5;
时域仿真波形
FFT(1)取样信号频谱
取样信号频谱放大:
FFT(2)重建信号频谱:
分析:当取样脉冲宽度较大时(此时为01),重建信号已经无法反映原信号,周期差别太大,幅度衰减明显。
五【思考题】
1. 分析分析两次实验中FFT(1)取样信号频谱有何区别,并解释其原因。
答:前者幅度不变,后者随着频率增大有衰减。因为随着取样脉冲宽度的增大,其频域Sa函数的衰减将越明显。
2. 观察第一次实验信号时域采样后,其对应的频谱周期延拓现象,其周期是多少? 答:由图可知:频域周期延拓后,周期为1/0.4=2.5
3. 观察并对照两组参数设置下出现的不同仿真现象,结合信号与系统相关理论分析不同采样函数占空比对信号频谱的影响。
答:占空比越大,频域Sa函数的衰减将会越明显,使得取样信号功率谱衰减更大。
实验七 SSB调制与解调(模块)
一【实验目的】
(1)了解产生 SSB 调制的基本原理。
(2)了解 SCICOS 中的超级模块。
(3)了解利用相干解调法解调幅度调制信号的方法。
二【实验原理】
SSB 调制
SSB AM 产生方法一:
SSB AM 产生方法二:
单边带调制信号表达式为:
SSB 解调
用相干解调或同步解调来还原幅度调制信号。其解调框图如下:
如图 上图5.3所示,载波应该提取自输入信号,通过平方环法或 COSTAS 环方法提取。由于这次实验是验证解调方法,假定已经获得了解调所用的载波的频率,所以直接使用调制端正弦波发生器产生的载波信号充当解调载波。 三【实验步骤】
SSB 调制
1、将正弦波发生器( sinusoid generator )、组合希尔伯特变换器(来自
Scicom_signalprocess 元件库)、组合移相器(来自 Scicom_signalprocess 元件库)、加法器
模块、乘法器模块、触发时钟(CLOCK_c)、示波器模块(MScope)、和频谱示波器(FFT)模块按下图连接。
2、按照实验教程甚至相关参数进行仿真。
SSB 解调
1、SSB 解调实验步骤与实验七中 DSB 解调步骤相同,只是将生成 DSB 超级模块(Modulator)换成生成 SSB 超级模块(Modulator)。
2、按照实验教程设置相关参数,进行仿真。
四【实验结果与分析】
SSB 调制
时域波形图
频域波形图
分析:频域波经放大后可知调制后的频率约为15Hz,这与题设基本一致,采取了下边带调制,原信号频率为5Hz,载波频率为20Hz,经下边带调制后为15Hz。
SSB 调制
时域波形图
频域波形图
分析:由图可知,信号经下边带解调后的频率约为3Hz,解调后的幅度为原幅度的0.5倍,这与解调原理刚好相吻合。
五【思考题】
1、SSB 信号的波形和频谱的特点是什么?
答:答:频谱只有上边带或者只有下边带,传输带宽较窄,信道利用率较高。
2、实验步骤 5 的参数之间有什么关系?为什么?改变参数值,配合实验加以解释。 答:两个参数乘积应为整数,使得希尔伯特变换较为准确
编程实验
一【实验要求】
假设基带信号为m(t)sin(2000t)2cos(1000t),载波频率为20kHz,仿真出SSB信号,观察已调信号的波形及频谱。
二【程序代码】
clear all
exec t2f.sci;
exec f2t.sci;
fs=800; //采样速率
T=200; //截短时间
N=T*fs; //采样点数
dt=1/fs; //时域采样间隔
t=[-T/2:dt:T/2-dt]; //时域采样点
df=1/T; //频域采样间隔
f=[-fs/2:df:fs/2-df]; //频域采样点数
fm1=1; //待观测正弦波频率,单位KHz,下同
fm2=0.5; //待观测余弦波频率
fc=20; //载波频率
//以上为初始化参数设置
m1=sin((2*%pi)*fm1*t); //待观测正弦波部分
M1=t2f(m1,fs); //傅里叶变换
MH1=-%i*sign(f).*M1; //希尔伯特变换
mh1=real(f2t(MH1,fs)); //希尔伯特反变换
m2=2*cos((2*%pi)*fm2*t); //待观测余弦波部分
M2=t2f(m2,fs); //傅里叶变换
MH2=-%i*sign(f).*M2; //希尔伯特变换
mh2=real(f2t(MH2,fs)); //希尔伯特反变换
s1=(1+(m1+m2)/abs(max(m1+m2))).*cos((2*%pi)*fc*t); //AM信号时域表达式
S1=t2f(s1,fs); //AM信号频域表达式
s2=(m1+m2).*cos((2*%pi)*fc*t); //DSB-SC信号时域表达式
S2=t2f(s2,fs); //DSB-SC信号频域表达式
s3=(m1+m2).*cos((2*%pi)*fc*t)-(mh1+mh2).*sin((2*%pi)*fc*t); //SSB信号时域表达式,以上边带为例
S3=t2f(s3,fs); //SSB信号上边带频域表达式
//以上是仿真计算部分
//以下为绘图部分
//AM信号
xset('window',1)
plot(f,abs(S1))
title('AM信号频谱')
xlabel('f')
ylabel('S(f)')
mtlb_axis([-25,25,0,max(abs(S1))]);
xset('window',2)
plot(t,s1)
title('AM信号波形')
xlabel('t')
ylabel('s(t)')
mtlb_axis([-3,3,-3,3]);
//DSB-SC信号
xset('window',3)
plot(f,abs(S2))
title('DSB-SC信号频谱')
xlabel('f')
ylabel('S(f)')
mtlb_axis([-25,25,0,max(abs(S2))]);
xset('window',4)
plot(t,s2)
title('DSB-SC信号波形')
xlabel('t')
ylabel('s(t)')
mtlb_axis([-1,4,-3,3]);
//SSB信号(以上边带为例)
xset('window',5)
plot(f,abs(S3))
title('SSB信号频谱')
xlabel('f')
ylabel('S(f)')
mtlb_axis([-25,25,0,max(abs(S3))]) xset('window',6)
plot(t,s3)
title('SSB信号波形')
xlabel('t')
ylabel('s(t)')
mtlb_axis([0,6,-3,3])
三【实验结果】
分析:对于AM信号,调制系数必须不大于1。观察时域波形即可发现,可观察出明显的载波信号,而且有包络信号,且为调幅信号,而观察频谱图即可发现符合代码中编写的AM频域特点。
对于DSB信号:波形会有明显的相位翻转的现象
对于SSB信号:仅取了DSB信号的半边,符合单边带调制的理论。
实验十二 ASK调制与解调(模块)
一【实验目的】
了解幅度键控(ASK)调制与解调的基本组成和原理。
二【实验原理】
用数字基带信号去控制正弦型载波的幅度称为振幅键控(ASK)。2ASK 是指二进制振幅键控又名 OOK,它以单极性不归零码序列来控制正弦载波的开启与关闭。其产生框图:
频谱分析图为:
解调:在加性高斯白噪声信道条件下,OOK信号的解调方法有相干解调和非相干解调。两种解调方法的原理框图如下
三【实验步骤】
调制
1、将正弦波发生器(sinusoid generator)、二进制随机数产生器(binary random generator, 来自 Scicom_sources)、乘法器模块、触发时钟(CLOCK_c)、示波器模块(MSCOPE)、按图
5.87 连接。
2、设置正弦波模块 Source Signal,产生频率为 1HZ 的信号。设置二进制序列产生器的
时钟频率,产生周期为 2s 的二进制序列。
3、设置相关参数进行仿真
解调
方法一:
1、将正弦波发生器(sinusoid generator)、二进制随机数产生器(binary random
generator)、乘法器模块、整流器(rectifier,修改自mathematical expressions模块)、模拟低通滤波器(analog low pass filter)、比较判断模块(switch,来自Branching元件库)、触发时钟(CLOCK_c)、示波器模块(MSCOPE)、按下图连接。其中整流器由mathematical expression 模块修改而来。
2、设置相关参数进行仿真
方法二:
1、将正弦波发生器(sinusoid generator)、二进制随机数产生器(binary random
generator)、乘法器模块、组合移相器器、低通滤波器(analog low pass filter)、触发时钟(CLOCK_c)、比较判断模块(switch)、示波器模块(MSCOPE)按图5.90连接。
2、由于是验证解调原理的实验,所以解调载波直接由调制载波充当。适当调节移相器的相
位,和低通滤波器的截止频率,使得解调结果正确。
3、设置相关参数进行仿真
四【实验结果与分析】
调制
解调
五【思考题】
1、从频域和时域分析,rectifier/analog low pass filter两个模块的作用。
答:rectifier为整流器,把双极性信号转化为单极性信号。Analog low pass filter负责包络检波,滤除高频成分。
2、MASK 调制解调模型如何构建?
答:MASK信号调制时首先将二进制信号变为M进制幅度序列,然后与二进制ASK信号相同;在解调时,信号乘上载波后,进行一次积分,然后进行判决即可。
设计类实验二 线路码型HDB3编码
一【实验目的】
1、了解几种常用线路码型及其编码规则。
2、掌握HDB3码的编码原理及其SCILAB实现。
3、学会使用HDB3码编码模块及其调试。
二【实验原理】
HDB3码的编码规则:
(1) 将消息代码变换成AMI码;
(2) 检查AMI码中的连0情况,当无4个以上的连0传时,则保持AMI的形式不变;若出现4个或4个以上连0时,则将1后的第4个0变为与前一非0符号(+1或-1)同极性的符号,用V表示(+1记为+V,-1记为-V
(3) 检查相邻V符号间的非0符号的个数是否为偶数,若为偶数,则再将当前的V符号的前一非0符号后的第1个0变为+B或-B符号,且B的极性与前一非0符号的极性相反,并使后面的非0符号从V符号开始再交替变化。
举例如下:
代码 1 0 1 0 1 1 0 0 0 0 0 1 1 0 0 0 0 1
HDB3
码 +10 -1 0 +1 -1 0 0 0 -1 0 +1 -1 +1 0 0 +1 -1 -V +B +V
HDB3码的特点如下:
(1) 基带信号无直流成分,且只有很小的低频成分;
(2) 连0串符号最多只有3个,利于定时信息的提取;
(3) 不受信源统计特性的影响。
虽然HDB3码的编码规则比较复杂,但译码却比较简单。HDB3是CCITT推荐使用的码型之一。
三【实验步骤】
1、设计过程中,由于HDB3是由AMI码衍生而来,故可以将模块分为AMI编码和B/V码编码两部分,然后再由超级模块将两部分进行组合即可。最后按照“模块性能测试”中提供的测试例程对模块进行测试。成功运行后,将结果与测试例程的结果进行对比,然后进行调试和优化等后续工作。
2、按照下图所示连接各模块
3、设置相关参数进行仿真
四【实验结果与分析】
编码图
分析:对照上图和HDB3的编码特点和检验规则,可知上图符合HDB3的编码规则。 五【思考题】
1、观察HDB3码的频谱图,分析其与原码频谱的不同。
答:原码为不归零码,其频谱宽度为Rb,HDB3码功率谱高频及低频分量很少,能量主要集中于二分之一码速附近。
2、HDB3码译码如何仿真实现,请给出SCICOS模块连接的工作原理图,以及相关波形图。
工作原理图
解码波形图
3、线路码在数字通信系统中起什么作用?除了HDB3码,还有哪些线路码?试列出各自的优缺点以及适用场合。
答:作用:去除直流分量,减少低频及高频分量,有利于正确判决,有利于观察误码情况,降低误码率。
除了HDB3码一歪,还有NRZ码,RZ码,差分码,AMI码,MANCHESTER码等。
数字基带系统的仿真(综合)
一【实验目的】
1、了解在理想限带及加性白高斯噪声信道条件下数字基带系统的基本原理和设计方法, 完成在仿真平台上的系统搭建与仿真。
2、进一步熟悉SCICOS 下的复杂系统设计。
3、运用工具库现有的通信工具模块搭建系统,利用其功能全面且封装性强的特点,针对数字基带系统进行横向功能分解,使系统设计更加精炼、准确。
4.、深入学习数字基带系统核心传输节点的性能,并掌握眼图示波器的使用方法,观察接收滤波器输出的眼图和功率谱密度,判断系统传输的正确性和精准性,调试以达到最佳传输效果。
二【实验原理】
若使得在接收端抽样时刻码间干扰为零,则系统的合成传递函数必须满足以下条件:
j2ft0G(f)C(f)G(f)Xfe T fWR升余
即要满足
TfCfRf2ft0
GTfCfGRfX升余f
f分别是发送滤波器、信道及接受滤波器的相频特性,t0式中的Tf、Cf、R
一时间延迟,t0是tTtCtR,其中tT、tC、tR分别是发送、信道、接收滤波器引入的时延,W为升余弦滤波器的截止频率。
在接收端抽样时刻无码间干扰条件下,引起误码的是加性噪声,此时,最佳接收的接受滤波器应匹配于所接收的确定信号,使接收端抽样时刻的信噪比最大。
设:限带信道是理想低通特性,并设信道不引入时延(tC0),
1fWfCf 0fW2W
则接收到的确定信号的频谱仅取决于发送滤波器的GT(f)特性,所以接收滤波器的GR(f)应与发送滤波器GT(f)共轭匹配,即
GRfGT*fej2ft0GRfej2ftR
这样,在理想限带信道情况下,既要使收端抽样时刻的抽样值无码间干扰,又要使得在抽样时刻抽样值的信噪比最大,则
2j2ft
X升余fGT(f)GR(f)GTfe
GTfX升余fej2ftT
综上所诉,数字PAM信号通过限带信道、并受到加性噪声干扰的情况下,在限带信道是理想低通条件下的最佳基带传输的发送及接收滤波的设计是:总的收发系统的传递函数要符合无码间干扰基带传输的升余弦特性;且又要考虑在抽样时刻信噪比最大的收、发滤波共轭匹配的条件。
可得无码间干扰的条件下,其系统框图如图1:
三【实验步骤】
方案一、编程实现理想基带传输系统的升余弦特性 1、按下图所示将各模块连接成数字基带传输系统。
用时钟控制产生二进制随机序列作为高斯白噪声源,同时分频后用于产生双极性二进制序列作为原信号。首先经过发送滤波器,然后经过一低通滤波器,即模拟限带信道,然后与噪声相加。在经过一接收滤波器,最后经过采样S/H和判决后输出解调信号。
2、 在同一参数被不止一个模块调用的情况下可以统一利用一个符号表示,不过需要在 Diagram菜单中的Context中对这个符号进行预先定义。
在本次系统设计中,需要在【Diagram】中的“Context”中进行内容设置。具体设计见讲义。 3、“sending filter”与“receiving filter”的实现 :
设置相关参数进行仿真
方案二、利用根升余弦滤波器模块实现基带传输系统的升余弦特性 1、 按照图下图搭建数字基带系统
2、【Diagram】中的“Context”中进行内容设置 3、设置相关参数进行仿真
四【实验结果与分析】 方案一
信号经过了取样、发送、信道噪声、接收滤波和抽样判决等步骤。仿真波形为:
眼图
方案二 仿真波形
采样信号频谱图:(无限宽)
接收滤波器输出的频谱图:
分析:接收滤波器输出的功率谱,为经过升余弦后的接收波形的频谱,可见在低频处有一波峰,为正弦信号的频谱。
原信号(黑)与接收匹配滤波器的输出信号(绿)波形
:
眼图仿真
关于眼图:
( 1 )最佳抽样时刻应 在 “眼睛” 张开最大的时刻。
( 2 )对定时误差的灵敏度可由眼图斜边的斜率决定。斜率越大,对定时误差就越灵敏。 ( 3 )在抽样时刻上,眼图上下两分支阴影区的垂直高度,表示最大信号畸变。 ( 4 )眼图中央的横轴位置应对应判决门限电平。
( 5 )在抽样时刻上,上下两分支离门限最近的一根线迹至门限的距离表示各相应电平的噪声容限,噪声瞬时值超过它就可能发生错误判决。
( 6 )对于利用信号过零点取平均来得到定时信息的接收系统,眼图倾斜分支与横轴相交的区域的大小,表示零点位置的变动范围,这个变动范围的大小对提取定时信息有重要的影响。
五【思考题】
1、 分析α的大小对系统带宽和误码性能的响应。
答:α越大,系统带宽越宽,误码性能越好。
2、 哪个方案更好,为什么?
答:方案二更好,方案二与方案一相比,抗噪能力更强。 附Context内容设置:
*************************t2f***********
function X=t2f(x) //傅式变换
H=fft(x);
X=[H(mtlb_imp(mtlb_a(N/2,1),N)),H(mtlb_imp(1,N/2))]*dt; endfunction
****************f2t**********************
function x=f2t(X) //傅式反变换
S=[X(mtlb_imp(mtlb_a(N/2,1),N)),X(mtlb_imp(1,N/2))]; x=ifft(S)/dt;
endfunction
*****************************************
dt=0.01; //时域采样间隔
L=32; //每个码元周期的抽样点数
M=16; //码元数
N=L*M; //总的抽样点数
Ts=L*dt; //码元间隔
Rb=1/Ts; //码元速率
df=1/(N*dt); //频域采样间隔
T=N*dt; //截断时间
Bs=N*df/2; //系统带宽
alpha=0.5; //滚降系数
t=linspace(-T/2,T/2,N); //频域横坐标
f=linspace(-Bs,Bs,N)+%eps; //时域横坐标
//升余弦信道
hr1=sin(%pi*t/Ts)./(%pi*t/Ts);
hr2=cos(((alpha*%pi)*t)/Ts)./(1-(((2*alpha)*t)/Ts).^2); hr=hr1.*hr2; //升余弦脉冲波形
HR=abs(t2f(hr)); //升余弦脉冲的傅式变换
GT=sqrt(HR);
GR=GT; //最佳系统的发送接收滤波器的傅式变换