ЛАБОРАТОРНАЯ РАБОТА № 2. Преобразование Лапласа. Нахождение оригинала функции по её изображению.. ЛР 2 МХТ 14. Лабораторная работа Преобразование Лапласа. Нахождение оригинала функции по её изображению
![]()
|
ЛАБОРАТОРНАЯ РАБОТА № 2. Преобразование Лапласа. Нахождение оригинала функции по её изображению. Цель работы: изучить преобразование Лапласа, его свойства, способы вычисления оригинала функции по её изображению, а также применение MATLAB к вычислению оригиналов функций по их изображениям. Закрепить полученные знания на практике. Приборы и оборудование: - Компьютер совместимый с IBM PC,128-512 Мб. ОЗУ; - Операционная система WINDOWS NT, XP, UNIX; - Математический пакет MATLAB Version 7.*. Форма отчётности студентов: индивидуальный отчёт с типовым титульным листом и результатами моделирования. Длительность работы: 4 академических часа. Защита работы: собеседование с преподавателем по контрольным вопросам, выполнение индивидуальных заданий. I. ТЕОРЕТИЧЕСКАЯ ЧАСТЬ. Функция, которая подвергается преобразованию Лапласа, должна обладать следующими свойствами: 1. функция должна быть определена и дифференцируема по всей положительной полуоси ; 2. функция должна быть тождественно равна 0 при ![]() ![]() ![]() 3. функция должна быть ограниченна, т.е. для функции ![]() ![]() ![]() ![]() ![]() Преобразование Лапласа это соотношение вида ![]() ![]() ![]() ![]() ![]() ![]() При этом ![]() ![]() Символическая запись преобразования Лапласа, а именно, ![]() ![]() Нахождение оригиналов функций по их изображениям Если изображение функции является дробно-рациональной функцией вида ![]() нахождение оригиналов по аналитическим формулам: 1. Корни простые, вещественные: ![]() 2 Корни простые, вещественные и один корень нулевой, т.е. ![]() ![]() 3. Корни комплексно-сопряженные: ![]() ![]() ![]() ![]() 4. Корни комплексно-сопряженные и один нулевой: ![]() II. ПРАКТИЧЕСКАЯ ЧАСТЬ. Определим порядок вычисления оригинала функции по его изображению: 1. Вычислить корни полинома ![]() ![]() Число корней ![]() ![]() ![]() 2. Выбрать формулу для расчёта оригинала. 3. Вычислить производную ![]() 4. Вычислить значения полиномов ![]() ![]() 5. Вычислить значения коэффициентов при функциях ![]() Пример 1: ![]() 1. ![]() 2. Корни вещественные, выбираем формулу (1): ![]() 3. Вычисляем ![]() ![]() Рассмотрим пример вычисления с использованием MATLAB: >> p=[1 5 4]; - коэффициенты полинома знаменателя; >> b=[3]; - коэффициенты полинома числителя; >>r=roots(p) – вычисление корней знаменателя; r = -4 -1 >>r1=r(1); r2=r(2) >> dp=polyder(p) – вычисление коэффициентов производной ![]() dp = 2 5 >> A1=polyval(dp,r1) – вычисление ![]() A1 = 3 >> A2=polyval(dp,r2) – вычисление ![]() A2 = -3 >>В1= polyval(b, r1) >>В1=3 >> C1=В1./A1 – вычисление коэффициента при первом корне; C1 = 1 >>В2= polyval(b, r2) >>В2= 3 >> C2=В2./A2 – вычисление коэффициента при втором корне; C2 = -1 >> t=[0:0.01:5] ;– интервал времени и шаг счёта; >> x=C1.*exp(r1*t)+C2.*exp(r2.*t); >> plot(t,x),grid on, xlabel('Time(sec)'), ylabel('x(t)') ![]() Пример 2: ![]() 1. ![]() 2. Корни вещественные, выбираем формулу (2): ![]() 3. Вычисляем ![]() ![]() Воспользуемся результатами расчёта из примера 1 и пересчитаем коэффициенты: >>C3=C1./(r1) – вычисление коэффициента при первом корне; C3 = -1 >> C4=C2./(r2) – вычисление коэффициента при втором корне; C4 = 0.2500 >>В0= polyval(b, 0)=3 >>А0=polyval(p,0)=4; >>С0=В0./А0 – вычисление коэффициента при нулевом корне; С0=0.75 >> x1=С0+C3.*exp(r1.*t)+C4.*exp(r2.*t); >> plot(t,x1),grid on, xlabel('Time(sec)'), ylabel('x1(t)') ![]() Пример 3: ![]() 1. ![]() 2. Корни комплексно-сопряжённые, выбираем (3) ![]() 3. Вычисляем ![]() 4. Расчёт ![]() ![]() ![]() Представим комплексный корень в тригонометрическом виде ![]() ![]() Алгоритм вычисления оригинала функции, если корни знаменателя комплексные, остаётся прежним. >> p=[1 2 5]; >> r=roots(p) r = -1.0000 + 2.0000i -1.0000 - 2.0000i >>r1=r(1)=-1+2.*i – выбираем один из корней и для него проводим вычисления; >> dp=polyder(p) – вычисление коэффициентов производной ![]() dp = 2 2 >> A1=polyval(dp, r1) – вычисление ![]() A1 =0 + 4.0000i >>В1= polyval(b, r1) >>B1=2; >> C1=B1./A1 – вычисление коэффициента при выбранном корне; C1 =0 - 0.5000i >> t =[0:0.01:5]; >> x= exp(-1.*t).*sin(2.*t); >> plot(t,x), grid on, xlabel('Time(sec)'), ylabel('x(t)'). ![]() Пример 4: ![]() 1. ![]() 2. Корни комплексно-сопряжённые и один нулевой, выбираем(4) ![]() 3. Вычисляем ![]() 4. Расчёт ![]() ![]() ![]() Воспользуемся результатами расчёта из примера 3 и пересчитаем коэффициенты >>C3=С1./r1 – вычисление коэффициента при выбранном корне; >>C3 =- 0.2000 +0.1000i >>A0=5; >>B0=2 >>C0=B0./A0– вычисление коэффициента при нулевом корне; |>>0.4000 >> x1=C0 – 0.4.*exp(-1.*t).*cos(2.*t) – 0.2.*exp(-1.*t).*sin(2.*t)); >> plot(t,x1), grid on, xlabel('Time(sec)'), ylabel('x1(t)') ![]() III. ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ Выбрать из таблицы 1 в соответствии с номером варианта изображения функции для выполнения заданий 1 и 2. Задание 1. 1. Вычислить самостоятельно оригинал х(t) по изображению ![]() 2. Определить начальное и конечное значение функции, используя теоремы о предельных значениях из свойств преобразования Лапласа. 3. Определить оригинал функции с применением вычислительных процедур в пакете MATLAB и построить график этой функции x(t). Задание 2. 1. Вычислить самостоятельно оригинал х(t) по изображению ![]() 2. Определить начальное и конечное значение функции, используя теоремы о предельных значениях из свойств преобразования Лапласа. 3. Определить оригинал функции с применением вычислительных процедур в пакете MATLAB и построить график этой функции x(t). Выбрать из таблицы 1 свой вариант изображения функции для вычисления оригинала функции для выполнения заданий 3 и 4. Задание 3. 1. Вычислить самостоятельно оригинал х(t) по изображению ![]() 2. Определить начальное и конечное значение функции, используя теоремы о предельных значениях из свойств преобразования Лапласа. 3. Определить оригинал функции с применением вычислительных процедур в пакете MATLAB и построить график этой функции x(t). Задание 4. 1. Вычислить самостоятельно оригинал х(t) по изображению ![]() 2. Определить начальное и конечное значение функции, используя теоремы о предельных значениях из свойств преобразования Лапласа. 3. Определить оригинал функции с применением вычислительных процедур в пакете MATLAB и построить график этой функции x(t). Варианты заданий
IV. СОДЕРЖАНИЕ ОТЧЁТА. Изображения заданных функций для каждого задания. Вычисление оригинала функции по её изображению с использованием Matlab. Начальные и конечные значения функций, определённые по предельным теоремам. Графики функций по каждому заданию. V. КОНТРОЛЬНЫЕ ВОПРОСЫ. Какими свойствами должна обладать функция, подвергаемая преобразованию Лапласа? Дайте определение преобразования Лапласа. Дайте определение передаточной функции в изображениях по Лапласу. Какие значения комплексной переменной называются полюсами? Какие значения комплексной переменной называются нулями? Укажите взаимосвязь и различие между операторной передаточной функцией и передаточной функцией в преобразованиях Лапласа. |