Labview решение системы линейных уравнений

Лабораторная работа 3. Решение алгебраических уравнений в среде LabVIEW

Пример 3.1. Определение токов в цепи с использованием формульного узла

1. Токи в цепи рис. 3.1 можно рассчитать по законам Ома и Кирхгофа

;

; (1)

.

Значения сопротивлений ,,и электродвижущей силы Е выбираются самостоятельно.

2. Расчет по этим формулам можно выполнить при помощи формульного узла Formula Node, который относится к элементам «Структуры» и вызывается правой клавишей мыши на панели блок-диаграмм по пути: All Functions => Structures => Formula Node. Появившаяся рамка формульного узла растяги­вается до нужного размера и в нее вписываются расчетные формулы (1). Не­известные записываются в левой части формул. Каждая формула пишется на отдельной строке и заканчивается точкой с запятой.

3. Затем в формулы нужно внести исходные данные и вывести результаты рас­чета. Для этого курсор устанавливается правой клавишей мыши на рамке фор­мульного узла и из всплывающего меню левой клавишей вызывается Add Input (добавить вход) для входных величин и Add Output (добавить выход) для выходных величин. В появившиеся рамки вписываются наименования этих величин.

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

5. При помощи управляющих элементов задаются исходные данные, после чего схема запускается на решение.

Вид формульного узла показан на рис. 3.2.

Примечание. Формульный узел не применяется для работы с комплексными числами.

Пример 3.2. Решение алгебраических уравнений в матричной форме

Расчет токов в цепи рис. 3.1 можно провести по линейным алгебраическим уравне­ниям, составленным по законам Кирхгофа:

(1)

Эти уравнения можно записать в матричной форме

. (2)

Для решения системы линейных алгебраических уравнений в среде LabVIEW существует элемент «Решение линейных уравне­ний», вызываемый по пути All Functions => Analyze => Mathe­matics => Linear Algebra => Solve Linear Equatons. Терминалы его можно раскрыть нажатием на иконку правой клавишей мыши и далее из всплыва­ющего меню Visible Items => Terminals. Вид элемента, исходный и с открытыми терминалами, изображен на рис. 3.3.

Для определения назначения терминалов можно снова нажать на иконку правой клавишей мыши и открыть переднюю панель Open Front Panel. Передняя панель имеет вид, представленный на рис. 3.4. Уравнение (2) вводится таким образом: мат­рица коэффициентов Input Matrix (первая матрица уравнения) — подается на левый верхний терминал, вектор заданных воздействий Known Vector (правая часть) — под­водится к левому среднему терминалу, а результирующий вектор Solution Vector (век­тор искомых токов) — снимается с правого верхнего терминала

Для того чтобы ввести матрицу, вызывается элемент Array (построение масси­ва). Элемент вызывается на лицевой панели по пути All Controls => Array&Cluster => Array. Появляется пустая ячейка, в которую вносится цифровой управляющий элемент (для матрицы коэффициентов и вектора управляющих воздействий) или цифровой индикатор (для вектора искомых токов). Затем элемент растягивается до нужной размерности матрицы инструментом «перемещение» (стрелка). Окна управляющих элементов и индикаторов имеют серый цвет, после внесения туда цифровых данных они становятся белыми. Одновременно с вызовом элемента Array на лицевой панели появляется его иконка на панели блок-диаграмм. Вид ли­цевой панели и панели блок-диаграмм, где выполнены необходимые соединения, представлен на рис. 3.5.

В результате выполнения операции получается вектор решений Output. В этом столбце представлены значения токов /, = ЗА, /2 = 2А, /3 = 1А.

Можно из этого столбца выделить значения отдельных неизвестных. Для этого используется функция извлечения элемента массива All Functions => Array => Index Array. Выходы функции подключаются к элементам индикации (рис. 3.6). После запуска программы индикаторы покажут величины, соответствующие зна­чениям строк выходного вектора.

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

Видео:Метод Крамера за 3 минуты. Решение системы линейных уравнений - bezbotvyСкачать

Реферат: Моделирование линейных непрерывных систем в среде LabVIEW

Видео:Решение системы линейных уравнений. Подстановка. С дробными выражениями.Скачать

Введение

LabVIEW (LaboratoryVirtualInstrumentEngineeringWorkbench) позволяет разрабатывать прикладное программное обеспечение для организации взаимодействия с измерительной и управляющей аппаратурой, сбора, обработки и отображения информации и результатов расчетов, а также моделирования как отдельных объектов, так и автоматизированных систем в целом. Разработчиком LabVIEW является американская компания National Instruments.

LabVIEW является открытой системой программирования и имеет встроенную поддержку всех применяемых в настоящее время программных интерфейсов, таких как Win32 DLL, COM.net, DDE, сетевых протоколов на базе IP, DataSocket и др. В состав LabVIEW входят библиотеки управления различными аппаратными средствами и интерфейсами, такими как PCI, CompactPCI/PXI, VME, VXI, GPIB (КОП), PLC, VISA, системами технического зрения и др. Программные продукты, созданные с использованием LabVIEW, могут быть дополнены фрагментами, азработанными на традиционных языках программирования, например C/С++, Pascal, Basic, FORTRAN. И наоборот можно использовать модули, разработанные в LabVIEW в проектах, создаваемых в других системах программирования. Таким образом, LabVIEW позволяет разрабатывать практически любые приложения, взаимодействующие с любыми видами аппаратных средств, поддерживаемых операционной системой компьютера.

среда программирование дифференциальное уравнение

Видео:Общее, частное, базисное решение системы линейных уравнений Метод ГауссаСкачать

Моделирование линейных непрерывных систем

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

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

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

Видео:Решение системы линейных уравнений с двумя переменными способом подстановки. 6 класс.Скачать

Численное решение дифференциальных уравнений

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

Здесь x= x (t ) — независимая функция (входной процесс), u= u (t ) — решение уравнения (выходной процесс).

Численное решение находится для дискретных значений аргумента t , отличающихся на шаг интегрирования Dt . В одношаговых разностных методах для нахождения следующего значения u _к = u (t _к ) требуется информация только об одном предыдущем шаге. Из одношаговых методов наибольшую известность получили методы Рунге-Кутта. В основу метода Рунге-Кутта первого порядка, называемого также явным или прямым методом Эйлера, положено разложение функции u (t ) в ряд Тейлора в окрестности точки A (t _k-1, , u _k-1 ):

В методах Эйлера (и Рунге-Кутта тоже) ограничиваются только двумя первыми членами разложения в ряд. Запишем значение u_k = u (t_k ), приняв в выражении (5.2) t= t_k и ограничившись двумя первыми членами ряда:

Это выражение является приближенным решением дифференциального уравнения (1) прямым методом Эйлера. Оно рекуррентное и позволяет найти значение выходного процесса u_k по значениям выходного и входного процессов в предыдущем такте.

На рис. 1 а ) проиллюстрировано решение прямым методом Эйлера.

а )	б )
Рис.1

Видим, что при использовании этого метода используется линейная экстраполяция и тангенс угла наклона экстраполирующей прямой равен производной функции u (t ) в точке А. Экстраполированное значение u_k отличается от точного на величину ошибки.

Неявный (обратный) метод Эйлера основан на разложении функцииu (t ) в ряд Тейлора в окрестности точки В (u _k, , t _k ) (см. рис.1 б ):

Приняв в этом выражении t= t_{k — 1} и ограничившись двумя первыми членами ряда, получим

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

Применим методы Эйлера для расчета переходной характеристики интегрирующей цепи. Передаточная функция интегрирующей цепи:

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

и в канонической форме:

Перепишем его в виде (1):

Запишем рекуррентную формулу для прямого метода Эйлера в соответствии с (5.3)

Формула для обратного метода Эйлера запишется в соответствии с (4)

Так как уравнение линейное, то значение y_k вычисляется в явной форме:

Методы Эйлера обладают низкой точностью. В более точных методах используются различные способы определения угла наклона экстраполирующей прямой, чтобы она прошла ближе к точному решению. Хорошей точностью обладает метод Рунге-Кутта четвертого порядка, который обычно и используется. Программы для численного решения дифференциальных уравнений имеются практически в любом пакете прикладных программ, в том числе и в LabVIEW.

Для вычислений по формулам (5.5) и (5.6) используем структуру FormulaNode. Внутри этой структуры запишем точное выражение для переходной характеристики:

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

и обратным методом Эйлера:

при нулевых начальных условиях: y (0) = 0, v (0) = 0.

В этих выражениях использованы различные обозначения для выходных переменных и принято x = 1 (t ) = 1, так как t > 0.

На рис.2 показана эта структура. В формулах Δt обозначена как dt .

Рис.2	Рис.3

Напомним, что для образования входных и выходных терминалов нужно щелкнуть ПКМ на границе структуры в предполагаемом месте терминала и в раскрывшемся меню выбрать AddInputили AddOutput.

Для формирования массивов выходных переменных структура FormulaNodeпомещается внутрь структуры ForLoop, при этом задержанные на интервал дискретизации отсчеты выходных переменных y ₁ и v ₁ получаются с помощью регистра сдвига (рис.3).

Прямой метод Эйлера при большом интервале дискретизации может дать неустойчивое решение. Это случится, если отклонение решения от входного процесса x_{k — 1} — y_{k — 1} (см формулу (5)) даст такое значение y_k . что отклонение на следующем шаге x_k — y_k будет той же величины, что и предыдущее, но обратным по знаку. Решение будет колебательным незатухающим.

В предыдущих лабораторных работах развертка графического индикатора Graphосуществлялась автоматически в соответствии с типом данных, подаваемых на вход графического индикатора. В этой работе мы сформируем данные так, чтобы по горизонтальной оси откладывалось время. Для этого надо сформировать кластер, куда кроме массива данных будет входить информация о времени. Используем ВП Bundle (Объединить), который находится в подпалитре Cluster (Кластер). На его входы elementподаются (см. рис.4): на верхний — время начала развертки — 0; на средний — интервал дискретизации — Δt; на нижний — массив данных

Видео:Решение системы уравнений методом ГауссаСкачать

Замена непрерывной передаточной функции дискретной

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

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

1

p 2

(1 +pT )

Видео:Базисные решения систем линейных уравнений (01)Скачать

Сопромат .in.ua

изучаем сопротивление материалов

Видео:15. Однородная система линейных уравнений / фундаментальная система решенийСкачать

Некоторые варианты решения системы уравнений

Видео:Базисные решения систем линейных уравнений (03)Скачать

Решение online

На странице “Решение системы линейных уравнений“ вы можете решить систему до 10 уравнений online.

Видео:Решение системы неравенств с двумя переменными. 9 класс.Скачать

Решение системы линейных алгебраических уравнений в OpenOffice.org

Из линейной алгебры нам известно, что для системы линейных уравнений A*x=b решением 1 есть x=A -1 *b, если детерминант матрицы A не равен нулю. Здесь A -1 – обратная матрица коэффициентов СЛАУ .

Решим систему из двух уравнений.

В ячейках (A1:B2) вводим данные.
Детерминант матрицы: MDETERM(A1:B2)

Затем получаем обратную мартицу A -1 :
=

Умножаем полученную обратную матрицу на столбик з правою частью систем уравнений:
=

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

Видео:Система линейных уравнений. Общее решение. Метод ГауссаСкачать

Пример решения системы линейных уравнений в Scilab

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

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

Более подробно с пакетом можно ознакомится на домашней странице Scilab или в википедии

Пример

Решим систему из двух уравнений вида A*x=b

Задаем матрицу коэффициентов: —>A=[1,2;3,4] A = 1. 2. 3. 4. Задаем вектор свободных коэффициентов: —>b=[2;1] b = 2. 1.

Определяем детерминант матрицы: —>det(A) ans = — 2 Так как определитель матрицы отличен от нуля, решение СЛАУ существует.

Решение СЛАУ методом обратной матрицы —>x=Ab x = — 3. 2.5 Запись Ab соответствует математической операции A -1 b. Записи x=Ab эквивалентна запись x=inv(A)*b

Эту систему уравнений также можно решить с помощью специальной функции linsolve (эта функция решает уравнение вида Ax+b=0, поэтому поменяем знак вектора b) x=linsolve(A,-b) x = — 3. 2.5 и получили ответ: x₀= -3; x₁= 2.5;

Проверяем A*x-b ans = 1.0D-14 * 0.2220446 0.3552714 результирующий вектор близок к нулю, то есть система решена верно.
Существуют и другие способы решения системы линейных уравнений в среде Scilab.

Видео:6 способов в одном видеоСкачать

Пример решения системы линейных уравнений в среде программирования R

R — язык программирования для статистической обработки данных и работы с графикой, а также свободная программная среда вычислений с открытым исходным кодом.

Пример

Решаем систему Ax=b (из двух уравнений)

A A [,1] [,2] [1,] 1 2 [2,] 3 4 > b b [1] 2 1 > x x [1] -3.0 2.5 и получили ответ: x₁= -3; x₂= 2.5;

Проверка: > A %*% x [,1] [1,] 2 [2,] 1

Определитель матрицы: > det(A) [1] -2

Видео:Неоднородная система линейных уравненийСкачать

Решение системы линейных уравнений Ax=b в Euler Math Toolbox

Домашняя страница Euler Math Toolbox . В пакет интегрирована система компьютерной алгебры Maxima. Система напоминает Matlab, но имеет собственный стиль, и немного иной синтаксис.
Домашняя страница автора René Grothmann

Пример

This is Euler Math Toolbox, Version 5.8. >A:=[1,2;3,4] 1 2 3 4 >b:=[2;1] 2 1 >x:=Ab -3 2.5 Получили ответ: x₀= -3; x₁= 2.5;

Определитель матрицы: >det(A) -2

Проверка: умножим A на x (операция умножения обозначается точкой «.») >A.x 2 1

Решить указанную систему можно также с помощью следующего выражения (inv(A) — обратная матрица) >x:=inv(A).b -3 2.5

1 Этот способ решения называется — решение СЛАУ методом обратной матрицы

📹 Видео

Системы уравнений с двумя переменными - 9 класс алгебраСкачать

Solving quadratic equations on LabVIEW / Решения квадратных уравнений в LabVIEW.Скачать

Неоднородные системы линейных уравненийСкачать

Решение систем линейных уравнений, урок 5/5. Итерационные методыСкачать

LLVM и GCC / Register allocation / IT-образование в вузах (часть 1)Скачать

Фундаментальная система решений системы линейных уравнений ФСР СЛАУСкачать

Решение системы линейных уравнений графическим методом. 7 класс.Скачать