matlab.rar_机械波模型_横波 模拟_横波模拟_流体运动_系统运动方程

% % yq_20158301044_3 clear;clc; % %声偶级（Dipole） theta=0:0.01:2*pi; D=abs(cos(theta)); subplot(4,3,1);polar(theta,D); title('Dipole'); % % 同向小球源辐射（Bipole） %l=1/2λ theta=0:0.01:2*pi; l=1;lambda=2; delta=l/2*sin(theta); k=2*pi/lambda; D=abs(sin(2*k*delta)./2./sin(k*delta)); subplot(4,3,4);polar(theta,D); title('l=1/2λ'); %l=λ theta=0:0.01:2*pi; l=2;lambda=2; delta=l/2*sin(theta); k=2*pi/lambda; D=abs(sin(2*k*delta)./2./sin(k*delta)); subplot(4,3,5);polar(theta,D); title('l=λ'); %l=2λ theta=0:0.01:2*pi; l=4;lambda=2; delta=l/2*sin(theta); k=2*pi/lambda; D=abs(sin(2*k*delta)./2./sin(k*delta)); subplot(4,3,6);polar(theta,D); title('l=2λ'); %%%程序最后有声柱大图代码； %%声柱(固定n=4) %l=1/2λ theta=0:pi/100:2*pi; k=1;lambda=2;n=4; l=1/2*lambda; delta=l/2*sin(theta); D=abs(sin(k*n*delta)./(n*sin(k*delta))); subplot(4,3,7);polar_dB(theta,D); title('l=1/2λ'); %l=λ theta=0:pi/100:2*pi; k=1;lambda=2;n=4; l=lambda; delta=l/2*sin(theta); D=abs(sin(k*n*delta)./(n*sin(k*delta))); subplot(4,3,8);polar_dB(theta,D); title('l=λ'); %l=2λ theta=0:pi/100:2*pi; k=1;lambda=2;n=4; l=2*lambda; delta=l/2*sin(theta); D=abs(sin(k*n*delta)./(n*sin(k*delta))); subplot(4,3,9);polar_dB(theta,D); title('l=2λ'); %%声柱（固定L=4λ） theta=0:0.01:2*pi; k=1;lambda=2;L=4*lambda; n1=2;n2=4;n3=8; l1=L/(n1-1);l2=L/(n2-1);l3=L/(n3-1); delta1=l1/2*sin(theta);delta2=l2/2*sin(theta);delta3=l3/2*sin(theta); D1=abs(sin(k*n1*delta1)./(n1*sin(k*delta1))); subplot(4,3,10);polar_dB(theta,D1);title('n=2'); D2=abs(sin(k*n2*delta2)./(n2*sin(k*delta2))); subplot(4,3,11);polar_dB(theta,D2);title('n=4'); D3=abs(sin(k*n3*delta3)./(n3*sin(k*delta3))); subplot(4,3,12);polar_dB(theta,D3);title('n=8'); %%声柱大图 %%声柱(固定n=4) theta=0:pi/100:2*pi; k=1;lambda=2;n=4; l1=1/2*lambda;l2=lambda;l3=lambda; delta1=l1/2*sin(theta);delta2=l2/2*sin(theta);delta3=l3/2*sin(theta); D1=abs(sin(k*n*delta1)./(n*sin(k*delta1))); subplot(2,3,1);polar_dB(theta,D1);title('l=1/2λ','position', [0.15, -1.3]); D2=abs(sin(k*n*delta2)./(n*sin(k*delta2))); subplot(2,3,2);polar_dB(theta,D2);title('l=λ','position',[0.15,-1.3]); D3=abs(sin(k*n*delta3)./(n*sin(k*delta3))); subplot(2,3,3);polar_dB(theta,D3);title('l=2λ','position',[0.15,-1.3]); %%声柱（固定L=4λ） theta=0:0.01:2*pi; k=1;lambda=2;L=4*lambda; n1=2;n2=4;n3=8; l1=L/(n1-1);l2=L/(n2-1);l3=L/(n3-1); delta1=l1/2*sin(theta);delta2=l2/2*sin(theta);delta3=l3/2*sin(theta); D1=abs(sin(k*n1*delta1)./(n1*sin(k*delta1))); subplot(2,3,4);polar_dB(theta,D1);title('n=2','position',[0.15,-1.3]); D2=abs(sin(k*n2*delta2)./(n2*sin(k*delta2))); subplot(2,3,5);polar_dB(theta,D2);title('n=4','position',[0.15,-1.3]); D3=abs(sin(k*n3*delta3)./(n3*sin(k*delta3))); subplot(2,3,6);polar_dB(theta,D3);title('n=8','position',[0.15,-1.3]);