声纳技术实验报告四.docVIP

1.实验目的 (1)、学习实际声纳信号的处理方法，加深对理论知识的理解。 (2)、练习MATLAB用于声纳信号处理的方法。 2.基本原理 由以前的实验可知，声纳信号处理的基本方法有谱分析和相关分析，分别为： 基于此，可以对声纳回波信号进行谱分析和互相关分析。 3.实验内容 1、FFT频率估计 CRLB下界 观测时间长度对频率估计性能影响 信噪比对频率估计性能影响 2、相关法时延估计 给出不同信噪比情况下时延估计精度性能Monte-calo分析 分析影响因素 4. 实验结果及数据分析 1、信号建模（sonar.m） clc; close all; clear; SNR=20; f0=15e3; % 单频信号的中心频率 Tw=0.006; % 信号的脉宽 fs=400e3; t=0:1/fs:(Tw-1/fs); Wn=2*[1000 2000]/fs; bb=fir1(128,Wn,bandpass); sigref=0.6*cos(2*pi.*f0.*t); sigref1=sigref sig1=[zeros(1,1000) sigref1 zeros(1,2000)];%发射信号 t1=((0:length(sig1)-1)./fs); noise=normrnd(0,1,1,length(sigref)); noise=filter(bb,1,noise) ; noise=noise/std(noise); noise=noise/10^(SNR/20); sigref2=sigref+noise; sig2=[zeros(1,2000) sigref2 zeros(1,1000)];%接收信号 t2=((0:length(sig2)-1)./fs); noise=normrnd(0,1,1,length(sigref)); noise=filter(bb,1,noise) ; noise=noise/std(noise); SNR=0; noise=noise/10^(SNR/20); sigref3=sigref+noise; sig3=[zeros(1,2000) sigref3 zeros(1,1000)];%接收信号 t3=((0:length(sig3)-1)./fs); noise=normrnd(0,1,1,length(sigref)); noise=filter(bb,1,noise) ; noise=noise/std(noise); SNR=-20; noise=noise/10^(SNR/20); sigref4=sigref+noise; sig4=[zeros(1,2000) sigref4 zeros(1,1000)];%接收信号 t4=((0:length(sig4)-1)./fs); noise=normrnd(0,1,1,length(sigref)); noise=filter(bb,1,noise) ; noise=noise/std(noise); SNR=-20; noise=noise/10^(SNR/20); sigref5=sigref+noise; sig5=[zeros(1,2000) sigref5 zeros(1,3000)];%接收信号 t5=((0:length(sig5)-1)./fs); figure plot(t1,sig1);title(发射信号); xlabel(时间/s) ylabel(幅值/V) grid on figure plot(t2,sig2) title(接收信号(信噪比20，观测时间0.014s)); xlabel(时间/s) ylabel(幅值/V) grid on figure plot(t3,sig3) title(接收信号(信噪比0，观测时间0.014s)); xlabel(时间/s) ylabel(幅值/V) grid on figure plot(t4,sig4) title(接收信号(信噪比-20，观测时间0.014s)); xlabel(时间/s) ylabel(幅值/V) grid on figure plot(t5,sig5) title(接收信号(信噪比-20，观测时间0.02s)); xlabel(时间/s) ylabel(幅值/V) grid on 结果： 2、频谱分析 sonar fftsig=fft(sig2,fs); %fftsig=fft(sig3,fs); %fftsig=