61257
Работа добавлена на сайт samzan.net:
Псковский государственный университет
Отчёт
по лабораторной работе №1
«Автоматизация расчета переходных процессов в схемах с управляемыми источниками»
Выполнил:
студент гр. 0024-04
Васильев А.С.
Проверил: преподаватель
Сохор Ю.Н.
Псков
2012г.
Цель работы: построить расчетную схему замещения с управляемым источником, составить алгоритм расчета в соответствии с методом контурных токов, построить график переходного процесса.
Задание:
R1 = 10 кОм
С1 = 1 мкФ
U = 1 В
Rн = 100 кОм
Rвх = 10 Мом
Rвых = 10 Ом
k = 104
Рис.1. Исходная схема
Порядок выполнения работы:
- Составить электрическую схему замещения устройства в соответствии с заданием.
- Представить все ветви в виде обобщенных ветвей.
- Составить алгоритм расчета переходного процесса, используя метод
контурных токов.
- Составить программную реализацию алгоритма на языке Scilab.
- Выполнить расчет переходного процесса.
1. Составим электрическую схему замещения устройства
Рис.2. Электрическая схема замещения устройства
Значения параметров схемы:
R1 = 104 Ом C1 = 10-6 Ф R0 = 0 Ом e1 = 1 В
R2 = 104 Ом C2 = 10-6 Ф Rin = 106 Ом e2 = k*URin В
R3 = 104 Ом C3 = 10-6 Ф Rout = 10 Ом k = 104
Rn = 105 Ом
Начальные условия:
Uc10 = 0
Uc20 = 0
Uc30 = 0
2. Представим все ветви в виде обобщённых ветвей и укажем токи
Рис.3. Преобразованная схема.
3. Составим алгоритм расчёта переходного процесса методом контурных токов
Расчетный шаг и время переходного процесса:
dt=0.00000005
Tsum=0.0000035
Выражения для расчёта начальных условий на ветвях с ёмкостями:
Uc10 = U10
Uc20 = U20
Uc30 = U30
Выражения для расчёта параметров обобщённых ветвей
а) неинерционные ветви:
z1 = R0 e1 = 1
z4 = R3 e4 = 0
z5 = R2 e5 = 0
z9 = Rn e9 = 0
z10 = R1 e10 = 0
б) инерционные ветви для метода Эйлера:
z2 = dt/C2 e2 = - U20
z3 = dt/C2 e3 = - U30
z6 = dt/C2 e6 = - U60
в) управляемый источник:
z7 = Rin; z8 = Rout; z87 = - k * z7; e8 = -k * e7;
Контурная топологическая матрица:
Ik1 Ik2 Ik3 Ik4
z1 1 0 0 0
z2 0 1 1 0
z3 -1 1 0 0
z4 0 0 -1 0
C = z5 0 1 0 0
z6 1 0 0 0
z7 1 0 -1 0
z8 0 0 1 -1
z9 0 0 0 1
z10 0 1 0 0
4. Программная реализация на языке SciLab
clear; // очистка данных
// Ввод параметров
U=1; R1=10000; R2=10000; R3=10000; k=10000; C1=0.000001; C2=0.000001; C3=0.000001; Rin=10000000; Rout=10; Rn=100000;
// Ввод начальных условий
Uc10=0; Uc20=0; Uc30=0;
// Ввод расчетного шага dt и времени переходного процесса
dt=0.00000005; Tsum=0.0000035;
Nt=floor(Tsum/dt); // расчет числа точек
// Ввод формул для присвоения начальных условий машинным переменным
U30=Uc10; U20=Uc20; U60=Uc30;
// Ввод выражений для расчета сопротивлений и источников неинерционных ветвей
z(1,1) = 0; e(1) = U;
z(4,4) = R3; e(4) = 0;
z(5,5) = R2; e(5) = 0;
z(7,7) = Rin; e(7) = 0;
z(9,9) = Rn; e(9) = 0;
z(10,10) = R1; e(10) = 0;
// Ввод выражений для расчета сопротивлений и источников инерционных ветвей
z(2,2) = dt/C2; e(2) = - U20;
z(3,3) = dt/C1; e(3) = - U30;
z(6,6) = dt/C3; e(6) = - U60;
// Ввод выражений для расчета сопротивлений и источников управляемых ветвей
z(7,7) = Rin; z(8,8) = Rout; z(8,7) = - k * z(7,7); e(8,8) = -k * e(7);
// Ввод контурной матрицы
C=[ 1 0, 0, 0;
0, 1, 1, 0;
-1, 1, 0, 0;
0, 0, -1, 0;
0, 1, 0, 0;
1, 0, 0, 0;
1, 0, -1, 0;
0, 0, 1, -1;
0, 0, 0, 1;
0, 1, 0, 0];
// Расчет контурных сопротивлений и проводимостей
Zk = C'*z*C; Yk=Zk^(-1);
// Ввод выражений для объявления массивов, выводимых на графики
Un(1)=0; In(1)=0; t(1)=0;
for n=2:Nt, //Задание цикла по n
//расчет контурных э.д.с.
ek = C'*e;
//расчет контурных токов
ik=Yk*ek;
//расчет токов и напряжений во всех ветвях
i = C*ik; U = z*i - e;
// Ввод выражений для расчета новых начальных условий
U20=U(2); U30=U(3); U60=U(6);
// Ввод выражений для расчета э.д.с., зависящих от начальных условий
e(2) = - U20; e(3) = - U30; e(6) = - U60;
// Ввод выражений для массивов тока и напряжения, выводимых на график
Un(n) = U(9); In(n) = i(9);
// расчет нового значения времени
t(n) = t(n-1) + dt;
end; // конец цикла по n
plot2d(t,[-1*Un, -10000*In ],leg='U@I') // построение графиков
xgrid(2) // задание сетки
5. После выполнения программы получаем графики изменения тока и напряжения на нагрузке Rn
2. тема- Музей Булгакова Нехорошая квартира и Патриаршие пруды ~ литературноисторические памятники Москвы.html
3. Собрание итогов сделанных форм готовых положений мир каков он есть для всякого кто абонирован в нем
4. Лекция 8 Права и свободы человека гражданские политические экономические
5. Вехи] Николай БердяевФИЛОСОФИЯ СВОБОДНОГО ДУХА ОГЛАВЛЕНИЕ Введение
6. реферату- Малокаліберна гвинтівкаРозділ- Військова справа ДПЮ Малокаліберна гвинтівка Малокаліберна гв
7. а Оглавление [1] История Узайра Ездры Эзры мир ему [2] Рассказ о Соломоне
8. Станислав Лем Межличностная коммуникация Я стараюсь максимально позитивно излагать материал и следу
9. тематического института Клэя в миллион долларов США
10. Остальцы Русской Древлеправославной веры
11. Утверждаю Зав
12. і Землеробство і скотарство не тільки були спроможні прогодувати населення країни а й виробляли продукти х
13. Облік і аудит в управлінні підприємницькою діяльністю в наступному складі-
14. Охрана Водных Экосистем
15. РЕФЕРАТ дисертації на здобуття наукового ступеня кандидата технічних наук Київ ~ 2002 Дисерт.
16. 45 лет Дмитрий 45 лет Наталья его соседка 45 лет Зорро её муж 45 лет Анатолий их сын 18 лет Отец Две
17. Национальный банк Украины и его функции
18. 069 мкм с твердым телом рубином получено в 1960 г
19. Лингвокультурное исследование фразеологизмов времени в русском и английском языках
20. за кордону надходить на внутрішнє ринок найбільш конкурентноздатна продукція вітчизняного виробника і фор