pinpu.m

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

zizhenyuan_shu=101;fs=40000;
Data_path='d:\bsdata\shuiping\20zheng\1.';
i=1
data_i=rename(i);
filename=strcat(Data_path,data_i);

[hp,hvx,hvy]=xindao2(filename,fs);
hp1=hp;
length1=size(hp1,2);
for n1=1:length1;
    if hp1(1,n1)~=0,break,end
end
% hp1=hp1(:,[n1:length1]);
hp1=hp1(:,[n1:(n1+3999)])
figure(1);
subplot(311)
plot(hp1);


fft_p1=fft(hp1);
subplot(312);
plot((0:length(fft_p1)-1)/length(fft_p1)*fs,abs(fft_p1));
xlabel('f/Hz');
ylabel('|H(f)|');

subplot(313);
plot((0:length(fft_p1)-1)/length(fft_p1)*fs,angle(fft_p1));
xlabel('f/Hz');
ylabel('arg{H(f)}');

Data_path='d:\bsdata\shuiping\20zheng\1.';
for j=5
    data_j=rename(j);
    filename=strcat(Data_path,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
    hp2=hp2(:,[n2:length2]);
    figure(2);
    subplot(311)
    plot(hp2);

    %信道二图
    fft_p2=fft(hp2);
    subplot(312);
    plot((0:length(fft_p2)-1)/length(fft_p2)*fs,abs(fft_p2));
    xlabel('f/Hz');
    ylabel('|H(f)|');

    subplot(313);
    plot((0:length(fft_p2)-1)/length(fft_p2)*fs,angle(fft_p2));
    xlabel('f/Hz');
    ylabel('arg{H(f)}');

    if size(hp2,2)>size(hp1,2);
        M=size(hp2,2)-size(hp1,2);
        b(j)=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);
        b(j)=max(abs(xcorr(hp1,[hp2,zeros(1,M)],'coeff')));
        figure(3);
        plot(abs(xcorr(hp1,[hp2,zeros(1,M)],'coeff')));
    end
end