【正文】 　Rx_decoded_phase(phase_negative) = rem(Rx_decoded_phase(phase_negative)+360,360)。　　%　% Extract phase differences (from the differential encoding)　　% the matlab diff( ) function is perfect for this operation　% again, normalize the result to be between 0 and 359 degrees　　%　Rx_decoded_phase = diff(Rx_phase)。　　phase_negative = find(Rx_phase 0)。)。　polar(Rx_phase_P, Rx_mag_P,39。%% PLOT EACH RECEIVED SYMBOL　　%figure (9)　　Rx_phase_P = angle(Rx_carriers)。OFDM Receive Spectrum, Phase39。FFT Bin39。Phase (degrees)39。b*39。b*39。go39。OFDM Receive Spectrum, Magnitude39。FFT Bin39。Magnitude39。b*39。%% Transform each symbol from time to frequency domain% take the fft of each column%　　Rx_spectrum = fft(Rx_Data_matrix)。　　Rx_Data = Tx_data + noise?！　oise_scale_factor = sqrt(noise_sigma)?！　?linear_SNR = 10^(SNR/10)。)　%　% ENDPLOT　%% Upconversion to RF 　　%　　% For this model, the baseband will be inserted directly into the channel　　% without conversion to RF frequencies.　%Tx_data = ofdm_modulation。)　　%title(39。)　　%xlabel(39。)%grid on%axis([0 40 max(average_fft_log)])　%ylabel(39?！?figure (6)　%plot((0:(avg_temp_time1))/avg_temp_time, average_fft_log)　%hold on　%plot(0:1/IFFT_bin_length:1, 35, 39?！? average_fft = average_fft + (subset_ofdm_f/averages)?！　?for a = 0:(averages1)% subset_ofdm = ofdm_modulation(((a*avg_temp_time)+1):((a+1)*avg_temp_time))。　　%averages = floor(temp_time/avg_temp_time)。)%　　% PLOT OFDM SIGNAL (spectrum)%symbols_per_average = ceil(symbols_per_carrier/5)。)　%title(39。)%xlabel(39。%figure (5)　%plot(0:temp_time1,ofdm_modulation)%grid on　%ylabel(39。, 1, IFFT_bin_length*(symbols_per_carrier+1))。reshape39?！? windowed_time_wave_matrix(i,:) = real(time_wave_matrix(i,:))。Separated Time Waveforms Carriers39。Time39。Amplitude39。)?！　? temp_bins(conjugate_carriers(f))=IFFT_modulation(2,conjugate_carriers(f))。　%for f = 1:carrier_count　% temp_bins(1:IFFT_bin_length)=0+0j。k39。b39。g39。r39。k39。b39。g39。r39。k39。b39。g39。r39。k39。b39。g39。r39。OFDM Time Signal, One Symbol Period39。Time39。Amplitude39。start39?！ime_wave_matrix = time_wave_matrix39。s spectrum (represented by a row of carriers) to the % time domain via IFFT%　　time_wave_matrix = ifft(IFFT_modulation39。OFDM Carrier Phase39。IFFT Bin39。Phase (degrees)39。b*39。b*39。go39。OFDM Carrier Frequency Magnitude39。IFFT Bin39。Magnitude39。b*39?！FFT_modulation(:,conjugate_carriers) = conj(plex_carrier_matrix)。%% Assign each carrier to its IFFT bin　　% each row of plex_carrier_matrix represents one symbol period, with% a symbol for each carrier% a matrix is generated to represent all IFFT bins (columns) and all % symbols (rows)　% the phase modulation for each carrier is then assigned to the % appropriate bin　　% the conjugate of the phase modulation is then assigned to the % appropriate bin% the phase modulation bins and their conjugates are symmetric about 　　% the Nyquist frequency in the IFFT bins% since the first bin is DC, the Nyquist Frequency is located 　% at (number of bins/2) + 1　% symmetric conjugates are generated so that when the signal is 　% transformed to the time domain, the time signal will be realvalued　　% example% 1024 IFFT bins　　% bin 513 is the center (symmetry point)　% bin 1 is DC% bin 514 is the plex conjugate of bin 512　% bin 515 is the plex conjugate of bin 511　　% ....　% bin 1024 is the plex conjugate of bin 2 (if all bins 　　% were used as carriers)　% So, bins 2512 map to bins 1024514　%　　IFFT_modulation = zeros(symbols_per_carrier + 1, IFFT_bin_length)。　　%% Convert the phase to a plex number% each symbol is given a magnitude of 1 to go along with its phase % (via the ones(r,c) function)　% it is then converted from polar to cartesian (plex) form　% the result is 2 matrices, X with the real values and Y with the imaginary% each X column has all the real values for a carrier, and each Y column 　% has the imaginary values for a carrier% a single plex matrix is then generated taking X for real and 　　% Y for imaginary　%[X,Y] = pol2cart(carrier_matrix, ones(size(carrier_matrix,1),size(carrier_matrix,2)))。　　for i = 2:(symbols_per_carrier + 1)　　 carrier_matrix(i,:) = rem(carrier_matrix(i,:)+carrier_matrix(i1,:),2^bits_per_symbol)。y] of a row of zeros with the carrier matrix, sweet!)　　% perform modulo N addition between symbol(n) and symbol(n1) to get the % coded value of symbol(n)　% for example:　　% bits_per_symbol = 2 (modulo 4)　　% symbol stream = 3 2 1 0 2 3　% start symbol = 0　　%% coded symbols = 0 + 3 = 3% 3 + 2 = 11 = 1　　% 1 + 1 = 2% 2 + 0 = 2　% 2 + 2 = 10 = 0　　% 0 + 3 = 3%　% coded stream = 0 3 1 2 2 0 3　%%carrier_matrix = [zeros(1,carrier_count)?！? end　　end　　%% Serial to Parallel Conversion　　% convert the serial modulo N stream into a matrix where each column 　% represents a carrier and each row represents a symbol　% for example:　　%　% serial input stream = a b c d e f g h i j k l m n o p　%% parallel carrier distribution =% C1/s1=a C2/s1=b C3/s1=c C4/s1=d　　% C1/s2=e C2/s2=f C3/s2=g C4/s2=h　% C1/s3=i C2/s3=j C3/s3=k C4/s3=l% . . . .　% . . . .%　carrier_matrix = reshape(modulo_baseband, carrier_count, symbols_per_carri