русс | укр

Языки программирования

ПаскальСиАссемблерJavaMatlabPhpHtmlJavaScriptCSSC#DelphiТурбо Пролог

Компьютерные сетиСистемное программное обеспечениеИнформационные технологииПрограммирование

Все о программировании


Linux Unix Алгоритмические языки Аналоговые и гибридные вычислительные устройства Архитектура микроконтроллеров Введение в разработку распределенных информационных систем Введение в численные методы Дискретная математика Информационное обслуживание пользователей Информация и моделирование в управлении производством Компьютерная графика Математическое и компьютерное моделирование Моделирование Нейрокомпьютеры Проектирование программ диагностики компьютерных систем и сетей Проектирование системных программ Системы счисления Теория статистики Теория оптимизации Уроки AutoCAD 3D Уроки базы данных Access Уроки Orcad Цифровые автоматы Шпаргалки по компьютеру Шпаргалки по программированию Экспертные системы Элементы теории информации

Разработка математической модели и её решение с использованием метода пространства состояний


Дата добавления: 2014-11-27; просмотров: 768; Нарушение авторских прав


 

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

1) Входные переменные, характеризующие внешнее воздействие на входы системы.

2) Переменные состояния - внутренние (промежуточные) переменные, совокупность которых полностью характеризует свойства системы.

3) Выходные переменные, представляющие реакцию системы на внешние воздействия и те состояния системы, которые представляют интерес для исследователя.

Собственно система, её входы и выходы - это три взаимосвязанных объекта, которые в каждом конкретном случае однозначно описывают систему. В зависимости от того, какой из объектов подлежит определению при остальных двух заданных различают три типа задач исследования проектирования: анализ, синтез и измерения. Решение любой из этих задач связано с исследованием состояний системы, множество которых образует пространство состояний.

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

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

Непрерывную детерминированную систему в каждый момент времени можно записать парой матричных уравнений:

1) Уравнение состояния системы:

 



 



2) Уравнение выхода:

 



 



Эти уравнения в зависимости от свойств системы принимают специфическую форму. Для линейных систем уравнения преобразуются так:

3) Уравнение состояния системы:

 



 



4) Уравнение выхода:

 



Получим уравнения состояния и её выходов для устройства, показанного на рис. 2.1.

 



 
 

 



 




Рисунок 2.1 - Структурная схема устройства

 



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

(2.1)

 



Из третьего и четвёртого уравнений системы (2.1) выражаем токи I1(p) и I2(p) через ток I3(p):

 



(2.2)

 



(2.3)

 



Подставляем найденные токи во второе уравнение системы:

 



(2.4)

 



Теперь можем получить передаточную функцию:

 



(2.5)

 



Перейдём от передаточной функции W(p) (2.5) к дифференциальному уравнению:

 



(2.6)

 



откуда получаем дифференциальное уравнение:

 



или

 



(2.7)

 



где

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

Уравнение (2.7) является дифференциальным уравнением третьего порядка. Приведём его к уравнениям первого порядка и решим их систему.

Выражаем дифференциальное уравнение (2.7) относительно старшей производной:

 



(2.8)

 



 



Осуществляем цепочку замен:

Пусть ,

тогда

 



I
Формируем структурную схему, где операцию интегрирования обозначим с помощью интегратора .

 



 
 

 




Рисунок 2.2 - Структурная схема алгоритма

 



 



Запишем матрицу коэффициентов переменных состояний:

 



(2.11)

 



На следующем этапе анализа системы составляем строки для подпрограммы, реализующей метод Рунге-Кутта, осуществляем запуск программы и получаем результат в виде числового и графического материала.

Для анализа системы зададимся в исходном случае следующими значениями сопротивления и ёмкости:

R=10Ом; С=0,02Ф.

 



Тогда коэффициенты в матрице (2.11) будут иметь следующие значения:

A0=1(Ом×Ф); A1=7(Ом×Ф); A2=4,96(Ом×Ф)2; A3=0,512(Ом×Ф).

 



Составляем строки для подпрограммы:

 



500 F(1)=H*Y(2)

510 F(2)=H*Y(3)

520 F(3)=H*(-1/A3*Y(1)-A1/A3*Y(2)-A2/A3*Y(3)+1/A3)

 



Осуществляем запуск программы RUNKUT.BAS (приложение 2), в режиме диалога вводим следующие значения:

 



МЕТОД РУНГЕ-КУТТА ДЛЯ N УРАВНЕНИЙ

НАЧ.И КОН.ЗНАЧЕНИЕ АРГУМЕНТА(X,XK)? 0,50

КОЛИЧЕСТВО ФУНКЦИЙ N? 3

ВВЕДИ КОЛИЧЕСТВО ТОЧЕК М? 800

ЧЕРЕЗ СКОЛЬКО ТОЧЕК ВЫВОДИТЬ НА ЭКРАН ?? 5

НАЧ.ЗНАЧЕНИЯ ФУНКЦИЙ

Y( 1 )=? 0

Y( 2 )=? 0

Y( 3 )=? 0

 



В результате получаем решение (приложение 3, а).

 



 



Изменим сопротивление R , ёмкость C:

R=1Ом; С=0,02Ф.

 



Тогда коэффициенты в матрице (2.11) будут иметь следующие значения:

A0=1(Ом×Ф); A1=0,7 (Ом×Ф); A2=0,0496(Ом×Ф)2; A3=0,064(Ом×Ф).

 



Составляем строки для подпрограммы:

 



500 F(1)=H*Y(2)

510 F(2)=H*Y(3)

520 F(3)=H*(-1/A3*Y(1)-A1/A3*Y(2)-A2/A3*Y(3)+1/A3)

 



Осуществляем запуск программы RUNKUT.BAS (приложение 2), в режиме диалога вводим следующие значения:

 



МЕТОД РУНГЕ-КУТТА ДЛЯ N УРАВНЕНИЙ

НАЧ.И КОН.ЗНАЧЕНИЕ АРГУМЕНТА(X,XK)? 0,50

КОЛИЧЕСТВО ФУНКЦИЙ N? 3

ВВЕДИ КОЛИЧЕСТВО ТОЧЕК М? 800

ЧЕРЕЗ СКОЛЬКО ТОЧЕК ВЫВОДИТЬ НА ЭКРАН ?? 5

НАЧ.ЗНАЧЕНИЯ ФУНКЦИЙ

Y( 1 )=? 0

Y( 2 )=? 0

Y( 3 )=? 0

 



В результате получаем решение (приложение 3, б).

 






<== предыдущая лекция | следующая лекция ==>
Что такое компьютер ? | Метод контурных токов


Карта сайта Карта сайта укр


Уроки php mysql Программирование

Онлайн система счисления Калькулятор онлайн обычный Инженерный калькулятор онлайн Замена русских букв на английские для вебмастеров Замена русских букв на английские

Аппаратное и программное обеспечение Графика и компьютерная сфера Интегрированная геоинформационная система Интернет Компьютер Комплектующие компьютера Лекции Методы и средства измерений неэлектрических величин Обслуживание компьютерных и периферийных устройств Операционные системы Параллельное программирование Проектирование электронных средств Периферийные устройства Полезные ресурсы для программистов Программы для программистов Статьи для программистов Cтруктура и организация данных


 


Не нашли то, что искали? Google вам в помощь!

 
 

© life-prog.ru При использовании материалов прямая ссылка на сайт обязательна.

Генерация страницы за: 1.047 сек.