Ѕаза знаний студента. –еферат, курсова€, контрольна€, диплом на заказ

курсовые,контрольные,дипломы,рефераты

»сследование движени€ центра масс межпланетных космических аппаратов —  осмонавтика

ѕосмотреть видео по теме ƒиплома

†TOC o "1-3" 1. ќглавление.................................................................................. 1

2. »сследовательска€ часть................................................... 3

2.1. ¬ведение.......................................................................................... 3

2.2.  раткие сведени€ об орбите........................................................... 4

2.2.1. ’арактеристика орбиты.............................................................. 4

2.2.2. —в€зь ћ ј с наземными пунктами управлени€........................ 5

2.2.3. ¬ыведение на рабочую орбиту................................................... 6

2.3. »сходные данные и цели работы................................................. 10

2.3.1. »сходные данные....................................................................... 10

2.3.2. ÷ели работы............................................................................... 12

2.4. ћоделирование движени€ центра масс ћ ј.............................. 13

2.4.1. ”равнени€ движени€ ћ ј........................................................ 13

2.4.2. ¬озмущающие ускорени€, действующие на ћ ј................... 15

2.4.3. –асчет параметров текущей орбиты ћ ј............................... 22

2.5.  оррекци€ траектории ћ ј....................................................... 24

2.5.1.  оррекци€ приведени€.............................................................. 24

2.5.2. –асчет потребного топлива....................................................... 26

2.5.3.  оррекци€ поддержани€........................................................... 27

2.6. ƒвижение ћ ј относительно центра масс................................. 28

2.6.1. ”равнени€ движени€ относительно центра масс ћ ј............ 28

2.6.2. —табилизаци€ углового положени€ при коррекции................ 28

3. ќ–√јЌ»«ј÷»ќЌЌќ-Ё ќЌќћ»„≈— јя „ј—“№................... 31

3.1. ќрганизаци€ и планирование выполнени€ темы........................ 31

3.2. ќпределение затрат труда............................................................ 31

3.3. –асчет сметы затрат на разработку программного продукта.... 35

4. ѕроћџЎЋ≈ЌЌјя Ё ќЋќ√»я » Ѕ≈«ќѕј—Ќќ—“№............... 39

4.1. ¬ведение........................................................................................ 39

4.2. јнализ вредных факторов........................................................... 39

4.3. “ребовани€ к видеотерминальным устройствам........................ 44

4.4. –асчет вредных излучений........................................................... 46

4.5. –ациональна€ организаци€ рабочего места................................ 46

4.6. –екомендации по снижению утомл€емости................................. 47

4.7. «ащита от напр€жени€ прикосновени€. «ануление..................... 48

4.8. ѕожарна€ безопасность................................................................ 49

5. —писок литературы............................................................... 53

6. ѕриложение. “екстџ ѕрограмм дл€ Borland C++ и Matlab 4.0 for windows........................................................................................... 54

6.1. ќсновной программный модуль main.cpp.................................. 54

6.2. ѕодпрограмма расчета возмущающих ускорений, параметров орбиты и коррекции sfun.cpp................................................................................................ 57

6.3. ‘айл начальной инициализации init.h......................................... 77

6.4. ‘айл описани€ переменных def.h................................................ 79

6.5. ‘айл sfun.h.................................................................................... 81

6.6. ‘айл rk5.h..................................................................................... 81

6.7. ѕрограмма построени€ временных диаграмм............................ 82


2. »——Ћ≈ƒќ¬ј“≈Ћ№— јя „ј—“№

2.1. ¬¬≈ƒ≈Ќ»≈

¬ данной работе проводитс€ исследование движени€ центра масс ћ ј под действием различных возмущающих ускорений (от не≠центральности гравитационного пол€ «емли, сопротивлени€ атмо≠сферы, прит€жени€ —олнца и Ћуны, из-за давлени€ солнечных лу≠чей) и создание математической модели движени€ ÷ћ ћ ј, по≠зво≠л€ющей учесть при интегрировании уравнений движени€ ÷ћ ћ ј эволюцию орбиты ћ ј.

¬ работе разрабатываетс€ алгоритм коррекции, ликвидирующий ошибки выведени€ ћ ј и рассчитываетс€ масса топлива, необхо≠дима€ дл€ проведени€ коррекции, необходимой из-за эволюции па≠раметров орбиты и из-за ошибок выведени€ ћ ј на рабочую ор≠биту.

“очность проведени€ коррекции зависит от точности направле≠ни€ корректирующего импульса, заданной в “«. Ѕыло проведено моде≠лирование системы коррекции в режиме стабилизации угло≠вого по≠ложени€ при работе корректирующей двигательной уста≠новки.

¬ работе привод€тс€ программы, реализующие интегрирование уравнений движени€ ÷ћ ћ ј, процесс осуществлени€ коррекции и расчет топлива дл€ коррекции.


2.2.  –ј“ »≈ —¬≈ƒ≈Ќ»я ќЅ ќ–Ѕ»“≈

ќсновными показател€ми эффективности космической группировки, €вл€ютс€:

- предельна€ производительность ћ ј в сутки на освещенной стороне «емли не менее 400-500 объектов.

- периодичность наблюдени€ районов съемки не реже одного раза в сутки.

–асположение плоскости орбиты по отношению к —олнцу выбрано таким образом, чтобы угол между линией узлов и следом терминатора на плоскости экватора «емли составл€л Dт = 30∞. ѕри† этом северный полувиток орбиты должен проходить над освещенной частью земной поверхности. ƒл€ определенности углу Dт приписываетс€ знак Ђ+ї в том случае, если восход€щий узел орбиты находитс€ над освещенной частью «емли, и знак Ђ-ї, если ¬” находитс€ над неосвещенной частью. ѕри выборе баллистического построени€ оперируют углом D, однозначно определ€ющимс€ пр€мым восхождением —олнца a0 и долготой восход€щего узла орбиты в абсолютном пространстве W: D = a0 - W. —оотношение между углом Dт и углом D: D º Dт - 90∞.

2.2.1. ’ј–ј “≈–»—“» ј ќ–Ѕ»“џ

ƒл€ решени€ задач наблюдени€ «емли из космоса с хорошим раз≠решением при жестких ограничени€х на массу  ј и минимизации затрат на выведение целесообразно использовать низкие круговые орбиты. ¬ этом классе орбит выдел€ют солнечно-синхронные ор≠биты со следующими свойствами:

- скорость прецессии плоскости орбиты в пространстве состав≠л€ет примерно 1∞ в сутки, что практически обеспечивает посто€н≠ство ориентации ее относительно терминатора «емли в течении всего срока активного существовани€  ј.

- близость наклонени€ плоскости орбиты к пол€рному, что обес≠печивает глобальность накрыти€ полюсами обзора поверхности «емли.

- возможность наблюдени€ районов на поверхности «емли при≠мерно в одно и то же местное врем€ при незначительном изменении углов места —олнца в точке наблюдени€.

¬сем этим услови€м удовлетвор€ют солнечно-синхронные ор≠биты с высотами от нескольких сот до полутора тыс€ч километров. Ќа больших высотах наклонение солнечно-синхронной орбиты от≠ли≠чаетс€ от пол€рного, и глобальность накрыти€ поверхности «емли не обеспечиваетс€. ƒл€ повышени€ эффективности наблюде≠ни€ це≠лесообразно выбрать орбиты с изомаршрутной трассой, у которых следы орбит ежесуточно проход€т на одними и теми же районами «емли, что позвол€ет обеспечивать периодичность на≠блюдени€ од≠ного и того же объекта, как минимум, раз в сутки с одного  ј.

ѕредварительные расчеты показали, что целесообразно исполь≠зовать орбиту с высотой Ќ = 574 км и наклонением плоскости ор≠биты к плоскости экватора «емли i = 97,6∞.

ћасса ћ ј может составить от 500 до 800 кг (что зависит от вида целевой аппаратуры, устанавливаемой на борту ћ ј). ƒл€ выведени€ ћ ј на орбиту используетс€ –Ќ ——-19 (Ђ–окотї) с разгонным блоком ЂЅризї.

2.2.3. —¬я«№ ћ ј — Ќј«≈ћЌџћ» ѕ”Ќ “јћ» ”ѕ–ј¬Ћ≈Ќ»я

”правление ћ ј осуществл€етс€ с наземных пунктов управлени€ на территории –оссии. »х количество и место расположени€ выбираетс€ таким образом, чтобы на любом витке можно было организовать сеанс св€зи с ћ ј хот€ бы с одного пункта управлени€. ”гол возывшени€ ћ ј над горизонтом наземного пункта управлени€ должен быть не менее 7∞, а дальность до ћ ј не должна превышать 2200 км.

¬ расчете зон св€зи были использованы следующие исходные данные:

- высота орбиты - 574 км.

- наклонение орбиты - 97,6∞.

- географическа€ долгота восход€щего узла первого витка - 4∞ в.д.

- минимальный угол возвышени€ ћ ј над местным горизонтом - 7∞.

»з рассматривавшихс€ возможных наземных пунктов управлени€ (ћосква, Ќовосибирск, ’абаровск, ћурманск,  алининград, ƒиксон,  омсомольск-на-јмуре, ѕетропавловск- амчатский), было выбрано три (ћосква, ƒиксон, ѕетропавловск- амчатский), обеспечивающие возможности св€зи с ћ ј на любом витке орбиты. ѕри этом зоны св€зи с ћ ј составл€ют от 3 до 9 минут на витке.

»нтергральные характеристики возможности св€зи с ћ ј:

- высота орбиты - 574 км.

- число витков, видимых из ћосквы, вит/сутки - 6.

- суммарное врем€ видимости из ћосквы, мин - 41.

- суммарное врем€ видимости с трех пунктов, мин - 153.

- максимальное врем€ видимости одного витка, мин - 9,1.

2.2.4. ¬џ¬≈ƒ≈Ќ»≈ ћ ј Ќј –јЅќ„”ё ќ–Ѕ»“”

¬ыведение ћ ј на орбиту с наклонением i = 97,6∞ и высотой Ќ = 574 км осуществл€етс€ ракетой-носителем Ђ–окотї с разгонным блоком ЂЅризї. ѕри выведении дл€ каждой отдел€ющейс€ части –Ќ (отработанна€ перва€ ступень, обтекатель, отработанна€ втора€ ступень) существует свой район падени€.

¬озможные варианты старта:

1. ѕолигон Ѕайконур.

»з-за отсутстви€ зон падени€ отдел€ющихс€ частей возможно сформировать опорную орбиту с наклонением i пор€дка 65∞. ƒл€ формировани€ опорной орбиты с наклонением близким пол€рному при использовании трассы с азимутом стрельбы более 180∞ (направление стрельбы на юг) - перва€ ступень падает в районе јшхабада, обтекатель сбрасываетс€ на высоте Ќ пор€дка 100 км, втора€ ступень падает за јравийским полуостровом. — точки зрени€ энергетики, выведение осуществл€етс€ не по оптимальной схеме, в результате чего на круговую орбиту высотой Ќ пор€дка 700 км выводитс€ ћ ј массой менее 600 кг.

2. ѕолигон Ћед€на€ (—вободный).

»з-за отсутстви€ зон падени€ отдел€ющихс€ частей возможно сформировать опорную орбиту с наклонением i пор€дка 54∞ и 65∞. ѕри северном запуске –Ќ перва€ ступень падает в районе заповедника в устье реки ќлейма (приток Ћены).

3.  осмодром ѕлесецк.

јзимуты пуска с космодрома ѕлесецк обеспечивают наклонени€ орбит i от 72∞ до 93∞. ‘ормирование требуемового наклонени€ i = 97,6∞ осуществл€етс€ с помощью разгонного блока ЂЅризї.

¬ результате работы двух ступеней –Ќ формируетс€ баллистическа€ траектори€ с наклонением i = 93∞. ¬ысота в момент окончани€ работы двигател€ второй ступени составл€ет Ќ = 190 км, наклонна€ дальность L = 300 км. ѕриблизительно через 1,2 секунды после прохождени€ команды на выключение двигател€ второй ступени проходит команда на запуск ƒ” –Ѕ. ѕосле выключени€ двигател€ второй ступени –Ќ происходит отделение от ракеты св€зки –Ѕ с  ј. ¬рем€ расцепки t = 318 секунд. јбсолютна€ скорость в момент отделени€ V = 5550 м/с. ќтдел€ема€ масса 6700 кг.

ƒвигательна€ установка –Ѕ ЂЅризї выполн€ет задачу доразгона  ј при формировании опорной орбиты.

’арактеристики двигател€ –Ѕ ЂЅризї:

- т€га R, кг - 2000.

- удельный импульс Rуд, сек - 324.

- количество включений, р - 7.

- интервал между включени€ми, сек - 20.

- врем€ функционировани€, час - 7.

¬ результате работы двигател€ –Ѕ ЂЅризї при первом включении происходит увеличение высоты баллистической траектории с Ќ = 190 км до Ќ = 270 км и к моменту окончани€ работы двигател€ (t = 905,5 сек) в точке с аргументом широты u = 104,1∞ формируетс€ опорна€ эллиптическа€ орбита с параметрами:

- высота в перигее Ќп, км - 190.

- высота в апогее Ќа, км - 574.

- больша€ полуось орбиты а, км - 6747.

- эксцентриситет e - 0,02548

- наклонение i, ∞ - 93,4.

- период обращени€ “, час - 1,53.

- аргумент периге€ w, ∞ - 128,38.

- долгота восход€щего узла в гринвичской — , фиксированной на момент старта Wг, ∞ - 48,37.

¬еличина импульса характеристической скорости, отрабатываемого при первом включении ƒ” –Ѕ dV1 = 2,36 км/с, врем€ работы пор€дка 600 сек.

–абота двигател€ при первом включении происходит вне зоны видимости Ќѕ” на территории –оссии. √еографические координаты, соответствующие этому моменту:

- широта j ї 76∞.

- долгота l ї 238∞.

¬ момент прохождени€ ћ ј периге€ опорной эллиптической орбиты (t = 1231 сек) географические координаты составл€ют:

- широта j ї 53∞.

- долгота l ї 227∞.

Ќа опорной эллиптической орбите ћ ј совершает пассивный полет до апоге€. ¬ районе апоге€ (t = 1,12 час) осуществл€етс€ второе включение ƒ” –Ѕ.

¬ результате приложени€ второго компланарного импульса характеристической скорости dV2 = 0,12 км/с, при втором включении (врем€ работы 20 сек) формируетс€ кругова€ орбита с параметрами:

- высота Ќ, км - 574.

- наклонение i, ∞ - 93,4.

- период обращени€ “, час - 1,6.

–абота двигател€ при втором включении происходит вне зоны видимости Ќѕ” на территории –оссии. √еографические координаты, соответствующие этому моменту:

- широта j ї 1,5∞.

- долгота l ї 35,8∞.

ƒл€ создани€ круговой, солнечно-синхронной орбиты необходимо изменить наклонение до i = 97,6∞. — этой целью осуществл€етс€ третье включение ƒ” –Ѕ при первом прохождении восход€щего узла орбиты (t = 1,3 час).

¬ результате приложени€ ортогонального импульса характеристической скорости dV3 = 0,62 км/с, при третьем включении (врем€ работы 90 сек) формируетс€ солнечно-синхронна€ кругова€ орбита с параметрами:

- высота Ќ, км - 574.

- наклонение i, ∞ - 97,6.

- период обращени€ “, час - 1,6.

- число оборотов в сутки N - 15.

–абота двигател€ при третьем включении происходит вне зоны видимости Ќѕ” на территории –оссии. √еографические координаты, соответствующие этому моменту:

- широта j ї 0∞.

- долгота l ї 28,1∞.

ѕосле выключени€ двигател€ при третьем запуске происходит отделение ћ ј от –Ѕ ЂЅризї.

 инематические параметры в гринвичской — , фиксированной на момент старта –Ќ и оскулирующие элементы орбиты на момент отделени€ от –Ѕ:

ѕараметр

«начение

t, сек

4946,5

X, м

4638800

Y, м

5120280

Z, м

689680

Vx, м/с

241,23

Vy, м/с

-1233

Vz, м/с

7473,5

l, ∞

28,1

T, c

5761,67

e

0,0009

i, ∞

97,595

Ra, м

6940000

Rп, м

6952000


2.3. »—’ќƒЌџ≈ ƒјЌЌџ≈ » ÷≈Ћ» –јЅќ“џ

2.3.1. »—’ќƒЌџ≈ ƒјЌЌџ≈

Ќоминальна€ орбита, необходима€ дл€ выполнени€ задач ћ ј, имеет следующие параметры:

- кругова€, e = 0.

- солнечно-синхронна€, скорость прецессии линии узлов орбиты W равна скорости обращени€ —олнца относительно «емли

W = 2p / 365,2422 = 0,0172 рад/сут = 0,98 ∞/сут.

- изомаршрутна€, за сутки ћ ј совершает целое количество оборотов (n = 15).

Ёто обеспечивает прохождение ћ ј над одними и теми же рай≠онами в одно и тоже местное врем€.

- период “ = 5765 с.

- высота орбиты Ќ = 574 км.

- наклонение орбиты i = 97,6∞.

- географическа€ долгота восход€щего узла орбиты lэ = 28,1∞.

ƒолгота восход€щего узла в геоцентрической экваториальной (абсолютной) системе координат OXYZ определ€етс€ как разность

lэ - s0,

где s0 - часовой угол, отсчитывающийс€ от гринвичского мери≠диана до оси X, направленной в точку весеннего равноденстви€.

„асовой угол зависит от даты старта и выбираетс€ из астроно≠ми≠ческого ежегодника. ¬ данной задаче дл€ моделировани€ вы≠бран часовой угол = 0.

—ледовательно долгота восход€щего узла орбиты W = lэ = 28,1∞.

»сход€ из “«, начальна€ точка выведени€ имеет следующие ко≠ор≠динаты в гринвичской системе координат, фиксированной на момент старта –Ќ:

ѕараметр

«начение

t, сек

4946.5

X, м

4638800

Y, м

5120280

Z, м

689506,95

Vx, м/с

241,23

Vy, м/с

-1233

Vz, м/с

7472,65

Ёлементы орбиты:

l, ∞

28,1

T, c

5761,67

e

0,0009

i, ∞

97,595

Ra, м

6940000

Rп, м

6952000

 инематические параметры в геоцентрической экваториальной системе координат:

t, сек

4946.5

X, м

6137262,9

Y, м

3171846,1

Z, м

689506,95

Vx, м/с

-201,3

Vy, м/с

-1247,03

Vz, м/с

7472,65

l, ∞

28,1

“очность выведени€:

- предельна€ ошибка по координате (3s) - 7 км.

- предельна€ ошибка по скорости (3s) - 5 м/с.

ѕересчитав ошибку по координате на ошибку по периоду выве≠дени€ орбиты получим предельную ошибку по периоду DT - 10 сек.

 оррел€ционна€ матрица ошибок выведени€ на момент выведе≠ни€ составл€ет:

„лены, сто€щие на главной диагонали представл€ют собой квад≠раты предельных ошибок - (3s)2.

K11 = K22 = K33 = (3s)2 = 72 = 49 км.

K44 = K55 = K66 = (3s)2 = 52 = 25 м/с.

ќстальные члены представл€ют собой вторые смешанные мо≠менты Kij = Kji = rijsisj или Kij = Kji = rjj(3si)(3sj), где rjj - коэффици≠енты св€зи величин i и j. ¬ данном случае вторые смешанные мо≠менты Kij = Kji = 0.

 инематические параметры в геоцентрической экваториальной системе координат на момент выведени€ с учетом ошибок выведе≠ни€:

t, сек

4946.5

X, м

6144262,9

Y, м

3178846,1

Z, м

696506,95

Vx, м/с

-206,3

Vy, м/с

-1252,03

Vz, м/с

7477,65

l, ∞

28,1

ѕараметры орбиты с учетом ошибок выведени€:

l, ∞

28,13

T, c

5795,7

W, ∞

28,13

p, км

6973,5

а, км

6973,6

e

0,00314

i, ∞

97,637

2.3.2. ÷≈Ћ» –јЅќ“џ

1) »сследование и моделирование движени€ ÷ћ ћ ј при воз≠действии на  ј возмущающих ускорений.

2) –азработка алгоритмов проведени€ коррекции траектории ћ ј, моделировани€ процесса, и расчет потребного топлива дл€ проведени€ коррекции траектории.

3) »сследование динамики системы коррекции траектории при стабилизации углового положени€ в процессе проведени€ коррек≠ции траектории ћ ј.


2.4. ћќƒ≈Ћ»–ќ¬јЌ»≈ ƒ¬»∆≈Ќ»я ÷≈Ќ“–ј ћј—— ћ ј

2.4.1.”–ј¬Ќ≈Ќ»≈ ƒ¬»∆≈Ќ»я  ј

–ассмотрим невозмущенное движение материальных точек ћ и m в некоторой инерциальной системе координат. ƒвижение совер≠ша≠етс€ под действием силы прит€жени€ Fz. —ила Fz дл€ материаль≠ной точки m определ€етс€ формулой:

где ¶ - посто€нна€ прит€жени€,

ro- единичный вектор, направленный от ћ к m,

где - радиус-вектор, проведенный из т.ћ до т.m.

r - относительное рассто€ние от ћ до m.

Ќа точку ћ действует сила Fz, равна€ по величине и направлен≠на€ в противоположную сторону.

Ќа основе второго закона Ќьютона уравнени€ движени€ матери≠альных точек ћ и m имеют вид:

(1), †(2)

или

(3), †(4)

где p1 - радиус-вектор, проведенный из начала инерциальной сис≠темы координат в точку m.

p2 - радиус-вектор, проведенный из начала инерциальной сис≠темы координат в точку ћ.

¬ычита€ из уравнени€ (3) уравнение (4), получим уравнение дви≠жени€ мате≠риальной точки m относительно прит€гивающего цен≠тра ћ:

“ак как m<<ћ, следовательно, можно пренебречь ускорением, которое  ј с массой m сообщает прит€гивающему центру ћ. “о≠гда можно совместить начало инерциальной системы координат с при≠т€гивающим центром ћ. —ледовательно,

“аким образом, уравнение невозмущенного движени€  ј отно≠сительно прит€гивающего центра ћ в инерциальной системе коор≠динат, центр которой находитс€ в ћ, имеет вид

где m = fM - гравитационна€ посто€нна€ «емли.

–ассмотрим возмущенное движение  ј в геоцентрической† эква≠ториальной (абсолютной) системе координат OXYZ:

- начало ќ - в центре масс «емли.

- ось X направлена в точку весеннего равноденстви€ g.

- ось Z совпадает с осью вращени€ «емли и направлена на —евер≠ный полюс «емли.

- ось Y дополн€ет систему до правой.

ƒвижение  ј в абсолютной системе координат OXYZ происхо≠дит под действием центральной силы прит€жени€ «емли Fz, а также под действием возмущающих сил Fв. ”равнение движени€ имеет вид

†или

где m = 597 кг - масса  ј.

¬ проекци€х на оси абсолютной системы координат OXYZ полу≠чим

†или

†или

†или

где axв, ayв, azв - возмущающиес€ ускорени€.

ќсновные возмущающиес€ ускорени€ вызываютс€ следующими причинами:

- нецентральностью пол€ прит€жени€ «емли.

- сопротивлением атмосферы «емли.

- вли€нием —олнца.

- вли€нием Ћуны.

- давлением солнечного света.

2.4.2. ¬ќ«ћ”ўјёў»≈ ”— ќ–≈Ќ»я, ƒ≈…—“¬”ёў»≈ Ќј ћ ј

1) ¬озмущающеес€ ускорение, вызванное нецентральностью гра≠витационного пол€ «емли.

–ассмотрим потенциал пол€ прит€жени€ «емли. ѕри точном рас≠чете параметров орбиты спутников, в качестве хорошего прибли≠же≠ни€ к действительной поверхности «емли принимают геоид. √еоид - это гипотетическа€ уровенна€ поверхность, совпадающа€ с поверх≠ностью спокойного океана и продолженна€ под материком.

»ногда в баллистике под геоидом понимают не поверхность, а тело, которое ограничено поверхностью мирового океана при не≠ко≠тором среднем уровне воды, свободной от возмущений. ¬о всех точках геоида потенциал прит€жени€ имеет одно и то же значение.

ѕотенциал прит€жени€ «емли можно представить в виде разло≠же≠ни€ по сферическим функци€м.

где mz = fMz - гравитационна€ посто€нна€ «емли.

r0 - средний экваториальный радиус «емли.

сnm, dnm - коэффициенты, определ€емые из гравиметрических дан≠ных, а также по наблюдени€м за движением »—«.

L - долгота прит€гивающей точки.

j - широта прит€гивающей точки.

Pnm(sinj) - присоединенные функции Ћежандра степени m и по≠р€дка n (при m ¹ 0).

Pnm(sinj) - многочлен Ћежандра пор€дка n (при m = 0).

—оставл€ющие типа (mz/r)(r0/r)ncn0Pn0(sinj) - называют зональ≠ными гармониками n-по≠р€дка. “.к. полином Ћежандра n-го по≠р€дка имеет n действительных корней, функци€ Pn0(sinj) будет ме≠н€ть знак на n широтах, сфера делитс€ на n+1 широтную зону, где эти составл€ю≠щие имеют попеременно Ђ+ї или Ђ-ї значени€. ѕо≠этому их называют зональными гармониками.

—оставл€ющие типа

(mz/r)(r0/r)ncnmcos(mL)Pnm(sinj) и (mz/r)(r0/r)ndnmsin(mL)Pnm(sinj)

- называют тессеральными гармониками n-пор€дка и степени m. ќни обращаютс€ в 0 на 2m меридианах, где cos(mL) = 0 и sin(mL) = 0 и на n-m параллел€х, где Pnm(sinj) = 0 или dmPnm(sinj)/d(sinj)m = 0, сфера делитс€ на n+m+1 трапецию, где эти составл€ющие сохра≠н€ют знак.

—оставл€ющие типа и

(mz/r)(r0/r)ncnncos(nL)Pnn(sinj) и (mz/r)(r0/r)ndnnsin(nL)Pnn(sinj)

- называют секториальными гармониками n-пор€дка и степени m. Ёти составл€ющие мен€ю знак только на меридианах, cos(nL) = 0 и sin(nL) = 0, на сфере выдел€ют 2n меридиональных секторов, где эти составл€ющие со≠хран€ют знак.

ћногочлен Ћежандра степени n находитс€ по следующей фор≠муле:

Pn0(z) = 1/(2nn!)´(dn(z2 - 1)n/dzn)

ѕрисоединенна€ функци€ Ћежандра пор€дка n и степени m нахо≠дитс€ по следующей формуле:

Pnm(z) = (1-z2)m/2´dmPn0(z)/dzm

¬озмущающа€ часть гравитационного потенциала «емли равна

Uв = UТ + DUТ = (U - mz/r) + DUТ

где DUТ - потенциал аномалий силы т€готени€ «емли.

UТ - часть потенциала «емли, котора€ учитывает несферичность «емли.

—ледовательно,

ѕерва€ зональна€ гармоника в разложении потенциала учиты≠вает пол€рное сжатие «емли.

«ональные гармоники нечетного пор€дка и тессеральные гармо≠ники, где n-m нечетное число - учитывают ассиметрию «емли отно≠сительно† плоскости экватора.

—екториальные и тессеральные гармоники - учитывают ассимет≠рию «емли относительно оси вращени€.

ѕерва€ зональна€ гармоника имеет пор€док 10-3, а все остальные - пор€док 10-6 и выше. ѕоэтому будем учитывать в разложении по≠тен≠циала прит€жени€ только зональную гармонику (n=2, m=0) и секторальную гармонику (n=2, m=2). “акже не будем учитывать потенциал аномалий силы т€го≠тени€ «емли DUТ.

“аким образом,

Uв = (mz/r)(r0/r)2[c20P20(sinj) + (c22cos(2L) + d22sin(2L))P22(sinj)],

где c20 = - 0,00109808,

c22 = 0,00000574,

d22 = - 0,00000158.

P20(x) = 1/222!´d2(x2 - 1)2/dx2.

—ледовательно P20(x) = (3x2 - 1)/2.

“ак как sinj = z/r, следовательно P20(sinj) = (3(z/r)2 - 1)/2.

P22(x) = (1 - x2)2/2´d2P20(x)/dx2 = 1/2´(1 - x2)´d2(3x2 - 1)/dx2

—ледовательно P22(x) = 3(1 - x2).

“ак как sinj = z/r, следовательно P22(sinj) = 3(1 - (z/r)2).

«начит

Uв по координатам X, Y, Z, причем r = Ö(x2 + y2 + z2).

—ледовательно,

2) ¬озмущающее ускорение, вызванное сопротивлением атмо≠сферы.

ѕри движении в атмосфере на  ј действует сила аэродинамиче≠ского ускорени€ Rx, направленна€ против вектора скорости  ј от≠но≠сительно† атмосферы:

где Cx = 2 - коэффициент аэродинамического сопротивлени€.

Sм = 2,5 м2 - площадь миделевого сечени€ - проекци€  ј на плос≠кость, пер≠пендикул€рную направлению скорости полета.

V - скорость  ј.

r - плотность атмосферы в рассматриваемой точке орбиты.

“ак как исследуема€ орбита - кругова€ с высотой Ќ = 574 км, бу≠дем считать, что плотность атмосферы одинакова во всех точках ор≠биты и равна плотности атмосферы на высоте 574 км. »з таб≠лицы стандартной атмосферы находим плотность наиболее близ≠кую к вы≠соте Ќ = 574 км. ƒл€ высоты Ќ = 580 км r = 5,098´10-13 кг/м3.

—ила аэродинамического ускорени€ создает возмущающее каса≠тельное ускорение aa:

Ќайдем проекции аэродинамического ускорени€ на оси абсолют≠ной системы координат axa, aya, aza:

aa направлено против скорости  ј, следовательно единичный век≠тор направлени€ имеет вид

ea = [Vx/|V|, Vy|V|, Vz/|V|], |V| = Ö(Vx2+Vy2 +Vz2)

“аким образом,

«начит

, ,

3) ¬озмущающее ускорение, вызванное давлением солнечного света.

ƒавление солнечного света учитываетс€ как добавок к посто€н≠ной т€готени€ —олнца - Dmc. Ёта величина вычисл€етс€ следующим об≠разом:

Dmc = pSмA2/m

где p = 4,64´10-6 Ќ/м2 - давление солнечного света на рассто€нии в одну астрономи≠ческую единицу ј.

A = 1,496´1011 м - 1 астрономическа€ единица.

m - масса  ј.

Sм = 8 м2 - площадь миделевого сечени€ - проекци€  ј на плос≠кость, пер≠пендикул€рную направлени€ солнечных лучей.

“аким образом,

Dmc = 1,39154´1015 м3/c2.

4) ¬озмущающее ускорение, возникающее из-за вли€ни€ —олнца.

”равнение движени€  ј в абсолютной системе координат OXYZ относительно «емли при воздействии —олнца:

где mz - посто€нна€ т€готени€ «емли.

mc - посто€нна€ т€готени€ —олнца.

r - радиус-вектор от «емли до  ј.

rc - радиус-вектор от «емли до —олнца.

“аким образом, возмущающее ускорение, возникающее из-за вли€ни€ —олнца:

.

«десь первое слагаемое есть ускорение, которое полу≠чил бы  ј, если он был неприт€гиваю≠щим, а «емл€ отсутствовала.

¬торое слагаемое есть ускорение, которое сообщает —олнце «емле, как неприт€гиваю≠щему телу.

—ледовательно, возмущающее ускорение, которое получает  ј при движении относительно «емли - это разность двух слагаемых.

“ак как rc>>r, то в первом слагаемом можно пренебречь r. —ледо≠ва≠тельно

| rc - r| = Ö((xc-x)2+(yc-y)2+(zc-z)2)

где xc, yc, zc - проекции радиуса-вектора —олнца на оси абсолют≠ной сис≠темы координат.

ћоделирование движени€ —олнца проводилось следующим об≠ра≠зом: за некоторый промежуток времени t —олнце относительно «емли сместитс€ на угол J = Jн + wct,

где Jн = W + (90 - D) - начальное положение —олнца в эклиптиче≠ской системе коор≠динат.

W = 28,1∞ - долгота восход€щего узла первого витка  ј.

D = 30∞ - угол между восход€щим узлом орбиты  ј и терминато≠ром.

wc - углова€ скорость —олнца относительно «емли.

wc = 2p/T = 2p/365,2422´24´3600 = 1,991´10-7 рад/c = 1,14´10-5 ∞/c

“аким образом, в эклиптической системе координат проекции† составл€ют:

xce = rccosJ

yce = rcsinJ

zce = 0

rc = 1,496´1011 м (1 астрономическа€ единица) - рассто€ние от «емли до —олнца

ѕлоскость эклиптики наклонена к плоскости экватора на угол e = 23,45∞, проекции rc на оси абсолютной системы координат можно найти как

xc = xce = rccosJ

yce = ycecose = rccosJcose

zce = rcsinJsine

“аким образом, проекции возмущающего ускорени€ на оси абсо≠лютной системы координат:

axc = - mcx/(Ö((xc-x)2+(yc-y)2+(zc-z)2))3

ayc = - mcy/(Ö((xc-x)2+(yc-y)2+(zc-z)2))3

azc = - mcz/(Ö((xc-x)2+(yc-y)2+(zc-z)2))3

— учетом солнечного давлени€

axc = - (mc-Dmc)x/(Ö((xc-x)2+(yc-y)2+(zc-z)2))3

ayc = - (mc-Dmc)y/(Ö((xc-x)2+(yc-y)2+(zc-z)2))3

azc = - (mc-Dmc)z/(Ö((xc-x)2+(yc-y)2+(zc-z)2))3

5) ¬озмущающее ускорение, возникающее из-за вли€ни€ Ћуны.

”равнение движени€  ј в абсолютной системе координат OXYZ относительно «емли при воздействии Ћуны:

где mл = 4,902´106 м3/c2- посто€нна€ т€готени€ Ћуны.

rл - радиус-вектор от «емли до Ћуны.

“аким образом, возмущающее ускорение, возникающее из-за вли€ни€ Ћуны:

“ак как rл>>r, то в первом слагаемом можно пренебречь r. —ледо≠ва≠тельно

|rл - r| = Ö((xл-x)2+(yл-y)2+(zл-z)2)

где xл, yл, zл - проекции радиуса-вектора Ћуны на оси абсолютной системы координат.

ƒвижение Ћуны учитываетс€ следующим образом: положение Ћуны в каждый момент времени рассчитываетс€ в соответствии с данными астрономического ежегодника. ¬се данные занос€тс€ в массив, и далее этот массив считаетс€ программой моделировани€ движени€  ј. ¬ первом приближении принимаетс€:

- орбита Ћуны - кругова€.

- угол наклона плоскости орбиты Ћуны к плоскости эклиптики i = 5,15∞.

- период обращени€ линии пересечени€ плоскостей лунной ор≠биты и эклиптики (по ходу часовой стрелки, если смотреть с север≠ного полюса) = 18,6 года.

”гол между плоскост€ми экватора «емли и орбиты Ћуны можно найти по формуле

cos(hл) = cos(e)cos(i) - sin(e)sin(i)cos(Wл)

где Wл - долгота восход€щего узла лунной орбиты, отсчитыва≠етс€ от направлени€ на точку весеннего равноденстви€.

e - угол между плоскост€ми эклиптики и экватора «емли.

¬еличина hл колеблетс€ с периодом 18,6 лет между минимумом при hл = e - i = 18∞18Т и максимумом при hл = e + i = 28∞36Т при W = 0.

ƒолгота восход€щего узла лунной орбиты Wл измен€етс€ с тече≠нием времени t на величину Wл = t´360/18,6´365,2422´24´3600.

ѕоложение Ћуны на орбите во врем€ t определ€етс€ углом

J л = t´360/27,32´24´3600.

ѕо формулам перехода найдем проекции вектора положени€ Ћуны на оси абсолютной системы координат:

xл = rл(cosJлcosWл - coshлsinJлsinWл)

yл = rл(cosJлsinWл + coshлsinJлcosWл)

zл = rлsinhлsinJл

rл = 3,844´108 м - среднее рассто€ние от «емли до Ћуны

“аким образом, проекции возмущающего ускорени€ на оси абсо≠лютной системы координат:

axл = - mлx/(Ö((xл!-x)2+(yл-y)2+(zл-z)2))3

ayл = - mлy/(Ö((xл!-x)2+(yл-y)2+(zл-z)2))3

azл = - mлz/(Ö((xл!-x)2+(yл-y)2+(zл-z)2))3

”равнени€ возмущенного движени€ при действии корректирую≠щего ускорени€ имеют вид:

или

d2x/dt2 = - (mz/r2)x + axu + axa + axc + axл + axк

d2y/dt2 = - (mz/r2)y + ayu + aya + ayc + ayл + ayк

d2z/dt2 = - (mz/r2)z + azu + aza + azc + azл + azк

2.4.3. –ј—„≈“ ѕј–јћ≈“–ќ¬ “≈ ”ў≈… ќ–Ѕ»“џ  ј

ѕолученна€ система уравнений движени€ ÷ћ  ј интегрируетс€ методом –унге- утта 5-го пор€дка с переменным шагом. Ќачаль≠ные услови€ x0, y0, z0, Vx0, Vy0, Vz0 - в абсолютной системе коорди≠нат, соответствуют началь≠ной точке вывода при учете ошибок вы≠ведени€. ѕосле интегриро≠вани€ мы получаем вектор состо€ни€  ј (x, y, z, Vx, Vy, Vz) в любой момент вре≠мени.

ѕо вектору состо€ни€ можно рассчитать параметры орбиты. со≠ответствующие этому вектору состо€ни€.

а) ‘окальный параметр - р.

р = C2/mz, где — - интеграл площадей.

C = r ´ V, |C| = C = Ö(Cx2+Cy2+Cz2)

Cx = yVz - zVy

Cy = zVx - xVz† - проекции на оси абсолютной — 

Cz = xVy - yVx

б) Ёксцентриситет - е.

e = f/mz, где f - вектор Ћапласа

f = V ´ C - mzr/r, |f| = f = Ö(fx2+fy2+fz2)

fx = VyCz - VzCy - mzx/r

fy = VzCx - VxCz - mzy/r - проекции на оси абсолютной — 

fz = VxCy - VyCx - mzz/r

в) Ѕольша€ полуось орбиты.

a = p/(1 - e2)

г) Ќаклонение орбиты - i.

Cx = Csin(i)sinW

Cy = - Csin(i)cosW

Cz = Ccos(i)

можно найти наклонение i = arccos(Cz/C)

д) ƒолгота восход€щего узла - W.

»з предыдущей системы можно найти

sinW = Cx/Csin(i)

cosW = - Cy/Csin(i)

“ак как наклонение орбиты измен€етс€ несильно в районе i = 97,6∞, мы имеем право делить на sin(i).

≈сли sinW => 0, W = arccos (-Cy/Csin(i))

≈сли sinW < 0, W = 360 - arccos (-Cy/Csin(i))

е) јргумент перицентра - w.

fx = f(coswcosW - sinwsinWcos(i))

fy = f(coswsinW + sinwcosWcos(i))

fz = fsinwsin(i)

ќтсюда найдем

cosw = fxcosW/f + fysinW/f

sinw = fz/fsin(i)

≈сли sinw > 0, w = arccos (fxcosW/f + fysinW/f)

≈сли sinw < 0, w = 360 - arccos (fxcosW/f + fysinW/f)

ж) ѕериод обращени€ - “.

T = 2pÖ(a3/mz)

√рафики изменени€ элементов орбиты при действии всех, рас≠смотренных выше, возмущающих ускорений в течение 2-х перио≠дов (“ = 5765 с) приведены на рис. 1-12.

√рафики изменени€ во времени возмущающих ускорений приве≠дены на рис. 13-18.


2.5. ѕ–ќ¬≈ƒ≈Ќ»≈  ќ––≈ ÷»» “–ј≈ “ќ–»» ћ ј

—уществующие ограничени€ на точки старта –Ќ и зоны падени€ отработавших ступеней –Ќ, а также ошибки выведени€ не позво≠л€ют сразу же после пуска реализовать рабочую орбиту.  роме того, эволюци€ параметров орбит под действием возмущающих ус≠корений в процессе полета ћ ј приводит к отклонению парамет≠ров орбиты  ј от требуемых значений. ƒл€ компенсации воздей≠ст≠ви€ указанных факторов осуществл€етс€ коррекци€ орбиты с по≠мощью корректирующей двигательной установки ( ƒ”), котора€ располагаетс€ на борту ћ ј.

¬ данной работе проведена разработка алгоритма коррекции, моделирование процесса коррекции и расчет топлива, необходи≠мого дл€ проведени€ коррекции.

»з-за различных причин возникновени€ отклонений элементов орбиты проводитс€:

- коррекци€ приведени€ - ликвидаци€ ошибок выведени€ и при≠ве≠дение фактической орбиты к номинальной с заданной точно≠стью.

- коррекци€ поддержани€ - ликвидаци€ отклонений параметров орбиты от номинальных, возникающих из-за действи€ возмущаю≠щих ускорений в процессе полета.

ƒл€ того, чтобы орбита отвечала заданным требовани€м, откло≠нени€ параметров задаютс€ следующим образом:

- максимальное отклонение наклонени€ орбиты Di = 0,1∞

- предельное суточное смещение  ј по долготе† Dl = 0,1∞

—ледовательно, максимальное отклонение периода орбиты DT = 1,6 сек.

јлгоритм коррекции следующий:

1)  оррекци€ приведени€.

2)  оррекци€ поддержани€.

2.5.1.  ќ––≈ ÷»я ѕ–»¬≈ƒ≈Ќ»я

ѕосле окончани€ процесса выведени€ ћ ј, провод€тс€ внешне-траекторные измерени€ (¬“»). Ёти измерени€ обеспечивают, по баллистическим расчетам, знание вектора состо€ни€ с требуемой точностью через 2 суток. ѕосле этого начинаетс€ коррекци€ приве≠дени€.

ѕредложена следующа€ схема проведени€ коррекции:

а)  оррекци€ периода.

б)  оррекци€ наклонени€.

 орректирующий импульс прикладываетс€ в апсидальных точ≠ках, либо на линии узлов в течение 20 сек и происходит исправле≠ние одного параметра орбиты. “аким образом используетс€ одно≠пара≠метрическа€, непрерывна€ коррекци€.

а)  оррекци€ периода.

ќсуществл€етс€ в два этапа:

- коррекци€ перицентра

- коррекци€ апоцентра

—начала осуществл€етс€ коррекци€ перицентра - приведение те≠кущего рассто€ни€ до перицентра rp к номинальному радиусу rн = 6952137 м. ѕо≠сле измерени€ вектора состо€ни€ рассчитываютс€ параметры ор≠биты. ƒалее определ€етс€ нужный корректирующий импульс† DVк. Ќа≠правление импульса (тормоз€щий или разгон€ю≠щий) зависит от взаимного расположени€ перицентра орбиты и радиуса номиналь≠ной орбиты. ƒл€ этого вычисл€етс€ Drp = rp - rн.

¬озможны ситуации:

1)

Drp< 0 - прикладываетс€ разгон€ющий импульс

2) Drp> 0 - прикладываетс€ тормоз€щий импульс

 ј долетает до апоцентра и в апоцентре прикладываетс€ кор≠рек≠тирующий импульс. ¬рем€ работы  ƒ” - 20 сек.

“ак как врем€ работы  ƒ” ограничено, а DVк может быть боль≠шим, следовательно, далее рассчитываетс€ максимальный импульс скоро≠сти DVmax за 20 сек работы двигател€:

DVmax = Pt/m = 25´20/597 = 0,8375 м/с

≈сли DVк > DVmax в апоцентре прикладываетс€ импульс DVк = DVmax. ¬ результате этого rpнемного корректируетс€. Ќа следую≠щем витке оп€ть рассчитыва≠етс€ DVк, и если на этот раз DVк < DVmax, в апоцентре прикладываетс€ импульс DVк.  ƒ” включаетс€ не на полную мощность P = (DVк/DVmax)Pmax.

¬рем€ включени€ = 20 сек.

Ёто происходит до тех пор, пока не приблизитс€ к rpс заданной точностью.

ѕосле того, как скорректирован перицентр, начинаетс€ коррек≠ци€ апоцентра. –ассчитываютс€ параметры орбиты и нужный кор≠ректирующий импульс, такой, чтобы ra = rн = 6952137 м. Ќаправле≠ние корректи≠рующего импульса также зависит от величин ra и rн.

¬ычисл€етс€† Dra = ra - rн.

¬озможна ситуаци€:

Dra> 0 - в перицентре прикладываетс€ тормоз€щий импульс.

 ј долетает до перицентра и в перицентре прикладываетс€† кор≠ректирующий импульс. ¬рем€ работы  ƒ” - 20 сек.

“ак как врем€ работы  ƒ” ограничено, а DVк может быть боль≠шим, следовательно, далее рассчитываетс€ максимальный импульс скоро≠сти DVmax за 20 сек работы двигател€:

DVmax = Pt/m = 25´20/597 = 0,8375 м/с

≈сли DVк > DVmax, в перицентре прикладываетс€ импульс DVк = DVmax. ¬ результате этого немного корректируетс€ ra. Ќа следую≠щем витке оп€ть рассчитыва≠етс€ DVк, и если на этот раз DVк < DVmax, в перицентре прикладываетс€ импульс DVк.  ƒ” включаетс€ не на полную мощность P = (DVк/DVmax)Pmax.

¬рем€ включени€ = 20 сек.

Ёто происходит до тех пор, пока ra не приблизитс€ к rн с задан≠ной точностью.

“аким образом осуществл€етс€ коррекци€ перехода.

б)  оррекци€ наклонени€.

ѕосле коррекции периода провод€тс€ внешне-траекторные изме≠рени€ и получают вектор состо€ни€  ј. ≈сли снова необходима коррекци€ периода ее провод€т еще раз и снова измер€ют вектор состо€ни€  ј.

ƒалее проводитс€ коррекци€ наклонени€ по такой же схеме.  ор≠рекци€ производитс€ в точке пересечени€ орбиты  ј с линией уз≠лов.

ѕосле того, как рассчитаны корректирующие импульсы скоро≠сти, по формулам перехода проекции вектора на оси абсолютной сис≠темы координат. ƒалее рассчитываетс€ корректирующее уско≠рение и подставл€етс€ в уравнени€ движени€ центра масс  ј. ѕо≠сле этого уравнени€ интегрируютс€ методом –унге- утта 5-го по≠р€дка с пе≠ременным шагом.

√рафики изменени€ элементов орбиты в процессе коррекции при≠ведени€ приведены на рис.19-30.

2.5.2. –ј—„≈“ ѕќ“–≈ЅЌќ√ќ “ќѕЋ»¬ј

ћасса топлива, необходимого дл€ проведени€ коррекции траек≠тории рассчитываетс€ по формуле ÷иолковского:

m = m0(1 - e-DVк/W)

m0 = 597 кг - начальна€ масса ћ ј (кг)

W = 2200 м/с - скорость истечени€ газов из сопла двигател€.

–езультаты проведени€ коррекции приведени€:

tн, с

tк, с

Dt, с

DVк, м/c

»мп.

m, кг

 оррекци€ периода

176242

262592

300

12,1

15

3,26

 оррекци€

наклонени€

273984

432298

580

24,11

29

6,48

2.5.3. ќ––≈ ÷»я ѕќƒƒ≈–∆јЌ»я

ќсновна€ задача ћ ј - проведение съемки определенных рай≠онов «емли по крайней мере один раз в сутки, т.е. трасса  ј должна проходить над заданным районом каждые сутки.

“ребовани€ дл€ проведени€ коррекции:

- предельное суточное смещение орбиты по долготе †Di = 0,1∞

- предельное отклонение наклонен舆 Dl = 0,1∞.

¬ пересчете отклонени€ †Dl на отклонение по периоду получим:

DT = 1,597 сек. - максимальное отклонение по периоду.

ѕри помощи программы моделировани€ было просчитано 3 ме≠с€ца и получено, что средний период изменилс€ на 3,2 сек, а накло≠нение - на 0,001∞.

“аким образом, коррекцию периода надо делать примерно 1 раз в 1,5 мес.

Ќужный импульс скорости - 1 м/с за врем€ активного существо≠вани€ - 5 лет - коррекцию периода надо провести 40 раз, DV = 40 м/с, масса топ≠лива = 10,8 кг.

«а 5 лет Di = 0,02∞ - коррекцию наклонени€ проводить не надо.

√рафики изменени€ элементов орбиты за 3 мес€ца приведены на рис.31-42.


2.6. ƒ¬»∆≈Ќ»≈ ћ ј ќ“Ќќ—»“≈Ћ№Ќќ ÷≈Ќ“–ј ћј——

2.6.1. ”–ј¬Ќ≈Ќ»я ƒ¬»∆≈Ќ»я ќ“Ќќ—»“≈Ћ№Ќќ ÷ћ  ј

ѕри рассмотрении движени€ относительно ÷ћ  ј используют уравнени€ Ёйлера:

Jxwx + (Jz-Jy)wywz = Mxy + Mxв

Jywy + (Jx-Jz)wxwz = Myy + Myв

Jzwz + (Jy-Jx)wywx = Mzy + Mzв

где Jx, Jy, Jz - главные моменты инерции,

My - управл€ющий момент,

Mв - возмущающий момент.

“ак как угловые скорости  ј малы, следовательно, можно пре≠небречь произведением угловых скоростей, значит, уравнени€ Ёй≠лера имеют вид:

Jxwx = Mxy + Mxв

Jywy = Myy + Myв

Jzwz = Mzy + Mzв

√лавные моменты инерции:

Jx = 532 кг´м2, Jy = 563 кг´м2, Jz = 697 кг´м2.

÷ентробежные моменты инерции принимаютс€ равными 0.

¬озмущающий момент Mв возникает из-за того, что двигатель коррекции расположен не в центре масс  ј, и реактивна€ т€га, ли≠ни€ действи€ которой находитс€ на удалении (плече) l от центра масс  ј, создает паразитный крут€щий момент Mв.

Mв = P´l,

где P = 25 H - т€га корректирующего двигател€,

l = 4 мм - плечо.

“аким образом, Mв = 25´0,0004 = 0,1 Ќм.

2.6.2. —“јЅ»Ћ»«ј÷»я ”√Ћќ¬ќ√ќ ѕќЋќ∆≈Ќ»я ѕ–»  ќ––≈ ÷»»

ќсновное требование, предъ€вл€емое в этом режиме:

- точность поддержани€ направлени€ импульса коррекции - не хуже 1 угл.мин.

÷елью данной главы €вл€етс€ исследование динамики системы при стабилизации углового положени€ при коррекции.

‘ункциональна€ схема ћ ј состоит из следующих эелементов:

1) ћ ј - малый космический аппарат.

ћ ј описываетс€ как абсолютно твердое тело.

2) ƒ”— - датчик угловой скорости.

¬ качестве ƒ”— используетс€ командный гироскопический при≠бор. ќн описываетс€ колебательным звеном с параметрами T = 1/30 c-1 и e = 0,7, а также нелинейным звеном с насыщением 2∞/сек.

3) ј÷ѕ - аналогово-цифровой преобразователь.

ѕреобразует аналоговый сигнал с ƒ”— в цифровой сигнал.

4) ÷јѕ - цифро-аналоговый преобразователь.

ѕреобразует цифровой сигнал с ÷¬ћ в аналоговый.

5) Ў»ћ - широтно-импульсный модул€тор.

ѕредназначен дл€ формировани€ скважности импульсов управ≠лени€ двигателем стабилизации, пропорциональной управл€ю≠щему напр€жению. ¬ этом случае мы имеем среднее значение управл€ющего момента, пропорциональное управл€ющему сиг≠налу.

“ак как динамика ÷јѕ, ј÷ѕ, Ў»ћ как электронных аналого≠вых приборов оказывает на систему незначительное вли€ние по сравнению с динамикой механических (ƒ”—, двигатели) динамиче≠ские звень€, описывающие эти элементы, можно заменить соответ≠ствующими коэффициентами усилени€. ¬ первом приближении значени€ коэффициентов не принципиально.

6) ƒвигатель стабилизации.

ƒвигатель описываетс€ нелинейностью с насыщением 0,127 Ќм и звеном запаздывани€ с “д = 0,05 сек.

“€га двигател€ 0,1 Ќ

7) ÷¬ћ.

¬ ÷¬ћ формируетс€ управление по углу и угловой скорости. «а≠кон управлени€ имеет вид:

e = K(K1j +K2j),   = 1,  1 = 550,  2 = 430.

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

—истема была промоделирована по каналу х. ƒл€ других каналов схемы моделировани€ будут аналогичными.

ƒл€ разомкнутой системы были по≠строены Ћј„’ и ‘„’. Ёти графики представлены на рис.43.

–езультаты моделировани€ замкнутой системы представ≠лены на рис.44-46.

“аким образом, в результате моделировани€ получено, что про≠цесс стабилизации углового положени€ происходит примерно за 15 сек., статическа€ точность поддержани€ углового положени€ - 0,62 угл.мин., что полностью удовлетвор€ет требовани€м технического задани€.


3. ќ–√јЌ»«ј÷»ќЌЌќ-Ё ќЌќћ»„≈— јя „ј—“№

3.1. ќ–√јЌ»«ј÷»я » ѕЋјЌ»–ќ¬јЌ»≈ ¬џѕќЋЌ≈Ќ»я “≈ћџ

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

ќрганизаци€ дипломной работы по любой теме складываетс€ из определЄнных этапов и подэтапов, каждый из которых хот€ и может иметь разное содержание, однако структурно занимает равное положение дл€ всех дипломных работ, выполн€емых в данной отрасли.

“аким образом, структура дипломной работы может быть сформирована по типовой схеме, упор€доченной в соответствии с конкретным видом исследовани€. —остав дипломной работы по заданной теме, а также потребные категории исследований по этапам и подэтапам представлены в табл.1.

Ётапы

—одержание

»сполнители

1.

“ехническое задание

ѕостановка задачи. ќпределение состава программного продукта.

–уководитель –азработчик

2.

Ёскизный проект

–азработка общего описани€ программного продукта.

–уководитель –азработчик

3.

“ехнический проект

–азработка структуры программного продукта.

–азработчик

4.

–абочий проект

ѕрограммировани€ и отладка программы. ѕроверка результатов и внесение корректив в программу.

–уководитель–азработчик

5.

¬недрение

ќформление необходимой документации.

–азработчик

3.2. ќѕ–≈ƒ≈Ћ≈Ќ»≈ «ј“–ј“ “–”ƒј

ѕервым шагом при определении себестоимости программного комплекса €вл€етс€ расчет трудоемкости создани€ и внедрени€. –асчет производитс€ по методике, приведенной в документе Ђ“иповые нормы времени на программирование задач дл€ Ё¬ћї. “иповые нормы времени предназначены дл€ определени€ затрат времени на разработку программных средств вычислительной техники (ѕ—¬“).

»сходными данными дл€ расчета трудоемкости, при разработке программы €вл€ютс€:

 оличество разновидностей форм входной информации - 2,

в том числе:

информации, получаемой от решени€ смежных задач - 1,

справочной, условно посто€нной информации (файл инициализации) - 1;

 оличество разновидностей форм выходной информации - 2,

в том числе:

печатных документов (временные диаграммы) - 1,

информации, наносимой на магнитные носители (файл инициализации) - 1;

—тепень новизны комплекса задач - √ (разработка программной продукции, основанной на прив€зке типовых проектных решений).

—ложность алгоритма - 3 (реализуютс€ стандартные методы решени€, не предусмотрено применение сложных численных и логических методов).

¬ид используемой информации:

количество разновидностей форм переменной информации (ѕ») - 1, в том числе: информации, получаемой от решени€ смежных задач - 1;

количество разновидностей форм нормативно-справочной информации (Ќ—») (файл инициализации) - 1;

язык программировани€ - Borland —++.

¬ид представлени€ исходной информации - группа 11 (требуетс€ учитывать взаимовли€ние различных показателей).

¬ид представлени€ выходной информации - группа 22 (вывод информационных массивов на машинные носители).

“рудоемкость разработки программного продукта tпп† может быть определена как сумма величин трудоемкостей выполнени€ отдельных стадий разработки программного продукта.

tпп = tтз + tэп + tтп + tрп + tв,†

где tтз - трудоемкость разработки технического задани€ на создание программного продукта,

tэп - трудоемкость разработки эскизного проекта программного продукта,

tтп - трудоемкость разработки технического проекта программного продукта,

tрп - трудоемкость разработки рабочего проекта программного продукта,

tв - трудоемкость внедрени€ программного продукта.

“рудоемкость разработки технического задани€ рассчитываетс€ по формуле

tтз = “зрз + “зрп,

где “зрз - затраты времени разработчика постановки задач на разработку “«, чел.-дни,

зрп - затраты времени разработчика программного обеспечени€ на разработку “«, чел.-дни.

«начени€ “зрз и “зрп рассчитываютс€ по формуле

зрз = tз зрз, “зрп = tз зрп,

где tз - норма времени на разработку “« дл€ программного продукта в зависимости от функционального назначени€ и степени новизны разрабатываемового программного продукта, чел.-дни (tз = 29),

 зрз - коэффициент, учитывающий удельный вес трудоемкости работ, выполн€емых разработчиком постановки задач на стадии “« ( зрз = 0,65),

 зрп - коэффициент, учитывающий удельный вес трудоемкости работ, выполн€емых разработчиком программного обеспечени€ на стадии “« ( зрп = 0,35).

зрз = 29*0,65 = 18,85 чел.-дней.

зрз = 29*0,35 = 10,15 чел.-дней.

†tтз = “зрз + “зрп =18,85 + 10,15 = 29 чел.-дней.

“рудоемкость разработки эскизного проекта рассчитываетс€ по формуле

tэп = “эрз + “эрп,

где “эрз - затраты времени разработчика постановки задач на разработку Ёѕ, чел.-дни,

эрп - затраты времени разработчика программного обеспечени€ на разработку Ёѕ, чел.-дни.

«начени€ “зрз и “зрп рассчитываютс€ по формуле

эрз = tэ эрз, “эрп = tэ эрп,

где tэ - норма времени на разработку Ёѕ дл€ программного продукта в зависимости от функционального назначени€ и степени новизны разрабатываемового программного продукта, чел.-дни (tэ = 41),

 эрз - коэффициент, учитывающий удельный вес трудоемкости работ, выполн€емых разработчиком постановки задач на стадии Ёѕ ( эрз = 0,7),

 эрп - коэффициент, учитывающий удельный вес трудоемкости работ, выполн€емых разработчиком программного обеспечени€ на стадии Ёѕ ( эрп = 0,3).

зрз = 41*0,7 = 28,7 чел.-дней.

зрз = 41*0,3 = 12,3 чел.-дней.

п = “зрз + “зрп = 28,7 + 12,3 = 41 чел.-дней.

“рудоемкость разработки технического проекта зависит от функционального назначени€ программного продукта, количества разновидностей входной и выходной информации и определ€етс€ как сумма времени, затраченного разработчиком постановки задач и разработчиком программного обеспечени€:

tтп = (tтрз + tтрпв р,

где tпрз, tпрп - норма времени на разработку “ѕ разработчиком постановки задач и разработчиком программного обеспечени€ соответственно, чел.-дни (tтрз = 9, tтрп = 8),

 в - коэффициент учета вида используемой информации,

 р - коэффициент учета режима обработки информации ( р = 1,1).

«начение коэффициента  в определ€етс€ по формуле

 в = ( пnп +  нсnнс +  бnб)/(nп + nнс + nб),

где  п,  нс,  б - значени€ коэффициентов учета вида используемой информации дл€ переменной, нормативно-справочной информации и баз данных соответственно ( п = 0,5,  нс = 0,43,  б = 1,25),

nп, nнс, nб - количество наборов данных переменной, нормативно-справочной информации и баз данных соответственно (nп = 1, nнс = 1, nб = 0).

 в = (0,5 + 0,43)/2 = 0,465

tтп = (8 + 9)*0,465*1,1 = 8,6955 чел.-дней.

“рудоемкость разработки рабочего проекта зависит от функционального назначени€ программного продукта, количества разновидностей входной и выходной информации, сложности алгоритма функционировани€, сложности контрол€ информации, степени использовани€ готовых программных модулей и рассчитываетс€ по формуле

tрп = (tррз + tррпк р  з иа,

где tррз, tррп - норма времени на разработку –ѕ разработчиком постановки задач и разработчиком программного обеспечени€ соответственно, чел.-дни (tррз = 5, tррп = 27),

 к - коэффициент учета сложности контрол€ информации ( к = 1,07),

 - коэффициент учета уровн€ используемового €зыка программировани€ (  = 1,0),

 з - коэффициент учета степени использовани€ готовых программных модулей ( з = 0,8),

 иа - коэффициент учета вида используемой информации, и сложности алгоритма программного продукта.

«начение коэффициента  иа определ€етс€ по формуле

 иа = ( Тпnп +  Тнсnнс +  Тбnб)/(nп + nнс + nб),

где  Тп,  Тнс,  Тб - значени€ коэффициентов учета сложности алгоритма программного продукта и вида используемой информации дл€ переменной, нормативно-справочной информации и баз данных соответственно ( Тп = 0,48,  Тнс = 0,29,  Тб = 0,24),

nп, nнс, nб - количество наборов данных переменной, нормативно-справочной информации и баз данных соответственно (nп = 1, nнс = 1, nб = 0).

 иа = (0,48 + 0,29)/2 = 0,385

tрп = (5 + 27)*1,07*1,1*0,8*0,385 = 11,6 чел.-дней.

“рудоемкость внедрени€ может быть рассчитана по формуле:

tв = (tврз + tврп)  к р з,

где tврз, tврп - норма времени на внедрение программного продукта разработчиком постановки задач и разработчиком программного обеспечени€ соответственно, чел.-дни (tврз = 8, tврп = 24).

tтп = (8 + 24)*0,8*1,07 = 27,392 чел.-дней.

tпп = 29 + 41 + 8,6955 + 11,6 + 27,392 = 117,6875 чел.-дней.

ѕродолжительность выполнени€ всех работ по этапам разработки программного продукта рассчитываетс€ по формуле

Ti = (ti + Q)/ni,

где ti - трудоемкость i-й работы, чел.-дни,

Q - трудоемкость дополнительных работ, выполн€емых исполнителем, чел.-дни,

ni - количество исполнителей, выполн€ющих i-ю работу, чел.

Tтз = tтз/2 = 29/2 = 14,5 (15) дней

Tэп = tэп/2 = 41/2 = 20,5 (21) дней

Tтп = tтп = 8,6955 (9) дней

Tрп = tрп/2 = 11,6/2 = 5,8 (6) дней

Tв = tв = 27,392 (28) дней

Tпп = S“i = 15 + 21 + 9 + 6 + 28 = 79 дней.

3.3. –ј—„≈“ —ћ≈“џ «ј“–ј“ Ќј –ј«–јЅќ“ ” ѕ–ќ√–јћћЌќ√ќ ѕ–ќƒ” “ј

—мета затрат на выполнение работ составл€етс€ по калькул€ционным стать€м. ¬ общем случае статьи, учитывающие расходы, следующие:

1. ћатериалы (суммарные затраты на материалы, приобретаемые дл€ разработки программного продукта).

2. —пециальные оборудование (суммарные расходы на аренду приборов, требуемых дл€ разработки программного продукта).

3. ќсновна€ заработна€† плата производственного персонала.

4. ƒополнительна€ заработна€† плата.

5. ќтчисление на социальное страхование.

6. Ќакладные расходы.

7. ѕроизводственные командировки.

8.  онтрагентские расходы.

ќднако затраты, св€занные с разработкой программного обеспечени€, нос€т специфический характер. –асходы по стать€м 7, 8 обычно крайне незначительны. —татьи 1, 2 св€заны с расходами на использование Ё¬ћ. Ёти расходы определ€ютс€, исход€ из затрат машинного времени и стоимости часа работ Ё¬ћ, а также стоимости необходимых материалов и покупных изделий, необходимых при работе на Ё¬ћ.

¬ результате можно определить следующие статьи расходов на разработку программного продукта:

1. —тоимость машинного времени, затраченного на разработку.

2. —тоимость материалов и покупных изделий.

3. ќсновна€ заработна€ плата исполнителей.

4. ƒополнительна€ заработна€† плата.

5. ќтчисление на социальное страхование.

6. Ќакладные расходы.

–асчет стоимости затраченного машинного времени.

эвм = tэвм эвми÷эвм эвмбд эвмэ,

где tэвм - врем€ использовани€ Ё¬ћ дл€ создани€ данного программного продукта, час (tэвм = 10).

 эвми - поправочный коэффициент учета времени использовани€ Ё¬ћ ( эвми= 0,6).

÷эвм - цена одного часа работы Ё¬ћ, руб (на €нварь 1997 года 8000 рублей).

 эвмбд - коэффициент учета степени использовани€ —”Ѕƒ ( эвмбд = 1,0).

 эвмэ - коэффициент учета быстодействи€ Ё¬ћ ( эвми= 1,2).

эвм = 10*0,6*1,2*8000 = 57600 рублей.

–асчет затрат на материалы.

¬ процессе разработки программного издели€ используютс€ следующие материалы:

бумага формата ј4 в количестве 500 листов общей стоимости 60000 рублей.

дискеты 3,5 дюйма в количестве 2 штук по цене 5000 рублей за штуку.

м = 60000 + 5000*2 = 70000 рублей.

“аблица 4.2. –асходные материалы.

ћатериал

÷ена, руб

 оличество, шт.

—тоимость, руб.

дискета 1,44 ћб

5000

2

10000

бумага формата ј4

500

60000

»того:

70000

–асчет основной заработной платы исполнител€.

зо = Sзiti/d,

где зi - средн€€ заработна€ плата i-го исполнител€ (300000 рублей в мес€ц),

ti - трудоемкость работ, выполн€емых i-м исполнителем (чел.-дни),

d - среднее количество рабочих дней в мес€це (d = 22).

зо = 79*300000/22 = 1077300 рублей.

–асчет дополнительной заработной платы.

¬ статье Удополнительна€ заработна€ платаФ учитываютс€ выплаты, предусмотренные законодательством о труде и коллективными договорами за непроработанное на производстве врем€.

ƒополнительна€ заработна€ плата определ€етс€ по установленному нормативу от основной заработной платы по формуле

зд = —зоa,

где a - коэффициент дополнительной заработной платы, a = 0,2.

зд = 1077300*0,2 = 215500 рублей.

–асчет отчислений на социальное страхование.

¬ статью Уотчислени€ на социальное страхованиеФ включено отчисление по единому установленному нормативу от суммы основной и дополнительной† заработной платы.

–азмер отчислений вычисл€етс€ по формуле

сс = (—зд + —зо)aсс,

где aсс - коэффициент, устанавливающий отчисление на социальное страхование и в фонд стабилизации, aсс = 0,4.

сс = (1077300 + 215500)*0,4 = 517120 рублей.

–асчет накладных расходов.

¬ статье Унакладные расходыФ учитываютс€ командировочные расходы, оплата подъемных при перемещени€х, арендна€ плата, оплата услуг сторонних организаций.

н =† —зоaн,

где aн - коэффициент накладных расходов, aн = 1,8.

н = 1,8*1077300 = 1939140 рублей

–асчет суммарных расходов.

— = —эвм + —м + —зо + —зд + —сс + —н =

=† 57,6 +70 + 1077,3 + 215,5 + 517,12 + 1939,14 =† 3876,66 тыс.рублей.

—мета затрат на разработку программного продукта приведена в таблице 4.3.

“аблица 4.3.

є

п/п

Ќаименование статей расходов

«атраты
(тыс.руб.)

”дельный вес, %

1

—тоимость машинного времени

57,6

1,4

2

ћатериалы

70

1,9

3

ќсновна€ заработна€ плата

1077,3

27,7

4

ƒополнительна€ заработна€ плата

215,5

5,5

5

ќтчислени€ на социальное страхование

517,12

13,4

6

Ќакладные расходы

1939,14

50,1

»того:

3876,66


4. ѕ–ќћџЎЋ≈ЌЌјя Ё ќЋќ√»я » Ѕ≈«ќѕј—Ќќ—“№

4.1. ¬¬≈ƒ≈Ќ»≈

¬ результате развити€ производственных сил общества возник≠ла проблема взаимодействи€ человека и машины. ќхрана труда и эргономика позвол€ют с научной точки зрени€ подойти к этой проблеме, способствуют изучению вли€ни€ окружающей среды на человека, который непосредственно контактирует с Ё¬ћ, опреде≠лению вредных и опасных производственных факторов, разраба≠тывают организационно-технические меропри€ти€, направленные на профилактику профессиональных заболеваний, создава€ здоровые и безопасные услови€ труда дл€ работающего.

ѕредметом исследовани€ эргономики в этой области стало согласование психо-физических возможностей человека со свойствами современных технических систем. “олько в этом случае можно рассчитывать на высокое качество и эффективность его труда. ќсобую актуальность эта проблема приобретает в св€зи с возросшим культурным уровнем современного персонала, предъ€вл€юшего повышенные требовани€ к содержанию и услови€м труда на рабочем месте (–ћ).

ѕод рабочим местом в эргатических системах (Ё—) согласно √ќ—“ 26387-84 понимаетс€ Ђчасть пространства в системе человек-машина (—„ћ), оснащенна€ средствами отображени€ информации, органами управлени€, вспомогательным оборудованием и предназначенна€ дл€ осуществлени€ де€тельности оператора —„ћї. —оответственно среда на –ћ определ€етс€ этим же √ќ—“ом как Ђсовокупность физических, химических и психологических факторов, воздействующих на оператора —„ћ, на его –ћ в ходе его де€тельностиї. 

4.2. јЌјЋ»« ¬–≈ƒЌџ’ ‘ј “ќ–ќ¬  

Ќормальна€ и безопасна€ работа инженера-программиста за экраном диспле€ во многом зависит от того, в какой мере услови€ его работы соответствуют оптимальным. ѕри этом под услови€ми работы подразумевают комплекс физических, химических, биологических и психофизических факторов, установленных стандартами по безопасности труда (——“Ѕ).

  физическим факторам относ€тс€:

- вибраци€ и шум из-за движущихс€ машин, механизмов и их элементов, запыленность и загазованность воздуха, температура поверхностей оборудовани€, материалов и воздуха;

- плотность воздуха, ее резкое изменение, подвижность и ионизаци€ воздуха;

- ионизирующие и электромагнитные излучени€, статические зар€ды и повышение напр€жени€ в цепи, электрические и магнитные пол€;

- отсутствие или недостаток естественного света, повышенна€ или пониженна€ освещенность, €ркость и контрастность, блесткость поверхности, пульсаци€ светового потока;

- ультрафиолетовое или инфракрасное излучение.

  химическим факторам относ€тс€:

- общетоксические, раздражающие, сенсибилизирующие, канцерогенные, мутагенные;

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

  биологическим факторам относ€тс€:

- микроорганизмы (бактерии, вирусы, грибы и т.д.);

- макроорганизмы (растени€ и животные). 

  психофизическим факторам относ€тс€ перегрузки:

- физические (статические, динамические, гиподинами€);

- нервно-психические (умственное перенапр€жение, монотонность труда, эмоциональные перегрузки). 

ѕри проектировании рабочего места инженера-программиста необходимо учитывать и нормировать все указанные группы факторов, поскольку при определенных услови€х они могут вызвать нежелательные функциональные сдвиги в организме оператора, снизить качество и эффективность его работы, оказать отрицательное вли€ние на его здоровье.  

Ќаиболее значительным фактором €вл€етс€ микроклимат, особенно температура и влажность воздуха. »сследовани€ показывают, что высока€ температура в сочетании с высокой влажностью воздуха оказывают большое вли€ние на работоспособность человека. –езко увеличиваетс€ врем€ сенсорных и моторных реакций, нарушаетс€ координаци€ движений, увеличиваетс€ количество ошибок. ¬ысока€ температура отрицательно сказываетс€ и на р€де психологических функций человека. ”меньшаетс€ объем оперативной пам€ти, резко суживаетс€ способность к ассоциаци€м. ѕри +110— начинаетс€ окоченение конечностей, така€ температура минимально допустима. Ќаиболее благопри€тный диапазон температур в летнее врем€ от +180— до +240—, в зимнее врем€ от +170 до +220—.

ƒвижение воздуха позвол€ет увеличить рабочий диапазон температур. “ак при скорости движени€ воздуха 0.1, 0.5, 0.9 м/с верхн€€ допустима€ граница рабочего диапазона сдвигаетс€ соответственно до +220, +240, +260— при интенсивном расходе энергии человеком пор€дка 1000 ƒж/ч.

јтмосферное давление в пределах 80-106 кѕа легко переносимо человеком. ѕри давлени€х, выход€щих за эти пределы, человеку требуетс€ предварительна€ акклиматизаци€. 

–езультаты работы инженера-программмиста в большой степени завис€т и от освещенности рабочего места. „тобы правильно спланировать рациональную систему освещени€, необходимо учитывать €ркость источников света, их расположение в помещении, €ркостной контраст между устройствами Ё¬ћ и фоном, блесткость поверхностей, качество и цвет светильников и поверхностей. ƒл€ малой и средней контрастности поверхностей Ё¬ћ при темном фоне наименьший уровень освещенности должен быть 150 лк. ƒл€ большой контрастности при светлом или темном фоне наименьший уровень освещенности 100 лк. 

¬ помещени€х, где эксплуатируют Ё¬ћ, необходимо предусматривать систему искусственного освещени€ из люминисцентных ламп дневного света или ламп† накаливани€. —уществуют пр€ма€, отраженна€ и диффузна€ системы искусственного освещени€. ѕри пр€мом освещении свет попадает на объект непосредственно от источников света. ѕри этом 90-100% мощности светильника направлено на рабочую поверхность, что вызывает €ркостные контрасты, резкие тени и блесткость (свойство €рко освещенной поверхности вызывать ослепление или дезадаптацию наблюдател€). ѕри освещении отраженным светом 90-100% света направл€етс€ на потолок и верхнюю часть стен, от которых свет более или менее равномерно отражаетс€ по всему помещению. ѕри этом достигаетс€ равна€ освещенность без теней и блесткости. ƒиффузное освещение обеспечивает рассе€нный свет, одинаково распределенный по всем направлени€м. “ака€ система освещени€ требует меньшей мощности, чем две предыдущие, но вызывает частичное образоование теней и блесткости.

 роме освещенности, большое вли€ние на де€тельность человека оказывает цвет окраски помещени€ и спектральные характеристики используемого цвета. –екомендуетс€, чтобы потолок отражал 80-90%, стены - 50-60%, панели - 15-20%, а пол - 15-30% падающего на них света.  роме того, цвет обладает некоторым психологическим и физиологическим действием. “ак, например, применение тонов теплой гаммы (красный, оранжевый, желтый) создает впечатление бодрости, возбуждени€ и замедленного течени€ времени. Ёти же цвета вызывают у человека ощущение тепла. 

Ѕольшое вли€ние на де€тельность инженера-программиста  оказывает и уровень акустического шума. Ўум резко снижает  производительность труда и увеличивает травматизм. ‘изиологически шум воздействует на органы зрени€ и слуха, повышает кров€ное давление, при† этом притупл€етс€ внимание.

Ўум оказывает также и эмоциональное воздействие: он €вл€етс€ причиной возникновени€ таких отрицательных эмоций, как досада, раздражение. ќсобенно непри€тны высокочастотные и прерывистые шумы. 

ќсновным из механических факторов производственной среды €вл€ютс€ вибрации. ќни не только вредно воздействуют на организм, но и мешают человеку выполн€ть как мыслительные так и двигательные операции. ѕод действием вибраций ухудшаетс€ зрительное воспри€тие, в осообенности на частотах между 25 и 40 √ц и между 60 и 90 √ц. Ќаиболее опасна вибраци€ с частотой 6-8 √ц, так как в этом диапазоне лежит собственна€ резонансна€ частота тела, головы и брюшной† полости человека.

  числу неблагопри€тных факторов относ€тс€ злектромагнитные пол€ (Ёћѕ) высоких частот. »х воздействие на человека может вызвать функциональные сдвиги в организме: быструю утомл€емость, головные боли, нарушение сна, раздражительность, утомление зрени€ и т.п. 

ѕредельно допустимые уровни Ёћѕ следующие: 

- в —¬„-диапазоне - мк¬т/см;

- в диапазоне до 300 ћ√ц по электрической составл€ющей - 5 ¬/м, по  магнитной составл€ющей - 5 ј/м. — учетом этого стандарта было исследовано свыше 150 мониторов различных типов.

Ќа жизнеде€тельность человека большое вли€ние оказывает газовый состав воздуха. «десь обычно исследуетс€ две группы факторов: изменение обычного† состава воздуха (кислорода и углекислого газа) и посторонние добавки к нему в результате работы техники. 

Ѕлагопри€тными услови€ми газового состава воздуха считаетс€ содержание кислорода 19-20%, углекислого газа около 1%; допустимые значени€, при† которых не происходит выраженного снижени€ работоспособности составл€ют: кислорода - 18-29%, углекислого газа - 1-2%. —нижение содержани€ кислорода ниже 16% и повышение содержани€ углекислого газа выше 3% €вл€ютс€ недопустимыми и могут привести к нежелательным последстви€м. ¬ажнейшим способом борьбы с неблагопри€тным воздействием на человека химических факторов €вл€етс€ соблюдение их предельно допустимых концентраций в производственных помещени€х. ѕредельно допустимыми считаютс€ такие максимальные концентрации вредных веществ, которые при ежедневной работе не могут вызывать у работающих заболевани€ или отклонени€ в состо€нии здоровь€. “акими концентраци€ми считаютс€, например, дл€ аммиака - 20 мг/м, анилина - 3 мг/м, ацетона - 200 мг/м, бензола - 5 мг/м, бензина - 100 мг/м, серной кислоты - 1 мг/м.

ѕри выполнении данной дипломной работы используютс€ следующие элементы вычислительной техники:

персональный компьютер IBM PC 486DX;

струйный принтер Canon Bubble Jet.

ѕерсональный компьютер питаетс€ напр€жением 220¬/50√ц, которое превышает безопасный предел 42 ¬. —ледовательно возникает опасность поражени€ электрическим током.

¬оздействие на человека электрического тока приводит к общим травмам (электроудары) и местным (ожоги, металлизаци€ кожи, электрические знаки, электроофтальми€, механические повреждени€).

¬озникновение рентгеновского излучени€ обусловлено наличием на† аноде электронно-лучевой трубки диспле€ напр€жени€ до 30 к¬ (а при напр€жении 3-500 к¬ присутствует рентгеновское излучение различной жесткости). ѕользователь попадает в зону м€гкого рентгеновского излучени€.

ѕри воздействии рентгеновского излучени€ на организм человека происходит:

образование чужеродных соединений молекул белка, обладающих даже токсическими свойствами;

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

ѕри работе за экраном электронно-лучевой трубки диспле€ пользователь попадает под воздействие ультрафиолетового излучени€ с длинами волн < 320 нм. “акже при образовании строчной и кадровой разверток диспл分 возникает излучение электромагнитных полей частотой до 100 к√ц. Ёто может €вл€тьс€ причиной возникновени€ следующих заболеваний:

обострение некоторых заболеваний кожи (угрева€ сыпь, себорроидна€ экзема, розовый лишай, рак кожи и др.);

нарушение в протекании беременности;

увеличение в 2 раза веро€тности выкидышей у беременных женщин;

нарушение репродуктивной функции† и† возникновение рака;

нарушение режима терморегул€ции организма;

изменени€ в нервной системе (потер€ порога чувствительности);

понижение/повышение артериального давлени€.

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

ѕри работе за персональным компьютером дл€ вывода информации на бумажный носитель примен€етс€ принтер. ѕринтер Canon Bubble Jet имеет уровень звука на рассто€нии 1 метр от корпуса 49 дЅ (используетс€ 1 час в течении смены), что соответствует норме. —ледовательно, вредного воздействи€ по звуку на пользовател€ не оказываетс€.

“аким образом пользователь, работающий с персональным компьютером подвергаетс€ воздействию следующих опасных и вредных факторов:

поражение электрическим током;

воздействие рентгеновского излучени€;

ультрафиолетовое излучение и излучение электромагнитных полей радиочастот;

воздействие статического электричества.

4.3. “–≈Ѕќ¬јЌ»я   ¬»ƒ≈ќ“≈–ћ»ЌјЋ№Ќџћ ”—“–ќ…—“¬јћ

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

¬идеотерминальное устройство должно соответствовать следующим требовани€м:

€ркость свечени€ экрана не менее 100 кд/м2;

минимальный размер свет€щейс€ точки не более 0,4 мм дл€ монохромного диспле€ и не более 0,6 мм дл€ цветного;

контрастность изображени€ знака не менее 0,8;

частота регенерации изображени€ при работе с позитивным контрастом в режиме обработки текста не менее 72 √ц;

количество точек на экране не менее 640;

экран должен иметь антибликовое покрытие;

размер экрана должен быть не менее 31 см по диагонали, а высота символов не менее 3,8 мм, при этом рассто€ние от экрана до глаз оператора должно быть 40Ц80 см.

ѕри работе с текстовой информацией наиболее предпочтительным €вл€етс€ предъ€вление чЄрных знаков на светлом (белом) фоне.

ћаксимальные значени€ напр€женности магнитного пол€, измеренные на рассто€нии 50 см от экранов наиболее распространЄнных мониторов.

ѕолоса частот

ћагнитное поле,
ј/м

Ќормы BGA

5 - 1000 √ц

0,2

160 - 0,8

10 - 150 к√ц

0,17

0,8 - 0,6

150 - 300 к√ц

-

0,6 - 0,42

0,3 - 30 ћгц

0,00000066

0,42 - 0,73

30 - 300 ћгц

0,00000066

0,73

ћаксимальна€ напр€женность электрического пол€, допускаема€ нормами BGA, равна 2,5 к¬/м. Ёто значение установлено из расчЄта того, чтобы при прикосновении к зар€женной провод€щей поверхности электрический разр€д не стал причиной шока.

ћаксимальные значени€ напр€женности электрического пол€, измеренные на рассто€нии 50 см от экранов наиболее распространЄнных мониторов.

ѕолоса частот

Ёлектрическое поле,
¬/м

Ќормы BGA

5 - 1000 √ц

4,8

2500 - 177

10 - 150 к√ц

4,8

87

150 - 300 к√ц

0,48

87

0,3 - 30 ћгц

0,0024

87 - 27,5

30 - 300 ћгц

0,0024

27,5

»змерени€ BGA показывают, что напр€женность электростатического пол€ около монитора может превысить 7 к¬/м. —огласно полученным SSI и SEMKO (Ўвеци€) данным, эти значени€ дл€ некоторых устройств достигают 50 к¬/м.

¬ –оссии нормирование электромагнитных полей осуществл€етс€ в соответствии с √ќ—“ 12.1.006-84 и санитарными нормами —Ќиѕ2963-84.

¬ зоне индукции нормируетс€ напр€женность электрического и магнитного пол€ в зависимости от частоты. ¬ зоне излучени€ нормируетс€ плотность потока энергии в зависимости от времени пребывани€.

Ќормир.

„астота f, ћ√ц

велич.

0.06-1.5

1.5-3.0

3.0-30

30-50

50-300

300-3*105

≈, ¬/м

50

50

20

10

5

нет

Ќ, ¬/м

5.0

Ц

Ц

0.3

Ц

нет

I, ¬б/м2

Ц

Ц

Ц

Ц

Ц

I0 = e/T

Ёлектромагнитные пол€ нормируютс€ следующим образом:

электрические: E = 6/ÖT; 1 £ T £ 9, где “- врем€ воздействи€;

магнитные: Hn £ 8 кј/м в течение рабочего дн€; e = 2 (¬т r/м2) Ц энергетическа€ нагрузка на организм.

4.4. –ј—„≈“ ¬–≈ƒЌџ’ »«Ћ”„≈Ќ»… ¬»ƒ≈ќƒ»—ѕЋ≈я

¬рем€ работы на персональном компьютере по санитарным нормам не должно превышать 4 часа.Ѕольшинство используемых в –оссии мониторов не соответствуют шведскому стандарту защита пользовател€ от излучений и имеют† на рассто€нии 5 см от экрана диспле€ имеют мощность дозы рентгеновского излучени€ 100 мк–/час. –ассчитаем, какую дозу рентгеновского излучени€ получит пользователь на различном рассто€нии от экрана диспле€.

Pr = P0e-mr, где

Pr - мощность дозы рентгеновского излучени€ на рассто€нии r, мк–/час;

P0 - уровень мощности дозы рентгеновского излучени€ на рассто€нии 5 см от экрана диспле€, мк–/ч.

m - линейный коэффициент ослаблени€ рентгеновского излучени€ воздухом, см-1;

r - рассто€ние от экрана диспле€, см;

¬озьмем m = 3.14*10-2 см-1.

r, см

5

10

20

30

40

50

60

70

80

90

100

P, мк–/ч

100

73

53

39

28

21

15

11

8

6

4

—реднестатистический пользователь располагаетс€ на рассто€нии 50 см от экрана диспле€. –ассчитаем дозу облучени€, которую получит пользователь за смену, за неделю, за год.

«а смену

4 часа

4*21

84 мк–/ч

«а неделю

5 дней

5*84

420 мк–/ч

«а год

44 рабочие недели

44*420

18480 мк–/ч

4.5. –ј÷»ќЌјЋ№Ќјя ќ–√јЌ»«ј÷»я –јЅќ„≈√ќ ћ≈—“ј

ƒл€ повышени€ производительности труда при работе за компьютером необходимо создать на рабочем месте наиболее благопри€тные услови€ с точки зрени€ эргономики и эстетики.

–азработка меропри€тий по рациональной организации рабочего места инженера-программиста и инженера-разработчика может идти в следующих направлени€х:

устранение неблагопри€тных факторов:

снижение шума в помещении;

правильный выбор источников освещени€;

устранение запылЄнности и загазованности.

оптимизаци€ условий труда на рабочем месте:

эргономические требовани€;

психологические требовани€.

создание комфортных условий отдыха в течение рабочего дн€.

ѕроизводственные помещени€ вычислительного центра должны проектироватьс€ в соответствии с требовани€ми —Ќиѕ 2.03.04-87 Ц Ујдминистративные и бытовые здани€ и помещени€ производственных предпри€тийФ.

ѕлощадь помещени€ следует принимать из расчЄта 6 м2 на одного работника. ѕри оснащении рабочих мест терминалами Ё¬ћ, печатающими устройствами и пр. площади помещени€ допускаетс€ увеличивать в соответствии с техническими услови€ми на эксплуатацию оборудовани€.  убатура должна быть не менее 19,5 м3 с учЄтом максимального числа одновременно работающих.

ћинимальна€ ширина проходов с передней стороны пультов и панелей управлени€ Ё¬ћ при однор€дном расположении должна быть не менее 1 м, при 2-х р€дном расположении не менее 1,2 м. ¬идеотерминалы должны располагатьс€ при однор€дном размещении на рассто€нии не менее 1 м от стен. –абочие места с диспле€ми должны располагатьс€ между собой на рассто€нии не менее 1,5 м.

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

4.6. –≈ ќћ≈Ќƒј÷»» ѕќ —Ќ»∆≈Ќ»ё ”“ќћЋя≈ћќ—“»

Ќеобходимо расположить экран диспле€ немного выше уровн€ глаз. Ёто создаст разгрузку тех групп окологлазных мышц, которые наиболее напр€жены при обычном направлении взгл€да - вниз или вперЄд.

ѕомещение, где наход€тс€ компьютеры и видеомониторы, должно быть достаточно просторным с посто€нным обновлением микроатмосферы. ћинимальна€ площадь на один видеомонитор - 9-10 м2.  райне нежелателен визуальный контакт работника с другими мониторами или телевизионными экранами. Ќеобходимо исключить наличие всевозможных бликов на экране монитора, часто возникающих на стекл€нных экранах. —ледует также избегать большой контрастности между €ркостью экрана и окружающего пространства - оптимальным считаетс€ выравнивание €ркости экрана и компьютера. «апрещаетс€ работа с компьютером в тЄмном или полутЄмном помещении.

¬ечернее освещение рабочего помещени€ желательно голубоватого цвета с €ркостью, примерно равной €ркости экрана. ¬ услови€х дневного освещени€ также рекомендуетс€ обеспечить вокруг монитора голубой фон - за счЄт окраски стен или хот€ бы наличи€ плакатов.

ƒл€ большего эргономического комфорта целесообразно расположить в кресле опору - в районе по€сничного изгиба позвоночника (в виде продолговатой подушечки или валика).

≈сли работник имеет те или иные рефракционные отклонени€ (близорукость, дальнозоркость и др.), то последние должны быть полностью коррегированы очками. ѕри более серьЄзных отклонени€х вопрос о возможности работы с видеотерминалами должен решатьс€ с участием врача-офтальмолога.

„ерез каждые 40-45 минут необходимо проводить физкультурную микропаузу: вращение глаз по часовой стрелке и обратно, лЄгкие гимнастические упражнени€ дл€ всего тела, например поднимание и опускание рук.

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

4.7. «јў»“ј ќ“ Ќјѕ–я∆≈Ќ»я ѕ–» ќ—Ќќ¬≈Ќ»я. «јЌ”Ћ≈Ќ»≈

«анулением называетс€ преднамеренное соединение нетоковедущих частей с нулевым защитным проводником (Ќ«ѕ). ќно примен€етс€ в трехфазных сет€х с глухозаземленной нейтралью в установках до 1000 вольт и €вл€етс€ основным средством обеспечени€ электробезопасности.

ѕри попадании напр€жени€ сети на корпус ѕЁ¬ћ возникает режим короткого замыкани€. ƒл€ защиты электрической сети от короткого замыкани€ и перегрузок примен€ютс€ автоматические выключатели или предохранители. ѕри проектировании защитного устройства необходимо рассчитать его номинальный ток срабатывани€ - Iном:

Ialarm ³ KIном, где

Iном = Ialarm/K

Iном - номинальный ток срабатывани€ защитного устройства, A;

K - коэффициент, учитывающий тип защитного устройства:

K = 3 - дл€ автомата с электромагнитным расцепителем,

K = 1.4 - дл€ теплового автомата,

Ialarm - ток короткого замыкани€, A.

–ассчитаем величину тока короткого замыкани€:

Ialarm = Uf/(Rn + Rm/3)

Rn = Rf + R1 + jx1

Uf = 220 ¬

Rm = 0,312W

Rf = 0,412W

jx1 = 0,6W

R1 = r/S

r - удельное сопротивление Ќ«ѕ, [Wmm2/m];

l - длина Ќ«ѕ, m;

rcu = 0,0175 W mm2 /m,

l = 50 m,

S = 1,5 mm2

R1 = 0,0175(50/15) = 0,58W

Rn = (0,412 + 0,58 + 0,6) = 1,59W

Ialarm = 130 A

Iном = 43 A

ƒл€ того, чтобы в случае короткого замыкани€ или других причин ѕЁ¬ћ отключалась от электрической сети необходимо в цепь питани€ поставить автомат с электромагнитным расцепителем с Iном = 43 A.

4.8. ѕќ∆ј–Ќјя Ѕ≈«ќѕј—Ќќ—“№

¬ помещени€х ¬÷ существуют все три основные фактора, необходимые дл€ возникновени€ пожара.

√орючими материалами на ¬÷ €вл€ютс€: строительные материалы дл€ акустической и эстетической отделки помещений, перегородки, двери, полы, изол€ци€ силовых и сигнальных кабелей, шкафы, жидкости дл€ очистки элементов и узлов Ё¬ћ и т.д.

ƒл€ отвода тепла от Ё¬ћ в производственных помещени€х ¬÷ посто€нно действует система кондиционировани€. ѕоэтому кислород, как окислитель процессов горени€, имеетс€ в любой точке помещений ¬÷.

»сточниками зажигани€ на ¬÷ могут оказатьс€ электронные схемы Ё¬ћ, приборы, приборы, примен€емые дл€ технического обслуживани€, устройства электропитани€, кондиционеры воздуха.

ѕо пожарной опасности ¬÷ относ€тс€ к категории У¬Ф (в производстве обращаютс€ твердые сгораемые вещества и материалы). »сход€ из этого ¬÷ проектируетс€ с II степенью огнеустойчивости.

ћинимальные пределы огнеустойчивости в часах:

Cтепень огнестойкости зда≠ний и сооружений

II

ќсновные строительные конструкции:

Ќесущие стены, стены лестничных клеток, колонны

2

Ћестничные площадки

1

Ќаружние стены из навесных панелей

0,25

¬нутренние не≠сущие стены, перего≠родки

0,25

Ќесу≠щие конст≠рук≠ции междуэтажных перекрытий

0,75

ѕлиты, настилы и др.

0,25

ƒл€ обнаружени€ начальной стадии загорани€ используют систему автоматической пожарной сигнализации (јѕ—). јѕ— состоит из пожарных извещателей, линий св€зи и приемных пультов (станций).

¬ помещени€х ¬÷ примен€т дымовые пожарные извещатели типа –»ƒ-1.

ѕринцип действи€ –»ƒ-1 основан на изменении величины электрического тока, протекающего через ионизационную камеру, при попадании в нее дыма.

“ехнические показатели дл€ –»ƒ-1:

чувствительный элемент

ионизационна€ камера

параметр срабатывани€

тлеющий фитиль

инерционность, сек

10

диапазон температур, —

-30 ...+50

относительна€ влажность, %

80

Ќорма расстановки пожарных извещателей в помещени€х с гладким полом:

“ип

«ащищаема€ площадь, i2

–ассто€ние между извещател€ми, м

максимальное

в узких коридорах

–»ƒ-1

100

12

15

Ћинии св€зи систем јѕ— с приемными станци€ми стро€тс€ по лучевому принципу. ѕриемные станции јѕ— устанавливаютс€ в помещении дежурного по ¬÷, где организуетс€ круглосуточное дежурство.

ѕриемные станции обеспечивают следующие функции:

прием сигналов от пожарных извещателей с индикацией номера луча;

непрерывный контроль состо€ни€ лучей по всей длине с автоматическим вы€влением характера повреждени€;

светова€ и звукова€ сигнализаци€ тревоги;

автоматическое переключение на резервный источник питани€ при сбо€х сети с включением соответствующей сигнализации.

Ќа ¬÷ используетс€ приемна€ станци€ –ќ”ѕ-1.

“ехнические характеристики устройства –ќ”ѕ-1:

извещателей –»ƒ-1, шт

до 300

шлейфов блокировки, компл.

до 30

напр€жение питани€, ¬

220±10

потребл€ема€ мощность, ¬т

не более 180

диапазон температур, —

+5 ... +50

относительна€ влажность, %

до 80

срок службы, лет

8

дополнительные функции

может управл€ть устройствами пожаротушени€

Ќа ¬÷ примен€ютс€ установки газового тушени€ пожара, действие которых основано на быстром заполнении помещени€ газом с низким содержанием кислорода. »спользуетс€ автоматическа€ установка газового пожаротушени€ (ј”√ѕ) с электрическим пуском.

“ехнические характеристики ј”√ѕ с электрическим пуском:

число пусковых баллонов, шт

2

число рабочих баллонов, шт

4

зар€д пускового баллона

сжатый воздух

зар€д рабочего баллона

фреон 114 ¬ч

вместимость пускового баллона, л

27

вместимость рабочего баллона, л

40

давление в пусковом баллоне, ћѕа

125

давление в рабочем баллоне, ћѕа

12,5

продолжительность пуска, с

65

масса батареи без зар€да, кг

480

ѕри использовании ј”√ѕ дл€ предотвращени€ отравлени€ персонала предусмотрена предупредительна€ звукова€ и светова€ сигнализаци€, срабатывающа€ при ручном, дистанционном и автоматическом включении за 30 секунд до начала выпуска газа.

–асчет необходимого количества баллонов с сжатым воздухом и огнегасительным средством:

 оличество огнегасительного вещества (фреона)

Gт = GвWпKу, где Gт - количество огнегасительного вещества,

Wп - расчетный объем защищаемого помещени€, м3,

Gв - огнегасительна€ концентраци€ газового состава, кг/м3,

Kу - коэффициент, учитывающий особенности процессов газообмена в защищаемом помещении.

ƒл€ ¬÷ Gв= 0,25 кг/м3, Kу = 1,2.

Wп = SH, где S - площадь помещени€, м2.

H - высота помещени€, м.

S = 100 м2. H = 3 м. Wп = 300 м3.

Gт = 0,25*300*1,2 = 90 кг.

Ќеобходимое количество баллонов

Nб = Gт/Vбra, где Vб - объем баллона, м3,

r - плотность, кг/л,

a - коэффициент наполнени€ баллона.

Vб = 40 л, r = 2,17 кг/л, a = 0,9.

Nб = 90/(40*2,17*0,9) = 2.

ќбъем воздушных баллонов

Wб = (–смин+1)(Wс+Wт)/(–макс-–бмин), где –смин и –бмин - конечное давление в воздушных баллонах и баллонах с огнегасительным средством, ћѕа,

макс - начальное давление воздуха в баллоне, ћѕа,

Wс иWт - объем баллонов с огнегасительным составом и трубопроводов, л.

смин = –бмин = 5 ћпа, Wс = 2*40 = 80 л, Wт = 20л, –макс = 125 ћѕа.

Wб = (5+1)(80+20)/(125-5) = 4,8 л.


5. —ѕ»—ќ  Ћ»“≈–ј“”–џ.

1. Ђќсновы теории полета космических аппаратовї / ѕод ред. √.—.Ќариманова, ћ. .“ихонравова. ћ., ћашиностроение, 1972.

2. ј.ѕ.–азыграев Ђќсновы управлени€ полетом космических аппаратовї. ћ., ћашиностроение, 1990.

3. √.√.Ѕебенин, Ѕ.—.—кребишевский, √.ј.—околов Ђ—истемы управлени€ полетом космических аппаратовї. ћ., ћашиностроение, 1978.

4.  .Ѕ.јлексеев, √.√.Ѕебенин Ђ”правление космическими летательными аппаратамиї. ћ., ћашиностроение, 1974.

5. ¬.¬.—олодовников, ¬.Ќ.ѕлотников,  .¬.яковлев Ђ“еори€ автоматического управлени€ технических системї. ћ., изд.ћ√“” им.Ѕаумана, 1993.

6. Ѕ.—трауструп Ђязык программировани€ —++ї. ћ., Ђ–адио и св€зьї, 1991.

7. ј.¬.Ѕошкин, ѕ.Ќ.ƒубнер Ђ–абота с —++ї. ћ., Ђёкисї, 1991.

8. ¬.¬.јрсеньев, Ѕ.ё.—ажин Ђћетодические указани€ к выполнению организационно-экономической части дипломных проектов по созданию программной продукцииї, ћ., изд. ћ√“” им.Ѕаумана, 1994.

9. √ќ—“ 2.103-68 Ќ»–. ћ.: »зд-во стандартов, 1968.

10. ¬. .«елинский ЂЌќ“ в проектно-конструкторской организацииї. ћ.: ЂЁкономикаї, 1969.

11. Ђ”правление трудовым коллективомї / √.ѕ.«айцев, Ё.¬.ћинько, Ќ.¬.јртамонова и др. —вердловск, »зд-во ”√”, 1989.

12. Ђ“иповые нормы времени на программирование задач дл€ Ё¬ћї, утвержденные постановлением √осударственного комитета ———– по труду и социальным вопросам и —екретариата ¬÷—ѕ— от 27 июл€ 1987 г. є454/22-70

13. ё.√.—ибиров Ђќхрана труда в ¬÷ї. ћ., Ђћашиностроениеї, 1985.

14. —ибиров ё.√., Ђќсновы инженерной психологииї / под ред. Ѕ.‘.Ћомова. ћ., Ђћашиностроениеї, 1986.

15. —Ќиѕ 2.09.04-87 Ђјдминистративные и бытовые здани€ и помещени€ производственных предпри€тийї.

16. Ђ«рениеї / под ред. Ќ.». удр€шовой, ћ., Ђћашиностроениеї, 1995.

17. Ђ¬ременные рекомендации труда операторов за диспле€миї. √ќ—“ 12.1.006-84.

18. —Ќиѕ2963-84 ЂЌормирование электромагнитных полейї.

19. Ђ—овременные нормы электростатического и электромагнитного излучени€ї, ЂComputer Worldї є7, 1995.3


6. ѕ–»Ћќ∆≈Ќ»≈. “≈ —“џ ѕ–ќ√–јћћ ƒЋя BORLAND C++ » MATHLAB 4.0 FOR WINDOWS

6.1. ќ—Ќќ¬Ќќ… ѕ–ќ√–јћћЌџ… ћќƒ”Ћ№ MAIN.CPP

#include

#include

#include

#include

#include "rk5.h"

#include "sfun.h"

#include "init.h"

#include

typedef long double real;

const float g_r = M_PI/180.;

const float r_g = 180./M_PI;

real t_beg;

real t_end;

real dt;

real toler;

int Np;

int Curp;

real dTp;

real mu_z;

real mu_s;

real mu_l;

real m;

real m_t;

real W;

real w_s;

real w_z;

real w_l;

real ww_l;

real xs,ys,zs;

real xl,yl,zl;

real Fz,Fs,Fl,Fa,U20;

real J1,J2,J3;

int nomin;

real par[8];

real parn[8];

real a_p,e_p,p_p,Om_p,i_p,om_p,Rp_p,Ra_p;

real y_main[6];

real prmt[5];

int Fl_u;

real u_last;

int Fl_ka;

int Fl_kp;

int Fl_ki;

int Fl_i;

int Fl_p;

int Fl_a;

int Fl_lu;

int Fl_pkT;

real dl;

real T_vd;

real dRa;

real dRp;

int Sig;

int Sig_a;

real Tkor;

real Tkore;

real Vkor[3];

real akor[3];

int Fl_l0;

int Fl_l1;

int Fl_pki;

real dV_ps;

real dV_as;

real dV_is;

real dV_ss;

ofstream m_y ("m_y.dat");

ofstream m_f ("m_f.dat");

ofstream m_s ("m_s.dat");

ofstream m_l ("m_l.dat");

ofstream m_par ("m_par.dat");

ofstream u_f ("u_f.dat");

ofstream u_par ("u_par.dat");

ofstream k_par ("k_par.dat");

void out_p(real,real *,real*,int,int,real*);

void main()

{

clrscr();

init_m();

real dery[]={ .167, .167, .167, .167, .166, .166};

int ihlf;

int ndim = 6;

Drkgs(prmt,y_main,dery,ndim,ihlf,fct,out_p);

clrscr();

if (ihlf<11)

† {

† cout << 'n' << "”спешное завершение моделировани€" << 'n';

† cout << " t0 = " << t_beg << " tk = " << t_end << " dt = " << dt;

† cout << 'n' << "„исло делений шага=" << ihlf;

† }

else

† {

† cout << 'n' << "Ќенормальное завершение моделировани€" << 'n';

† cout << " t0 = " << t_beg << " tk = " << t_end << " dt = " << dt;

† cout << 'n' << "„исло делений шага=" << ihlf;

† }

getch();

m_y.close();

m_f.close();

m_s.close();

m_l.close();

m_par.close();

u_f.close();

u_par.close();

k_par.close();

}

void out_p(real x,real *y,real*,int,int,real*)

{

if (x >= (dTp*Curp))

† {

† Curp++;

† gotoxy(1,20);

† cout << "ѕроцесс выполнени€:" << float(Curp)*100./Np << " % " << 'n';

† cout.precision(7);

† m_y << x << 't' << y[0] << 't' << y[1] << 't' << y[2] << 't'

††††† << y[3] << 't' << y[4] << 't' << y[5] << 'n';

† m_f << x << 't' << Fz << 't' << Fs << 't' << Fl << 't' << Fa

††††† << 't' << U20 << 'n';

† m_s << x << 't' << xs << 't' << ys << 't' << zs << 'n';

† m_l << x << 't' << xl << 't' << yl << 't' << zl << 'n';

† m_par << x << 't' << par[0] << 't' << par[1] << 't' << par[2]

††††††††† << 't' << par[3] << 't' << par[4] << 't' << par[5]

††††††††† << 't' << par[6] << 't' << par[7] << 'n';

† }

if (Fl_u && (par[7] > parn[7]))

† {

† Fl_u = 0;

† dl = -(w_z-w_s)*(par[6]-parn[6]);

† u_par << x << 't' << par[0] << 't' << par[1] << 't' << par[2]

††††††††† << 't' << par[3] << 't' << par[4] << 't' << par[5]

††††††††† << 't' << par[6] << 't' << par[7] << 'n';

† u_f << x << 't' << Fz << 't' << Fs << 't' << Fl

††††† << 't' << Fa << 't' << U20 << 'n';

† }

if ((x > 79000) && (x < 81000))

† {

† k_par << x << 't' << par[5] << 't' << par[7] << 'n';

† }

}

6.2. ѕќƒѕ–ќ√–јћћј –ј—„≈“ј ¬ќ«ћ”ўјёў»’ ”— ќ–≈Ќ»…, ѕј–јћ≈“–ќ¬ ќ–Ѕ»“џ »  ќ––≈ ÷»» SFUN.CPP

#include "sfun.h"

const real p = 4.64e-6;

const real sm_s = 8.;

const real A = 1.496e11;

const real Cx = 2.;

const real sm_a = 2.5;

const real ro = 5.098e-13;

void korr (real& t, real *f, real *dery);

void fct(real& t, real *f, real *dery)

{

real x = f[0];

real y = f[1];

real z = f[2];

real Vx = f[3];

real Vy = f[4];

real Vz = f[5];

real Tet_s = (28.1+60*g_r)+w_s*t;

real e_0 = 23.45*g_r;

xs = A*cos(Tet_s);

ys = A*sin(Tet_s)*cos(e_0);

zs = A*sin(Tet_s)*sin(e_0);

real Tet_l = 0+w_l*t;

real Om_l = 0-ww_l*t;

real i_l = acos(cos(e_0)*cos(5.15*g_r)-sin(e_0)*sin(5.15*g_r)*cos(Om_l));

real rsr_l = 3.8448e8;

xl = rsr_l*(cos(Tet_l)*cos(Om_l)-cos(i_l)*sin(Tet_l)*sin(Om_l));

yl = rsr_l*(cos(Tet_l)*sin(Om_l)+cos(i_l)*sin(Tet_l)*cos(Om_l));

zl = rsr_l*sin(i_l)*sin(Tet_l);

real R_ka = sqrt(x*x+y*y+z*z);

real Fz_x = -mu_z*x/pow(R_ka,3.);

real Fz_y = -mu_z*y/pow(R_ka,3.);

real Fz_z = -mu_z*z/pow(R_ka,3.);

real mu_sd = p*sm_s*A*A/m;

real R_s = sqrt((x-xs)*(x-xs)+(y-ys)*(y-ys)+(z-zs)*(z-zs));

real Fs_x = -(mu_s-mu_sd)*x/pow(R_s,3.);

real Fs_y = -(mu_s-mu_sd)*y/pow(R_s,3.);

real Fs_z = -(mu_s-mu_sd)*z/pow(R_s,3.);

real R_l = sqrt((x-xl)*(x-xl)+(y-yl)*(y-yl)+(z-zl)*(z-zl));

real Fl_x = -mu_l*x/pow(R_l,3.);

real Fl_y = -mu_l*y/pow(R_l,3.);

real Fl_z = -mu_l*z/pow(R_l,3.);

real V_ka = sqrt(Vx*Vx+Vy*Vy+Vz*Vz);

real Fa_x = (-Cx*sm_a/(2*m))*ro*V_ka*Vx;

real Fa_y = (-Cx*sm_a/(2*m))*ro*V_ka*Vy;

real Fa_z = (-Cx*sm_a/(2*m))*ro*V_ka*Vz;

const real c20 = -1.09808e-3;

const real c22 = 5.74e-6;

const real d22 = -1.58e-6;

const real r_e = 6378137.;

real cr = mu_z*r_e*r_e/pow(R_ka,5);

real lr = 2*atan(y/x);

real mr = 3*(c22*cos(lr)+d22*sin(lr));

real U20_x = cr*x*(c20*(1.5-7.5*z*z/pow(R_ka,2))+mr*(5*z*z/pow(R_ka,2)-3));

real U20_y = cr*y*(c20*(1.5-7.5*z*z/pow(R_ka,2))+mr*(5*z*z/pow(R_ka,2)-3));

real U20_z = cr*z*(c20*(4.5-7.5*z*z/pow(R_ka,2))+5*mr*(z*z/pow(R_ka,2)-1));

dery[0] = Vx;

dery[1] = Vy;

dery[2] = Vz;

dery[3] = (Fz_x+U20_x+Fs_x+Fl_x+Fa_x+akor[0]);

dery[4] = (Fz_y+U20_y+Fs_y+Fl_y+Fa_y+akor[1]);

dery[5] = (Fz_z+U20_z+Fs_z+Fl_z+Fa_z+akor[2]);

Fz = sqrt(Fz_x*Fz_x+Fz_y*Fz_y+Fz_z*Fz_z);

Fs = sqrt(Fs_x*Fs_x+Fs_y*Fs_y+Fs_z*Fs_z);

Fl = sqrt(Fl_x*Fl_x+Fl_y*Fl_y+Fl_z*Fl_z);

Fa = sqrt(Fa_x*Fa_x+Fa_y*Fa_y+Fa_z*Fa_z);

U20 = sqrt(U20_x*U20_x+U20_y*U20_y+U20_z*U20_z);

parn[3] = parn[3]+w_s*t;

par_or(f,par);

korr(t,f,dery);

if ((u_last-par[7]) > 300*g_r)

† Fl_u = 1;

u_last = par[7];

}

void korr(real& t, real *f, real *)

{

if (t > (Tkor+172800.))

† {

† if ((fabs(dl) > 0.1*g_r) && (!Fl_ka) && (!Fl_kp) && (!Fl_ki))

††† {

††† Fl_kp = 1;

††† Fl_ka = 0;

††† Fl_ki = 0;

††† cout << "–езультат измерений накоплен" << 'n';

††† cout << "Ќеобходима коррекци€ периода. dl=" << dl*r_g << "град." << 'n';

††† cout << "ѕериод ном.=" << parn [6] << "ѕериод тек.=" << par[6] << 'n';

††† cout << "ѕараметры орбиты" << 'n';

††† cout << " Rp = " << par[2]*(1-par[1]) << 'n';

††† cout << " Ra = " << par[2]*(1+par[1]) << 'n';

††† cout << " p = " << par[0] << 'n';

††† cout† << " a = " << par[2] << " e = " << par[1] << "n T = "

††††††††† †<< par[6] << " w = " << par[5]*r_g << " u = " << par[7]*r_g

††††††††† †<< 'n';

††† clrscr();

††† }

† }

Fl_a = 0;

Fl_p = 0;

Fl_lu = 0;

real da;

if (par[5] > par[7])

† da = fabs(par[5]-par[7]-M_PI);

else

† da = fabs(par[5]-par[7]+M_PI);

if (da < .1*g_r)

† {

† Fl_a = 1;

† }

if (fabs(par[5] - par[7]) < .1*g_r)

† {

† Fl_p = 1;

† }

if (par[7] < .1*g_r )

† {

† Fl_lu = 1;

† }

real Vk;

if (T_vd)

if (t >= (T_vd +20))

† {

† T_vd = 0;

† akor[0] = 0;

† akor[1] = 0;

† akor[2] = 0;

† cout << "¬ыкл.дв. n t = " << t;

† }

if (((Fl_kp && Fl_a) || (Fl_ka && Fl_p) || (Fl_ki && Fl_lu)) && (!T_vd))

† {

† cout << " n  оррекци€ n";

† cout << "n Ќачало t=" << t << "сек n";

† int sim;

† if ((t-Tkor) < 2500)

††† {

††† cout << "Ќе корректировать!";

††† return;

††† }

† Tkor = t;

† real R_t = sqrt(f[0]*f[0]+f[1]*f[1]+f[2]*f[2]);

† real V_t = sqrt(f[3]*f[3]+f[4]*f[4]+f[5]*f[5]);

† real R_n = parn[0];

† if (Fl_a)

††† {

††† dRa = R_t-R_n;

††† dRp = par[2]*(1-par[1])-R_n;

††† cout << "јпоцентр dRp:" << dRp << "м n";

††† cout << "dRa:" << dRa << "м n";

††† cout << "w=" << par[5]*r_g << "u=" << par[7]*r_g << 'n';

††† real l,ln;

††† l = -(w_z-w_s)*par[6];

††† ln = -(w_z-w_s)*parn[6];

††† dl = -(w_z-w_s)*(par[6]-parn[6]);

††† cout << "T=" << par[6] << "“ном=" << parn[6] << " T-Tном="

††††††††† †<< par[6]-parn[6] << 'n' << "l=" << l*r_g << "lном="

††††††††† †<< ln*r_g << "l-lном=" << (l-ln)*r_g << "dl=" << dl

††††††††† †<< 'n';

††† if (dRp > 0)

††††† Sig_a = -1;

††† else

††††† Sig_a = 1;

††† cout << "«нак ускорени€:" << Sig_a << 'n';

††† clrscr();

††† real Rp = par[2]*(1-par[1]);

††† real Ra_p = par[2]*(1+par[1]);

††† real Rp_p2 = Rp;

††† real Ra_p2 = R_t;

††† cout << "Rp=" << Rp_p2 << "Ra=" << Ra_p2 << 'n';

††† cout << "Ra_p=" << Ra_p << "n Rt=" << R_t << 'n';

††† if (fabs(Rp - R_n) < 500)

††††† {

††††† Fl_kp = 0;

††††† Fl_ka = 1;

††††† cout << "«акончить коррекцию в апоцентре n" << "dRp=" << Rp-R_n

††††††††† †† << "dRa=" << dRa << "t=" << t << 'n';

††††† cout << "ѕараметры орбиты: n" << "Rp=" << par[2]*(1-par[1])

††††††††† †† << "Ra=" << par[2]*(1+par[1]) << "n p=" << par[0]

††††††††† †† << "a=" << par[2] << "e=" << par[1] << "n T="

††††††††† †† << par[6] << "w=" << par[5]*r_g << "u=" << par[7]*r_g

††††††††† †† << 'n';

††††† cout << "—уммарный импульс дл€ коррекции перицентра=" << dV_ps << 'n';

††††† clrscr();

††††† }

††† else

††††† {

††††† if (R_t > R_n)

††††††††† {

††††††††† Rp_p = R_n;

††††††††† Ra_p = R_t;

††††††††† a_p = (Ra_p+Rp_p)/2.;

††††††††† e_p = 1-Rp_p/a_p;

††††††††† p_p = a_p*(1-e_p*e_p);

††††††††† Vk = sqrt(mu_z/p_p)*(1-e_p);

††††††††† }

††††† else

††††††††† {

††††††††† Rp_p = R_t;

††††††††† Ra_p = R_n;

††††††††† a_p = (Ra_p+Rp_p)/2.;

††††††††† e_p = 1-Rp_p/a_p;

††††††††† p_p = a_p*(1-e_p*e_p);

††††††††† Vk = sqrt(mu_z/p_p)*(1+e_p);

††††††††† }

††††† real dV = Vk-V_t;

††††† real dVmax = 20*25./m;

††††† cout << "n dVтреб=" << dV << "dVmax за 20 сек=" << dVmax;

††††† if (fabs(dV) > dVmax)

††††††††† {

††††††††† akor[0] = Sig_a*(25./m)*f[3]/V_t;

††††††††† akor[1] = Sig_a*(25./m)*f[4]/V_t;

††††††††† akor[2] = Sig_a*(25./m)*f[5]/V_t;

††††††††† cout << "n dV=" << dV << "dVmax=" << dVmax;

††††††††† cout << "n  орректирующее ускорение:" << akor[0] << 't' << akor[1]

††††††††† †††† << 't' << akor[2] << 't' <<

††††††††† †††† sqrt(akor[0]*akor[0]+akor[1]*akor[1]+akor[2]*akor[2]) << 'n';

††††††††† dV_ps = dV_ps+dVmax;

††††††††† cout << "—уммарный импульс=" << dV_ps << 'n';

††††††††† }

††††† else

††††††††† {

††††††††† akor[0] = Sig_a*(fabs(dV)/dVmax)*(25./m)*f[3]/V_t;

††††††††† akor[1] = Sig_a*(fabs(dV)/dVmax)*(25./m)*f[4]/V_t;

††††††††† akor[2] = Sig_a*(fabs(dV)/dVmax)*(25./m)*f[5]/V_t;

††††††††† cout << "n dV=" << dV << "dVmax=" << dVmax;

††††††††† cout << "n  орректирующее ускорение:" << akor[0] << 't' << akor[1]

††††††††† †††† << 't' << akor[2] << 't' <<

††††††††† †††† sqrt(akor[0]*akor[0]+akor[1]*akor[1]+akor[2]*akor[2]) << 'n';

††††††††† dV_ps = dV_ps+fabs(dV);

††††††††† cout << "—уммарный импульс=" << dV_ps << 'n';

††††††††† }

††††† if (dVmax > fabs(dV))

††††††††† {

††††††††† dVmax = fabs(dV);

††††††††† real Vk_r = Sig_a*dVmax+V_t;

††††††††† real Ra_r = R_t;

††††††††† real e_r = -(Vk_r*Vk_r*Ra_r/mu_z)+1;

††††††††† real a_r = Ra_r/(1+e_r);

††††††††† real p_r = a_r*(1-e_r*e_r);

††††††††† real Rp_r = a_r*(1-e_r);

††††††††† cout << "ѕараметры орбиты: n" << " Rp_r = " << Rp_r

††††††††† †††† << " Ra_r = " << Ra_r << "n p_r = " << p_r << " a_r = "

††††††††† †††† << a_r << " e_r = " << e_r << 'n';

††††††††† }

††††† else

††††††††† {

††††††††† real Vk_r = Sig_a*dVmax+V_t;

††††††††† real Ra_r = R_t;

††††††††† real e_r = -(Vk_r*Vk_r*Ra_r/mu_z)+1;

††††††††† real a_r = Ra_r/(1+e_r);

††††††††† real p_r = a_r*(1-e_r*e_r);

††††††††† real Rp_r = a_r*(1-e_r);

††††††††† cout << "ѕараметры орбиты: n" << " Rp_r = " << Rp_r

††††††††† †††† << " Ra_r = " << Ra_r << "n p_r = " << p_r << " a_r = "

††††††††† †††† << a_r << " e_r = " << e_r << 'n';

††††††††† }

††††† T_vd = t;

††††† cout << "¬кл.дв. t=" << T_vd << 'n';

††††† }

††† }

† if (Fl_p)

††† {

††† dRp = R_t-R_n;

††† dRa = par[2]*(1+par[1])-R_n;

††† cout << " ѕерицентра - dRp:" << dRp << "м n";

††† cout << "dRa:" << dRa << "м. n";

††† cout << "w=" << par[5]*r_g << "u=" << par[7]*r_g << 'n';

††† real l,ln;

††† l = -(w_z-w_s)*par[6];

††† ln = -(w_z-w_s)*parn[6];

††† dl = -(w_z-w_s)*(par[6]-parn[6]);

††† cout << "T=" << par[6] << "Tном=" << parn[6] << "T-Tном="

††††††††† †<< par[6]-parn[6] << 'n' << "l=" << l*r_g << "lном="

††††††††† †<< ln*r_g << "l-lном=" << (l-ln)*r_g << "dl=" << dl << 'n';

††† if (dRa > 0)

††††† Sig_a = -1;

††† else

††††† Sig_a = 1;

††† cout << "«нак ускорени€:" << Sig_a << 'n';

††† clrscr();

††† real Ra = par[2]*(1+par[1]);

††† real Rp_p1 = R_t;

††† real Ra_p1 = Ra;

††† cout << "Rp=" << Rp_p1 << "Ra=" << Ra_p1 << 'n';

††† if ((fabs(Ra-R_n) < 500) || (fabs(dl*r_g) < .0001))

††††† {

††††† cout << "«акончить коррекцию в перицентре n" << "dRa="

††††††††† †† << Ra-R_n << "dRp=" << dRp << "t=" << t << 'n';

††††† cout << "ѕараметры орбиты: n " << "Rp="

††††††††† †† << par[2]*(1-par[1]) << "Ra=" << par[2]*(1+par[1])

††††††††† †† << " n p=" << par[0] << "a=" << par[2] << "e="

††††††††† †† << par[1] << " n T=" << par[6] << "w=" << par[5]*r_g

††††††††† †† << "u=" << par[7]*r_g << 'n';

††††† cout << "—уммарный импульс дл€ коррекции перицентра=" << dV_as << 'n';

††††† clrscr();

††††† Fl_ka = 0;

††††† Fl_ki = 1;

††††† }

††† else

††††† {

††††† if (R_t > R_n)

††††††††† {

††††††††† Rp_p = R_n;

††††††††† Ra_p = R_t;

††††††††† a_p = (Ra_p+Rp_p)/2.;

††††††††† e_p = 1-Rp_p/a_p;

††††††††† p_p = a_p*(1-e_p*e_p);

††††††††† Vk = sqrt(mu_z/p_p)*(1-e_p);

††††††††† }

††††† else

††††††††† {

††††††††† Rp_p = R_t;

††††††††† Ra_p = R_n;

††††††††† a_p = (Ra_p+Rp_p)/2.;

††††††††† e_p = 1-Rp_p/a_p;

††††††††† p_p = a_p*(1-e_p*e_p);

††††††††† Vk = sqrt(mu_z/p_p)*(1+e_p);

††††††††† }

††††† real dV = Vk-V_t;

††††† real dVmax = 20*25./m;

††††† cout << "n dVнадо=" << dV << " dVmax за 20 сек=" << dVmax;

††††† if (fabs(dV) > dVmax)

††††††††† {

††††††††† akor[0] = Sig_a*(25./m)*f[3]/V_t;

††††††††† akor[1] = Sig_a*(25./m)*f[4]/V_t;

††††††††† akor[2] = Sig_a*(25./m)*f[5]/V_t;

††††††††† cout << "n dV=" << dV << "dVmax=" << dVmax;

††††††††† cout << "n  орректирующее ускорение:" << akor[0] << 't' << akor[1]

††††††††† †††† << 't' << akor[2] << 't' <<

††††††††† †††† sqrt(akor[0]*akor[0]+akor[1]*akor[1]+akor[2]*akor[2]) << 'n';

††††††††† dV_as = dV_as+dVmax;

††††††††† cout << "—уммарный импульс=" << dV_as << 'n';

††††††††† }

††††† else

††††††††† {

††††††††† akor[0] = Sig_a*(fabs(dV)/dVmax)*(25./m)*f[3]/V_t;

††††††††† akor[1] = Sig_a*(fabs(dV)/dVmax)*(25./m)*f[4]/V_t;

††††††††† akor[2] = Sig_a*(fabs(dV)/dVmax)*(25./m)*f[5]/V_t;

††††††††† cout << "n dV=" << dV << " dVmax=" << dVmax;

††††††††† cout << "n  орректирующее ускорение:" << akor[0] << 't' << akor[1]

††††††††† †††† << 't' << akor[2] << 't' <<

††††††††† †††† sqrt(akor[0]*akor[0]+akor[1]*akor[1]+akor[2]*akor[2]) << 'n';

††††††††† dV_as = dV_as+fabs(dV);

††††††††† cout << "—уммарный импульс=" << dV_as << 'n';

††††††††† }

††††† if (dVmax > fabs(dV))

††††††††† {

††††††††† dVmax = fabs(dV);

††††††††† real Vk_r = Sig_a*dVmax+V_t;

††††††††† real Ra_r = R_t;

††††††††† real e_r = -(Vk_r*Vk_r*Ra_r/mu_z)+1;

††††††††† real a_r = Ra_r/(1+e_r);

††††††††† real p_r = a_r*(1-e_r*e_r);

††††††††† real Rp_r = a_r*(1-e_r);

††††††††† cout << "ѕараметры орбиты: n" << "Rp_r=" << Rp_r

††††††††† †††† << "Ra_r=" << Ra_r << "n p_r=" << p_r << "a_r="

††††††††† †††† << a_r << "e_r=" << e_r << 'n';

††††††††† }

††††† else

††††††††† {

††††††††† real Vk_r = Sig_a*dVmax+V_t;

††††††††† real Ra_r = R_t;

††††††††† real e_r = -(Vk_r*Vk_r*Ra_r/mu_z)+1;

††††††††† real a_r = Ra_r/(1+e_r);

††††††††† real p_r = a_r*(1-e_r*e_r);

††††††††† real Rp_r = a_r*(1-e_r);

††††††††† cout << "ѕараметры орбиты: n" << "Rp_r=" << Rp_r

††††††††† †††† << "Ra_r=" << Ra_r << "n p_r=" << p_r << "a_r="

††††††††† †††† << a_r << "e_r=" << e_r << 'n';

††††††††† }

††††† T_vd = t;

††††† cout << "¬кл.дв. t=" << T_vd << 'n';

††††† }

††† }

† if (Fl_lu)

††† {

††† real di = par[4]-parn[4];

††† cout << "Ћини€ узлов - di: " << di*r_g << "град n";

††† cout << "w=" << par[5]*r_g << "u=" << par[7]*r_g << 'n';

††† real l,ln;

††† l = -(w_z-w_s)*par[6];

††† ln = -(w_z-w_s)*parn[6];

††† dl = -(w_z-w_s)*(par[6]-parn[6]);

††† cout << "T=" << par[6] << "Tном=" << parn[6] << "T-Tном="

††††††††† †<< par[6]-parn[6] << 'n' << "l=" << l*r_g << "lном="

††††††††† †<< ln*r_g << "l-lном=" << (l-ln)*r_g << "dl=" << dl

††††††††† †<< "n i=" << par[4]*r_g << "iном=" << parn[4]*r_g << 'n';

††† cout << "ѕараметры орбиты: n " << "Rp="

††††††††† †<< par[2]*(1-par[1]) << "Ra=" << par[2]*(1+par[1])

††††††††† †<< " n p=" << par[0] << "a=" << par[2] << "e="

††††††††† †<< par[1] << " n T=" << par[6] << "w=" << par[5]*r_g

††††††††† †<< "u=" << par[7]*r_g << " n i=" << par[4]*r_g << 'n';

††† clrscr();

††† real Vk_x,Vk_y,Vk_z;

††† if (fabs(di) < .0001*g_r)

††††† {

††††† Fl_ki = 0;

††††† cout << "«акончить коррекцию наклонени€ n "

††††††††† †† << "di=" << (par[4]-parn[4])*r_g << "t=" << t << 'n';

††††† cout << "ѕараметры орбиты: n " << "Rp="

††††††††† †† << par[2]*(1-par[1]) << "Ra=" << par[2]*(1+par[1])

††††††††† †† << " n p=" << par[0] << "a=" << par[2] << "e="

††††††††† †† << par[1] << " n T=" << par[6] << "w=" << par[5]*r_g

††††††††† †† << "u=" << par[7]*r_g << " n i=" << par[4]*r_g << 'n';

††††† cout << "—уммарный импульс=" << dV_is

††††††††† †† << 'n';

††††† clrscr();

††††† }

††† else

††††† {

††††† real teta;

††† ††if (par[7] > par[5])

††††††††† teta = 2*M_PI+par[7]-par[5];

††††† else

††††††††† teta = par[7]-par[5];

††††† real Vr_i = sqrt(mu_z/par[0])*par[1]*sin(teta);

††††† real Vn_i = sqrt(mu_z/par[0])*(1+par[1]*cos(teta));

††††† V_t = sqrt(f[3]*f[3]+f[4]*f[4]+f[5]*f[5]);

††††† Vk_x = -Vn_i*cos(parn[4])*sin(par[3])+Vr_i*cos(par[3]);

††††† Vk_y = Vn_i*cos(parn[4])*cos(par[3])+Vr_i*sin(par[3]);

††††† Vk_z = Vn_i*sin(parn[4]);

††††† Vk = sqrt(Vk_x*Vk_x+Vk_y*Vk_y+Vk_z*Vk_z);

††††† real dV_x = Vk_x-f[3];

††††† real dV_y = Vk_y-f[4];

††††† real dV_z = Vk_z-f[5];

††††† real dV = sqrt(dV_x*dV_x+dV_y*dV_y+dV_z*dV_z);

††††† real dVmax = 20*25./m;

††††† cout << "Vнач=" << V_t << "Vк=" << Vk << "teta=" << teta*r_g

††††††††† †† << 'n';

††††† cout << "dV=" << dV << "dVmax за 20 сек=" << dVmax;

††††† if (dV > dVmax)

††††††††† {

††††††††† akor[0] = (25./m)*dV_x/dV;

††††††††† akor[1] = (25./m)*dV_y/dV;

††††††††† akor[2] = (25./m)*dV_z/dV;

††††††††† cout << "n  орректирующее ускорение:" << akor[0] << 't' << akor[1] <<

††††††††† ††††† 't' << akor[2] << 't' <<

††††††††† ††††† sqrt(akor[0]*akor[0]+akor[1]*akor[1]+akor[2]*akor[2]) << 'n';

††††††††† dV_is = dV_is+dVmax;

††††††††† cout << "—уммарный импульс=" << dV_is << 'n';

††††††††† }

††††† else

††††††††† {

††††††††† akor[0] = (fabs(dV)/dVmax)*(25./m)*dV_x/dV;

††††††††† akor[1] = (fabs(dV)/dVmax)*(25./m)*dV_y/dV;

††††††††† akor[2] = (fabs(dV)/dVmax)*(25./m)*dV_z/dV;

††††††††† cout << "n  орректирующее ускорение:" << akor[0] << 't' << akor[1]

††††††††† †††† << 't' << akor[2] << 't'<<

††††††††† †††† sqrt(akor[0]*akor[0]+akor[1]*akor[1]+akor[2]*akor[2]) << 'n';

††††††††† dV_is = dV_is+fabs(dV);

††††††††† cout << "—уммарный импульс=" << dV_is << 'n';

††††††††† }

††††† T_vd = t;

††††† cout << "¬кл.дв. t=" << T_vd << 'n';

††††† }

††† }

† if ((!Fl_ka) && (!Fl_kp) && (!Fl_ki))

††† {

††† cout << " оррекци€ окончена!" << 'n';

††† real m_t;

††† dV_ss = dV_ps+dV_as+dV_is;

††† m_t = m*(1-exp(-dV_ss/W));

††† cout << "ѕотребный импульс: n - перицентра dV_ps="

††††††††† †<< dV_ps << "n апоцентра dV_as=" << dV_as

††††††††† †<< "n —уммарный импульс=" << dV_ss << "ћасса топлива=" << m_t

††††††††† †<< 'n';

††† dV_ps = 0;

††† dV_as = 0;

††† dV_is = 0;

††† dV_ss = 0;

††† m_t = 0;

††† }

† }

}

void par_or(real *f, real *par)

{

real x = f[0];

real y = f[1];

real z = f[2];

real Vx = f[3];

real Vy = f[4];

real Vz = f[5];

real c1 = (y*Vz-z*Vy);

real c2 = (z*Vx-x*Vz);

real c3 = (x*Vy-y*Vx);

real C = sqrt(c1*c1+c2*c2+c3*c3);

par[0] = (C/mu_z)*C;

real R_ka = sqrt(x*x+y*y+z*z);

real V_ka = sqrt(Vx*Vx+Vy*Vy+Vz*Vz);

real f1 = (Vy*c3-Vz*c2)-(mu_z*x/R_ka);

real f2 = (Vz*c1-Vx*c3)-(mu_z*y/R_ka);

real f3 = (Vx*c2-Vy*c1)-(mu_z*z/R_ka);

real F = sqrt(f1*f1+f2*f2+f3*f3);

real h = V_ka*V_ka-(2*mu_z/R_ka);

if ((1+h*C*C/(mu_z*mu_z)) < 0)

† {

† cout << " ќшибка! n";

† getch();

† }

par[1] = F/mu_z;

if ((1-par[1]*par[1]) < 0)

† {

† cout << " (1-e*e) < 0 ќшибка! n";

† getch();

† }

par[2] = par[0]/(1-par[1]*par[1]);

par[4] = acos(c3/C);

real s_Om = c1/(C*sin(par[4]));

real c_Om = -c2/(C*sin(par[4]));

if (s_Om >= 0)

† par[3] = acos(c_Om);

else

† par[3] = 2*M_PI-acos(c_Om);

real c_om = (f1*cos(par[3])+f2*sin(par[3]))/F;

real s_om = f3/(F*sin(par[4]));

if (s_om > 0)

† par[5] = acos(c_om);

else

† par[5] = 2*M_PI - acos(c_om);

if (par[2] < 0)

† {

† cout << " ќшибка! n";

† getch();

† }

par[6] = 2*M_PI*sqrt((par[2]/mu_z)*par[2]*par[2]);

real c_u = (x*cos(par[3])+y*sin(par[3]))/R_ka;

real s_u = z/(R_ka*sin(par[4]));

if (s_u > 0)

† par[7] = acos(c_u);

else

† par[7] = 2*M_PI - acos(c_u);

}

#include "rk5.h"

#include

void Drkgs(real *prmt,real *y,real *dery,int ndim,int& ihlf,

††† void (*fct)(real &,real*,real*),

††† void (*out_p)(real,real*,real*,int,int,real*))

{

static real a[] = { 0.5, 0.292893218811345248, 1.70710678118665475,

†††††††††††††††††† ††† 0.16666666666666667 };

static real b[] = { 2.0, 1.0, 1.0, 2.0 };

static real c[] = { 0.5, 0.292893218811345248, 1.70710678118665475, 0.5 };

real *aux[8];

int i,j,imod,itest,irec,istep,iend;

real delt,aj,bj,cj,r,r1,r2,x,xend,h;

for (i=0; i<8; i++) aux[i] = new real[ndim];

for (i=0; i

x = prmt[0];

xend = prmt[1];

h = prmt[2];

prmt[4] = 0.0;

fct(x,y,dery);

r = h*(xend-x);

if (r <= 0.0)

† {

† ihlf = 13;

† if (r == 0.0) ihlf = 12;

† goto l39;

† }

for(i=0; i

† {

† aux[0][i] = y[i];

† aux[1][i] = dery[i];

† aux[2][i] = 0.0;

† aux[5][i] = 0.0;

† }

irec = 0;

h = h+h;

ihlf = -1;

istep = 0;

iend = 0;

l4: r = (x+h-xend)*h;

if (r >= 0.0)

† {

† iend = 1;

† if (r > 0.0) h = xend-x;

† }

out_p(x,y,dery,irec,ndim,prmt);

if (prmt[4] != 0.0) goto l40;

itest = 0;

l9: istep++;

j = 0;

l10: aj = a[j];

bj =b[j];

cj = c[j];

for (i=0; i

† {

† r1 = h*dery[i];

† r2 = aj*(r1-bj*aux[5][i]);

† y[i] = y[i]+r2;

† r2 = r2+r2+r2;

† aux[5][i] += r2-cj*r1;

† }

if (j-3 < 0)

† {

† j++;

† if (j-2 != 0) x = x+0.5*h;

† fct(x,y,dery);

† goto l10;

† }

if (itest <= 0)

† {

† for (i=0; i

† itest = 1;

† istep = istep+istep-2;

l18: ihlf++;

† x = x-h;

† h = 0.5*h;

† for (i=0; i

††† {

††† y[i] = aux[0][i];

††† dery[i] = aux[1][i];

††† aux[5][i] = aux[2][i];

††† }

† goto l9;

† }

imod = istep/2;

if (istep-imod-imod != 0)

† {

† fct(x,y,dery);

† for (i=0; i

††† {

††† aux[4][i] = y[i];

††† aux[6][i] = dery[i];

††† }

† goto l9;

† }

delt = 0.0;

for (i=0; i

delt += aux[7][i]*fabs(aux[3][i]-y[i]);

if (delt-prmt[3] > 0.0)

† {

† if (ihlf-10 >= 0)

††† {

††† ihlf = 11;

††† fct(x,y,dery);

††† goto l39;

††† }

† for (i=0; i

† istep = istep+istep-4;

† x = x-h;

† iend = 0;

† goto l18;

† }

fct(x,y,dery);

for (i=0; i

† {

† aux[0][i] = y[i];

† aux[1][i] = dery[i];

† aux[2][i] = aux[5][i];

† y[i] = aux[4][i];

† dery[i] = aux[6][i];

† }

out_p(x-h,y,dery,ihlf,ndim,prmt);

if (prmt[4] != 0) goto l40;

for (i=0; i

† {

† y[i] = aux[0][i];

† dery[i] = aux[1][i];

† }

irec = ihlf;

if (iend > 0) goto l39;

ihlf--;

istep = istep/2;

h = h+h;

if (ihlf < 0) goto l4;

imod = istep/2;

if ((istep-2*imod != 0) || (delt-0.02*prmt[3] > 0.0)) goto l4;

ihlf--;

istep = istep/2;

h = h+h;

goto l4;

l39: out_p(x,y,dery,ihlf,ndim,prmt);

l40: for (i=0; i

return;

}

6.3. ‘ј…Ћ Ќј„јЋ№Ќќ… »Ќ»÷»јЋ»«ј÷»» INIT.H

ifndef _INIT

#define _INIT

#include "def.h"

#include

#include

ifstream if_init;

void nex_ln (void);

void init_m()

{

Np = 150;

t_beg = 0;

t_end = 8000000;

dt = 2;

toler = .05;

dTp = (t_end-t_beg)/float(Np);

Curp = 0;

J1 = 532;

J2 = 563;

J3 = 697;

m = 597.;

W = 2200;

mu_z = 3.9858e14;

mu_s = 1.3249e20;

mu_l = 4.9027e12;

w_s = 2*M_PI/(365.2422*24*3600);

w_z = 2*M_PI/(24*3600);

w_l = 2*M_PI/(27.32*24*3600);

ww_l = 2*M_PI/(18.6*365.2422*24*3600);

parn[0] = 6952137.;

parn[1] = 0;

parn[2] = 6952137;

parn[3] = 28.1*g_r;

parn[4] = 97.6*g_r;

parn[5] = 63.1968*g_r;

parn[6] = 5769.;

parn[7] = 5.751*g_r;

Fl_u = 1;

u_last = parn[7];

Fl_ka = 0;

Fl_kp = 0;

Fl_ki = 0;

Fl_p = 0;

Fl_a = 0;

Fl_i = 0;

Fl_pkT = 0;

Tkor = 0;

T_vd = 0;

akor[0] = 0;

akor[1] = 0;

akor[2] = 0;

dV_ps = 0;

dV_as = 0;

dV_is = 0;

dV_ss = 0;

Fl_l0 = 0;

Fl_l1 = 0;

Fl_pki = 0;

real x0 = 6137262.9+7000;

real y0 = 3171846.1+7000;

real z0 = 689506.95+7000;

real Vx0 = -201.288+5;

real Vy0 = -1247.027+5;

real Vz0 = 7472.65+5;

prmt[0] = t_beg;

prmt[1] = t_end;

prmt[2] = dt;

prmt[3] = toler;

prmt[4] = 0.0;

y_main[0] = x0;

y_main[1] = y0;

y_main[2] = z0;

y_main[3] = Vx0;

y_main[4] = Vy0;

y_main[5] = Vz0;

}

void nex_ln (void)

{

char ch;

if_init.get(ch);

while (ch != 'n')

† if_init.get(ch);

}

#endif

6.4 ‘ј…Ћ ќѕ»—јЌ»я ѕ≈–≈ћ≈ЌЌџ’ DEF.H

#ifndef _DEFH

#define _DEFH

#include

typedef long double real;

extern const float g_r;

extern const float r_g;

extern int Np;

extern int Curp;

extern real dTp;

extern real t_beg;

extern real t_end;

extern real dt;

extern real toler;

extern real J1,J2,J3;

extern real mu_z;

extern real mu_s;

extern real mu_l;

extern real m;

extern real m_t;

extern real W;

extern real w_s;

extern real w_z;

extern real w_l;

extern real ww_l;

extern real xs,ys,zs;

extern real xl,yl,zl;

extern real Fz,Fs,Fl,Fa,U20;

extern int nomin;

extern real par[8];

extern real parn[8];

extern real a_p,e_p,p_p,Om_p,i_p,om_p,Rp_p,Ra_p;

extern real y_main[6];

extern real prmt[5];

extern int Fl_u;

extern real u_last;

extern int Fl_ka;

extern int Fl_kp;

extern int Fl_ki;

extern int Fl_i;

extern int Fl_p;

extern int Fl_a;

extern int Fl_lu;

extern int Fl_pkT;

extern real dl;

extern real T_vd;

extern real dRa;

extern real dRp;

extern int Sig;

extern int Sig_a;

extern real Vkor[3];

extern real akor[3];

extern real Tkor;

extern real Tkore;

extern real dV_ps;

extern real dV_as;

extern real dV_is;

extern real dV_ss;

extern int Fl_l0;

extern int Fl_l1;

extern int Fl_pki;

#endif

6.5 ‘ј…Ћ SFUN.H

#ifndef _SFUN

#define _SFUN

#include "def.h"

#include

#include

#include

void out_p(real x,real *y,real*,int,int,real *);

real interpl(real*,real*,int,real);

void fct(real& ,real *y,real *dery);

void par_or(real *,real *);

#endif

6.5 ‘ј…Ћ RK5.H

#ifndef _RK5

#define _RK5

#include "def.h"

#include

#include

#include "sfun.h"

void Drkgs(real *prmt,real *y,real *dery,int ndim,int& ihlf,

† void (*fct)(real&,real*,real*),

† void (*out_p)(real,real*,real*,int,int,real*));

#endif

6.6 ѕ–ќ√–јћћј ѕќ—“–ќ≈Ќ»я ¬–≈ћ≈ЌЌџ’ ƒ»ј√–јћћ

clc

g_r = pi/180;

r_g = 180/pi;

load m_y.dat

t = m_y(:,1);

x = m_y(:,2);

y = m_y(:,3);

z = m_y(:,4);

Vx = m_y(:,5);

Vy = m_y(:,6);

Vz = m_y(:,7);

clear m_y;

s_tmp = size(t);

s_m = s_tmp(1);

clear s_tmp;

load m_f.dat

Fz = m_f(:,2);

Fs = m_f(:,3);

Fl = m_f(:,4);

Fa = m_f(:,5);

U20 = m_f(:,6);

clear m_f;

load m_s.dat

xs = m_s(:,2);

ys = m_s(:,3);

zs = m_s(:,4);

clear m_s;

load m_par.dat

p = m_par(:,2);

e = m_par(:,3);

a = m_par(:,4);

Om = m_par(:,5);

i = m_par(:,6);

omg = m_par(:,7);

T = m_par(:,8);

u = m_par(:,9);

clear m_par;

p_n = 6952137.;

e_n = 0;

a_n = 6952137.;

Om_n0 = 28.1*g_r;

i_n = 97.6*g_r;

omg_n = 346.725*g_r;

T_n = 5765;

ws = 2*pi/(365.2422*24*3600);

for j = 1:s_m, tmp(j) = Om_n0+ws*t(j);

end

Om_n = tmp';

clear tmp;

map = [1,1,1];

colormap(map);

plot(t,p,'y-',[min(t) max(t)],[p_n p_n],'r-'), title (' ‘окальный параметр '), grid on;

print -dwin;

pause;

plot(t,p-p_n,'y-'), title (' dp '), grid on;

print -dwin;

pause;

plot(t,e,'y-',[min(t) max(t)],[e_n e_n],'r-'), title (' Ёксцентриситет '), grid on;

print -dwin;

pause;

plot(t,e-e_n,'y-'), title (' de '), grid on;

print -dwin;

pause;

plot(t,a,'y-',[min(t) max(t)],[a_n a_n],'r-'), title (' Ѕольша€ полуось орбиты '), grid on;

print -dwin;

pause;

plot(t,a-a_n,'y-'), title (' da '), grid on;

print -dwin;

pause;

plot(t,Om*r_g,'y-',t,Om_n*r_g,'r-'), title (' ƒолгота восход€щего узла '), grid on;

print -dwin;

pause;

plot(t,Om*r_g-Om_n*r_g,'y-'), title (' dOm '), grid on;

print -dwin;

pause;

plot(t,i*r_g,'y-',[min(t) max(t)],[i_n*r_g i_n*r_g],'r-'), title (' Ќаклонение '), grid on;

print -dwin;

pause;

plot(t,i*r_g-i_n*r_g,'y-'), title (' di '), grid on;

print -dwin;

pause;

plot(t,T,'y-',[min(t) max(t)],[T_n T_n], 'r-'), title (' ѕериод '), grid on;

print -dwin;

pause;

plot(t,T-T_n,'y-'), title (' dT '), grid on;

print -dwin;

pause;

plot3(x,y,z,'b')

axis([min(x) max(x) min(y) max(y) min(z) max(z)])

set(gca,'box','on')

title (' ѕоложение ћ ј ')

hold on

plt = plot3(0,0,0,'.','erasemode','xor','markersize',24);

dk = ceil(length(y)/2500);

for k = 1:dk:length(y)

set(plt,'xdata',x(k),'ydata',y(k),'zdata',z(k))

drawnow

end

hold off

pause;

plot(t,Fz,'y-'), title (' √равитаци€ «емли ' ), grid on;

print -dwin;

pause;

plot(t,Fs,'y-'), title (' √равитаци€ —олнца и солнечное давление '), grid on;

print -dwin;

pause;

plot(t,Fl,'y-'), title (' √равитаци€ Ћуны '), grid on;

print -dwin;

pause;

plot(t,Fa,'y-'), title (' —опротивление атмосферы '), grid on;

print -dwin;

pause;

plot(t,U20,'y-'), title (' Ќецентральность гравитационного пол€ «емли '), grid on;

print -dwin;

pause;

plot(t,Fz+Fs+Fl+Fa+U20,'y-'), title (' —уммарное возмущающее ускорение '), grid on;

print -dwin;

pause;

clear all

clc

g_r = pi/180;

r_g = 180/pi;

p_n = 6952137.;

e_n = 0;

a_n = 6952137.;

Om_n0 = 28.1*g_r;

i_n = 97.6*g_r;

omg_n = 346.725*g_r;

T_n = 5765;

load u_par.dat

t_u = u_par(:,1);

p_u = u_par(:,2);

e_u = u_par(:,3);

a_u = u_par(:,4);

Om_u = u_par(:,5);

i_u = u_par(:,6);

omg_u = u_par(:,7);

T_u = u_par(:,8);

u_u = u_par(:,9);

clear u_par;

load u_f.dat;

Fz_u = u_f(:,2);

Fs_u = u_f(:,3);

Fl_u = u_f(:,4);

Fa_u = u_f(:,5);

U20_u = u_f(:,6);

clear u_f;

s_tmp = size(t_u);

s_m_u = s_tmp(1);

clear s_tmp;

ws = 2*pi/(365.2422*24*3600);

for j = 1:s_m_u, tmp(j) = Om_n0+ws*t_u(j);

end

Om_n_u = tmp';

clear tmp;

plot(t_u,p_u,'y-',[min(t_u) max(t_u)],[p_n p_n],'r-'), title (' ‘окальный параметр '), grid on;

print -dwin;

pause;

plot(t_u,p_u-p_n,'y-'), title (' dp '), grid on;

print -dwin;

pause;

plot(t_u,e_u,'y-',[min(t_u) max(t_u)],[e_n e_n],'r-'), title (' Ёксцентриситет '), grid on;

print -dwin;

pause;

plot(t_u,e_u-e_n,'y-'), title (' de '), grid on;

print -dwin;

pause;

plot(t_u,a_u,'y-',[min(t_u) max(t_u)],[a_n a_n],'r-'), title (' Ѕольша€ полуось орбиты '), grid on;

print -dwin;

pause;

plot(t_u,a_u-a_n,'y-'), title (' da '), grid on;

print -dwin;

pause;

plot(t_u,Om_u*r_g,'y-',t_u,Om_n_u*r_g,'r-'), title (' ƒолгота восход€щего узла '), grid on;

print -dwin;

pause;

plot(t_u,Om_u*r_g-Om_n_u*r_g,'y-'), title (' dOm '), grid on;

print -dwin;

pause;

plot(t_u,i_u*r_g,'y-',[min(t_u) max(t_u)],[i_n*r_g i_n*r_g],'r-'), title (' Ќаклонение '), grid on;

print -dwin;

pause;

plot(t_u,i_u*r_g-i_n*r_g,'y-'), title (' di '), grid on;

print -dwin;

pause;

plot(t_u,T_u,'y-',[min(t_u) max(t_u)],[T_n T_n], 'r-'), title (' ѕериод '), grid on;

print -dwin;

pause;

plot(t_u,T_u-T_n,'y-'), title (' dT '), grid on;

print -dwin;

pause;

clear all

†TOC o "1-3" 1. ќглавление.................................................................................. 1 2. »сследовательска€ часть................................................... 3 2.1. ¬ведение............................................

 

 

 

¬нимание! ѕредставленный ƒиплом находитс€ в открытом доступе в сети »нтернет, и уже неоднократно сдавалс€, возможно, даже в твоем учебном заведении.
—оветуем не рисковать. ”знай, сколько стоит абсолютно уникальный ƒиплом по твоей теме:

Ќовости образовани€ и науки

«аказать уникальную работу

—вои сданные студенческие работы

присылайте нам на e-mail

Client@Stud-Baza.ru