equation.asv

就是自适应的一写程序

%equation_dfe-rls
close all
clear all
clc
signal_length=100;                                                 
%信息序列的长度
training_sequence_length=25;
%训练序列的长度
time=1;                                                    
%蒙特卡洛仿真的次数
K1=5;
%前馈滤波器的阶数为K1+1
K2=5;
%后馈滤波器的阶数为K2
N=K1+K2+1;
%均衡器的阶数
channel=[0.26,
         0.93,
         0.26]; 
%针对p496页的信道情况
path_number=length(channel);
%信道多径的径数
snr=20;
%信噪比
err_result=zeros(training_sequence_length,1); 
%求均方误差
for k=1:time
    signal_in=randint(signal_length,2);  
    %产生信息序列
    signal_out_all=qpsk_mapping(signal_in);
    %进行QPSK映射
    signal_receive_all=channelpass(signal_out_all,channel,snr);
    %映射后的信号通过信道
    signal_out_training=[signal_out_all(1:1:training_sequence_length)];
    [signal_equation,signal_decision,filter_coefficient]=dferls_equation(K1,K2,signal_receive_all,signal_out_training);
    signal_demapping=demapping(signal_decision);
    bit_error=xor(signal_demapping,signal_in([(training_sequence_length+1):(length(signal_in)+path_number-K1-2)],:));
    bit_error_rate=sum(sum(bit_error))/(length(bit_error)*2);
end
%求均方误差
scatterplot(signal_out_all),title('调制之后的星座图')
%调制之后的星座图
scatterplot(signal_receive_all) ,title('通过信道之后的星座图')
%通过信道之后的星座图
scatterplot(signal_equation) ,title('均衡之后的星座图')
%均衡之后的星座图