xiangguanxishu.m

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

zizhenyuan_shu=101;fs=40000;
Data_path='e:\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='e:\206\80m\'
for k=1:1:41;
    data_k=int2str(k);
    for j=1:1:101
        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(2)
        plot((0:3999)*1000/fs,hp2);
        xlabel('time(ms)');
        ylabel('幅度')
        a(j,k)=max(abs(xcorr(hp1,hp2,'coeff')))
        figure(3)
        plot(abs(xcorr(hp1,hp2,'coeff')));
        %                 if size(hp2,2)>size(hp1,2)
        %                     M=size(hp2,2)-size(hp1,2)
        %                     a(j,k)=max(abs(xcorr([hp1,zeros(1,M)],hp2,'coeff')))
        %                     figure(3)
        %                     plot(abs(xcorr([hp1 ,zeros(1,M)],hp2,'coeff')));
        %                 else size(hp2,2)<size(hp1,2)
        %                     M=size(hp1,2)-size(hp2,2)
        %                     a(j,k)=max(abs(xcorr(hp1,[hp2,zeros(1,M)],'coeff')))
        %                     figure(3)
        %                     plot(abs(xcorr(hp1,[hp2,zeros(1,M)],'coeff')));
        %                 end
    end
    fid=fopen('xishu80m.dat','w');
    fprintf(fid,'%6.4f\n',a);
    fclose(fid)
end