基于BEM源码的规则

通信系统时变信道建模，各类BEM模型建模及仿真 对BEM模型的CE-BEM、PBEM、DPS-BEM、DKL-BEM进行仿真，用到LS估计技术

clear all
clc

load（‘channel_data10.mat‘）;
M=9;
N=30*（M+1）;
N_end=20000;
number=1;
delta=20;

;%% 构建基函数

for m=1:M
    for n=1:delta:N
        B（mn）=exp（（2*pi*j*（m-0.5*M）*n）/（N_end））;
    end
    B1=B（:1:delta:N）;
   
    for n=1:N
        B2（mn）=exp（（2*pi*j*（m-0.5*M）*n）/（N_end））;
    end
end

B1=B1‘;
[RC]=size（B1）;
% H1=（h（1:））‘;

for t=1:N_end
    if（rem（t-1N）==0）
        H1=（h（numbert:delta:t+N-1））‘;
        %%
%         for m=1:M
%             for n=1:N
%                 B（mn）=exp（（2*pi*j*（m-0.5*M）*（t+n））/（N_end））;
%             end
%         end
%         B=B‘;
%         [RC]=size（B）;
        
        
        %% 基系数的计算
        b1=B1\H1;
        if（rank（H1‘）==R）
            b1=inv（B1）*H1
        end
        
        %% NMSE的计算
        
        H_estimate（1t:t+N-1）=B2‘*b1;%为什么波形相似了，但是幅度不相等
    end
end
H_e=abs（h（number:））-abs（H_estimate）;
% NMSE= sum（abs（（h（number:）-H_estimate））.^2）/sum（abs（h（number:））.^2）;
NMSE= sum（（H_e）.^2）/sum（abs（h（number:））.^2）

figure;
plot（abs（H_estimate））
figure;
% plot（abs（H1））
plot（abs（h（number:）））;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1151  2014-09-20 21:10  BEM源码\111\BEM_CE_1.m

     文件       1310  2014-09-21 14:47  BEM源码\111\BEM_CE_d.m

     文件       1291  2014-09-21 10:18  BEM源码\111\BEM_CE_M.m

     文件       1344  2014-09-21 18:36  BEM源码\111\BEM_CE_N.m

     文件       1403  2014-09-21 15:22  BEM源码\111\BEM_d.mat

     文件       1092  2014-09-21 11:15  BEM源码\111\BEM_DKL.m

     文件       1900  2014-09-21 15:20  BEM源码\111\BEM_DKL_d.m

     文件       1790  2014-09-21 11:29  BEM源码\111\BEM_DKL_M.m

     文件       1422  2014-09-21 18:36  BEM源码\111\BEM_DKL_N.m

     文件       1082  2014-09-21 10:24  BEM源码\111\BEM_DPS.m

     文件        412  2014-09-21 12:13  BEM源码\111\BEM_M.mat

     文件       1451  2014-09-21 19:55  BEM源码\111\BEM_N.mat

     文件       1201  2014-09-21 15:02  BEM源码\111\BEM_P_d.m

     文件       1193  2014-09-21 11:02  BEM源码\111\BEM_P_M.m

     文件       1387  2014-09-21 18:36  BEM源码\111\BEM_P_N.m

     文件    1102344  2014-09-16 00:43  BEM源码\111\channel_data_90_10.mat

     文件    1315350  2014-09-16 00:43  BEM源码\111\channel_data_90_10.txt

     文件      17215  2014-09-22 18:05  BEM源码\111\D.fig

     文件        133  2014-09-20 11:12  BEM源码\111\data_10.m

     文件        283  2014-09-20 10:56  BEM源码\111\data_90_10.m

     文件        150  2014-09-21 15:22  BEM源码\111\d_figure.m

     文件      16945  2014-09-22 18:06  BEM源码\111\M.fig

     文件        157  2014-09-21 17:26  BEM源码\111\M_figure.m

     文件      16598  2014-09-23 08:16  BEM源码\111\N.fig

     文件      10840  2014-09-23 08:16  BEM源码\111\N.png

     文件        200  2014-09-21 18:17  BEM源码\111\N_figure.m

     文件       1079  2014-09-21 08:58  BEM源码\222\t2_1\BEM_DKL.m

     文件       2305  2014-09-21 21:03  BEM源码\222\t2_1\BEM_DKL_d.m

     文件       2173  2014-09-21 21:13  BEM源码\222\t2_1\BEM_DKL_M.m

     文件       1871  2014-09-21 21:30  BEM源码\222\t2_1\BEM_DKL_N.m

............此处省略55个文件信息