Нелинейности с симметричными характеристиками, их моделирование и фиксация на входе и выходе каждого звена средствами пакета Matlab. Изучение процессов в нелинейной системе с нелинейным элементом. Исследование систем методом гармонического баланса.
Аннотация к работе
Исследование нелинейных систем методом фазовой плоскости 19 3. Работы выполняются с использованием пакета Matlab. Характеристики нелинейных элементов: а - идеальное двухпозиционное реле; б - усилитель с ограничением и зоной нечувствительности; в-трехпозиционное реле; г - двухпозиционное реле с гистерезисом; д - люфт Цель работы - моделирование указанных нелинейностей и фиксация процессов на входе и выходе каждого нелинейного звена средствами пакета Matlab (c использованием его расширения - пакета моделирования динамических систем Simulink). Для проведения моделирования создать 2 файла: файл-сценарий TN_prog.m и файл модели TN_mod.mdl. Изучение типичных нелинейностей Используемые обозначениЯ: НЭ - нелинейный элемент. clear all % очистка памЯти close all % закрытие всех предыдущих рисунков Задание значениЯ переменной, определЯющей положение переключателЯ %конфигурации в файле TN_mod.mdl %1 - НЭ-идеальное двухпозиционное реле %2 - НЭ-усилитель с ограничением и зоной нечувствительности %3 - НЭ-трехпозиционное реле без гистерезиса %4 - НЭ-двухпозиционное реле с гистерезисом %5 - НЭ-люфт config = 5; %Определение значениЯ строковой переменной nlin switch config case 1, nlin = ид. 2-х поз. реле; case 2, nlin = ус-ль с огр. и зоной нечувст.; case 3, nlin = 3-х поз. реле без гист.; case 4, nlin = 2-х поз. реле с гист.; case 5, nlin = люфт; end %Константы, описывающие нелинейные элементы c = 6; b = 1; b2 = 7; %Параметры моделированиЯ t_end = 10; %времЯ моделированиЯ, с step = 1e-3;%шаг моделированиЯ, с %Описание генератора затухающих колебаний N = 15; %количество колебаний за времЯ моделированиЯ A = 10; %амплитуда в начальный момент времени lambda = 2*pi*N/t_end; %угловаЯ частота затухающих колебаний ksi = 1.5/sqrt((pi*N)^2 2.25); %декремент затуханиЯ %ksi = 0; omega = (2*pi*N)/(t_end*sqrt (1-ksi^2));%угловаЯ чатота гармонических колебаний gamma = ksi*omega; %коэффициент затуханиЯ %вызов модели open_system (TN_mod.mdl); %запуск моделированиЯ sim (TN_mod); %Построение процессов во времени (рис. Структурная схема модели Рис. 1.3. %Используемые обозначениЯ: НЭ - нелинейный элемент, ЛЧ - линейнаЯ часть.