BYR Achieve · 镜像论坛

[求助]关于ofdm符号定时的Schmidl算法

2008/1/15镜像同步0 回复

下边是我自己写的一段关于符号定时的代码，总觉得生成的M(d)的图形怪怪的，请达人们给指点一下，不胜感激 %子载波个数为256 N=256; %已知一个长度为256的伪随机序列 txtmp=[1 1 1 1 1 1 1 1 1 0 0 0 0 1 1 1 1 0 1 1 1 0 0 0 0 1 0 1 1 0 0 1 1 0 1 1 0 1 1 1 1 0 1 0 0 0 0 1 1 1 0 0 1 1 0 0 0 0 1 0 0 1 0 0 0 1 0 1 0 1 1 1 0 1 0 1 1 1 1 0 0 1 0 0 1 0 1 1 1 0 0 1 1 1 0 0 0 0 0 0 1 1 1 0 1 1 1 0 1 0 0 1 1 1 1 0 1 0 1 0 0 1 0 1 0 0 0 0 0 0 1 0 1 0 1 0 1 0 1 1 1 1 1 0 1 0 1 1 0 1 0 0 0 0 0 1 1 0 1 1 1 0 1 1 0 1 1 0 1 0 1 1 0 0 0 0 0 1 0 1 1 1 0 1 1 1 1 1 0 0 0 1 1 1 1 0 0 1 1 0 1 0 0 1 1 0 1 0 1 1 1 0 0 0 1 1 0 1 0 0 0 1 0 1 1 1 1 1 1 1 0 1 0 0 1 0 1 1 0 0 0 1 0 1 0 0 1 1 0 0 0 1 1 0 0 0 ]; %将二进制转化为四进制0、1、2、3，这一步不能保证生成的四进制也为伪随机序列 txtmptmp=reshape(txtmp,N/2,2); txtmptmptmp=bi2de(txtmptmp,'left-msb'); %************利用查表法生成复随机序列********************** QAMTable=[7+7i,-7+7i,-7-7i,7-7i]; tx=QAMTable(txtmptmptmp+1); %*************在奇数子载波的位置插入零********************* x=zeros(1,N); index = 1; for n=1:2:N x(n)=tx(index); index=index+1; end; %**************利用IFFT变换生成Schmidl训练符号*************** y=ifft(x,N); %y的长度为256 %*****************添加循环前缀****************************** prefix=y(1,N-N/8+1:N); ycp=zeros(1,N+N/8); ycp=[prefix';y']'; %*****************计算符号定时***************************** %*****序列P、R、M的长度均为2L即N**************************** L=N/2; P=zeros(1,N); for d=1:1:N for m=1:1:L if (d+m < N+N/8+1) index1 =d+m; else index1 =mod(d+m,N+N/8); end; if (d+m+L < N+N/8+1) index2 =d+m+L; else index2 =mod(d+m+L,N+N/8); end; ptemp = conj(ycp(index1))*ycp(index2); P(d)=P(d)+ptemp; end; end; R=zeros(1,N); for d=1:1:N for m=1:1:L if (d+m+L < N+N/8+1) index =d+m+L; else index =mod(d+m+L,N+N/8); end; rtemp = conj(ycp(index))*ycp(index); R(d)=R(d)+rtemp; end; end; %******************计算阈值M(d)***************************** M=zeros(1,N); for d=1:1:N M(d)=conj(P(d))*P(d)/abs(R(d))/abs(R(d)); end; %M=power(abs(P),2)./power(abs(R),2); plot(M); %问题就在这里了，M(d)的图形怎么理解？好像跟理论上的不一样。。。。

订阅后，新回复会通过你的通知中心匿名送达。