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個文件信息