clear?all;
close?all;
IFFT_bin_length=1024;
carrier_count=200;
bits_per_symbol=2;
symbols_per_carrier=50;
SNR=10;
baseband_out_length=carrier_count*symbols_per_carrier*bits_per_symbol;%總比特數
carriers=（1:carrier_count）+（floor（IFFT_bin_length/4）-floor（carrier_count/2））;
conjugate_carriers=IFFT_bin_length-carriers+2;
%發送端
%產生隨即二進制數據：
baseband_out=round（rand（1baseband_out_length））;
convert_matrix=reshape（baseband_outbits_per_symbollength（baseband_out）/bits_per_symbol）;
for?k=1:（length（baseband_out）/bits_per_symbol）
???modulo_baseband（k）=0;
???for?i=1:bits_per_symbol
???????modulo_baseband（k）=modulo_baseband（k）+convert_matrix（ik）*2^（bits_per_symbol-i）;
???end
end
%串并轉換
carrier_matrix=reshape（modulo_basebandcarrier_countsymbols_per_carrier）‘;
%對每個載波的符號進行差分編碼
carrier_matrix=[zeros（1carrier_count）;carrier_matrix];
for?i=2:（symbols_per_carrier+1）
????carrier_matrix（i:）=rem（carrier_matrix（i:）+carrier_matrix（i-1:）2^bits_per_symbol）;
end
%把差分符號代碼轉換成相位
carrier_matrix=carrier_matrix*（（2*pi）/（2^bits_per_symbol））;
%把相位轉換成復數
[XY]=pol2cart（carrier_matrixones（size（carrier_matrix1）size（carrier_matrix2）））;
complex_carrier_matrix=complex（XY）;
%分配載波到指定的IFFT位置
IFFT_modulation=zeros（symbols_per_carrier+1IFFT_bin_length）;
IFFT_modulation（:carriers）=complex_carrier_matrix;
IFFT_modulation（:conjugate_carriers）=conj（complex_carrier_matrix）;
%畫出頻域中的OFDM信號代表
figure（1）
stem（0:IFFT_bin_length-1abs（IFFT_modulation（21:IFFT_bin_length））‘b*-‘）
grid?on
axis（[0?IFFT_bin_length?-0.5?1.5]）
ylabel（‘Maginitude‘）
Xlabel（‘IFFT?Bin‘）
title（‘OFDM?Carrier?Frequenxy?Magnitude‘）

figure（2）
plot（0:IFFT_bin_length-1（180/pi）*angle（IFFT_modulation（21:IFFT_bin_length））‘go‘）
hold?on
stem（carriers-1（180/pi）*angle（IFFT_modulation（2carriers））‘b*-‘）
stem（conjugate_carriers-1（180/pi）*angle（IFFT_modulation（2conjugate_carriers））‘b*-‘）
axis（[0?IFFT_bin_length?-200?+200]）
grid?on
ylabel（‘Phase（degrees）‘）
Xlabel（‘IFFT?Bin‘）
title（‘OFDM?Carrier?Phase‘）
%通過IFFT將頻域轉化為時域信號，得到時域信號
time_wave_matrix=ifft（IFFT_modulation‘）;
time_wave_matrix=time_wave_matrix‘;

figure（3）
plot（0:IFFT_bin_length-1time_wave_matrix（2:））
grid?on
ylabel（‘Amplitude‘）
xlabel（‘Time‘）
title（‘OFDM?Time?SignalOne?Symbol?Period‘）

%畫出每個載波對應的時域信號（分離的OFDM信號）3131313131313131
for?f=1:carrier_count
????temp_bins（1:IFFT_bin_length）=0+0j;
????temp_bins（carriers（f））=IFFT_modulation（2carriers（f））;
????temp_bins（conjugate_carriers（f））=IFFT_modulation（2conjugate_carriers（f））;
????temp_time=ifft（temp_bins‘）;

end
?figure（4）
????plot（0:IFFT_bin_length-1temp_time）
????hold?on

grid?on
ylabel（‘Amplitude‘）
xlabel（‘Time‘）
title（‘Separated?Time?Waveforms?Carriers‘）
for?i=1:symbols_per_carrier+1
?????windowed_time_wave_matrix（i:）=real（time_wave_matrix（i:））.*hamming（IFFT_bin_length）‘;
????windowed_time_wave_matrix（i:）=real（time_wave_matrix（i:））;
end
%串并轉換
ofdm_modulation=reshape（windowed_time_wave_matrix‘1IFFT_bin_length*（symbols_per_carrier