球体的声散射 - 声波在目标上的反射与散射仿真程序.docxVIP

程序 1： %刚性球体散射声场形态函数、回波波形仿真 %对应水声学原理第五章 5.5 节的部分内容 clear all close all clc Cw=1500;%水中声速 a=2;%%%%Radius of a sphere 球体半径 p0=1;%%%%Amplitude of incident wave 入射波振幅 r0=2000;%%%%Range between target and observer 目标与观测者之间的距离 f=2:1400;%%%%Sound frequency 声波频率Omega=2*pi*f;%%%%Angular frequency 角频率k=Omega/Cw;%%%%Wave number in water 水中波数x=k*a;%%%%ka N=80; n=(0:N); ConstVector=ones(size(n)); nn=length(f); Vector=ones(size(f)); %jnx=ConstVector*sqrt(pi./(2*x)).*(besselj(n+0.5,x)).;%%第一类球 bessel 函数 jnpx=sqrt(pi./2)*((-n-1)*(x.^(-1.5)).*(besselj(n+0.5,x)).+ConstVector*(x.^(-0.5)).*(besselj(n- 0.5,x)).);%%