tic
clear all;
clc;close all;
f=10; % %%频率f0
w=2*pi*f;
H=50; % %%海水深度为200m
c0=1500; %%初始声速
rou=1024; %%海水密度
dr=1; %%步长?????
rmax=1000; % %%最大水平距离?????
k=w/c0; %%k波束
nmax=floor((f*2*H/c0)+0.5);
N=nmax+10;
khankel=2; %%第二类汉克尔函数
rio=0.8; % %%半径比率
% o=[-2.159,2.513,0.38375];
o=[0,0,-25];
load 'D:\XM\等效源波叠加法\自编\球壳\ndir.mat'; %结构表面法向方向
load 'D:\XM\等效源波叠加法\自编\球壳\uloc.mat'; %结构表面坐标
load 'D:\XM\等效源波叠加法\自编\球壳\U.mat'; %结构表面振速
qloc(:,1)=(uloc(:,1)-o(1))*rio+o(1);
qloc(:,2)=(uloc(:,2)-o(2))*rio+o(2); %等效源坐标
qloc(:,3)=(uloc(:,3)-o(3))*rio+o(3);
% scatter3(uloc(:,1),uloc(:,2),uloc(:,3))
% hold on
% scatter3(qloc(:,1),qloc(:,2),qloc(:,3))
% axis([-4,4,-4,4]
D=JZB(uloc,qloc,ndir,H,c0,khankel,f,N);
ND=(D)^(-1);
det(D)
ND*D;
[su,ss,sv]=csvd(D);
lambda=l_curve(su,ss,U);
Q=tikhonov(su,ss,sv,U,lambda);
jzb (2)_水声简正波_
版权申诉
4星 · 超过85%的资源 143 浏览量
2021-10-03
00:12:48
上传
评论
收藏 1KB ZIP 举报
海四
- 粉丝: 54
- 资源: 4715