%DVB-T?2K?Transmission
%The?available?bandwidth?is?8?MHz
%2K?is?intended?for?mobile?services
clear?all;
close?all;
%DVB-T?Parameters
Tu=224e-6;?%useful?OFDM?symbol?period
T=Tu/2048;?%baseband?elementary?period
G=1/4;?%choice?of?1/4?1/8?1/16?and?1/32
delta=G*Tu;?%guard?band?duration
Ts=delta+Tu;?%total?OFDM?symbol?period
Kmax=1705;?%number?of?subcarriers
Kmin=0;
FS=4096;?%IFFT/FFT?length
q=10;?%carrier?period?to?elementary?period?ratio
fc=q*1/T;?%carrier?frequency
Rs=4*fc;?%simulation?period
t=0:1/Rs:Tu;
%Data?generator?（A）
M=Kmax+1;
rand（‘state‘0）;
a=-1+2*round（rand（M1））.‘+i*（-1+2*round（rand（M1）））.‘;
A=length（a）;
info=zeros（FS1）;
info（1:（A/2））?=?[?a（1:（A/2））.‘];?%Zero?padding
info（（FS-（（A/2）-1））:FS）?=?[?a（（（A/2）+1）:A）.‘];
%Subcarriers?generation?（B）
carriers=FS.*ifft（infoFS）;
tt=0:T/2:Tu;
figure（1）;
subplot（211）;
stem（tt（1:20）real（carriers（1:20）））;
subplot（212）;
stem（tt（1:20）imag（carriers（1:20）））;
figure（2）;
f=（2/T）*（1:（FS））/（FS）;
subplot（211）;
plot（fabs（fft（carriersFS））/FS）;
[Pxxf]=pwelch（carriers[][][]2/T）;
plot_Pxx1=10*log10（Pxx）;
subplot（212）;
plot（fplot_Pxx1）;
%?D/A?simulation
L?=?length（carriers）;
chips?=?[?carriers.‘;zeros（（2*q）-1L）];
p=1/Rs:1/Rs:T/2;
g=ones（length（p）1）;?%pulse?shape
figure（3）;
stem（pg）;
dummy=conv（gchips（:））;
u=[dummy（1:length（t））];?%?（C）
figure（4）;
subplot（211）;
plot（t（1:400）real（u（1:400）））;
subplot（212）;
plot（t（1:400）imag（u（1:400）））;
figure（5）;
ff=（Rs）*（1:（q*FS））/（q*FS）;
subplot（211）;
plot（ffabs（fft（uq*FS））/FS）;
[Pxxf]=pwelch（u[][][]Rs）;
plot_Pxx2=10*log10（Pxx）;
subplot（212）;
plot（fplot_Pxx2）;
[ba]?=?butter（131/20）;?%reconstruction?filter
[HF]?=?FREQZ（baFSRs）;
figure（6）;
plot（F20*log10（abs（H）））;
uoft?=?filter（bau）;?%baseband?signal?（D）
figure（7）;
subplot（211）;
plot（t（80:480）real（uoft（80:480）））;
subplot（212）;
plot（t（80:480）imag（uoft（80:480）））;
figure（8）;
subplot（211）;
plot（ffabs（fft（uoftq*FS））/FS）;
[Pxxf]=pwelch（uoft[][][]Rs）;
plot_Pxx3=10*log10（Pxx）;
subplot（212）;
plot（fplot_Pxx3）;
%Upconverter
s_tilde=（uoft.‘）.*exp（1i*2*pi*fc*t）;
s=real（s_tilde）;?%passband?signal?（E）
figure（9）;
plot（t（80:480）s（80:480））;
figure（10）;
subplot（211）;
plot（ffabs（fft（sq*FS））/FS）;
[Pxxf]=pwelch（s[][][]Rs）;
plot_Pxx4=10*log10（Pxx）;
subplot（212）;
plot（fplot_Pxx4）;
figure（11）;
subplot（211）;
plot（ffabs（fft（（（real（uoft）.‘）.*cos（2*pi*fc*t））q*FS））/FS）;
[Pxxf]=pwelch（（（real（uoft）.‘）.*cos（2*pi*fc*t））[][][]Rs）;
plot_Pxx5=10*log10（Pxx）;
subplot（212）;
plot（fplot_Pxx5）;
figure（12）;
subplot（211）;
plot（ffabs（fft（（（imag（uoft）.‘）.*sin（2*pi*fc*t））q*FS））/FS）;
[Pxxf]=pwelch（（（imag（uoft）.‘）.*sin（2*pi*fc*t））[][][]Rs）;
plot_Pxx6=10*log10（Pxx）;
subplot（212）;
plot（fplot_Pxx6）;

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

?????文件???????2783??2007-06-01?13:04??OFDM?Simulation?Using?Matlab\ofdm_4.m

?????文件?????246272??2007-06-01?13:28??OFDM?Simulation?Using?Matlab\OFDM仿真利用MATLAB.doc

?????目錄??????????0??2008-07-05?17:02??OFDM?Simulation?Using?Matlab

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

???????????????249055????????????????????3