3. Определение параметров колебательного звена по временным характеристикам

disp('Лабораторная работа №1')

%

clear all;

close all;

disp('Определение параметров коллебательного звена по временным характеристикам ')

for i = 1:100

%вводим начальные данные

R = input('enter resistance, kOm: R = ');

L = input('enter inductance, Gn: L = ');

C = input('enter capacity, mkF: C = ');

%считаем постоянную времени и показатель коллебательности

T=(L*C*0.001)^1/2;

ksi=R*1000/2*(C*0.001/L)^1/2;

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

W = tf([1],[T^2 2*T*ksi 1])

%открываем файл в Simulink

open kolleb.mdl

%строим график переходного процесса

figure(1);

stepplot(W);

grid on;

hold all;

title('Переходные процессы коллебательных звеньев с разными T и ksi');

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

hust = dcgain(W);

[y ~] = step(W);

hmax = max(y);

sigma = (hmax - hust)/hust*100

%выводим легенду

r(i,1) = R;

c(i,1) = C;

l(i,1) = L;

L1(i,:)=['R ='];

L2(i,:)=['C ='];

L3(i,:)=['L ='];

Legen = [L1, num2str(r),repmat(';', size(L1,1), 1),L2, num2str(c),repmat(';', size(L1,1), 1),L3, num2str(l)];

legend(Legen);

%запрос о продолжении

k = input('Хотите продолжить? Y/N:', 's');

if (k == 'n') || (k == 'N'), break; end;

end

ПРИЛОЖЕНИЕ 2

ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ДИНАМИЧЕСКИХ ЗВЕНЬЕВ ПО ИХ ЧАСТОТНЫМ ХАРАКТЕРИСТИКАМ

1. Определение параметров апериодического звена по частотным характеристикам

disp('Лабораторная работа №2')

clear all;

close all;

disp('Определение параметров апериодического звена по частотным характеристикам ')

for i = 1:100

%вводим начальные данные

R = input('enter resistance, kOm: R = ');

C = input('enter capacity, mkF: C = ');

%считаем постоянную времени

T=R*C;

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

W = tf([1],[T 1])

%открываем файл в Simulink

open aperiod.mdl

%строим АФХ

figure(1);

nyquist(W);

grid on;

hold all;

title('АФХ апериодических звеньев с разными постоянными времени');

%выводим легенду

r(i,1) = R;

c(i,1) = C;

L1(i,:)=['R ='];

L2(i,:)=['C ='];

Legen = [L1, num2str(r),repmat(';', size(L1,1), 1),L2, num2str(c)];

legend(Legen);

%строим ЛАЧХ и ФЧХ

figure(2);

bodeplot(W);

grid on;

hold all;

title('ЛАЧХ и ФЧХ апериодических звеньев с разными постоянными времени');

%выводим легенду

r(i,1) = R;

c(i,1) = C;

L1(i,:)=['R ='];

L2(i,:)=['C ='];

Legen = [L1, num2str(r),repmat(';', size(L1,1), 1),L2, num2str(c)];

legend(Legen);

%запрос о продолжении

k = input('Хотите продолжить? Y/N:', 's');

if (k == 'n') || (k == 'N'), break; end;

end

2. Определение параметров дифференциального звена по частотным характеристикам

disp('Лабораторная работа №2')

clear all;

close all;

disp('Определение параметров дифференциального звена по частотным характеристикам ')

for i = 1:100

%вводим начальные данные

R = input('enter resistance, kOm: R = ');

C = input('enter capacity, mkF: C = ');

%считаем постоянную времени

T=R*C;

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

W = tf([T 0],[T 1])

%открываем файл в Simulink

open differenc.mdl

%строим АФХ

figure(1);

nyquist(W);

grid on;

hold all;

title('АФХ дифференциальных звеньев с разными постоянными времени');

%выводим легенду

r(i,1) = R;

c(i,1) = C;

L1(i,:)=['R ='];

L2(i,:)=['C ='];

Legen = [L1, num2str(r),repmat(';', size(L1,1), 1),L2, num2str(c)];

legend(Legen);

%строим ЛАЧХ и ФЧХ

figure(2);

bodeplot(W);

grid on;

hold all;

title('ЛАЧХ и ФЧХ дифференциальных звеньев с разными постоянными времени');

%выводим легенду

r(i,1) = R;

c(i,1) = C;

L1(i,:)=['R ='];

L2(i,:)=['C ='];

Legen = [L1, num2str(r),repmat(';', size(L1,1), 1),L2, num2str(c)];

legend(Legen);

%запрос о продолжении

k = input('Хотите продолжить? Y/N:', 's');

if (k == 'n') || (k == 'N'), break; end;

end