%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?
%?%?%?%?%?%???調整xt?yt信號的時間軸，相等，以便提取ht信道沖激響應???port_4??????把步長調整為相等得步長???????%?%?%?%?%?%?%?
%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%
%?adjust_double_multipulse3＋clean1124+reversal_signal2
clear?all
close?all
clc
tic
format?long
n=10000;
load?xt.txt
load?o6.txt
cst_out=o6;
%?xt=xxxt;
%?c1=o6;
%?n1=length（c1）;
%?cst_out（1:n11）=c1（1:n11）;
%?cst_out（1:n12）=c1（1:n12）;
%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?
%?%?%?%?%?%????????????port_4??????把步長調整為相等得步長???????%?%?%?%?%?%?%?
%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?
%?figure（1）
%?plot（cst_out（:1）cst_out（:2））
cst_x4=cst_out（:1）;
cst_y4=cst_out（:2）;
x_line4=linspace（cst_out（11）cst_out（end1）n）;
y_line4=spline（cst_x4cst_y4x_line4）;
figure（1）
plot（x_line4y_line4）
%?x_temp4=linspace（x_line4（1）x_line4（end）n2）;
%?y_temp4=spline（
port44（:1）=x_line4*1e-0;
port44（:2）=y_line4;
save???yt1sametime.txt???port44?-ascii;
step_port44=port44（21）-port44（11）
%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?
%?%?%?%?%?%???????port_source????1e0??%?%?%?%?%?%?%?
%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?%?
clear?cst_out;
cst_out（:1）=xt（:1）*1e9;
cst_out（:2）=xt（:2）;
len_time=（cst_out（end1）-cst_out（11））*1e-9;????%time?all?length
step_time=floor（len_time/step_port44）;???????????%time?step?new
x_line_1=0:step_port44:step_time*step_port44;
x_line_1（step_time+1）=cst_out（end1）*1e-9;
y_line=spline（cst_out（:1）*1e-9cst_out（:2）x_line_1）;
figure（2）
plot（cst_out（:1）cst_out（:2））
x_line_2=（x_line_1）;
x_line_2（length（x_line_1）+1）=（step_time+2）*step_port44;
y_line（length（x_line_1）+1）=0;
aa_line=0:step_port44:（step_time+2）*step_port44;
y_line_2=spline（x_line_2y_lineaa_line）;
figure（3）
plot（aa_liney_line_2）
%?x_time=x_line（2）-x_line（1）;
x_len=length（aa_line）
for?i=x_len+1:n
????aa_line（i）=step_port44*i;
????y_line_2（i）=0;
end
port_s（:1）=aa_line;
port_s（:2）=y_line_2;
save???xtsametime.txt???port_s?-ascii;
figure（4）
plot（aa_liney_line_2）
%?%?%?%?%?----------------------------------------------------

%?clc
clear?all
close?all
clc
load?xtsametime.txt?%發射二階高斯信號
%?load?tem_obs.txt?%載入模板信號（距離源1m）
load?yt1sametime.txt?%載入接收信號
%?xt=xt2;
%?xt1（:1）=xt（1:12001）;xt1（:2）=xt（1:12002）;
xt=xtsametime;
yt1=yt1sametime;
%?yt1（:1）=yt1（:1）*1e-9;
figure（3）
plot（yt1（:1）yt1（:2））
tt=yt1（:1）;
dt=tt（2）-tt（1）;
x=xt（:2）;
s=xt（:2）;?%s模板信號

y=yt1（:2）;?%預處理的接收信號
L?=?length（y）;?%L=1001
xxx?=?zeros（3*length（s）1）;
xxx（1）?=?1;
y=conv（yxxx）;?%在y中預留3*length（s）
c=zeros（2*length（y）-11）;?%?產生clean?map，預留一倍空間
Rss=xcorr（s）;?%自相關，Rss

[at]=max（abs（Rss））;?%Rss的絕對值峰值a，出現位置t
%??Rss=Rss./max（abs（Rss））;
Rsy=xcorr（ys）;?%模板信號與接收信號的互相關，Rsy?
d=Rsy;?
[AT]=max（abs（Rsy））;?%Rsy的絕對值峰值A，出現位置T
%?Thre=A;
flag=1;
threshold?=45;?
while?flag
????[aktk]=max（abs（Rsy））

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????8379??2019-07-14?19:47??CLEAN算法實現超寬帶信道估計代碼\clean算法的整體代碼，如果只需要會clean算法，只看main就可以.m

?????文件?????340000??2019-07-14?20:19??CLEAN算法實現超寬帶信道估計代碼\ht1.txt

?????文件?????172654??2012-07-19?19:22??CLEAN算法實現超寬帶信道估計代碼\o6.txt

?????文件??????15295??2012-07-19?19:22??CLEAN算法實現超寬帶信道估計代碼\xt.txt

?????文件?????340102??2019-07-14?19:43??CLEAN算法實現超寬帶信道估計代碼\xtsametime.txt

?????文件?????340000??2019-07-14?19:43??CLEAN算法實現超寬帶信道估計代碼\yt1sametime.txt

?????文件?????367166??2019-07-14?20:19??CLEAN算法實現超寬帶信道估計代碼\yt1_check.txt

?????文件?????679966??2019-07-14?19:45??CLEAN算法實現超寬帶信道估計代碼\yt_reversal.txt

?????文件????????272??2019-07-15?10:04??CLEAN算法實現超寬帶信道估計代碼\先讀這個.txt

?????目錄??????????0??2019-07-15?10:05??CLEAN算法實現超寬帶信道估計代碼

-----------?---------??----------?-----??----

??????????????2263834????????????????????10