| 查看: 11316 | 回復(fù): 389 | |||
[交流]
關(guān)鍵論文接受,散金幣300
|
|||
| 關(guān)鍵論文接受,散金幣300。祝福各位蟲子多發(fā)牛文。。。。。。。 |
» 搶金幣啦!回帖就可以得到:
+3/164
+1/88
+1/88
+1/82
+3/71
+2/58
+1/49
+1/40
+3/26
+1/16
+1/9
+1/6
+1/6
+1/5
+2/4
+1/4
+1/3
+1/3
+1/2
+1/2
|
% 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)]); [ 發(fā)自手機版 http://www.gaoyang168.com/3g ] |







| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 310分求調(diào)劑 +4 | 成功上岸wang 2026-04-01 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 318求調(diào)劑,計算材料方向 +5 | 吸喵有害笙命 2026-04-01 | 5/250 |
|
|
[考研] 332求調(diào)劑 +8 | Lyy930824@ 2026-03-29 | 8/400 |
|
|
[考研] 349求調(diào)劑 +6 | 吃的不少 2026-04-01 | 6/300 |
|
|
[考研] 總分343,求生物學(xué)調(diào)劑 +6 | 深刻的凱撒 2026-03-26 | 6/300 |
|
|
[考研] 調(diào)劑 +5 | 好好讀書。 2026-03-28 | 7/350 |
|
|
[考研] 279求調(diào)劑 +6 | 莫xiao 2026-04-01 | 6/300 |
|
|
[考研] 一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +10 | @taotao 2026-03-31 | 11/550 |
|
|
[考研] 333求調(diào)劑 +4 | 阿科逸 2026-03-31 | 4/200 |
|
|
[考研] 初試301,代碼085701環(huán)境工程,本碩一致,四六級已過,有二區(qū)一作,共發(fā)表5篇論文 +3 | axibli 2026-04-01 | 3/150 |
|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +12 | Koxui 2026-03-28 | 12/600 |
|
|
[考研] 352分-085602-一志愿985 +6 | 海納百川Ly 2026-03-29 | 6/300 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-29 | 10/500 |
|
|
[考研] 262求調(diào)劑 +7 | ZZ..000 2026-03-30 | 8/400 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +5 | 小張zxy張 2026-03-26 | 10/500 |
|
|
[考研] 298求調(diào)劑 +3 | 什么是胖頭魚 2026-03-30 | 5/250 |
|
|
[考研] 283求調(diào)劑(080500) +14 | A child 2026-03-27 | 14/700 |
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 中國科學(xué)院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|