Форум программистов, компьютерный форум, киберфорум
Наши страницы

Кусовая на С++ "Системы управления движением космических аппаратов" - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Подсчитать сумму элементов, кратных 9 http://www.cyberforum.ru/cpp-beginners/thread1162562.html
Подсчитать сумму элементов, кратных 9
C++ Запись в конец файла Создать текстовый файл message.txt (например в редакторе «Блокнот»), содержащий текстовое сообщение, состоящее из нескольких строк текста. Необходимо программно обратиться к существующему файлу и... http://www.cyberforum.ru/cpp-beginners/thread1162542.html
Указать параметром шаблона статическую матрицу C++
int matrix; List<int**> list; list << matrix; no match <int**> and <int(*)> Как правильно указать параметр шаблона тогда?
C++ Прием MPI_Bcast
При передаче сообщений с помощью функции MPI_Send, принять процессом-получателем можно с помощью функции MPI_Recv, а как принимать сообщения, отправленные MPI_Bcast?
C++ Написать три функции, которые будут вызываться из функции main() http://www.cyberforum.ru/cpp-beginners/thread1162451.html
Доброго времени суток... уважаемые форумчане может кто поможет с моей бедой. а беда вроде проста но мне она не понятна... а вот и сама беда Во всех вариантах необходимо написать три функции, которые...
C++ Напечатать все простые числа, не превосходящие заданное число M Люди добрые помогите разобраться вот с этим.. может кто поможет с решением столь не легкой задачи. Напечатать все простые числа, не превосходящие заданное число M. Для ускорения вычислений полезно... подробнее

Показать сообщение отдельно
valentina123423
0 / 0 / 0
Регистрация: 11.11.2013
Сообщений: 4
03.05.2014, 17:21  [ТС]
В архиве это и есть все задание, полностью.

Добавлено через 4 минуты
Исходные данные
Режим управления:

Система демпфирования остаточных угловых скоростей КА, возникающих после отделения от РБ.

Начальные условия моделирования:


=


Тензор инерции:



Исполнительные органы: реактивные двигатели малой тяги.

Измерительные приборы: магнитометры.

Функция GetMagDat выглядит следующим образом:

void GetMagData(double t, double *Q, double *b);

Входными переменными функции являются:
t – время, прошедшее от момента отделения КА,
Q – кватернион, характеризующий положение КА относительно J2000, определяемый из уравнения (2).
Выходной переменной является b – вектор магнитной индукции, в проекции на оси связанной системы координат КА, измеряемый в нанотеслах (1 нТл = .
В КР используется DLL в среде разработки Borland Builder v 6.0.
.
При моделировании работы полагается, что магнитометр измеряет магнитное поле без ошибок.
Функция GetMagData:
1. Осуществляет решение уравнений орбитального движения КА методом конических сечений для Кеплеровой невозмущенной орбиты, без учета нецентральности гравитационного поля Земли и влияния планет.
2. Осуществляет вычисление магнитного поля Земли согласно международной модели IGRF-11.

Постановка задачи

Смоделируем управляемое движение КА применяя режим демпфирования остаточных угловых скоростей КА возникающих после отделения от РБ

Для этого выполним следующую последовательность действий:

1) Проинтегрируем уравнение динамики Эйлера

(1)
Где:
– вектор угловой скорости вращения КА в осях связанной системы координат (ССК),
где


– тензор инерции космического аппарата,



- суммарный момент сил, как внешних возмущающих, так и управляющих. Где М=0.

Уравнение (1) справедливо для случая, когда исполнительными органами КА не являются гироскопические исполнительные органы (двигатели-маховики, гиродины и проч.), а органы управления, работа которых не влияет на кинетический момент КА как замкнутой системы – магнитные исполнительные органы, реактивные двигатели стабилизации и проч.

2) Вычислим

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


Где

В результате вычисления получим

3) Проинтегрируем кинематическое соотношение Пуассона

(2)
Где,
– кватернион ориентации КА, имеющий форму:
, – единичный вектор ори поворота.
– угол поворота, – вектор угловой скорости вращения КА, определяемый из уравнений (1). Оператор в данном случае можно интерпретировать следующим образом, переписав уравнение (2):


Где =

Интегрирование (1) вместе с уравнением (2) позволяют полностью описать вращательное движение КА как твердого тела относительно некоторой инерциальной СК.

В результате вычисления получим

4) Вычислим единичные вектора вектора магнитной индукции в проекции на оси ССК КА в моменты времени и , используя GetMagData.

Полагая, что достаточно мало для того, что бы изменение проекции вектора магнитного поля на оси ССК за счет собственного вращения КА было существенно больше, чем его изменение за счет орбитального движения, воспользуемся соотношениями:


Где
– приращение угла поворота за счет вращения КА по информации магнитометров,
– орт вектора поворота.

В процессе вычисления получим .

5) В результате интегрирования получим графики зависимостей:
, а так же и .





Формулы посмотрите в файле МОСКОВСКИЙ АВИАЦИОНЫЙ ИНСТИТУТ в архиве
0
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru