216849 (588966), страница 5
Текст из файла (страница 5)
сигнала со сдвигом
Исходный сигнал во временной Восстановленный сигнал
области во временной области
Виды частотных окон Ханна ИПХ фильтров Ханна
для трех диапазонов для трех диапазонов
Для первого диапазона (1-2.5) кГц
Для второго диапазона (2-5) кГц
Для третьего диапазона (4-8) кГц.
Заключение
В данной работе была разработана структура тракта прослушивания гидроакустических сигналов на выходе сформированного пространственного канала (канала наблюдения) в тракте шумопеленгования с использованием многоэлементной антенной решетки для трех частотных диапазонов.
Разработан программный макет тракта прослушивания, удовлетворяющий всем требования для моделирования сигналов и процессов их обработки в тракте ШП.
Установлена взаимосвязь основных параметров тракта прослушивания с базовыми параметрами тракта ШП.
Требования задания полностью выполнены и подтверждены графиками, полученными в результате моделирования.
Список литературы
-
Липатов В.В. Электромагнитные поля в морской воде [Книга]. - Ленинград : ГМТУ, 1990.
-
Рогожников К. И. Морские информационные системы [Книга]. - Санкт-Петербург : АМУР-ПРЕСС, 2002. - стр. 106.
-
Лоскутова Г.В., Полканов К.И. Пространственно-частотные и частотно – волновые методы описания и обработки гидроакустических полей. [Книга]. - Санкт-Петербург .: Наука, 2007.-239с.
-
Ю.А.Корякин, С.А. Смирнов, Г.В.Яковлев. Корабельная гидроакустическая техника: состояние и актуальные проблемы.- СПб.:Наука.-410с.
-
Марпл-мл.С.Л. Цифровой спектральный анализ и его приложения. М.: Мир, 1990.
-
Смарышев М.Д. Направленность гидроакустических антенн. Л. Судостроение, 1973.
-
Гусев В.Г. Системы пространственно-временной обработки гидроакустической информации. Л.: Судостроение, 1988 г.
-
Дьяконов В.П., MATLAB 6 универсальная интегрированная система компьютерной математики. С-Пб: Питер, 2001 г.
-
Лазарев Ю.А. Моделирование процессов и систем в MATLAB. СПб.: Питер, 2005 г
-
Поршнев С.В. MATLAB 7. Основы работы и программирования. М.: Бином-Пресс, 2006 г.
-
Сергиенко А.Б. Цифровая обработка сигналов. СПб.: Питер, 2003 г.
Приложение 1. Текст программы в Matlab
fd=24000;
dt=1/fd;
t=(0:dt:2);
d=0.1;
c=1500;
n=512;
M=30;
jk=sqrt(-1);
df=fd/n;
fk=df*(0:n-1);
f=input ('введите частоту сигнала (1-2.5)kHz-I чд,(2-5)kHz-II чд,(4-8)kHz-III чд =');
fk0=floor(f/df+0.5);
nd=input('Введите номер диапазона=');
k1=[21 42 85];
k2=[53 107 171];
k0_=[7 10 7];
kl_=[10 8 6];
kH=k1(nd);
kB=k2(nd);
kl=kl_(nd);
k0=k0_(nd);
fm=5;
am=0.0;
al=-pi/4;
s=cos(2*pi*f*t);
s=cos(2*pi*f*t).*(1+am*cos(2*pi*fm*t));
tau0=d/c*sin(al);
tau=(0:M-1)'*tau0;
Ttau=repmat(tau,1,size(t,2));
%sound (s,fd)
T=repmat(t,M,1);
x=cos(2*pi*f*(T-Ttau)).*(1+am*cos(2*pi*fm*(T-Ttau)));
fk1=[fk(1:n/2+1) fk(n/2+2:n)-fd];
kolf=exp(jk*2*pi*tau*fk1);
fo=zeros(1,n);
fo(kH-kl:kB+kl)= hann (kB-kH+2*kl+1);
Ro=repmat(fo,M,1);
Per=1*n/2;
n_per=n-Per;
K=Per/2+1;
for i=1:100;
In=(i-1)*n_per+1;
Ik=In+n-1;
Y0=fft(x(:,In:Ik),n,2);
Y=Y0.*Ro;
z=sum(Y.*kolf,1);
V=zeros(1,n);
V(k0:k0+kB-kH+2*kl)=z(kH-kl:kB+kl);
U=ifft(z,n);
Ik1=i*n_per; In1=(i-1)*n_per+1;
w(In1:Ik1)=real(U(K:K+n_per-1));
U=ifft(V,n);
w1(In1:Ik1)=real(U(K:K+n_per-1));
end
nBCE=size(w,2);
figure,plot(fd/nBCE*(0:nBCE-1),20*log10(abs(fft(w))));
xlabel('Гц')
ylabel('дБ')
figure, plot (w)
xlabel('T/dt')
sound (w,fd);
figure,plot(fd/nBCE*(0:nBCE-1),20*log10(abs(fft(w1))));
xlabel('Гц')
ylabel('дБ')
figure, plot (w1)
xlabel('T/dt')
sound (w1,fd);
figure,plot(abs(ifft(fo,512))),grid on;
figure,plot (fo);