Matla峰均功率比设计.doc

设计题目 峰均功率比设计。 实验内容及要求 生成四个长度为64的序列,分别为s1,s2,s3,s4,对每个序列乘上[1,-1,j,-j]中的任意一个因子后求和.求所有可能的和序列中峰均功率比最低的序列及相应的因子。 三、实验过程(详细设计) 2.相关知识： 峰均比，或称峰值因数(crest factor)，简称PAR（peak-to-average ratio），或叫峰均功率比(简称PARR,peak-to-average power ratio)。 先说定义：峰均比是一种对波形的测量参数,等于波形的振幅除以有效值(RMS)所得到的一个比值峰值功率既是只以某种概率出现的肩峰的瞬时功率。通常概率取为0.01%。 平均功率是系统输出的实际功率。 在某个概率下峰值功率跟平均功率的比就称为某个概率下的峰均比，比如PAR=9.1@0.1%,各种概率的峰均比就形成了CCDF曲线（互补累积分布函数）。 在概率为0.01%处的PAR,一般称为CREST因子。 3.峰均功率比详细设计： 1）利用相关函数生成s1，s2，s3.s4序列； 2）利用四重循环，再加判断语句完成因子的24种不同的组合； 3）求每一种因子的功率最大值，功率平均值，以及每一种因子对应的峰均功率比； 4）求所有可能中的最小峰均功率比，以及对应的因子和序列。 4.MATLAB源代码： clc;clear; %clc，清屏；clear，清楚变量，以免在程序运行中出现错误。 s=rand(4,64); %生成s1,s2,s3,s4。 a=[1,-1,j,-j]; %原始因子组合。 b=zeros(1,4); %生成一个1*4的零矩阵，用于存放一种因子。 d=zeros(24,4); %生成一个24*4的零矩阵，用于存放24种不同的排列的因子。 x=1; %标记。 c=zeros(24,64); %用于存放S与因子相乘之后的和。 p=zeros(24,64); %用于存放每一个元素的功率。 for m1=1:4 %用四重循环，再加判断语句完成因子的24种不同的组合。 b(1)=a(m1); for m2=1:4 if m2~=m1 b(2)=a(m2); else continue; end for m3=1:4 if m3~=m2 m3~=m1 b(3)=a(m3); else continue; end for m4=1:4 if m4~=m3 m4~=m2 m4~=m1 b(4)=a(m4); else continue; end end d(x,:)=b;x=x+1; %将每一种因子的组合放在d数组中，x加1。 end end end for m6=1:24 c(m6,:)=d(m6,:)*s; p(m6,:)=abs(c(m6,:)).^2; %每一行求模,再求功率。 pmax(m6)=max(p(m6,:)); %求每一行功率中的最大值，即每一种因子对应的功率的最大值。 avg(m6)=mean(p(m6,:)); %求每一种因子对应的功率平均值。 pm(m6)=pmax(m6)/avg(m6); %求每一种因子对应的峰均功率比。 end pmin=min(pm); %求所有可能中的最小峰均功率比。 y=1; %标记y。 for n=1:24 if pm(n)==pmin yinzi(y,:)=d(n,:); %求最小峰均功率比对应的因子。 xuli