声场防真文件代码.rar_find87z_声场_声场仿真_探头声场_探头声场仿真

InitiateConst; [X,Y,Z,deltaS]=TransDevide(TransR,OpenR,deltaR); Zsize=(ZMax-ZMin)/ZStep+1; Xsize=2*XMax/XStep+1; Ysize=2*YMax/YStep+1; Ssize=size(deltaS,2); H=zeros(Xsize*Ysize*Zsize,Ssize); Pnow=0; for ztmp=ZMin:ZStep:ZMax for ytmp=-YMin:YStep:YMax for xtmp=XMin:XStep:XMax Pnow=Pnow+1; r=sqrt((X-xtmp).^2+(Y-ytmp).^2+(Z-ztmp).^2); Ind=-(j*K+Alpha)*r; H(Pnow,:)=j*Rou*Frequency*exp(Ind).*u/r; end; end; end; SoundP=H*deltaS'; SoundP=reshape(SoundP,Xsize,Ysize,Zsize); SoundI=abs(SoundP).^2/(2*Rou*c); save('intensity','SoundI'); XSound=-XMax:XStep:XMax; YSound=-YMax:YStep:YMax; ZSound=ZMin:ZStep:ZMax; Soundxy=reshape(SoundI(:,:,-ZMin/ZStep+1),2*XMax/XStep+1,2*YMax/YStep+1); figure(1); mesh(XSound,YSound,log10(Soundxy)); figure(2); mesh(XSound,YSound,Soundxy); figure(3); contourf(XSound,YSound,log10(Soundxy)); Soundxz=reshape(SoundI(:,YMax/YStep+1,:),2*XMax/XStep+1,(ZMax-ZMin)/ZStep+1); figure(4); mesh(ZSound,XSound,log10(Soundxz)); figure(5); mesh(ZSound,XSound,Soundxz); figure(6); contourf(ZSound,XSound,log10(Soundxz)); figure(7); Soundz=reshape(SoundI(XMax/XStep+1,YMax/YStep+1,:),1,(ZMax-ZMin)/ZStep+1); plot(ZSound,log10(Soundz)); figure(8); Soundx=reshape(SoundI(:,YMax/YStep+1,-ZMin/ZStep+1),1,2*XMax/XStep+1); plot(XSound,log10(Soundx));