CLEAN算法实现超宽带信道估计代码.rar

CLEAN算法实现超宽带信道估计代码.rar

% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % 
% % % % % %   调整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