注水定理和矩阵分解.docxVIP

随机信号、非周期信号、准周期信号和周期信号，由于其时间是无限的，总为功率信号能量信号只可能是有限时间内存在的确定性信号公式 1AWGN信道离散事件加性高斯白噪声信道： ，信道容量由香农公式得出：接收信噪比定义为x[i]的功率除以n[i]的功率，即用公式（1）可计算N次观测下出信号和噪声的功率，但由于白噪声单边功率谱密度为N0/2，所以噪声功率可用N0B表示。信道容量和注水定理信道容量就是下面这个公式：C=log(1+P/N)其中C是信道容量，P是信号功率，N是噪声功率。 这几个字母，是香农老人家最牛逼的成果，指引着通信技术的发展方向。如果信号的功率为P=P1+P2，可以验证一个公式：C=log(1+(P1+P2)/N)=log(1+P1/N)+log(1+P2/(P1+N))=C1+C2。C1=log(1+P1/N)可理解为信号功率P1，噪声为N时的容量，C2=log(1+P2/(P1+N))可理解为信号功率P2，噪声为(P1+N)时的容量。有人把这个叫superposition coding。也就是说，这两份功率，第一份功率产生了一个容量，同时等效成了对第二份功率的噪声。了解了这个结论，我们来想，如果有N个并行的信道，这N个信道上已经有一些噪声或者信号的功率。如果我再有一小份功率，把它分配到哪个信道上能获得最大的信道容量呢？从上面的分析我们知道，已经存在的功率，无论是噪声还是信号，对后来的信号来说都是噪声，因此，这一小份功率分配到累计功率最小的信道上获得的容量最大。这样分配的结果就是N个信道上的功率是相同的，就是注水定理。注水定理可以看作是信道容量的动态过程。matlab中的矩阵分解?矩阵分解是指根据一定的原理用某种算法将一个矩阵分解成若干个矩阵的乘积。常见的矩阵分解有LU分解（三角分解）、QR分解（正交变换）、Cholesky分解，以及Schur分解、Hessenberg分解、奇异分解等。?(1) LU分解（三角分解）矩阵的LU分解就是将一个矩阵表示为一个交换下三角矩阵和一个上三角矩阵的乘积形式。线性代数中已经证明，只要方阵A是非奇异（即行列式不等于0）的，LU分解总是可以进行的。MATLAB提供的lu函数用于对矩阵进行LU分解，其调用格式为：[L,U]=lu(X)：产生一个上三角阵U和一个变换形式的下三角阵L(行交换)，使之满足X=LU。注意，这里的矩阵X必须是方阵。[L,U,P]=lu(X)：产生一个上三角阵U和一个下三角阵L以及一个置换矩阵P，使之满足PX=LU。当然矩阵X同样必须是方阵。（设P?是一个?m×n?的?(0,1)?矩阵，如?m≤n且?P*P′=E，则称?P为一个?m×n的置换矩阵。）实现LU分解后，线性方程组Ax=b的解x=U\(L\b)或x=U\(L\Pb)，这样可以大大提高运算速度。?例7-2??用LU分解求解例7-1中的线性方程组。命令如下：A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4];b=[13,-9,6,0];[L,U]=lu(A);x=U\(L\b)或采用LU分解的第2种格式，命令如下：[L,U ,P]=lu(A);x=U\(L\P*b)?(2) QR分解（正交变换）对矩阵X进行QR分解，就是把X分解为一个正交矩阵Q和一个上三角矩阵R的乘积形式。QR分解只能对方阵进行。MATLAB的函数qr可用于对矩阵进行QR分解，其调用格式为：[Q,R]=qr(X)：产生一个一个正交矩阵Q和一个上三角矩阵R，使之满足X=QR。[Q,R,E]=qr(X)：产生一个一个正交矩阵Q、一个上三角矩阵R以及一个置换矩阵E，使之满足XE=QR。实现QR分解后，线性方程组Ax=b的解x=R\(Q\b)或x=E(R\(Q\b))。?例7-3??用QR分解求解例7-1中的线性方程组。命令如下：A=[2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4];b=[13,-9,6,0];[Q,R]=qr(A);x=R\(Q\b)或采用QR分解的第2种格式，命令如下：[Q,R,E]=qr(A);x=E*(R\(Q\b))?(3) Cholesky分解如果矩阵X是对称正定的，则Cholesky分解将矩阵X分解成一个下三角矩阵和上三角矩阵的乘积。设上三角矩阵为R，则下三角矩阵为其转置，即X=RR。MATLAB函数chol(X)用于对矩阵X进行Cholesky分解，其调用格式为：R=chol(X)：产生一个上三角阵R，使RR=X。若X为非对称正定，则输出一个出错信息。[R,p]=chol(X)：这个命令格式将不输出出错信息。当X为对称正定的，则p=0，R与上述格式得到的结果相同；否则p为一个正整数。如果X为满秩矩阵，则R为一个阶数为q=p-1的上三角阵，且满足RR=X(1:q,1: