%function?a1

%author：liuqingwei?May?12th?2006?9:40?PM.
%OFDM?Channel?Estimation?based?on?Comb?Pilot
%IFFT_bin_length:?IFFT和FFT的點數
%carrier_count:?子載波個數
%bits_per_symbol:?每符號上的比特數
%symbols_per_carrier:?每楨的OFDM符號數
%Ｘ：欲發送的二進制比特流

clear?all;
clc;
IFFT_bin_length=512;
carrier_count=270;
bits_per_symbol=2;
symbols_per_carrier=12;
LI=7?;?%導頻之間的間隔
Np=ceil（carrier_count/LI）+1;%導頻數?%加1的原因：使最后一列也是導頻

N_number=carrier_count*symbols_per_carrier*bits_per_symbol;
carriers=1:carrier_count+Np;

GI=8;?????????????%?guard?interval?length
N_snr=20;??????????%?每比特信噪比
snr=8;?????????????%信噪比間隔
%------------------------------------------------------------
%?vector?initialization


X=zeros（1N_number）;
X1=[];
X2=[];
X3=[];
X4=[];
X5=[];
X6=[];
X7=[];
Y1=[];
Y2=[];
Y3=[];
Y4=[];
Y5=[];
Y6=[];
Y7=[];
XX=zeros（1N_number）;
dif_bit=zeros（1N_number）;
dif_bit1=zeros（1N_number）;
dif_bit2=zeros（1N_number）;
dif_bit3=zeros（1N_number）;
X=randint（1N_number）;%產生二進制隨即序列（非0即1）
%--------------------------------------------------------
%QPSK調制：（1?1）->pi/4;（0?1）->3*pi/4;（0?0）->-3*pi/4;（10）->-pi/4;

s=（X.*2-1）/sqrt（2）;
sreal=s（1:2:N_number）;
simage=s（2:2:N_number）;
X1=sreal+j.*simage;

%---------------------------------------------------------
%產生隨機導頻信號
%--------------------------------------------------------
train_sym=randint（12*symbols_per_carrier）;
t=（train_sym.*2-1）/sqrt（2）;
treal=t（1:2:2*symbols_per_carrier）;
timage=t（2:2:2*symbols_per_carrier）;
training_symbols1=treal+j.*timage;
training_symbols2=training_symbols1.‘;
training_symbols=repmat（training_symbols21Np）;

%disp（training_symbols）
pilot=1:LI+1:carrier_count+Np;
if?length（pilot）~=Np
????pilot=[pilotcarrier_count+Np];
end
%--------------------------------------------------------
%串并轉換
X2=reshape（X1carrier_countsymbols_per_carrier）.‘;
%---------------------------------------------------------
%插入導頻
signal=1:carrier_count+Np;
signal（pilot）=[];
X3（:pilot）=training_symbols;
X3（:signal）=X2;
%X3=cat（1training_symbolsX2）;
IFFT_modulation=zeros（symbols_per_carrierIFFT_bin_length）;
IFFT_modulation（:carriers）=X3;
%IFFT_modulation（:conjugate_carriers）=conj（X3）;
X4=ifft（IFFT_modulationIFFT_bin_length2）;
%X5=X4.‘;
%加保護間隔（循環前綴）
for?k=1:symbols_per_carrier;
???for?i=1:IFFT_bin_length;
??????X6（ki+GI）=X4（ki）;
???end
???for?i=1:GI;
??????X6（ki）=X4（ki+IFFT_bin_length-GI）;????
???end
end
%---------------------------------------------------------
%并串轉換
X7=reshape（X6.‘1symbols_per_carrier*（IFFT_bin_length+GI））;

%---------------------------------------------------------
%信道模型:帶多普勒頻移的瑞利衰落信道

fd=100;?%多普勒頻移
r=6;???%多徑數
a=[0.123?0.3?0.4?0.5?0.7?0.8];?%多徑的幅度
d=[2?3?4?5?9?13];?%各徑的延遲
T=1;??%系統采樣周期
th=[90?0?72?144?216?288]*pi./180;%相移
h=zeros（1carrier_count）;
hh=[];
????for?k=1:r
????????%deta=[zeros（1d（k）-1）1zeros（1carrier_count-d（k））];
????????h1=a（k）*exp（j*（（2*pi*T*fd*d（k）/carrier_count）））;
????????%h1=a（k）*exp（j*（（2*pi*T*fd*d（k）/carrier_count）））;
????????hh=[hhh1];
????end
????h（d+1）=hh;
%noise=randn（1length（X7））+j

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????7897??2007-04-29?09:57??OFDM\a1.m

?????文件???????1668??2006-05-15?11:10??OFDM\OFDM_BPSK_BER\BPSKmodanddemod.m

?????文件???????1833??2006-05-15?11:10??OFDM\OFDM_BPSK_BER\increaseber.m

?????文件???????4707??2006-05-15?11:10??OFDM\OFDM_BPSK_BER\ofdm_ber.m

????..AD...?????????0??2007-04-26?14:52??OFDM\OFDM_BPSK_BER

?????目錄??????????0??2007-04-28?11:27??OFDM

-----------?---------??----------?-----??----

????????????????16323????????????????????7