Блок измерителей угловой скорости на основе тетрады ВОГ для космического аппарата

Математическая модель тетрады чувствительных элементов прибора БИУС-ВО. Принцип действия чувствительного элемента прибора БИУС-ВО – волоконно–оптического гироскопа. Разработка методики оценки шумовых составляющих канала измерения угловой скорости.

Рубрика Коммуникации, связь, цифровые приборы и радиоэлектроника
Вид дипломная работа
Язык русский
Дата добавления 24.09.2012
Размер файла 1,7 M

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Rш1, Rш2 - шунтирующий резистор в цепи обратной связи интегратора U3 при разомкнутом и замкнутом ключе (Vкч=0,Vкч=1);

d - минимальная величина входного резистора R3;

С1 - емкость конденсатора в обратной связи интегратора И3;

V01, V02- порог срабатывания и отпускания компаратора К4;

Тк0 - постоянная времени RC-цепи на выходе компараторов К1 и К2.

1.6.5 Математическая модель ЯЛ-121

Схема ячейки логической ЯЛ-121 представлена на Рис. 8. Логическая ячейка построена в виде унифицированного самонастраивающегося устройства, содержащего RS-триггеры (ячейки памяти «», Я6, Я7) Д триггеры (делитель частоты Я1 и формирователи импульсов заданной длительности Я3, Я4, Я5, Я5) и двухвходовые логические элементы «ИЛИ-НЕ» (микросхемы 765ЛЕ5В-1«Н») и «И-НЕ» (микросхемы 765ЛА9В-1«Н»).

Ячейки памяти «» и Я6 выполнены на двух элементах «ИЛИ-НЕ» первые входы которых охвачены перекрестными связями, а вторые входы являются входами ячеек «» и Я6. Ячейка памяти Я7 выполнена на базе Д-триггера (микросхема 765ТМ2В-1«Н») управляющими входами соединенными с общей шиной. ЯЛ-121 имеет два основных выхода q+ и q-, а также один дополнительный выход ИУ. В соответствии со знаком скорости формируется, с одного из основных выходов, последовательность импульсов длительностью 3,6 мкс, частота которых по окончанию переходного процесса пропорциональна измеряемой скорости. С дополнительного выхода формируется импульс управления ИУ типа «меандр» скважностью 2, частота которого в 2 раза меньше частоты импульсов с основных выходов ЯЛ-121. При этом самонастройка ЯЛ-121 осуществляется следующим образом. С момента появления скорости состояние выходов Л1 и Л1 делителя частоты Я1 может не соответствовать (при смене знака скорости всегда не соответствует) направлению скорости. В этом случае, после формирования «ЛОГ. 1» с одного из выходов (Vк1 или Vк2) преобразователя ПС-186 состояние выходов ячейки «» не меняется, поэтому формируется «ЛОГ. 1» с выхода Vк4 преобразователя ПС-186. Состояние выходов Л1 и Л1 делителя частоты Я1 изменится и, тем самым, с дополнительного выхода ЯЛ-121 изменится фаза импульса ИУ и, одновременно, подается «ЛОГ. 0» на одну из ячеек «ИЛИ-НЕ» на основном выходе ЯЛ-121 q+ и q-.

Размещено на http://www.allbest.ru/

Поэтому вновь сформированная “ЛОГ. 1” с одного из выходов (Vк1 или Vк2) преобразователя ПС-186 изменяет состояние выходов ячейки «», и, тем самым, с дополнительного выхода ЯЛ-121 вновь изменится фаза импульса ИУ и, одновременно, с выходов формирователей Я3 или Я4 появляются короткие импульсы, подготавливающие к приходу коротких импульсов Л5 и Л5 соответствующую ячейку памяти (Я6 и Я7), а также формирующие короткий импульс на одном из основных выходов ЯЛ-121 (q+ и q-). При этом частота поочередно следующих коротких импульсов с Я3 и Я4, соответствуют частоте следования импульсов Vк1 и Vк2, соответственно. Частота следования коротких импульсов Л5 и Л5 соответствует частоте следования переднего (заднего) фронта импульса Vк3, формируемого с выхода преобразователя ПС-186 и, тем самым, частота коротких импульсов (длительность 3,6 мкс) на одном из основных выходов ЯЛ-121 удваивается и, тем самым, обеспечивается цена импульса не более 1”, а именно 0,96”.

Математическая модель ЯЛ-121, описывающая формирование импульса управления ИУ и последовательность импульсов частотой fq с выхода, соответствующего направлению скорости, имеет следующий вид.

Импульс управления определяется выражением:

1, если Л2=0, Л2=0

0, в остальных случаях,

где 1, если a=1 (=0), Л1=0

0, в остальных случаях

1, если а=0, Л1=1 (Л1=0);

0, в остальных случаях

а(0), при (tj-1+аб) t<(tj+аб)

(0) , при (tj+аб) t<(tj+1+аб),

1, если (0)=0

0, если (0)=1

t = tj , j= 1,2,3..., если Vki(tj - h) = 0 Vki(tj)=1, i =1,2

Л1(0), при (tк-1+я1) t<(tк+я1)

Л1(0), при (tк+я1) t<(tк+1+я1)

1, если Л1(0)=0

0, если Л1(0)=1

t = tк , к= 1,2,3..., если Vk4(tк - h) = 0, Vk4(tк)=1.

Последовательность импульсов с выходов ЯЛ-121 описывается следующим образом

1, если Л1=1, Л12=0

q+= (2.6.64)

0, в остальных случаях

1, если Л1=0, Л12=0

q-= (2.6.65)

0, в остальных случаях,

Где 1, если Л10=0, Л11=0

Л12= (2.6.66)

0, в остальных случаях

0, если Л3=0, Л4=0

1, в остальных случаях

0, если Л8=1, Л9=1

Л10= (2.6.68)

1, в остальных случаях

0, если Л5=1, Л7=1

1, в остальных случаях

0, если Л5=1, Л6=1

1, в остальных случаях

1, при (tа+я3+я7) t<(tк31+я51+я7)

0, в остальных случаях

t = tа , если а(tа - h) = 0 а(tа)=1

t = tк31 , если Vk3(tк31 - h) = 1 Vk3(tк31) = 0

1, при(t+я4+я6) t<(tк3+я5+я6)

Л6= (2.6.74)

0, в остальных случаях

t = t , если а(t - h) = 1, а(t)=0

t = tк3 , если Vk3(tк3 - h) = 0, Vk3(tк3) = 1

1, при(tк3+я5) t<(tк3+я5+л5),

0, в остальных случаях

1, при( tк31+ я51) t<(tк31+ я51+л51),

0, в остальных случаях

1, при(t+я4) t<(t+я4+л4)

0, в остальных случаях

1, при(tа+я3) t<(tа+я3+л3)

0, в остальных случаях

Текущее значение частоты импульсов с выходов q+ и q- представлено в виде:

fq+ = 1/ (t - t-1)

fq- = 1/ (t 1- t1-1)
t = t , = 1,2,3... при Л1=1, если Л12(t - h) = 1, Л12(t)=0;
t = t1 , 1= 1,2,3... при Л1=0, если Л12(t1 - h)=1, Л12(t1)=0;
Среднее значение частоты импульсов с выходов q+ и q- описывается следующим образом:
max {}
tтек - tнач
max {1}
tтек - tнач
Расчетное значение частоты импульсов с выходов q+ и q- будет иметь вид:
Fqрасч = КВОГ
Где Кс1 Кял
57,3 20
Расчетное значение угла определяется из выражения:
расч = 0t + (1t2)/2 + c (1-cos2fct)/(2fc)

Текущее значение угла, определяемое ВОГом, описывается соотношением

тек = ( - 1)/КВОГ;

Погрешность измерения угла вычисляется следующим образом

= расч - тек

Где а(0) - начальный уровень сигнала “а”;
(0) - начальный уровень сигнала “”;
Л1(0) - начальный уровень сигнала Л1;
Л1(0) - начальный уровень сигнала Л1;
а - запаздывание установки уровня сигналов “а” и“”;
я1 - запаздывание установки уровня сигналов Л1 и Л1;
я3 , л3- запаздывание и длительность высокого уровня сигнала Л3;
я4 , л4- запаздывание и длительность высокого уровня сигнала Л4;
я5 , л5- запаздывание и длительность высокого уровня сигнала Л5;
я51 ,л51- запаздывание и длительность высокого уровня сигнала Л5;
я6- запаздывание установки уровня сигнала Л6;
я7- запаздывание установки уровня сигнала Л7.

1.6.6 Математическая модель УП-118

Схема усилителя - преобразователя УП-118 представлена на Рис. 9.

Размещено на http://www.allbest.ru/

УП-118 [8] формирует с основного выхода знакопеременный прямоугольной формы сигнал Vв(t) и с дополнительных выходов - опорные напряжения +Vоп и -Vоп. УП-118 содержит цепь из суммирующего усилителя УС, инвертора с выхода и входа которого формируются опорные напряжения Vоп, двух, управляемых синхроимпульсом Ф6, последовательно-параллельных коммутатора К1 и К2, выходы которых соединены с входами, управляемого синхроимпульсом Ф7, третьего коммутатора К3, выход которого является основным выходом УП-118. При этом, коммутаторы К1 и К2 на одном из входов содержат, соответственно, повторители П1 и П2 с делителями напряжения на входе, соединенными со входом и выходом инвертора. На инвертирующем входе усилителя УС суммируются сигналы V0 от источника ВИП-019 и сигнал Vф2 со второго выхода ФЧВ-022, тем самым, во время трехкратного увеличения (синхроимпульс Ф6=0) амплитуды сигнала Vв(t) изменится уровень Uоп при отклонении фазовых сдвигов (сигнал Vф2 0). Математическая модель УП-118, описывающая формирование сигнала Vв(t) с периодически в Коп раз изменяемой амплитудой, для слежения за крутизной фазовых модуляторов имеет следующий вид.

Величина опорного напряжения, формируемого с выхода УП-118, определяется соотношением:

V -оп= - (КупV0 + Куп1Vф2(t))

Амплитуда напряжения вспомогательной модуляции Vв управляемая синхроимпульсом Ф6 формируется следующим образом:

- V -оп, если Ф6=0

- V -оп /Коп, если Ф6=1

Выходное напряжение с УП-118 будет иметь вид:

Vв, если Ф7=1

Vв(t) = (2.6.93)

-Vв, если Ф7=0

Где Ф7 - синхронизирующий импульс типа “меандр” частотой вспомогательной модуляции, равной 24 Гц.

Ф6 - синхроимпульс длительностью прерывания, равной двум периодам вспомогательной модуляции и частотой следования, равной 3 кГц (кратной частоте вспомогательной модуляции 24 кГц).

V0 - номинальная амплитуда сигнала Vв(t), приведенная к одному из входов суммирующего усилителя в УП-118;

Куп , Куп1 - коэффициенты передачи суммирующего усилителя;

Коп - коэффициент передачи делителя напряжения в УП-118.

Vф2(t) - напряжение, формируемое со второго выхода ФЧВ-022, поступающее на суммирующий усилитель УП-118 и предназначенное для стабилизации масштабного коэффициента ВОГ.

1.6.7 Математическая модель ФСИ-022

Схема формирователя синхронизирующих импульсов ФСИ-022 представлена на Рис. 10. ФСИ-022 [8] содержит кварцевый генератор Г, многоканальный делитель частоты, реализованный на базе двух микросхем 765ИЕ10В-1 «Н», и логическое устройство реализованное на базе микросхем 765ТМ2В-1 «Н», 765ЛА9В-1 «Н», 765ЛН2В-1 «Н» и 765ЛЕ5В-1 «Н».

ФСИ-022 формирует синхроимпульсы Ф1 - Ф6 и сигналы Ф7, Ф8 типа «меандр» скважностью 2, временные диаграмме которых описываются математической моделью вида:

0, если(t1+1/f48) t( t2+1/f48), где t1=(0.5+1)/f384, t2=1/ f384, 1=0,1,2...

1, в остальных случаях

1, если(t3+2/f3) t( t4+2/f3), где t3=91/ f384, t4= t3+2/f384, 2=0,1,2...

0, в остальных случаях

1, если(t5+3/f3) t( t6+3/f3), где t5=99/ f384, t6= t5+2/f384, 3=0,1,2...

0, в остальных случаях

1, если(t7+4/f3) t( t8+4/f3), где t7=(115+4)/ f384, t8=(115+41)/f384, 4=0,1,2...

0, в остальных случаях

1, если(t9+5/f3) t( t10+5/f3), где t9=(123+5)/f384, t10=(123+51)/f384, 5=0,1,2...

0, в остальных случаях

0, если(t11+6/f3) t( t12+6/f3), где t11=(107+6)/f384, t12=(107+61)/f384, 6=0,1,2...

1, в остальных случаях

0, если М24=0

1, если М24=1

1, если М48=0

0, если М48=1

М24=0.5[1+sign(sin2f24(t-24))], где 24=7/f384

М48=0.5[1+sign(sin2f48(t-48))], где 48=3/f384

Где fj (j=3,24,48,384) - частота типа “меандр", используемая в математической модели для формирования синхроимпульсов Ф1-Ф8.

1,4, 41, 5, 51, 6, 61 - временные задержки;

i i =1-6 - счетчики количества импульсов Фi, где i - соответствует номеру синхроимпульса.

Вторичный источник питания ВИП-019, выполненный по схеме однотактного обратноходового преобразователя, содержит силовой трансформатор, выпрямительные диоды, фильтры разнополярных выходных каналов, блок запуска и схему управления, на управляющий вход которой подается сигнал Ф8 с ФСИ-022.

оптический гироскоп канал измерение скорость

Размещено на http://www.allbest.ru/

Таблица 2.2

Шаг дискретизации

Т0, с

0,5 10-7

Ки1

3,8*106

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

P0, Вт

10-7

Ки2

3,3*104

p

0 (0,25)

Кф

0,2

fp, Гц

24000

Кд3

0,6

p, рад

(0)

Кд4

0,6

dp

1500

Кд5

0,6

np

150

Кд6

0,6

kp

2

Кт

1

pp

0(0,5)

К

0,6

Fpш

1 (0)

Т, с

1,2*10-3

Kc, с

3,56 (-3,56)

Ти1, с

104

Ксном, с

3,56

Ти2, с

104

tc, с

0 (0,2)

Ти11,с

2,65*10-4

tм, с

0 (0,2)

Ти12, с

0,56*10-3

Км1 рад/В

1

Тф,с

0,03

Км2, рад/В

1 (1,2)

Тд3, с

0,5*10-5

0, с

5,28*10-6

Тд4, с

0,5*10-5

Размещено на http://www.allbest.ru/

V0, В

0,375

Тд5, с

0,5*10-5

Куп

1

Тд6, с

0,5*10-5

Куп1

1

Vим0, В

1 (5)

Коп

3

Vим1, В

0 (1;5)

Размещено на http://www.allbest.ru/

Кут1

4

Vим2, В

0 (1;5)

Кут2

3

fим, Гц

0 (10)

Uнас, В

14

Uнас3, В

14

Rш1, Ом

109

Uнас4, В

14

Rш2, Ом

4

Кфд, а/Вт

0,5

R2, Ом

5,62*103

Кпм, Ом

1,3*106

R3, Ом

56*103

Кду1

30

R4, Ом

56*103

Кд

0,6

d, Ом

0

Кду2

17

Vд1, В

1

Тфд, с

10-6

Vд2, В

1

Тпм, с

4*10-2

д, /с

3

Тду, с

6*10-2

С1, Ф

1,54*10-8

Тд, с

10-6

V01, В

0,5

Тду2, с

0,48

V02,В

0,5

Тко, с

2,2*10-5

Uнас0, В

14

Размещено на http://www.allbest.ru/

а(0)

0 (1)

Uнас1, В

14

(0)

1 (0)

Uнас2, В

14

Л1

0 (1)

F3, Гц

3000

Л1

1 (0)

F24, Гц

24000

а

00,2*10-6

F48, Гц

48000

я1

00,2*10-6

F384, Гц

384000

я3

00,2*10-6

1

0

я4

00,2*10-6

4

0

я5

00,2*10-6

41

2

я_5

00,2*10-6

5

0

я6

00,2*10-6

51

2

я7

00,2*10-6

6

0

л3

3,6*10-6

61

16

л4

3,6*10-6

0, /с

0 (1;5)

л5

3,6*10-6

1,/с2

0 (1)

л_51

3,6*10-6

c/с

0 (1;5)

Кял

2

fc, Гц

0 (10)

1.7 Основные результаты математического моделирования

В таблице 2.2 приведены исходные данные для функциональных узлов ВОГ. При этом для обеспечения устойчивости и полосы пропускания ВОГ, близкой к расчетной, в схеме ФЧВ-022 пассивный фильтр содержит конденсатор 0,1 мкф, на входах суммирующего интегратора И1 установлено сопротивление 30 кОм. В МИОЭ крутизна фазовых модуляторов принята равной 1 рад/В. Фильтрующая цепь на неинвентирующем входе интегрирующего сумматора И2 содержит делитель напряжения с коэффициентом передачи 0,2. Номинальный уровень опорного напряжения Uоп, формируемого в УП-118, составляет 1,18 В.

Результаты моделирования ВОГ представлены в Приложении к дипломной работе на рис. 14п-31п в виде графиков номинальной частоты fqрасч и частоты импульсов fq+ и fq- с выходов q+ и q-, соответственно, а также графиков расчетного угла расч и текущего угла тек, определяемого по числу импульсов с выходов q+ и q-.

На рисунках приложения 3 приведены характеристики переходного процесса [3,10] в I и II диапазоне, соответственно. В I диапазоне процесс устанавливается за время 0,12 с, во II диапазоне процесс устанавливается за время менее 0,01 с. С момента времени 0,2 с в 1,2 раза увеличена крутизна фазовых модуляторов относительно номинальной крутизны фазовых модуляторов (1 рад/В), тем самым имитируется дрейф фазовых сдвигов излучений и проверяется работа канала стабилизации масштабного коэффициента ВОГа (КВОГном=3745,5Гц//с).

На рис. приложения 3 приведены графики частоты импульсов fq (рис. 14п) и угла (рис. 15п) при скачке скорости = +1/с.

На рис. приложения 3 приведены графики частоты импульсов fq (рис. 16п) и угла (рис. 17п) при скачке скорости = +5/с.

На рис. приложения 3 приведены графики частоты импульсов fq и угла при переключении диапазонов измерения при достижении скорости = +3/с.

На рис. приложения 3 приведены графики частоты импульсов fq (рис. 18п) и угла (рис. 19п) при переключении из I-го диапазона измерения во II-ой диапазон (скорость 0 = +2,8/с и ускорение 1 = +1/с2).

На рис. приложение 3 приведены графики частоты импульсов fq (рис. 20п ) и угла (рис. 21п) при переключении из II-го диапазона измерения в I-ый диапазон (скорость 0 = +3,2/с и ускорение 1 = -1/с2 ).

На рис. приложения 3 приведены графики частоты импульсов fq (рис. 22п) и угла (рис. 23п) при синусоидальном изменении скорости с амплитудой 3/с и частотой 5,0 Гц.

На рис. приложения 3 приведены графики частоты импульсов fq при скачке напряжения имитации Vим = 1В, при номинальной (рис. 24п) и уменьшенной в 2 раза (рис. 25п) мощности излучения Р0 = 0,510-7Вт.

На рис. приложения 3 приведены графики частоты импульсов fq (рис. 26п) и угла (рис. 27п) при амплитудной модуляции приращения мощности излучения до уровня 1,5Р0 с кратными частотами 1,5 кГц, 3кГц, 6кГц, 12кГц, 24кГц, 48кГц, 96кГц, 192кГц, 384кГц через интервал 0,1с. Таким образом оценивается влияние кратных гармоник мощности излучения на масштабный коэффициент ВОГ.

На рис. приведены графики частоты импульсов fq (рис. 28п) и угла (рис. 29п) при амплитудной модуляции приращения мощности излучения до уровня 1,5Р0 с частотой 24 кГц и фазовым сдвигом рад.

На рис. приведены графики частоты импульсов fq (рис. 30п) и угла (рис. 31п) при амплитудной модуляции приращения мощности излучения до уровня 1,1Р0 с частотой 24 кГц и фазовым сдвигом рад.

Из анализа результатов моделирования следует:

Погрешность измерения угла во время переходного процесса при увеличении скорости уменьшается, а именно, с - 0,02 при скачке скорости +1/с (рис. 1) до - 0,004 при скачке скорости +5/с. При этом в установившемся состоянии частота импульсов fq с выхода ВОГ пропорциональна скорости и определяется соотношением:

fq = fqрасч/(1-)

где fqрасч = КВОГ,

h fqрасч

КВОГ = 3745,5 Гц/(/с) - масштабный коэффициент ВОГ,

h = 0,5 10-7 - шаг дискретизации при математическом моделировании.

Методическая погрешность измерения угла (в установившемся состоянии) определяется следующим выражением:

1/ КВОГ

Выражение (2.5.2) получено из разности между углом, определяемым по числу импульсов с выхода ВОГ и расчетным углом, а именно:

= тек - расч

тек = /КВОГ

расч = Тизм

где интервал измерения Тизм описывается соотношением:

Тизм = Трасч + t

где t Трасч, приращение к сумме расчетных интервалов (Трасч = 1/fрасч) при наличии импульсов в интервале Тизм.

Преобразуя выражение (2.7.3) через (2.7.5) и (2.7.4) с учетом (2.7.6) получим выражение (2.7.2), из которого следует, что методическая погрешность измерения зависит от масштабного коэффициента ВОГ и составляет не более 0,96''.

Дополнительная погрешность измерения угла при наличии дрейфа фазовых сдвигов излучений, имитируемых скачкообразным увеличением в 1,2 раза крутизны фазовых модуляторов, составляет +0,004 и +0,008, при измерении скорости +1/с (рис. 14п-15п) и +5/с (рис. 16п-17п), соответственно. При этом в установившемся состоянии частота импульсов fq с выхода ВОГа пропорциональна скорости и определяется соотношением (2.7.1).

Дополнительная погрешность измерения угла при переключении диапазона измеряемой скорости составляет +0,11 и -0,05 при переключении из I во II диапазон (рис. 18п-19п), и из II в I диапазон (рис. 20п-21п), соответственно.

Амплитудно-фазовые искажения угла при измерении знакопеременной скорости (рис. 22п-23п) характеризуют наличие в 1 диапазоне полосы пропускания не более 10 Гц.

Модель позволяет задавать входное воздействие в виде скорости движения объекта. Для дополнительного исследования интенсивности излучения и параметров УФТ-023 на вход интегратора первого канала ФЧВ-022 подается напряжение имитации (рис. 24п-25п), которое позволяет определить отклонения этих параметров. При штатной работе прибора напряжение имитации должно отсутствовать.

Дополнительная погрешность измерения угла при наличии амплитудной модуляции мощности излучения (рис. 26п-27п) зависит от уровня и фазы амплитудной модуляции на частоте вспомогательной фазовой модуляции, а именно, при уровне амплитудной модуляции мощности излучения 1,5Р0 (рис. 28п-29п), 1,1Р0 (рис. 30п-31п) погрешность измерения угла через интервал 0,1 с. составляет не более 10”, 2.5”, 0.6”, соответственно.

2. Расчетно-конструкторская часть

2.1 Расчет основных параметров оптической цепи прибора.

Основными параметрами оптической цепи волоконно-оптического датчика угловой скорости считаются: длинна измерительного контура (волновода) и диаметр катушки волоконно-оптического контура.

Для определения оптимальных значений этих параметров необходимо рассчитать амплитуду и фазу вспомогательной модуляции.

2.1.1 Определение амплитуды вспомогательной модуляции.

Параметры волоконно-оптического контура (L,D) рассчитываются путем математического анализа схемы обработки сигнала [16].

Зависимость мощнсти излучения, приходящего на фотодиод в идеальном случае имеет вид

Pфд = ?Pист ( 1 +cosДцс )

где Рфд- мощнсть излучения, приходящего на фотодиод,(Вт);

Рист-мощность излучения источника,(Вт);

-фаза Саньяка,(рад).

Рфд=?·2·10-3·(1+cos0,2)=1,9·10-3 (Вт).

Данная зависимость имеет ряд недостатков, таких как: четность функции cos, из чего следует невозможность определения направления движения; производная по Дцс есть sin, это говорит о низкой чувствительности прибора; нелинейность функции. По этому необходим сдвиг фаз на р/2 ( рад).

Статистический сдвиг фаз практически не реализуем. В данном случае применяется динамический сдвиг фаз на цmsinщt.

После динамического воздействия зависимость примет вид

Рфд=?Рист[1+cos(цc + цmsinцt)]

где цm - амплитуда сигнала вспомогательной модуляции,(рад);

щ- частота сигнала вспомогательной модуляции (с-1).

Проведя косинусное преобразование получим

cos(Дце + цmsinщt)=cosДцеcos(цmsinщt)-sinДцеsin(цmsinщt)

Рфд=?Рист[1 + cosДцеcos(цmsinщt)-sinДцеsin(цmsinщt)

Функции sin(цm sinщt) и cos(цmsinщt) можно представить как Функция Бесселя неотрицательного целого порядка. Они могут быть получены как коэффициенты разложения рядов Фурье.

cos(цmsintщ)=є0(цm)+2 є2k (цm) cos2ktщ

sin(цmsintщ)= 2є2k-1цmsin(2k-1)tщ

є, при (|argцm|<р)

где Г'-гамма функция.

После подстановки (3.1..5) в (3.1.4) мы получим следующие выражение

Рфд=?Рист{1+cosДцc[(цm)+2є2(цm)cos2щt+…]

-sinДцc[2є1(цm)sinщt+2є3(ц)sin3щt+…]}

Полученное выражение представляет собой набор гармоник, которые в сумме составляют сигнал, приходящий на УФТ.

Если избавится от вторых гармоник, содержащих чётные функции cos, для чего на выходе фотодиода ставиться полосовой усилитель, настроенный на 1-ю гармонику, то получим мощность пропорциональное

Рфд~Ристє1(цm) sinщt sinДцc

После синхронного детектирования на входе усилителя фототока получаем мощность издучения

Рфд~Ристє1(цm) sinДцc

Наибольшая чувствительность прибора в данном случае будет достигаться при максимальном значении є1(цm).

Анализируя график функции Бесселя первого порядка можно сказать, что его максимальное значение будет при цm=1,84(рад). Подставляя полученное значение в (3.1.7) определим

Рфд~2·10-3·0,6·sin0,2=2,4·10-3 (Вт).

При подаче на модулятор напряжения равного

U=a sinщt,

где a- амплитуда напряжения, (рад);

циклическая частота, (с-1):

щ=2рf

f- частота, (Гц);

фаза каждого луча изменится на соответствующие величины:

для первого луча

ц1=a sinщt

для второго луча

ц1=a sinщ(t+ф)

цm=ц1-ц2=a·[sinщt - sinщ(t+ф)]

Представим разность синусов как произведение , т.е.

sinщt - sinщ(t+ф)=Ѕ sin·cos

Подставляя (3.1.13) в (3.1.12) и произведя некоторые преобразования, получим

цm=

Анализируя полученное выражение (3.1.14) можно заменить, что при

щф=р

мы получим максимальное значение разности динамического воздействия цm. Тогда после преобразования (3.1.14), применив формулы приведения, мы получим

цm==·sinщt.

Пренебрегая переменной составляющей в (3.1.16) видно, что

a=2цm=3,68

2.1.2 Определение диаметра ВОК

Из формулы (2.2.8) выразим произведение диаметра и длины волоконно-оптического контура

В первом приближении функцию sinДцc можно считать прямолтнейной при значении Дцc=±0,2(рад). В ВОГ для удобства и простоты съема сигнала на выходе должна быть прямолинейна зависимость.

Подставляя максимально-допустимые значения Дцс и максимальную измеряемую скорость вращения Щ=3(град./сек)=0,052(рад/сек) в формулу можно получить произведение диаметра и длинны волокна

(м2).

Исходя из технического задания и конструктивных соображений диаметр можно принять равным D=0,15 (м). Тогда допустимая длинна волокна будет равняться

(м).

2.1.3 Выбор частоты вспомогательной модуляции

Зная что оптимальная длинна волоконного контура равна L=1016 (м), можно рассчитать необходимое время обхода контура одним лучом [16]

где L- длинна волоконного контура, (м);

n- коэффициент преломления света в кварце, n=1,46;

с- скорость света в вакууме, с=3·108 (м).

Тогда подставляя значения получим

= 4,94·10-6.

Подставив значение ф в (2.15) получим

=635949,9 (с-1).

Выразим из (2.9а) f

Подставим в полученную формулу щ, получим

(Гц).

Наиболее стабильная частота генерирования имеется у кварцевого генератора, которая равна 384 кГц. Разделив данную частоту на 4 мы получим 96 кГц, наиболее приближенную к найденой. Примем f =96·103 Гц.

2.1.4 Определение длины волоконно-оптического контура.

Так как частота генерирования импульсов вспомогательной модуляции равна 96 кГц, по формуле(3.1.9а) найдём [16]

Щ=2р·96·103= 603185,8 (с-1).

Тогда время обхода контура одним лучом будет равно

ф=р/603185,8=5,2·10-6 (с).

Подставляя полученные значения в (3.1.19) найдём длину волоконно-оптического контура

1070,2 (м).

В приведённых расчётах были определены следующие основные параметры оптической цепи прибора:

длинна волоконно-оптического контура- L= 1070,2 (м);

диаметр волоконно-оптического контура- D= 0,15 (м);

амплитуда вспомогательной модуляции- a= 3,68 (рад);

частота вспомогательной модуляции- f= 96 (кГц).

2.2 Расчёт усилителя фототока.

2.2.1 Переход к прямоугольной вспомогательной модуляции.

Вспомогательное синусоидальное динамическое воздействие имеет ряд недостатков, таких как:

необходимость применения полосового усилителя, что усложняет схему усилителя фототока (далее по тексту УФТ);

остаточное присутствие второй гармоники на выходе УФТ, несмотря на применение полосовых усилителей;

необходимость съема сигнала на линейных участках синусоиды.

Как видно на рис.1(приложение 1) при появление фазы Саньяка (то есть наличие вращательного движения вокруг чувствительной оси прибора) на выходе интерферометра появляется первая гармоника, на которой сказывается влияние второй, в виде дополнительных парных пиков на максимумах. Хотя на чистоту первой настроены полосовые усилители, полностью отсечь влияние не желательных воздействий не удается.

В альтернативу синусоидальной вспомогательной модуляции можно сопоставить прямоугольную форму модулирующего напряжения, которая имеет ряд дополнительных преимуществ.

Применение прямоугольной формы модуляции позволяет избавится от вредного воздействия второй гармоники. Что несомненно упрощает электрическую схему усилителя фототока, отсутствием потребности и необходимости полосовых усилителей.

При переходе на прямоугольную форму модулирующего напряжения можно уйти в область более низких частот. В данном приборе реализована частота 24 кГц. В приложении приведены графики , анализ которых даст наглядное представление о преимуществах применения прямоугольного динамического воздействия.

При использовании прямоугольной модуляции фаза Саньяка Дцс будет иметь вид функции типа «меандр». Что значительно упрощает усиление и съем сигнала.

2.2.2 Расчет усилителя фототока.

Расчет усилителя фототока произведем исходя из требуемого коэффициента усиления сигнала, который пропорционален измеряемой угловой скорости вращения Щ.

Излучение, испускаемое источником, приходит на фотодиод со значительным ослаблением мощности. Это связанно с оптическими потерями в волоконно-оптическом тракте. Основные потери оптического тракта можно свести в таблицу.

Таблица 3.1

№ п/п

Наименование элемента оптического тракта.

Количество проходов через элемент.

Величина потерь, дб

1

Деполяризатор.

1

0.5.

2

Ответвитель.

2

1

3

МИОЭ.

2

7

4

ВК.

1

4

5

Ответвитель (деление излучения)

2

3

6

Сварка волокна.

9

0.3

Определим сумму потерь.

= 0,5 + 1·2.+ 7·2+4+3·2+0,3·9=29,2 (дб).

Рассчитаем коэффициент ослабления излучения испускаемого источника при прохождении оптического тракта логарифмическим преобразованием

=10 lg косл

где косл- коэффициент ослабления излучения оптическим трактом.

Выразив косл через получим

косл =10/10=102,9=794,3.

Зная мощность источника излучения Рист= 2мВт и коэффициент ослабления можно определить мощность принимаемую фотодиодом. Так же необходимо заметить, что фотодиод принимает не весь диапазон мощности, а только его часть (1/4 max). Наглядно в этом убедиться можно проанализировав рис.16. Это связано с тем, что измеряемая фаза Саньяка находится в пределах. Мощность излучения на фотодиоде определим по формуле [16]

Рфд=

Подставляя значения получим

Рфд=

Ток выделяемый фотодиодом равен произведению мощности излучения и коэффициента передачи фотодиода (Кфд= 0,8 А/Вт)

Подставляя значения получим

Нагрузочное сопротивление фотодиода равно Rф = 43·103 (Ом).

Зная сопротивление и силу тока по закону Ома можно определить

Vmax = J·Rф

Подставляя значения получим

Vmax = 0,5·10-6·43·103 =21,5·103 (B)

Максимальное напряжение сигнала используемое для последующих преобразований должно равняться Vmax = 15 В.

Таким образом зная необходимое напряжение и напряжение на выходе фотодиода можно определить необходимый коэффициент усиления УФТ.

Куфт=

где V max-необходимое напяжение, В;

V фд max- напряжение на выходе фотодиода.

Куфт=

Приняв количество каскадов равное трем, выберем коэффициент усиления каждого каскада Кус=29 для двух и один каскад с регулируемым коэффициентом усиления. Это допущение необходимо потому, что в зависимости от комплектации прибора коэффициент ослабления мощности излучения может меняться.

Выбрав схему усилителя, приведенную на рис.15, найдем номиналы сопротивления. Коэффициент усиления при таком включении равен

Кус=

Номинал сопротивления R2 выберем равным 43 кОм.

Величина этого сопротивления влияет на ток обратной связи. Таким образом, выразив R1 из (2.20) получим

R1=

Подставляя значения

R1= (кОм)

На выходе каждого операционного усилителя поставим фильтр (рис.16), настроный на частоты от 49с-1 и выше.

Постоянная времени равна

Т=

Подставляя значение получим

Т=

Таким образом

RC=T

Для того, чтобы не покупать резисторы разных номиналов выберем

R=43кОм.

Тогда из (3.2.10) выразим

С= Т/R

С =

Для надежности работы УФТ емкость продублируем.

Второй каскад усилителя выберем с такой же схемой и номиналами.

После второго каскада установим синхронный детектор работающий по сигналу с синхроимпульсного блока. На входе детектора установим разделяющее сопротивление с номиналом 1500 Ом.

На выходе детектора установим «ключ», который будет удерживать потенциал получения сигнала. Это необходимо для того, чтобы далее схема обработки сигнала работала со стробами длительностью TRC.

Подставляя значения получим

T=

Таким образом номиналы элементов необходимо подбирать исходя из достаточной величины постоянной времени RC.

Так как время обхода контура одним лучом равно ф=5 мкс, то

TRC=2,5·10-6

Зададимся большим сопротивлением R=510 кОм и вычислим величину емкости по

С=

После каждого каскада на опорное напряжение установим фильтр,для того, чтобы работа одного каскада никак не сказывалось на работе другого.

Принятые значение элементов указаны в спецификации на электрическую схему УФТ.

2.3 Расчет фазового модулятора.

Основным параметром фазового модулятора является коэффициент передачи

Кфн= /V

Произведем расчет коэффициента передачи исходя из максимальной измеряемой скорости [16]. Максимально измеряемая скорость в данном приборе равна

= 3 /с=10800 дуг.с.

В связи с тем, что масштабный коэффициент электронной части равен км1 (дуг.с/имп.),то частота следования импульсов на выходе прибора должна быть равна 10800 Гц.

Частота пилы компенсирующей модуляции в 2 раза меньше. Тогда период модуляции будет равен

Подставляя значения

Зная амплитуду моделирующего напряжения (V=5В) и время обхода контура лучом света можно составить пропорцию

Выразим из этого соотношения V

Подставляя значения

V

Подставляя полученные значения в (3.3.1) найдем

Кфм=0,191/0,141=1,36(рад/В).

В предложенной схеме обработки сигнала фазовый модулятор является одним из составляющих элементов МИОЭ. Поэтому рассчитанный коэффициент будет являться необходимым требованием к МИОЭ.

2.4 Анализ погрешностей

2.4.1 Разработка методики оценки шумовых составляющих канала измерения угловой скорости.

Модели, дискретная и непрерывная, полученные выше, являются идеализированными, т.к. учитывают динамические свойства каждого первичного измерителя и не учитывают его помех, привносимых в штатную информацию в процессе реальной работы.

Волоконно-оптический гироскоп представляет собой довольно сложную оптико-электронную систему. При конструировании реального прибора оптические элементы и электронные устройства должны выбираться и компоноваться так, чтобы минимизировать влияние внешних возмущений (температурных градиентов, механических и акустических вибраций, магнитных полей и др.). В самом приборе, кроме этого, существует ряд внутренних источников шумов и нестабильностей.

Условно эти шумы и нестабильности можно разделить на быстрые и медленные возмущения. Быстрые возмущения оказывают случайное кратковременное усредненное влияние на чувствительность ВОГ; они отчетливо проявляются при нулевой скорости вращения.

Медленные возмущения вызывают медленный дрейф сигнала, приводящий к долговременным уходам в считывании показаний ВОГ.

Для ВОГ принципиально неустранимыми шумами являются так называемые квантовые или фотонные шумы (их также называют дробовыми шумами). Эти шумы появляются лишь в присутствии полезного оптического сигнала на входе фотодетектора и обусловленны случайным распределением скорости прихода фотонов на фотодетектор, что приводит к случайным флуктуациям тока фотодетектора. Фотонные шумы являются следствием квантовой природы светового излучения. Уровень фотонных шумов зависит от интенсивности оптического излучения.

В практике экспериментальных исследований было обнаружено, что одним из главных источников шума в системе ВОГ является обратное рэлеевское рассеяние в волокне, а в некоторых системах еще и отражение от дискретных оптических элементов, используемых для ввода излучения в систему. Физически эти шумы появляются из-за рассеяния светового излучения прямого луча на микрочастицах и неоднородностях среды распространения.

Температурные градиенты, изменяющиеся во времени в волоконном контуре, приводят к появлению сигнала, эквивалентного некоторому значению скорости вращения, т. е. вызывают ошибку в считывании угловой скорости. Весьма сушественна в ВОГ стабильность масштабного коэффициента Квог. Согласно (2), (10) постоянство масштабного коэффициента определяется стабильностью площади витка контура и длины волны . Площадь витка является функцией температуры и материала катушки, на которую наматывают контур. Стабильность длины волны зависит от излучения оптического источника. Следует отметить, что шумы, связанные с изменением длинны волны излучения источника ВОГ, незначительны в большинстве систем регистрации фазы.

Рассмотрим теперь шумы, появляющиеся в ВОГ из-за нелинейного характера взаимодействия излучения со средой, в которой оно распространяется.

Несмотря на очень низкие уровни излучения, распространяющегося в ВОГ, нелинейные эффекты могут быть весьма значительными. Нелинейный электрооптический эффект носит название Керра и состоит в изменении фазового набега световой волны, распространяющейся в среде под действием интенсивности излучения, при этом фаза изменяется в зависимости от квадрата амплитуды излучения. При исследованиях ВОГ было замечено, что эффект Керра вносит значительный вклад в паразитный дрейф прибора.

Нестабильность характеристик ВОГ, приводящая к дрейфу в приборе, может быть обусловлена влиянием внешнего магнитного поля (эффект Фарадея).

При механическом несовершенстве конструкции ВОГ серьезным источником шумов могут быть акустические поля, механические вибрации и ускорения.

Для полноты описания возмущений ВОГ следует сказать о таких возмущениях как низкочастотный шум фотодетектора, спонтанные и стимулированные шумы лазерного источника излучения, рассеяние на фотонах и неоднородностях в среде. Однако, практически, уровень интенсивности этих шумов невысок.

Такова общая картина источников шумов и нестабильностей ВОГ. В зависимости от варианта конструкции ВОГ те или иные источники шумов и нестабильностей могут играть большую или меньшую роль. Основными источниками являются шумы обратного рэлеевского рассеяния, нелинейный электрооптический эффект, температурные градиенты, внешнее магнитное поле, а также нестабильность интенсивности и длины волны источника излучения. Принципиально неустранимым шумом является фотонный шум полезного сигнала, появляющийся в системе регистрации и определяющий фундаментальный предел чувствительности (точности) ВОГ.

Как следует из вышеизложеного, полное математическое описание и моделирование шумов в ВОГ представляет собой самостоятельную и достаточно сложную задачу. Для синтеза математической модели прибора, работающей в масштабе времени, близком к реальному, данная задача практически не реализуема. Для практики имеет значение комплексная оценка шумов ВОГ и наложение этих шумов на реальный выходной сигнал прибора.

Как правило, помехи носят характер случайных процессов во времени. Для разработки математической модели шумовой составляющей реального прибора существует два подхода.

1. Модель составляется для источника шума внутри каждого измерительного канала, исследуется его прохождение через передаточную функцию канала, приведенную к источнику шума, а на выходе полученные случайные величины суммируются и оценивается функция плотности распределения вероятности суммарной помехи, по которой можно восстановить в реальном масштабе времени случайную функцию, соответствующую совокупной шумовой составляющей прибора.

2. Производится непосредственный замер выходного сигнала изделия, выделяется систематическая и шумовая составляющие, оценивается среднеквадратическое отклонение, соответствующее случайной составляющей в запуске, и полученные значения накладываются на выходную информацию в математической модели прибора.

Оба эти метода имеют как свои достоинства, так и недостатки. К достоинствам первого метода следует отнести то, что он может дать реальную картину и позволяет оценить природу и “вес” каждой составляющей помехи. К недостаткам - то, что исследователь вынужден доверять паспортным характеристикам элементов, привносящих шумовую составляющую в измерительный канал, проверить которую не всегда представляется возможным. К достоинствам второго метода следует отнести то, что полученный результат - есть реальный результат, снятый с конкретного прибора. Однако, недостатком является тот факт, что число экспериментов ограничено и повторяемость результата всегда может быть поставлена под сомнение. Кроме того, можно оценить лишь совокупный шум, так как разделить помехи не представляется возможным. Очевидно, что наиболее оптимальным является вариант исследования двумя этими методами с последующим сравнением результатов.

В настоящей работе используется комплексный подход к вышеозначенной проблеме, заключающийся в следующем.

С одного из образцов реального волоконно-оптического гироскопа через плату АЦП снимается шумовой сигнал, носящий характер случайного процесса во времени.

При помощи разработанного программно-математического обеспечения строится функция плотности распределения вероятности и определяется характер распределения.

Разрабатывается математическая модель входного комплексного сигнала помехи с некоторым толерантным интервалом распределения.

При помощи метода импульсных переходных функций определяется реакция каждого измерительного канала на входное возмущающее воздействие при помощи разработанного программного обеспечения.

Полученных выходной сигнал помехи накладывается в каждый момент опроса математической модели ИУС ВО на выходной сигнал, тем самым достигается наибольшая достоверность выходной информации прибора, наиболее приближенная к реальности.

Оценка шумовой составляющей волоконно-оптического гироскопа.

В рамках решаемой задачи с целью формирования реальной шумовой составляющей были проведены замеры с аналогового выхода ВОГ посредством АЦП ЛАИ-24 с частотой опроса 30 Гц и временем опроса - 15 мин. В результате было получено 27000 мгновенных значений шумовой составляющей. Полученные значения были записанны в файл “ shum.dat “.

На Рис.13 приведен график шумовой составляющей ВОГ, построенный с помошью Excel по значениям, сформированным выше указанным образом.

С помошью программы “funcrasp.cpp” по массиву случайных величин из файла “ shum.dat “ ( логическая схема и листинг программы приведены ниже в соответствующих разделах данного диплома) была сформирована функция плотности распределения вероятости. На Рис. 15 приведен график функции плотности распределения вероятности шума, из котороговидно, что случайный процесс распределен по нормальному закону в толерантном интервале [0.003652 0.006654] с математическим ожиданием 0.005171 и среднеквадратическим отклонением 0.000374 . Так как выходным сигналом являлось напряжение, то размерностью данных величин являются вольты. Масштабный коэффициент, связывающий угловую скорость и напряжение равен 0.6, так как угловой скорости 3/c соответствует выходное напряжение 5 В. Отсюда, среднеквадратическое отклонение шумовой составляющей ВОГ ВОГ=0.8078/ч.

Шумовая составляющая ВОГ Рис. 13

Рис.14

GL=0.003652В М=0.005171В =0.000374В GP=0.006654В

Таким образом, на базе анализа экспериментальных данных получено, что шумы ВОГ распределены по нормальному закону.

Однако данные значения взяты с конкретных образцов приборов и от испытания к испытанию могут изменяться. Кроме того, в случае использования математической модели может возникнуть необходимость в составе математической модели изделия произвольно формировать шумы определенного уровня, а не использовать значения конкретного эксперимента. С этой целью предлагается алгоритм и программа формирования случайного процесса по заданному толерантному интервалу с нормальной функцией плотности распределения вероятности.

Алгоритм задания помех в виде случайной функции с некоторым законом распределения плотности вероятности.

Предлагаемый алгоритм состоит из следующих шагов:

Задаются значения математических ожиданий М и среднеквадратических отклонений , а также число разбиений N интервала распределения.

C шагом

рассчитывается по формуле (32) массив значений для функции плотности распределения вероятностей (i) ; (i=1, N+1).

Производится нормирование полученного массива таким образом, что значения пронормированного массива * (i) располагались в интервале [0;1].

Нормирование производится согласно соотношениям:

!Ошибка в формуле

Здесь в числителе стоит частичная сумма элементов массива (i) на каждом шаге изменения i. В знаменателе - полная сумма этого массива. В результате выполнения этой процедуры значения пронормированного массива * (i) будут лежать в интервале [0;1] и соответствовать значениям требуемой функции плотности распределения вероятностей каждой из помех.

В алгоритмическом языке СИ [6] имеется функция генерации случайных чисел random (i,j), при помощи которой можно получить некоторый массив случайных чисел, распределенных по равномерному закону в интервале [0;1]. Присваивая некоторому числу X значение random (i,j), получаем случайное число, лежащее в интервале [0;1]. Т.к. задачей является получение случайной функции (t), распределенной по нормальному закону, то задача состоит в преобразовании равномерного закона распределения случайной функции X (t) на интервале [0;1] в нормальный закон распределения функции (t) на интервале [М - 3 ; М +3] и получение значений функции (t).

Организуется цикл получения случайных значений X= random (i,j), причем число этих значений соответствует числу точек рассматриваемого временного интервала - N1. На каждом шаге производится сравнение числа X с соответствующим значением * (i) и определяется значение k - порядковый номер значения * (i).

Из соотношения i = kH+ (М -3)+XH определяется значение функции i в момент времени ti. Здесь (М -3) определяет левую границу интервала. KH определяет левую границу попадания случайной величины в соответствующий ей интервал разбиения отрезка [М -3; М +3], а величина XH определяет место случайной величины в данном интервале.

На Рис.15 показан график выходного сигнала ВОГ, снятый с 1 канала тетрады при подаче по оси Х постоянно действующего сигнала, угловой скорости, величиной 3/c, в 1 диапазоне работы каждого измерительного канала. Воздействия по осям Y и Z принимались равными нулю. Помеха накладывалась из экспериментального файла данных, соответствующего каждому измерительному каналу. Кривая 1 на каждом рисунке соответствует идеальному сигналу прибора, кривая 2 - сигналу прибора с наложением шумовой составляющей. На Рис.16 помеха накладывалась методом моделирования с заданными математическим ожиданием и СКО, близкими по своим значениям к экспериментальным. Из графиков видна практическая идентичность результатов, что говорит о корректности предлагаемого метода математического моделирования помехи.

Выходной сигнал ВОГ. Рис. 15

кривая 1 - идеальный сигнал

кривая 2 - сигнал с наложенной шумовой составляющей

Выходной сигнал ВОГ Рис. 16

кривая 1 - идеальный сигнал

кривая 2 - сигнал с наложенной шумовой составляющей, сгенерированной математическим методом.

Алгоритм оценки выходного сигнала ВОГ в случае случайного характера помех на входе канала измерения угловой скорости.

Для определения вероятностных характеристик погрешностей любой системы автоматического регулирования, к которой в нашем случае относится ВОГ, необходимо задать вероятностные характеристики входных сигналов, а также динамические свойства системы, определяемые ее передаточными функциями.

Если на вход системы поступает стационарная случайная функция X(t), а с выхода снимается случайная функция Z(t), то их взаимосвязь можно выразить через импульсную переходную функцию системы выражением

Z(t)=

где К() выражается через передаточную функцию Ф(j) системы формулой

K(t)=

Передаточную функцию Ф(j) можно представить в виде:

Ф(j)=R()+jS(), a ejt=cost + jsint, тогда вместо варажения (3.4.6) получим

K(t)=

Второй интеграл равен нулю, так как функция времени К(t) является действительной, следовательно,

K(t)=

Таким образом, поставленная задача сводится к вычислению двух неопределенных интегралов (3.4.6) и (3.4.8). Однако специфика решения настоящей задачи заключается в том, что динамические свойства системы определяются ее полосой пропускания, что говорит об ограниченности частотного диапазона (в нашем случае от 0 до 20 Гц). Поэтому пределы интегрирования в выражении (3.4.8) будут: нижний 0, а верхний 220.

Вычисление интегралов (3.4.6) и (3.4.8) не может ограничиться вычислениями значений подынтегральных выражений в узлах [5]. Необходимо аппроксимировать значение подынтегральной функции. В качестве метода аппроксимации был выбран метод сплайн - аппроксимации, как наиболее оптимальный для решения задач подобного класса [2]. Кратко изложим суть данного метода.

Под сплайном SDm(f, x) [2] будем понимать функцию, являющуюся многочленом степени m на каждом из отрезков [xi-1, xi]:

SDm(f, x) = Pim(x) = ai0 + …+aimxm при xi-1<x<xi,

и удовлетворяющую условиям непрерывности производных до порядка m-1 в точках x1, … xi-1:

(xi) = Pi+1

При k=0, 1, …m-1; i=1, 2, …N-1.

Предположим, что отрезок, на котором определена исходная функция плотности распределения вероятности (ф.п.р.в.) [3], разбит на N интервалов, т.е. i изменяется от 1 до N+1, где N+1- число узлов.

На каждом интервале xi-xi-1 будем аппроксимировать функцию кубическим сплайном, соблюдая условие (3.4.9), т.е. в данном конкретном случае, соблюдение условий непрерывности в узлах по 1-й и 2-й производным.

Допустим, что S(xi) = fi, а SІ(xi) = Mi. Т.к. кубический сплайн - кривая 3-го порядка, то её 2-я производная меняется по линейному закону

SІ(x) = Mi + Mi-1

Определим функцию S(x) как функцию вида:

S(x) = Mi-1 + Mi + c1x + c2

Внутри интервала [xi-1, xi], используя условия непрерывности функции S(x) на концах данного интервала S(xi) = fi и S(xi-1) = fi-1. Подставив в (3.4.11) вместо x соответственно значения xi и xi-1, получим:

S(xi) = Mi-1 + Mi + c1x + c2;

S(xi-1)= Mi-1 + Mi + c1x + c2.

Откуда имеем:

fi - Mi = c1x + c2

fi-1 - Mi-1 = c1xi-1 + c2

Полученные выражения можно рассматривать как требования для некоторой линейной функции c1x + c2, чтобы она в точке xi проходила через ординату

fi - Mi , а в точке xi-1 - через ординату fi-1 - Mi-1.

Очевидно, что такой функцией является функция вида:

c1x + c2 = +

S(x) = Mi-1 + Mi + + для xi-1 Ј x Ј xi.

Предполагая, что шаг интерполяции постоянен, и, обозначив xi - xi-1 = hi, получим

S(x) = Mi-1 + Mi + + для

xi-1 Ј x Ј xi

Таким образом, из уравнения (3.4.14) мы можем определить значение функции S(x), аппроксимированной сплайном. Необходимо определить значения коэффициентов Mi (второй производной). Используя уравнение (3.4.9) для кубического сплайна, получим

S'(x) = - Mi-1 + Mi + ·

Подставив в (3.4.15) xi, и, решив его относительно xi и xi-1, и, приравняв эти два значения, получим

Mi + + Mi-1 - Mi = -Mi + + Mi + Mi+1;

Полагая hi+1 = xi+1 - xi, получим

Mi-1 + Mi + Mi+1 =

- краевые условия

Таким образом, решение задачи аппроксимации некоторой заданной функции сплайном сводится к решению системы линейных алгебраических уравнений (3.4.16) при условии равенства нулю значений вторых производных на границах интервала [2]. Полученную систему можно записать в виде

CM = d

Где d = (d1, …dN-1)Т, где di вычисляются по формуле:

di =

Элементы матрицы C = (cij)i,j=1 задаются соотношениями

cij =

Так как шаг интерполяции считается постоянным, т.е. hi = hi+1 = h, получим

di =

cij =

Уравнение (3.4.18) примет вид

Умножив обе части на , получим

Mi-1 + 4Mi + Mi+1 =

С учётом выражения (3.4.22) можно записать новый вид

di = 6

cij =

Систему уравнений (3.4.17) можно решить методом прогонки [2]; в этом случае, используя (3.4.23) и (3.4.24), перепишем систему (3.4.17) в виде:

Mi-1 + 4Mi + Mi+1 = di

Введём новые обозначения:

Xi = Ai

f(xi) = Bi

Прогоночные коэффициенты определим из соотношения:

Mi-1 = MiAi + Bi

MiAi + Bi + 4Mi + Mi+1 = di

Mi(Ai+4) + Mi+1 + Bi = di

Откуда

Mi = = Mi+1Ai+1 + Bi+1

Заметим, что при i=1 A1 = 0 и B1 = 0, тогда A2 = - ; B2 =

Применение метода сплайн - аппроксимации обеспечивает необходимую точность и быстродействие для решения поставленной задачи [2] и определения массивов значений погрешности прибора, накладываемых на выходную информацию в каждый момент опроса математической модели ИУС ВО. Следует отметить тот факт, что данная задача решается автономно, в математической модели ИУС ВО используются только конечные результаты в виде файлов данных, к которым обращается головная задача. Это имеет свои преимущества, позволяя исследователю варьировать вариантами задаваемых помех, не усложняя головной задачи и не замедляя тем самым ход ее решения.

2.4.2 Расчет надежности при отказах

В данном подпункте рассматривается сравнительный анализ применения двух вариантов избыточных систем, а именно триады ортогонально расположенных датчиков первичной информации и тетрады. Тетрада образуется ось одного из датчика совпадает с одной из осью объекта, а три других расположены по образующей конуса с углом полураствора 54,7?.

Гироскопические датчики первичной информации (ГДПИ) являются достаточно сложными электромеханическими приборами и оптическими приборами, и уровень их надежности ниже требований, предъявляемых к уровню надежности системы и ее элементов в целом. Эффективным способом преодоления этого недостатка ГДПИ является применение структурного резервирования, позволяющего повысить такие показатели эффективности, как точность и надежность. Правда, это увеличивает массу и габариты прибора, поэтому кратность резервирования в приборе не превышает трех, т.е. K?3.

Надежность измерения параметров движения КА при структурном резервировании ГДПИ может быть увеличена принципиально двумя путями.

Первый путь заключается в резервировании самих ГДПИ в пределах одного измерительного базиса. Второй путь - в резервировании измерительных базисов.

Рассмотрим бесплатформенный инерциальный блок (БИБ) с произвольным числом одноосных измерителей N. Для оценки ортогональных составляющих измеряемого вектора в базовой системе координат будем использовать метод наименьших квадратов. При равноточных измерителях и аддитивных ошибках измерения с нулевым математическим ожиданием и дисперсией, равной 2, ковариационная матрица ошибок оценки ортогональных составляющих измеряемого вектора имеет вид:


Подобные документы

  • Конструкция и принцип действия поплавкового датчика угловой скорости КХ79-060. Расчет потребляемой мощности, коэффициента демпфирования и момента инерции поплавкового гидроузла. Математическая модель ДУС с цифровой обратной связью. Анализ погрешностей.

    дипломная работа [1,6 M], добавлен 23.01.2012

  • Противоречивые требования, предъявляемые к системе стабилизации линии визирования. Задача эффективного преобразования сигнала угловой скорости гироскопа в цифровую форму. Выбор элементной базы для аппаратной реализации на основе поставленных требований.

    дипломная работа [2,5 M], добавлен 01.12.2014

  • Модернизация поплавкового датчика угловой скорости (ДУС) путем введения цифровой обратной связи, разработка его структурной схемы с процессором. Математическая модель ДУС с цифровым регулятором. Расчет основных параметров. Анализ погрешностей датчика.

    дипломная работа [3,2 M], добавлен 30.01.2012

  • Частотный метод измерения высоты и составляющих скорости. Канал оценки составляющих скорости. Вычислительные требования к блоку измерителя и модуляции. Разработка схемы электрической принципиальной. Математическое моделирование усилителя ограничителя.

    дипломная работа [861,7 K], добавлен 24.03.2014

  • Электромагнитные тахометры угловой скорости. Тахометрический генератор постоянного тока. Тахометрические генераторы на переменном токе. Электромагнитные тахометры линейной скорости. Импульсные тахометры угловой скорости. Гирометры.

    курсовая работа [2,7 M], добавлен 08.10.2006

  • Разработка конструкторского расчета по техническому проектированию измерителя угловых скоростей на основе гексоды датчиков угловой скорости для космического корабля. Параметры троек неортогонально ориентированных ДУСов с электрическими обратными связями.

    дипломная работа [2,1 M], добавлен 23.01.2012

  • Патентно-аналитический обзор по датчикам измерения скорости, основания их классификации. Принцип действия и технические характеристики электромагнитных датчиков скорости. Использование эффекта Холла для конструирования датчика скорости автомобиля.

    курсовая работа [607,5 K], добавлен 13.01.2015

  • Реализация датчика угловой скорости вращения электродвигателя программным способом, анализируя количество опросов порта в течении периода импульсов, поступающих в заданный порт. оценка возможности уменьшения погрешности. Разработка и описание алгоритма.

    контрольная работа [70,2 K], добавлен 27.11.2012

  • Разработка электрической принципиальной схемы прибора. Описание ее элементов. Расчет усилителя, конденсатора для сглаживания пульсаций, напряжения на вторичной обмотке трансформатора. Выбор микросхемы стабилизатора напряжения и диодного выпрямителя.

    курсовая работа [1,9 M], добавлен 28.12.2014

  • Необходимость измерения скорости и направления кровотока. Доплеровские методы и аппараты. Доплеровские системы с двухмерной визуализацией. Разработка электрической принципиальной схемы и конструкции ультразвукового датчика прибора для измерения кровотока.

    дипломная работа [611,7 K], добавлен 07.05.2010

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.