tisi1.m

水声信道传输特性进行仿真的程序

% 功能
% 变量
% 调用的函数
% 运行的结果
clear all;
clc;
close all;


fs=40000;
Data_path='D:\bsdata\shuiping\80zheng\1.';
% Tisi=zeros(1,91);
for i=6;

    data_i=rename(i);
    filename=strcat(Data_path,data_i);

    [hp,hvx,hvy]=xindao2(filename,fs);
    hp1=hp;
    length0=size(hp1,2);
    for n1=1:length0;
        if hp1(1,n1)~=0,break,end
    end
    length1=4000-(length0-n1+1);
    if length1>0;
        hp1=[hp1,zeros(1,length1)];
    end            %n1后的数不足3999个，就补零
    hp1=hp1(:,[n1:(n1+3999)]);
    figure(1);
    plot((0:3999)*1000/fs,hp1);
    xlabel('time(ms)')
    ylabel('幅度')


    [A_max,A_n]=max(abs(hp1));
    A_limit=0.3*A_max;
    A_tisi=find(abs(hp1)>=A_limit);
    Tisi(1,i)=A_tisi(end)*1000/fs;
end
% R=1:1:100;
% figure(2)
% plot(R,Tisi);
% ylabel('Tisi(ms)')
% xlabel('R/Km')

% for n2=(A_n+1):1:3999;
%     if hp1(1,n2)>=A_limit;
%         for n3=(n2+1):1:3999;
%             if hp1(1,n3)>0,break,end
%         end
%         if hp1(1,n3)<=A_limit;
%              
%        break, end
%     end
% end
% tisi1=(n2-A_n)*1000/fs
% for k2=(n3+1):1:3999;
% if hp1(1,k2)>=A_limit;
%         for k3=(k2+1):1:3999;
%             if hp1(1,k3)>0,break,end
%         end
%         if hp1(1,k3)<=A_limit;
%              
%        break, end
%     end
% end
% tisi2=(k2-A_n)*1000/fs
% tisi=(tisi1+tisi2)/2