Пример 5. Расчет частотных характеристик цепи с резонансом напряжений.
% Программа для исследования резонанса напряжений
% Задание вектора угловой частоты
w=50:5:1000;
% Задание параметров цепи
U=10;
R=300;
L=1;
C=0.000003; % РАСЧЕТ
% Вычисление модуля полного входного сопротивления цепи
Z=abs(R+(w.*L-1./(w.*C))*j);
% Определение тока, напряжения на L и С,
% коэффициента мощности cosf
I=U./Z;
UL=I.*w.*L;
UC=I./(w.*C);
cosf=R./Z;
%
% Построение графика полного сопротивления цепи subplot(2,2,1), plot(w,Z)
% Установка русского шрифта и его размера set(gca,'FontName','Arial Unicode MS','FontSize',10) % Название осей
xlabel('Угловая частота - w, рад/с') ylabel('Полное сопротивление - Z, Ом') % Название графика title('Зависимость Z=f(w)') % Включение отображения сетки grid on
%
% Построение графика входного тока subplot(2,2,2), plot(w,I)
% Установка русского шрифта и его размера set(gca,'FontName','Arial Unicode MS','FontSize',10) % Название осей
xlabel('Угловая частота - w, рад/с')
ylabel('Входной ток - I, А') % введение название графика title('Зависимость I=f(w)') % Включение отображения сетки grid on
%
% Построение графика напряжений на L и C subplot(2,2,3), plot(w,UL,w,UC) % Установка русского шрифта и его размера set(gca,'FontName','Arial Unicode MS','FontSize',10) % Название осей
xlabel('Угловая частота - w, рад/с') ylabel('Напряжения UL и UC, В') % введение название графика title('Зависимость Ul=f(w) и UC=f(w)') % Включение отображения сетки grid on
%
% Построение графика коэффициента мощности subplot(2,2,4), plot(w,cosf) % Установка русского шрифта и его размера set(gca,'FontName','Arial Unicode MS','FontSize',10) % Название осей
xlabel('Угловая частота - w, рад/с') ylabel('Коэффициент мощности') % введение название графика title('Зависимость Cosf=f(w)') % Включение отображения сетки grid on