function [Outputdata,h]=BELLHOPchannel(Pos,delay,amp,Inputdata,Fs)
t=size(amp,1);
i=size(amp,2);
u=cell(t,1);
for counter1=1:t
u{counter1}=zeros(3,i);
end
for counter1=1:t %取出时延和幅度值
u{counter1}(1,1)=Pos.r.range(counter1);
u{counter1}(2,:)=delay(counter1,:);
u{counter1}(3,:)=abs(amp(counter1,:));
end
row=find(u{t}(3,:)>=3e-4);%找出门限值以上的数(去除幅度较小的值)
Gain1=u{t}(3,row);
Delay1=u{t}(2,row);
RA=angle(amp(t,:));
ra1=RA(row);%自多途时每一径的角度,单位为rad
o=length(Delay1);
%去除自多途
for i=1:o
if(i>length(Delay1))
break;
end
tt=find(abs(Delay1-Delay1(i))<=0.001);
p=find(Gain1(tt)==max(Gain1(tt)));
tt(p(1:end))=[];
Gain1(tt(1:end))=[];
Delay1(tt(1:end))=[];
ra1(tt(1:end))=[];
end
ra2=ra1;%处理后每一径的角度,单位为rad
%求相对时延
g=find(Delay1==min(Delay1));
g=Delay1(g);
for pp=1:length(Delay1)
Delay1(pp)=Delay1(pp)-g;
end
%幅度归一化
g=find(Gain1==max(Gain1));
g=Gain1(g);
for pp=1:length(Gain1)
Gain1(pp)=Gain1(pp)/g;
end
%Delay1
%Gain1
stem(Delay1,Gain1,'.')
Gain=[Gain1];
Delay=[round(Delay1.*Fs)];
h=zeros(1,length(Delay));
h(Delay+1)= Gain;
stem(h);
Outputdata = conv(Inputdata,h);
% Len_in = length(Inputdata);
% Len_out = Len_in +1000;
% Outputdata = zeros(1,Len_out);
%
% for path = 1:length(Delay)
% CHANNEL = [zeros(1,Delay(path)) Inputdata zeros(1,Len_out-Delay(path)-Len_in)];
% Outputdata = Outputdata+Gain(path)*CHANNEL;
% end
%for counter1=1:t
%z=max(r{counter1}(3,:));
%for u=1:i
% if r{counter1}(3,u)~=0
% r{counter1}(3,u)=r{counter1}(3,u)-z;
% end
%end
%end