OFDM.rar_QPSK 信道估计_ofdm 插值_ofdm插值_信道估计 qpsk_信道插值

%%%%%%%%%% 加载数据 %%%%%%%%%% load('hw5.mat') load('signal.mat') %%%%%%%%%% 符号映射 %%%%%%%%%% diff = zeros(25600,4); diff(:,1) = abs(signal_origin_vector - (-1-1i)); diff(:,2) = abs(signal_origin_vector - (-1+1i)); diff(:,3) = abs(signal_origin_vector - (1-1i)); diff(:,4) = abs(signal_origin_vector - (1+1i)); min_diff = min(diff')'; for k = 1:4 diff(:,k) = diff(:,k) == min_diff; end flag = 1 * diff(:,1) + 2 * diff(:,2) + 3 * diff(:,3) + 4 * diff(:,4); %%%%%%%%%% 恢复二进制比特流 %%%%%%%%%% bit = zeros(51200,1); %储存恢复的二进制比特流 for k = 1:25600 switch flag(k) case 1 bit(2*k-1) = 0; bit(2*k) = 0; case 2 bit(2*k-1) = 0; bit(2*k) = 1; case 3 bit(2*k-1) = 1; bit(2*k) = 0; case 4 bit(2*k-1) = 1; bit(2*k) = 1; end end %%%%%%%%%% 计算误比特率 %%%%%%%%%% compare = (bit ~= data_std); SBR = sum(compare) / 51200; disp(SBR); %%%%%%%%%% 计算误码率 %%%%%%%%%% compare2 = (compare(1:2:51200) + compare(2:2:51200)) == 2; SER = (sum(compare) - sum(compare2)) / 25600; disp(SER);