c0=1500;z0=30;H=100;f=25;w=2*pi*f;N=3;r=1:0.1:10000;n=3;
x1=0;y1=0;TL1=0;TL2=0;
x2=0;y2=0;z=10;
kzn=(n-1/2)*pi/H;
qn=((w/c0)^2-kzn^2)^(1/2);
for n=1:N
x1=2/H*sin(kzn*z)^2*sin(kzn*z0)^2*besselh(0,2,r*qn).*besselh(0,2,r*qn);
y1=x1+y1;
end
TL1=-10*log10(y1);
for m=1:N
for n=1:N
if m~=n
kzm=(m-1/2)*pi/H;
qm=((w/c0)^2-kzm^2)^(1/2);
x2=4/H^2*sin(kzn*z)*sin(kzn*z0)*sin(kzm*z)*sin(kzm*z0)*besselh(0,2,r*(qm+qn));
y2=x2+y2;
end
end
end
TL2=-10*log10(y2);
TL=TL1+TL2;
figure(1);
plot(r,TL)
title('传播损失')
xlabel('Rang(m)');
ylabel('TL(dB)');
c0=1500;z0=30;H=100;f=25;w=2*pi*f;N=3;r=1:0.1:10000;
x1=0;y1=0;TL1=0;TL2=0;
x2=0;y2=0;z=90;
kzn=(n-1/2)*pi/H;
qn=((w/c0)^2-kzn^2)^(1/2);
for n=1:N
x1=2/H*sin(kzn*z)^2*sin(kzn*z0)^2*besselh(0,2,r*qn).*besselh(0,2,r*qn);
y1=x1+y1;
end
TL1=-10*log10(y1);
for m=1:N
for n=1:N
if m~=n
kzm=(m-1/2)*pi/H;
qm=((w/c0)^2-kzm^2)^(1/2);
x2=4/H^2*sin(kzn*z)*sin(kzn*z0)*sin(kzm*z)*sin(kzm*z0)*besselh(0,2,r*(qm+qn));
y2=x2+y2;
end
end
end
TL2=-10*log10(y2);
TL=TL1+TL2;
hold on
plot(r,TL)
legend('10m','90m')
title('传播损失')
xlabel('Rang(m)');
ylabel('TL(dB)');
11111.zip_传播损失图_正弦信号_水声学_海底 声场_简正
版权申诉
5星 · 超过95%的资源 117 浏览量
2022-07-15
01:18:59
上传
评论 3
收藏 2KB ZIP 举报
钱亚锋
- 粉丝: 88
- 资源: 1万+
最新资源
- baseuavAntColonyOptimization-master.zip
- 碳排放权交易明细数据(2024年5月更新).xlsx
- 特殊文件属性命令chattr和lsattr
- HTML、CSS 和 JavaScript动态、交互式的网页 .txt
- b0cd8f9b23d4e5e381b6a8fd8ee0e907.JPG
- ff45d61c5900e45634cf4cac6cff61a1.JPG
- springboot.springboot.springboot.springboot.txt
- linux-进程与服务管理
- 毕业设计基于Django+MySQL+Redis实现简单的天气预报系统python源码.zip
- 基于Streamlit的口罩人脸识别系统python源码+模型+使用说明.zip
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
评论14