actup2.2水声信道仿真模型

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