149614 (731939), страница 3

Файл №731939 149614 (Математическое моделирование физических задач на ЭВМ) 3 страница149614 (731939) страница 32016-08-01СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 3)

(3.1)

Так как каждая ветвь соединяет два узла—выходит из одного узла и входит в другой, то столбец матрицы состоит из двух ненулевых элементов +1, —1 (их сумма равна нулю), так что достаточно заполнить таблицу для ny-1 узлов, которая является редуцированной матрицей соединений А. Эту незави­симую матрицу можно получить из полной матрицы Аa вычеркиванием строки, соответствующей выбранному базисно­му узлу.

Приняв в качестве базисного узел 4 и соответственно вычеркивая четвертую строку в (3.1), получим редуцированную матрицу соединений:

(3.2)

Строка матрицы А показывает, какие ветви выходят из каждого независимого узла графа цепи (и входят в него), а столбец – к каким узлам присоединена ветвь.

В отличие от полной матрицы Аа у реду­цированной матрицы соединений связного графа множест­во всех строк линейно независимо. Отсюда можно сделать вывод о том, что система уравнений равновесия токов в ny-1 узлах цепи линейно независима. Если ввести вектор токов пв, ветвей:

i=[i1, i2, … inв]Т, (3.3)

то систему независимых уравнений в nу-1 узлах по ЗТК в соответствии со смыслом матрицы А можно записать в виде:

Ai=0, (3.4)

где 0=[0 0 … 0]T - нулевой вектор размерности nу-1.

Для графа цепи рис. 3.1, б с матрицей соединений (3.2) имеем:

Транспонированная матрица соединений имеет вид:

(3.5)

Строка этой матрицы показывает, между какими узлами присоединена каждая ветвь.

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

2. Уравнения контурных токов

Метод контурных токов применим к цепям с планарным графом (рис. 3.3, а). В качестве переменных принимают замкну­тые контурные токи, проходящие по ветвям, образующим все внутренние ячейки графа.

Если намечать контуры – периметры ячеек – по порядку, начиная с одного края цепи, то легко убедиться, что в каждый последующий контур вносится новая ветвь, не вошедшая в предыдущие контуры. Отсюда следует, что уравнения равно­весия напряжений в таких контурах будут линейно независимы. Число внутренних ячеек равно nx=nв-ny+1.

Каждой ячейке при­писывается один кон­турный ток, замыкаю­щийся по ветвям, об­разующим ячейку. Об­щее число перемен­ных – контурных токов – равно числу ячеек. Направления всех контурных токов принимают одинаковыми – по часовой стрел­ке. Как видно из рис. 3.3, a, по каждой ветви цепи, за исключением периферийных ветвей, замыкаются два контурных тока, направленные в противоположные стороны. Запишем уравнения соединений.

  1. Приравнивая нулю суммы напряжений ветвей всех неза­висимых контуров (ячеек), имеем nх уравнений по ЗНК .

  2. Выражая ток каждой ветви через разность двух (в общем случае) замыкающихся по смежным контурам контурных токов, получим nв уравнений по ЗТК:

.

Как видно, токи всех ветвей, т. е. поведение всей цепи, полностью определяются nх контурными токами, число кото­рых меньше числа ветвей.

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

.

Для получения уравнений относительно выбранных перемен­ных необходимо:

  1. с помощью уравнений ветвей в уравнениях равновесия напряжений заменить напряжения всех ветвей токами;

  2. токи ветвей в получившейся системе заменить, согласно , контурными токами.

Получим уравнения для одной из ячеек, например первой (рис. 3.3, в), образованной тремя ветвями. Основным уравнени­ем равновесия напряжений в первом контуре будет:

u1+u2+u3=0 (*)

Токи ветвей ячейки:

.(**)

Уравнения ветвей:

(***)

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

,

После группировки имеем:

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

Аналогичные уравнения получим для остальных контуров. Если число контуров равно п, то предположив для общности число ветвей каждой ячейки также равным п, можно записать систему уравнений контурных токов:

Коэффициент Rkkсобственное сопротивление контура, рав­ное сумме сопротивлений всех ветвей ячейки, а коэффициент Rik=Rki (i<>k) – взаимное сопротивление контуров, равное сопротивлению общей для контуров i и k ветви, взятому с отрицательным знаком, которым учитываются встречные на­правления контурных токов в рассматриваемой ветви.

Каждое уравнение системы выражает условие равно­весия напряжений ветвей контура – резистивных и источников напряжения (в правой части). Слагаемое на главной диагонали дает напряжение всех резистивных ветвей только от собственного контурного тока, а слагаемое Рkjij=ukj - напря­жение на взаимном сопротивлении контуров только от тока в j-м контуре.

Составление уравнений сводится к записи симметричной матрицы параметров контурных токов:

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

При введении вектора искомых контурных токов уравнения (3.10) в матричной форме можно записать в виде:

.

3. Алгоритм формирования узловых уравнений

Для ввода графа или соединений цепи производят последовательную нумерацию:

  1. всех узлов от 1 до п=пy-1 (опорному узлу присваивают нулевой номер);

  2. всех ветвей от 1 до nв. Как указывалось, всю информацию о структуре графа содержит матрица соединений. Но вводить в память эту матрицу в виде двумерного массива нерационально из-за большого числа нулевых элементов. Поэтому структуру графа вводят с помощью таблицы соединений – одномерного массива троек целых чисел (k, i, j), где kномер ветви; i—номер узла, откуда ветвь выходит; j—номер узла, куда ветвь входит. Тройки чисел дают ненулевые элементы aik=1 и аjk=-1 матрицы А. По заданной таблице соединений можно получить все необходимые для формирования уравнений матрицы.

Формировать узловые уравнения на ЭВМ можно перемножив произведение первых двух матриц па транспонированную матрицу соединений, получим матрицу узловых проводимостей; перемножение транспониро­ванной матрицы АТ на вектор токов источников дает вектор узловых токов. Действия над матрицами (транспонирование, сум­мирование, перемножение) легко программируется. Упомяну­тые матрицы являются разреженными, т. е. содержат много нуле­вых элементов. Поэтому приведенный алгоритм, включающий мно­го действий умножения на нулевой элемент, применять невыгодно.

Более рационально формировать матрицу Gy и вектор iy непосредственно - по мере поступления данных составных ветвей, исходя из смысла собственной и взаимной проводи­мостей. Вначале матрицу узловых проводимостей и вектор узловых токов принимают равными нулю: Gy=0 и iy=0, затем к ним добавляют элементы, вносимые каждой составной ветвью. Данные ветви (k, i, j) с проводимостью Gk и током источника тока iok войдут в собственные проводимости Gii, Gjj, узлов i, j как добавки Gk, во взаимные проводимости Gij, Gji как добавки Gk и в элементы вектора узловых токов iyi, iyi - как добавки ±i0k. Добавки ветви в матрицу проводимости Gy и вектор тока iy можно представить в виде:

Если ветвь присоединена к базисному узлу (j=0) то она внесет добавку Gk, только в собственную проводимость Gii и добавку i0k в составляющую iiy вектора узловых токов. Учет данных последней ветви завершает формирование узловых проводимостей и вектора узловых токов.

Заключение

Таким образом в настоящей дипломной работе:

  • проведен литературный обзор по моделированию процессов в радиотехнических цепях, методов и средств расчетов токов и направлений в них;

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

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

  • приведены методические пояснения при работе с разработанной программой, а также расчетов и анализа электрических цепей;

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

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

В заключение хочу выразить благодарность своим научным руководителям: Цыпишка Дмитрию Ивановичу, Брагарь Луке Федеровичу и заведующиему кафедрой кафедрой общей физики и методики преподавания физики Стамову Ивану Григорьевичу, за методическую помощь, оказанную при написании дипломной работы.

Использованная литература

  1. К.С. Демирчан, П.А. Бутырин. Моделирование и машинный расчет электрических цепей. М., «Высшая школа», 1988г.

  2. В. Нерретер. Расчет электрических цепей на ПЭВМ. М., «Энергоатомиздат», 1991г.

  3. Пантюшин В.С. Сборник задач по электротехнике и основам электронники. М., «Высшая школа», 1979г.

  4. П.Н. Махтанов. Основы анализа электрических цепей. Линейные цепи. М., «Высшая школа», 1990г.

  5. «Электротехника». Под редакцией проф. В.С. Пантюшина. М., «Высшая школа», 1976г.

  6. В.Г. Абрамов, Н.П. Трифонов, Г.Н. Трифонова. «Введение в язык Паскаль». М., «Наука», 1988г.

  7. Ж.Джонс, К. Харроу. «Решение задач в системе Турбо Паскаль». М., «Финансы и статистика», 1991г.

  8. К. Боон. «Паскаль для всех». М., «Энергоиздат», 1988г.

  9. Д. Прайс. «Программирование на языке Паскаль». Практическое руководство. М., «Мир», 1987г.

Приложение

I. Рисунки с видами экрана при работе с программой 20

Рисунок №1. Общий вид экрана. 20

Характеристики

Тип файла
Документ
Размер
409,5 Kb
Тип материала
Предмет
Учебное заведение
Неизвестно

Список файлов реферата

Свежие статьи
Популярно сейчас
Зачем заказывать выполнение своего задания, если оно уже было выполнено много много раз? Его можно просто купить или даже скачать бесплатно на СтудИзбе. Найдите нужный учебный материал у нас!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6644
Авторов
на СтудИзбе
294
Средний доход
с одного платного файла
Обучение Подробнее