小灰灰的FPGA的博客

RS码，通过选择一个合适的生成多项式g(x)，并使得对每个信息段计算得到的码字多项式都是g(x)的倍式，也就是码字多项式除以g(x)的余式为0。decoded=rsdec(code,n,k,genpoly)表示对应code=rsenc(msg,n,k,genpoly)的译码。RS(n,k)码可以由m、n和k三个参数表示，其中m表示码元符号取自域GF(2^m)，n表示码字长度，k表示信息段长度。decoded=rsdec(code,n,k)表示对应code=rsenc(msg,n,k)的译码。

将接收信号矢量r=[r1,r2]投影到M个可能发送的信号矢量上，再根据最大的投影分量判决发送信号。从图中对比，可以看出，仿真得到的误符号率在不同信噪比情况下与理论值均比较吻合，而误比特率在低信噪比时，与理论值误差较大，信噪比较高时，误差明显缩小。上式中nc(t)和ns(t)是加性噪声的两个正交分量，是零均值且互不相关的高斯随机过程，方差为N0/2。θm=2πm/M表示载波的M个可能的相位，用于传送发送信息。若采用Gray编码是，误比特率为Pb=PM/k，k=log2(M)，每个符号传输的比特数。

基带信号sm(t)的频谱被搬移到载波频率fc上，带通信号是一个双边带抑制载波DSBSC的AM信号。n(t)是带通噪声过程，n(t)=nc(t)×cos2πfct-ns(t)×sin2πfct。数字PAM又称为幅移键控ASK，在数字基带PAM中，信号波形即sm(t)=Am×g(t)载波PAM信号的产生，基带信号波形sm(t)与正弦载波cos2πfct相乘。已调信号的频谱为：S(f)=(Am/2)×[G(f+fc)+G(f-fc)]接收信号r(t)=Am×g(t)×cos2πfct+n(t)。

信道信号形成器的作用是把原始基带信号变换成适合于信道传输的基带信号，主要通过码型变换和波形变换来实现，目的是与信道匹配，便于传输，减小码间串扰，利于同步提取和抽样判决。传输信号通过加性高斯白噪声信道awgn，叠加了噪声n(t)，n(t)是功率谱密度为(N0/2)的白色高斯随机过程的一个样本函数。数字基带传输——在某些具有低通特性的有效信道中，特别是传输距离不太远的情况下，数字基带信号可以直接传输，称为数字基带传输。信号相关器将接收到的信号r(t)与两个可能的发送信号s0(t)和s1(t)做互相关。

在一条公路上一端安装固定的基站，另一端完全反射电磁波的墙面，基站到反射墙的距离为d。移动台接收两路径信号，一路为基准发射的信号，另一路是反射墙反射过来的信号。多径的存在是因为发射机和接收机之间建筑物和其他物体的反射、绕射、散射等引起的。无线信道的特性是不确定的，随机变化的。若信道没有频率选择性，则最大的时延扩展Tmax要远远小于符合周期Ts,此时，所有的延迟路径分量到达的时段仅为一个符合时间的一小部分。多径衰落信道的建模和仿真，两个重要的参数是多径扩展和多普勒带宽。E(t)—时刻t移动台接收到的合成信号。

decoded=vitdec(code,treills,tblen,opmode,dectype)表示对码字code进行Viterbi译码，treills表示产生码字的卷积编码器，tblen表示回溯的深度，opmode表示译码器的操作模式，dectype表示译码器判决的类型，比如软判决和硬判决。序列译码法，译码先从码树的起始节点开始，把接收到的第一个子码的n个码元与自始节点出发的两条分支按照最小汉明距离比较，验证差异最小的分支走向第二个节点。概率译码有维特比译码和序列译码，是卷积码最主要的译码方法。

若错误过多，超过纠错码的纠错能力，但能检测出来，则接收端通过反馈信道，要求发送端重新传送有错的消息。发送端发出能够检测错误的码，接收端收到通过信道传来的码后，译码器根据该码的编码规则，判决收到的码序列中有无错误产生。分组码把信息序列以k个码元分组，通过编译器将每组的k元信息，按照一定规律产生r个多余码，作为监督码或校验码，从而输出长为n=k+r的一个码组。通俗来说，信道编码就是在发送端的信息码元序列中，以某种特定的编码规则，加入监督码元，在接收端再利用该规则进行检查识别，从而发现错误，纠正错误。

循环冗余校验码是利用线性编码理论，在发送端根据要传送的k位二进制码序列，以一定的规则产生一个校验用的r位监督码，也就是CRC码，并附加在信息位后面，构成一个新的共n=k+r位的二进制码序列。对于一个给定的(n,k)码，根据G(x)可以生成k位信息的校验码，G(x)称为这个CRC校验码的生成多项式。系统码的生成矩阵G=[Ik,P]，Ik是k×k的单位矩阵，P为一个k×(n-k)的单位矩阵。在发送信息时，用信息多项式G(x)表示，将G(x)左移r位，表示为C(x)×2。输入的信息位X=[X1,X2,……

由于序列较短的伪随机数之间的相关特性较大，对于实时性要求高、信息帧较短的通信系统，矩阵分组式交织性能优于伪随机和半伪随机式交织器。但由于信息帧长度增加，交织长度也增加，采用矩阵分组交织器，交织前后信息序列的不动点增多，伪随机数产生更加均匀，交织前后的序列相关性减小，所以采用随机交织器，半伪随机交织方式作为折中方案。送入交织器后，从存储器里按行输出，送入突发差错的信道，信道输出再送入反交织器，完成交织器的相反变换，即按行写入按列读出。反交织器的输出，阵列存储器中按列读出的信息，其差错规律就变成独立差错。

16QAM调制器，输入二进制数据经过串/并转化和2/4电平变换后得到两路码元宽度增大4倍的双极性四电平码，分别进行正交调制，合成后的信号即为16QAM信号。上式中nc(t)和ns(t)是加性噪声的两个正交分量，是零均值且互不相关的高斯随机过程，方差为N0/2。在矩形信号星座中，M=2^k,其中k为整数，QAM信号等价于在正交载波上的两个PAM信号。QAM信号的波形可以表示成啷个标准正交信号波形φ1(t)和φ2(t)的线性组合。上式中Amc和Ams是承载信息的正交载波的信号幅度；g(t)是信号脉冲。

差分编码的相位调制和解调过程，检测器的接收信号被解调成M个可能发送的相位中的一个。在检测器之后是一个相位比较器，比较相邻信号间隔上已调信号的相位，以便提取信息。在二进制PSK中，信息比特1通过载波相位对前一载波相位180°相移发送，信息比特0则是通过与前一载波相位相同发送。在四相PSK中，相继区间之间的相对相移0°,90°,180°和270°，分别对应信息比特00,01,11和10。在实际中，载波相位是从接收信号通过某些非线性运算提取的，会引入相位模糊。微信公众号上线，搜索公众号。

接收到的信号幅度Ai能够取±d和±3d，最佳幅度检测器要将输入r与4种可能的传输电平比较，并选择欧氏距离上最接近r的幅度电平，并选取最小距离的幅度。传输信号通过加性高斯白噪声信道awgn，叠加了噪声n(t)，n(t)是功率谱密度为(N0/2)的白色高斯随机过程的一个样本函数。其中Am为第m个波形的幅度，信号幅度取4种等间隔值，即Am=(2m-3)d，m=0,1,2,3。信号波形——sm(t)=Am×g(t),0≤t≤T,m=0,1,2,3。g(t)为矩形脉冲，g(t)为1/T的开方，0≤t≤T。

信号在频率表示为V(f)C(f),若信道带宽限于WHz内，当|f|＞W时，C(f)=0,也就是V(f)中频率的绝对值高于W的任何频率分量都不能通过该信道。当α=0时，升余弦频率响应特性，变成一个理想的，带宽为1/(2T)的理想低通滤波器，频率为1/(2T)，也就是所说的Nyquist频率。无失真信道或理想信道——对于所有的|f|≤W,A(f)为常数，且θ(f)是频率的线性函数，即包络延时特性τ(f)是一个常数。对于理想信道，当|f|≤W时，C(f)=1,脉冲x(t)具有频谱特性X(f)=|G(f)|^2。

通过调制，不仅可以进行频谱搬移，把调制信号的频谱搬移到所希望的位置上，从而将调制信号转换成适合于信道传输的已调信号，而且它对系统的传输有效性和可靠性有很多的影响。AM信号的频谱由载频分量和上、下两个边带组成，上边带的频谱结构与原调制信号的频谱结构相同，下边带是上边带的镜像。h(t)=δ(t),即滤波器H(ω)=1为全通网络，调制信号m(t)叠加直流A0后与载波相乘，形成AM信号。AM信号是带有载波的双边带信号，带宽是基带信号带宽fH的两倍，B(AM)=2fH。若调制信号没有直流分量，即m(t)的均值为0。

在DSBSC信号中，已调信号由m(t)和cos(2πfct)相乘(混频)给出。DSB信号的包络不再与调制信号变化规律一致，需采用相干解调，即同步检波。DSBSC信号的带宽是基带信号带宽fH的两倍，即B(DSB)=2fH。DSBSC信号的带宽是基带信号带宽fH的两倍，即B(DSB)=2fH。当y(t)通过带宽为W的低通滤波器时，高频分量被滤除，低频分量被解调。DSBSC信号功率为P(DSB)=1/2Ps,Ps为消息信号功率。DSBSC信号功率为P(DSB)=1/2Ps,Ps为消息信号功率。

DSB信号的上、下两个边带是完全对称的，都携带了调制信号的全部信息，因此仅传输其中一个边带即可，这种只传输一个边带的通信方式称为单边带通信。SSB信号的带宽是DSBSC信号和AM信号带宽的1/2，和消息信号带宽相同，即Bssb=fH，SSB信号的功率为Pssb=Ps/4，Ps为消息信号功率。产生SSB信号最直观的方法是让双边带信号通过一个边带滤波器，保留所需要的一个边带，滤除不要的边带。SSB信号是抑制载波的已调信号，它的包络不能直接反映调制信号的变化，需要采用相干解调。

调频信号的瞬时频率正比于调制信号的幅度，调频信号的解调器必须能产生正比于输入频率的输出电压。已调信号频谱不再是原调制信号频谱的线性搬移，而是频谱的非线性变换，会产生与频谱搬移不同的新的频率成分，故又称为非线性调制。幅度、频率均含调制信息的调幅调频信号，包络检波器将其幅度变化取出，并滤去直流后得到mo(t),其中Kd称为鉴频器灵敏度。频率调制FM，使高频载波的频率按调制信号的规律变化而振幅保持恒定的调制方式，可以看作是载波角度的变化。带宽为Bfm=2(βf+1)W=2×(3+1)×50=400。

随机数是根据一定算法产生的伪随机数，把产生的伪随机数的初始条件(即随机数种子)设为一致，也就是前后两次调用state相同，则产生的加性高斯白噪声结果也一致。具体表现为信号围绕平均值的一种随机波动过程，均值为0，方差是噪声功率的大小。函数awgn(x,snr,sigpower)，把加性高斯白噪声叠加到输入信号x中，snr以dBW的形式制定噪声的功率,输入信号的功率为sigpower,单位dBW。函数awgn(x,snr)，把加性高斯白噪声叠加到输入信号x中，snr以dB的形式制定噪声的功率。

它是时间t的函数，但在任意时刻上观察到的值是不确定的，是一个随机变量。当然，对通信系统的热噪声进行建模时候，通常假设这样的噪声是白色高斯随机过程，即功率谱密度S(f)对所有的f是一个常数。若随机过程的功率谱在某中心频率±f0附加的一个频段内有值，而在频率之外功率谱密度为0，称这个随机过程是带通过程，若通带B<<f0，称为一个窄带过程。宽平稳过程，平稳随机过程的数学期望及方差与t无关，自相关函数只与时间间隔τ有关，R(t1,t1+τ)=R(τ)平稳随机过程的功率谱密度S(f)定义为自相关函数的傅里叶变换。

称为带通信号x(t)的两个正交分量，xc(t)为同相分量，xs(t)为正交分量。Hilbert（希尔伯特）变换可以将实信号转换为复信号，常用于提取信号的相位信息，使得我们能够更全面地分析和处理信号，这在雷达通信领域中都很有用，尤其是在带通信号处理中。一个带通信号x(t)，x(t)的正频域表示式：X+(f)=2u(f)X(f)其中X(f)为x(t)的傅里叶变换，u(f)为单位阶跃函数。信号a(t)为x(t)的包络，θ(t)为s(t)的相位。信号x+(t)称为解析信号或x(t)的预包络。

如果x(n)是M点序列，h(n)是L点序列，y(n)=x(n)*h(n),即y(n)使x(n)和h(n)的线性卷积，y(n)序列的点数为(M+L-1)X(JΩ)=|X(JΩ)|e^(Jφ(Ω)),其中|X(JΩ)|表示幅频特性，φ(Ω)表示相频特性。f=ifourier(F)表示函数F的逆变换，默认的独立变量是ω，返回的是关于x的函数。一个单位脉冲响应为h(n)的系统对输入序列x(n)的输出为y(n)=x(n)*h(n)设序列x(n),h(n)都是N点序列，其DFT分别为X(k),H(k),Y(k)

当Δf时1/(2T)的倍数时候，信号是正交的，也就是要保证信号的正交性，必须满足相继频率之间的最小频率间隔是1/(2T)。从图中可以看出载波频率4Hz的两边，分布了对应符号0,1,3,2的频率2.5Hz,3.5Hz,4.5Hz,5.5Hz。最佳检测器，又称为平方律检测器，当传输信号波形等概率时，通过计算平方包络(rk)^2= (rkc)^2+ (rks)θm表示第m个信号由于传输延时而产生的相移，n(t)表示加性带通噪声。Δf是相继两个频率之间的频率间隔，Δf=f(m)-f(m-1)2}的最大值的信号。

离散信号x(t),t表示在时间轴上的离散点上的取值，可以表示为nTs，其中Ts表示相邻两个点之间的时间间隔，也就是所谓的抽样周期，n取整数。当序列n=0时，序列h(-m)是序列h(m)的时序取反的结果，时序取反使得h(m)以纵轴为中心翻转180°,相乘后求和的计算称为卷积和。一个离散时间系统，抽象看作一种变换或者一种映射，把输入序列x(n)变换为输出序列y(n),即y(n)=T[x(n)].一般情况下，当x(n)及h(n)的离散序列长度分别为N1、N2时，卷积y(n)的长度则为N1+N2-1。

经典滤波器，假设输入信号x(n)中的有效信号和噪声信号分布在不同频带，当x(n)通过一个线性滤波系统后，可以有效地减少或去除噪声信号成分。经典滤波器主要有低通滤波器(Low Psss Filter,LPF)、高通滤波器(High Pass Filter,HPF)、带通滤波器(Band Pass Filter,BPF)、带阻滤波器(Band Stop Filter,BSF)、全通滤波器(All Pass Filte,APF)。通带定义为|ω|≤ ωp，1-α1≤|H(e^(jω))|≤1。

subplot(m,n,p)命令将当前图形窗口分成m×n个绘图区，即每行n个，共m行，区号按行优先编号，且选定第p个区为当前绘图区。使用图形保持命令hold，在使用plot函数画图时，可以在已存在图形窗口中继续使用plot命令添加新的图形内容。loglog(x,y)——双对数坐标，横坐标和纵坐标都采用x,y的对数。semilogx(x,y)——单对数x坐标，横坐标采用x的对数。semilogy(x,y)——单对数y坐标，横坐标采用y的对数。句柄显示在图形窗口的标题栏中，即图形窗口标题。

square(T,DUTY)——产生指定占空比的方波，DUTY指定信号为正值的区域在一个周期T内所占的比例，取值为0-100，当DUTY为50时，产生方波信号，与square(T)函数相同。sawtooth(T，WIDTH)——对时间变量T产生三角波，WIDTH指定三角波的尺度值，取值为0-1，当WIDTH为0.5时，产生对称的三角波信号，当WIDTH为1时，产生锯齿波信号。freqz有两种用法：freqz(b,a,n,Fs)及[h,f]=freqz(b,a,n,Fs)

函数运算主要是调用Matlab系统内置的运算函数来求取矩阵的行列式(det(A))，求秩(rank(A)),求特征值和特征矢量([V,D]=eig(A)),求Jordan标准型(jordan(A))和矩阵分解。Matlab中的变量命名以字母开头，后接字母、数字或下划线的字符序列，最多63个字符，区分字母的大小写。m的随机矩阵，矩阵中的每一行，每一列元素都服从均值为0，方差为1的高斯分布。m的随机矩阵，矩阵中的每一行，每一列元素都服从[0,1]上的均匀分布。randn(m,n)或randn(m)产生m。