значение
Класс: tsSignature
Целевая сила под заданным углом и частотой
Синтаксис
tsval = value(tssig,az,el,freq)Описание
Входные параметры
Выходные аргументы
Примеры
Целевая сила твердого цилиндра
Задайте целевую силу (TS) твердого цилиндра 5 м длиной, погруженного в воду, и постройте значения TS вдоль сокращения азимута. Примите приближение короткой длины волны. Цилиндрический радиус составляет 2 м. Скорость звука составляет 1 520 м/с.
L = 5; a = 2;
Создайте массив целевых сильных мест в двух длинах волн. Во-первых, задайте область значений азимута и углов повышения, по которым задан TS. Затем используйте аналитическую модель, чтобы вычислить целевую силу. Создайте изображение TS.
lambda = [0.12, .1]; c = 1520.0; az = [-20:0.1:20]; el = [-10:0.1:10]; ts1 = ts_cylinder(L,a,az,el,lambda(1)); ts2 = ts_cylinder(L,a,az,el,lambda(2)); tsdb1 = 10*log10(ts1); tsdb2 = 10*log10(ts2); imagesc(az,el,tsdb1) title('Target Strength') xlabel('Azimuth (deg)') ylabel('Elevation (deg)') colorbar
Создайте объект tsSignature и постройте сокращение повышения в
азимуте.
tsdb(:,:,1) = tsdb1; tsdb(:,:,2) = tsdb2; freq = c./lambda; tssig = tsSignature('Pattern',tsdb,'Azimuth',az,'Elevation',el,'Frequency',freq); ts = value(tssig,30,el,freq(1)); plot(el,tsdb1) grid title('Elevation Profile of Target Strength') xlabel('Elevation (deg)') ylabel('TS (dBsm)') function ts = ts_cylinder(L,a,az,el,lambda) k = 2*pi/lambda; beta = k*L*sind(el')*ones(size(az)); gamma = cosd(el')*ones(size(az)); ts = a*L^2*(sinc(beta).^2).*gamma.^2/2/lambda; ts = max(ts,10^(-5)); end function s = sinc(theta) s = ones(size(theta)); idx = (abs(theta) <= 1e-2); s(idx) = 1 - 1/6*(theta(idx)).^2; s(~idx) = sin(theta(~idx))./theta(~idx); end
