Инструкция по выполнению работы

Основная часть команд вводится в командном окне среды Matlab. Команды, которые надо применять в других окнах, обозначены иконками соответствующих программ.

Этап выполнения задания Команды Matlab
1. Очистите рабочее пространство Matlab (память). clear all
2. Очистите окно Matlab. clc
3. Посмотрите краткую справку по команде tf. help tf
4. Определите адрес файла, который выполняет эту команду. which('tf')
5. Введите передаточную функцию[6] как объект tf. n = [n2 n1 n0] d = [1 d2 d1 d0] f = tf ( n, d )
6. Проверьте, как извлечь из этого объекта числитель и знаменатель передаточной функции. [n1,d1] = tfdata ( f, 'v' )
7. Найдите нули и полюса передаточной функции. z = zero ( f ) p = pole ( f )
8. Найдите коэффициент усиления звена в установившемся режиме. k = dcgain ( f )
9. Определите полосу пропускания системы (наименьшую частоту, на которой АЧХ становится меньше, чем дБ). b = bandwidth ( f )
10. Постройте модель системы в пространстве состояния. f_ss = ss ( f )
11. Сделайте так, чтобы коэффициент прямой передачи звена был равен 1. f_ss.d = 1
12. Найдите новый коэффициент усиления звена в установившемся режиме. k1 = dcgain ( f_ss )
13. Как связаны коэффициенты и ? Почему?  
14. Постройте модель исходной системы в форме «нули-полюса». f_zp = zpk ( f )
15. Проверьте, какие переменные есть в рабочем пространстве. who или whos (в чем разница?)
16. Постройте на графике расположение нулей и полюсов системы. pzmap ( f )
17. Определите коэффициенты демпфирования и собственные частоты для всех элементарных звеньев (первого и второго порядка). [wc,ksi,p] = damp ( f )
18. Запустите модуль LTIViewer. ltiview
19. Загрузите модель f. File – Import
20. Постройте импульсную характеристику (весовую функцию) этой системы. ПКМ – Plot Types - Impulse
21. Загрузите модель f_ss. File – Import
22. Проверьте, построена ли импульсная характеристика второй системы? ПКМ – Systems
23. Отключите систему f. Почему одинаковы построенные импульсные характеристики разных систем? ПКМ – Systems
24. Подключите обе системы. ПКМ – Systems
25. Постройте переходные характеристики систем. ПКМ – Plot Types – Step
26. Сделайте, чтобы на графике для каждой функции были отмечены: · максимум · время переходного процесса[7] · время нарастания (от 10% до 90% установившегося значения) · установившееся значение ПКМ – Characteristics: · Peak Response · Settling Time · Rise Time · Steady State  
27. Щелкая мышью по меткам-кружкам, выведите на экран рамки с численными значениями этих параметров и расположите их так, чтобы все числа были видны.  
28. Экспортируйте построенный график в отдельное окно. File – Print to Figure
29. Скопируйте график в буфер обмена в формате векторного метафайла. print -dmeta
30. Вставьте график из буфера обмена в отчет (Microsoft Word). ПКМ - Вставить
31. Закройте окно LTIViewer.  
32. Создайте массив частот для построения частотной характеристики[8] (100 точек в интервале от до с равномерным распределением на логарифмической шкале). w = logspace(-1, 2, 100);
33. Рассчитайте частотную характеристику исходной системы [9]… r = freqresp ( f, w ); r = r(:);
34. … и постройте ее на осях с логарифмическим масштабом по оси абсцисс. semilogx ( w, abs(r) )
35. Скопируйте график в буфер обмена в формате векторного метафайла. print -dmeta
36. Вставьте график из буфера обмена в отчет (Microsoft Word). Объясните, где на графике можно найти коэффициент усиления в статическом режиме и как определить полосу пропускания системы. ПКМ – Вставить
37. Закройте все лишние окна, кроме командного окна Matlab.  
38. Постройте сигнал, имитирующий прямоугольные импульсы единичной амплитуды с периодом 4 секунды (всего 5 импульсов). [u,t] = gensig('square',4);
39. Выполните моделирование и постройте на графике сигнал выхода системы f при данном входе. lsim (f, u, t)
40. Скопируйте график в буфер обмена в формате векторного метафайла. print -dmeta
41. Вставьте график из буфера обмена в отчет (Microsoft Word). ПКМ – Вставить

 

Таблица коэффициентов

Вариант
1. 1.0 1.10 0.100 3.0000 3.1600 1.2000
2. 1.1 1.54 0.495 2.8000 2.9200 1.2000
3. 1.2 1.08 0.096 2.3727 2.2264 0.9091
4. 1.3 1.04 0.091 2.1909 2.0264 0.9091
5. 1.4 -1.54 0.252 1.8333 1.5278 0.6944
6. 1.5 -0.90 -0.240 1.6667 1.3611 0.6944
7. 1.6 0.80 -0.224 1.3286 0.8959 0.4592
8. 1.7 1.36 0.204 1.1857 0.7673 0.4592
9. 1.8 -1.98 0.432 1.2000 0.7644 0.3556
10. 1.9 -0.76 -0.399 1.3333 0.8711 0.3556
11. 2.0 0.60 -0.360 1.2000 0.7406 0.2734
12. 2.1 1.68 0.315 1.3250 0.8281 0.2734
13. 2.2 -2.42 0.616 1.3059 0.7696 0.2076
14. 2.3 -0.46 -0.552 1.4235 0.8401 0.2076
15. 2.4 0.24 -0.480 1.3889 0.7531 0.1543
16. 2.5 2.25 0.500 1.5000 0.8086 0.1543
17. 2.6 0.26 -0.780 1.2421 0.6139 0.1108
18. 2.7 -0.27 -0.810 1.1368 0.5717 0.1108
19. 2.8 0.28 -0.840 0.8000 0.3700 0.0500
20. 2.0 0.50 -0.380 1.2600 0.7804 0.2844
21. 1.2 1.10 0.100 2.4234 2.1864 0.9091
22. 2.4 0.22 -0.520 1.4024 0.7835 0.1543
23. 1.8 -1.96 0.4269 1.2200 0.7673 0.3656
24. 2.9 3.19 0.870 0.7000 0.3500 0.0500

 

Контрольные вопросы к защите

1. Что такое

· передаточная функция

· нули и полюса передаточной функции

· импульсная характеристика (весовая функция)

· переходная функция

· частотная характеристика

· модель в пространстве состояний

· модель вида «нули-полюса»

· коэффициент усиления в статическом режиме

· полоса пропускания системы

· время переходного процесса

· частота среза системы

· собственная частота колебательного звена

· коэффициент демпфирования колебательного звена

2. В каких единицах измеряются

· коэффициент усиления в статическом режиме

· полоса пропускания системы

· время переходного процесса

· частота среза системы

· собственная частота колебательного звена

· коэффициент демпфирования колебательного звена

3. Как связана собственная частота с постоянной времени колебательного звена?

4. Может ли четверка матриц

быть моделью системы в пространстве состояний? Почему? Какие соотношения между матрицами должны выполняться в общем случае?

5. Как получить краткую справку по какой-либо команде Matlab?

6. В чем разница между командами Matlab

who и whos clear all и clc

7. Как ввести передаточную функцию ?

8. Как влияет изменение коэффициента прямой передачи (матрицы в модели в пространстве состояний) на статический коэффициент усиления?

9. Какие возможности предоставляет модуль LTIViewer?

10. Что можно сказать об импульсной характеристике системы f_ss? Почему она не была построена верно?

11. Как найти

· коэффициент усиления в установившемся режиме по АЧХ

· полосу пропускания системы по АЧХ

12. Как скопировать график из окна Matlab в другую программу?

13. Как построить массив из 200 значений в интервале от до с равномерным распределением на логарифмической шкале?

14. Какие величины откладываются по осям на графике АЧХ?