迭代注水算法

ofdm中的迭代注水，从而实现功率的最小化或速率的最大化。

%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