xiangguanpicture.m

来自「水声信道传输特性进行仿真的程序」· M 代码 · 共 75 行

zizhenyuan_shu=101;fs=40000;
Data_path='d:\bsdata\206\80m\21.';
i=51
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('幅度')

Data_path='d:\bsdata\206\80m\'
n3=1;
for k=21;
    data_k=int2str(k);

    for j=47:1:50
        data_j=rename(j);
        filename=strcat(Data_path,data_k,'.',data_j);

        [hp,hvx,hvy]=xindao2(filename,fs);
        hp2=hp;
        length2=size(hp2,2);
        for n2=1:length2;
            if hp2(1,n2)~=0,break,end
        end
        length3=4000-(length2-n2+1)
        if length3>0;
            hp2=[hp2,zeros(1,length3)];
        end            %n2后的数不足3999个，就补零
        hp2=hp2(:,[n2:(n2+3999)])

        figure(3)
        subplot(4,2,n3)
        plot(abs(xcorr(hp1,hp2,'coeff')));

        n3=n3+1;

    end
        for j=52:1:55
            data_j=rename(j);
            filename=strcat(Data_path,data_k,'.',data_j);
    
            [hp,hvx,hvy]=xindao2(filename,fs);
            hp2=hp;
            length2=size(hp2,2);
            for n2=1:length2;
                if hp2(1,n2)~=0,break,end
            end
            length3=4000-(length2-n2+1)
            if length3>0;
                hp2=[hp2,zeros(1,length3)];
            end            %n2后的数不足3999个，就补零
            hp2=hp2(:,[n2:(n2+3999)])
    
            figure(3)
            subplot(4,2,n3)
            plot(abs(xcorr(hp1,hp2,'coeff')));
            n3=n3+1;
        end
end