Целевая прочность при заданных угле и частоте
tsval = value(tssig,az,el,freq)tsval, силы цели, заданной объектом сигнатуры силы цели, tssig, вычисленный в азимутальном az, вертикальные el, и частотные freq. Если заданные азимут и повышение находятся вне области, в которой задана целевая сигнатура прочности, целевое значение прочности, tsval, возвращается следующим -Inf в дБсм.
perturbations | Возмущение, заданное для объекта |
perturb | Применить возмущения к объекту |
Задайте целевую прочность (TS) жесткого цилиндра длиной 5 м, погруженного в воду, и постройте графики значений TS вдоль азимутального разреза. Предположим, что приближение коротковолновой длины. Радиус цилиндра составляет 2м. Скорость звука - 1520 м/с.
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
