xindao_x_zhenyuan_y_yuanjianju_pinyu.m

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

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

zizhenyuan_shu=101;
fs=10000;
%Data_path='e:\yao\002.';
Data_path='e:\yao\002.';


i=1
    data_i=rename(i);
    filename=strcat(Data_path,data_i);
    [hp,hvx,hvy]=xindao2(filename,fs);
    hp1=hp;
figure(1);
plot(hp1);


i=2
    data_i=rename(i);
    filename=strcat(Data_path,data_i);
    [hp,hvx,hvy]=xindao2(filename,fs);
    hp2=hp;
figure(2);
plot(hp2);

size(hp1)
a=max(abs(xcorr([hp1 0 0 0 0],hp2,'coeff')));

Hp1=fft(hp1);
Hp1abs=abs(Hp1);
Hp1phase=angle(Hp1);

figure;
subplot(211);
plot(Hp1abs);
subplot(212);
plot(Hp1phase);