24小时热门版块排行榜

返回列表

末页

pjstang

木虫 (著名写手)

PC-EPI: 1
应助: 200 (高中生)
金币: 2177.2
帖子: 1555
在线: 292.1小时
虫号: 826940

[交流] 关键论文接受，散金币300

关键论文接受，散金币300。祝福各位虫子多发牛文。。。。。。。

回复此楼

» 本帖已获得的红花（最新10朵）

lichuang5577

» 猜你喜欢

职称评审没过，求安慰已经有21人回复
垃圾破二本职称评审标准已经有15人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有20人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
聘U V热熔胶研究人员已经有10人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复
三无产品还有机会吗已经有6人回复

» 抢金币啦！回帖就可以得到:
查看全部散金贴

1楼2014-05-22 23:12:46

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

回帖支持 ( 显示支持度最高的前 50 名 )

laolua

铁杆木虫 (小有名气)

应助: 0 (幼儿园)
金币: 7300.5
帖子: 193
在线: 64.5小时
虫号: 2639632

★
pjstang(金币+1): 谢谢参与

祝福楼主，也预祝自己的论文接收，

[ 发自小木虫客户端 ]

赞一下(1人)

回复此楼

7楼2014-05-22 23:16:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

mvpyqz

至尊木虫 (知名作家)

应助: 236 (大学生)
金币: 11502
帖子: 6086
在线: 394.1小时
虫号: 1105587

★
pjstang(金币+1): 谢谢参与

% CRB_final

%CRB CALCULATION

clear all
nFFT       = 64; % fft size
MSym       = 64; % lenght of each OFDM block
nb    = 150; % number of symbols
montcarlo=10;
EsN0dB    = [-10:5:25]; % bit to noise ratio

% Discrete Fourier Transform Matrix F
for iii=1:nFFT;
for jjj=1:nFFT;
      F(iii,jjj)=(1/sqrt(nFFT))*exp(j*(2*pi*(iii-1)*(jjj-1))/nFFT);
end;
end;

% multipath channel
nTap = 3;
k=0:2;
multipath=exp(-(k)/10);
multipath=multipath/sum(multipath);

for k=1:3
ht(:,k)=(randn(1,1)+ sqrt(-1)*randn(1,1))*sqrt(0.5)*sqrt(multipath(k)); %% rayleigh fading
end

% computing and storing the frequency response of the channel, for use at recevier
hF =sqrt(nFFT)*F(:,1:nTap)*(ht.');
hF_diag=diag(hF);

% producing matrix Hh , circulant channel matrix
htt=[ht(1,

zeros(1,nFFT-nTap)];
for i=1:nFFT
for j=1:nFFT
      H_cir(i,j)=htt(mod((i-j),nFFT)+1);
end
end

nt = (1/sqrt(2))*[randn(MSym,nb)+ sqrt(-1)*randn(MSym,nb)];


%Precoding

t=nFFT/4;
for m=1:nFFT
for q=1:nFFT
      if (m==q)&(m==t)
         W(m,q)=2;
      elseif m==q
            W(m,q)=1;
      elseif q==t
            W(m,q)=1;
      else
            W(m,q)=0;
      end
   end
end

P=W*(W');
% noise production

  b =(1/sqrt(2))*[randn(1,MSym)+ sqrt(-1)*randn(1,MSym)];
      moy=mean(b);
      bc = b - moy;
      varb = (1/(MSym))*trace((diag(bc)*(diag(bc))'));
%%%%%%%%%%%%%%%%%%%%%%

Rs=1;
for ii = 1:length(EsN0dB)

      noise = sqrt(((0.5*10^(-EsN0dB(ii)/10))/varb))* conj(bc');
      Var = std(noise)^2;

Rx=Rs* H_cir*(F')*P*F*( H_cir')+Var*eye(nFFT);
for l=1:nTap
      for i=1:nFFT
         for j=1:nFFT
            if (mod((i-j),nFFT))==(l-1)
                  Hl(i,j)=1;
            else
                  Hl(i,j)=0;
            end
         end
      end
      Al=Rs*(Rx^(-.5))*Hl*(F')*P*F*(hF_diag')*(Rx^(-.5));
      G1(:,l)=reshape((Al+(Al')),[],1);
      G2(:,l)=reshape((sqrt(-1)*Al-sqrt(-1)*(Al')),[],1);

end

v=reshape(((Rx^(-.5))*H_cir*(F')*P*F*(H_cir')*(Rx^(-.5))),[],1);
u=reshape((Rx^(-1)),[],1);
G=[G1 G2];
delta=[u v];
P_delta=eye(nFFT*nFFT)-delta*(((delta')*delta)^(-1))*(delta');
CRB=(1/150)*pinv((G')*P_delta*G);
CRB_plot(ii)=(trace(CRB));

end

plot(EsN0dB,CRB_plot)
axis([-10 25 10^(-5) 10^(-1)]);

[ 发自手机版 http://muchong.com/3g ]

赞一下(1人)

回复此楼

27楼2014-05-22 23:23:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

普通回帖