Исследование методов численного интегрирования таблично заданной функции

Министерство образования и науки Российской Федерации

Федеральное государственное бюджетное образовательное учреждение

высшего профессионального образования

«Владимирский государственный университет

имени Александра Григорьевича и Николая Григорьевича Столетовых»

Институт архитектуры, строительства и энергетики

Кафедра "Электротехника и электроэнергетика"

 Расчетно-графическая работа

по дисциплине «Математические задачи электроэнергетики»

на тему «Исследование методов численного интегрирования таблично заданной функции»

Вариант №3

Выполнил:

студент 2 курса

группы ЗЭЭу-115

Богомолов А.Б.

Проверил(а)

доц. Шмелёв В.Е.

Владимир 2017

Задача.  Пусть некоторая однофазная нагрузка питается от электронного преобразователя частоты. Формы напряжения и тока искажены по отношению к синусоиде:

t=linspace(0,50,2001); % Моменты времени, мс (здесь один период напряжения)

deu=[randn(l,40) 0]; deu(end)=deu(l);

tde=linspace(0,50,41);

u=100*sin(2*pi/50*t)+10*interp1 (tde,deu,t,'cubic');

 % Мгновенное напряжение, В

dei=[randn(l,40) 0]; deu(end)=deu(l);

ii=10*sin(2*pi/50*t-pi/6)+interpl (tde,dei,t,'cubic'); % Мгновенный ток, A

1.        Построить графики мгновенных значений напряжения, тока и потребляемой мощности.

2.        Определить постоянные составляющие, средние модули, действующие значения напряжения и тока.

3.        Определить активную и полную мощность, потребляемую нагрузкой.

Определить коэффициент мощности.

t=linspace(0,50,2001); % Моменты времени, мс (здесь один период напряжения)

deu=[randn(l,40) 0]; deu(end)=deu(l);

tde-linspace(0,50,41);

u=100*sin(2*pi/50*t)+10*interp1(tde,deu,t,'cubic'); % Мгновенное напряжение, В

dei=[randn(l,40) 0]; deu(end)=deu(l);

ii= 10*sin(2*pi/50*t-pi/6)+interp 1 (tde,dei,t,'cubic'); % Мгновенный ток, A

plot(t,u,'k-Vlinewidth',2) % Осциллограмма напряжения

х1аЬе1('Время, мс')

у1аЬе1('напряжение, В')

grid on

figure

plot(t,ii,'b-','linevvidth',2) % Осциллограмма тока

 х1аЬе1('Время, мс')

 у1аЬе1(Ток, А')

grid on

p=u.*ii; % Мгновенная мощность, Вт

 figure

plot(t,p,'r-',1inewidth',2) % Осциллограмма мгновенной мощности

 х1аЬе1('Время, мс')

 у1аЬе1('Мгновенная мощность, Вт')

grid on

dt=diff(t); % Конечная разность массива значений времени

T=t(end)-t(l); % Период изменения тока и напряжения

% Постоянная составляющая напряжения, В:

UO=dot(mean([u( 1 :end-1 );u(2:end)]),dt)/T;

% Постоянная составляющая тока, А:

IO=dot(mean([ii( 1 :end- l);ii(2:end)]),dt)/T;

% Среднее значение модуля напряжения, В:

Usr=dot(mean(abs([u( 1 :end-1 );u(2:end)])),dt)/T;

% Среднее значение модуля тока, А:

Isr=dot(mean(abs([ii(l rend-1 );ii(2:end)])),dt)/T;

% Действующее значение напряжения, В: