Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Применение MathCAD в решении задач электротехники. Ч.1 Линейные электрические цепи

.pdf
Скачиваний:
0
Добавлен:
29.11.2025
Размер:
2.36 Mб
Скачать

 

 

 

 

 

 

 

 

 

ЗАДАЧА 25.

 

 

 

 

 

 

 

 

ГАРМОНИЧЕСКИЙ АНАЛИЗ ПЕРИОДИЧЕСКОЙ ФУНКЦИИ U = F(T)

Задана несинусоидальная функция напряжения u(t)

в виде графиче-

ской диаграммы в интервале полного периода T.

Требуется выполнить гар-

монический анализ заданной функции, т. е. представить ее в виде гармони-

ческого ряда Фурье. Решение задачи выполняется в 2 этапа.

 

 

 

180

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

160

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

140

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

120

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

80

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

40

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u1(t) 20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

0

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

40

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

80

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

120

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

140

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1. Аппроксимация функции

 

 

 

 

 

 

 

Заданную несинусоидальную функцию разбивают на 20-30 участков и по графической диаграмме определяют координаты точек стыка отдельных участков, результаты оформляют в виде матриц. Затем функцию аппроксимируют кубическими сплайнами по стандартной программе.

tn := ( 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 )T ×10− 3

un:=(78 148 161 125 81 50 32 21 25 43 52 12

-59 -120 -127 -93 -68 -67 -62 -10 78)T

cs := cspline (tn , un)

u(t)

:= interp (cs , tn ,un ,t)

101

2. Гармонический анализ функции

Гармонический анализ функции выполняется по классическим формулам математики, при этом постоянная составляющая определяется отдельно, а амплитуды гармоник в матричной форме:

T := .02

ω := 2

π

j :=

 

k := 1 .. 5

-1

T

 

 

 

 

 

2.1. Постоянная составляющая:

1 T

Uo := × u(t) dt = 11.102

T 0

2.2. Амплитуды и начальные фазы отдельных гармоник (до 5-ой включительно):

2 T

Umk := × u( T 0

Umk =

100.303

49.564

30.196

20.804

0.107

 

 

T

 

 

2

 

t) ×sin (k ×ω ×t) dt + j ×

×

u(t) ×cos(k ×ω ×t) dt

 

T

0

 

arg(Umk) =

19.707 ×deg

59.755

-42.163

29.178

102.331

102

 

2.3. Гармонический состав функции:

 

 

 

 

 

5

(

 

 

×sin (k ×ω ×t + arg(Umk)))

 

ug(t) := Uo +

Umk

 

 

 

 

k = 1

 

 

 

 

 

Для проверки качества выполненного гармонического анализа строят

в одной системе координат графические диаграммы исходной u(t) и рас-

четной ug(t) функций.

 

 

 

 

 

 

 

 

 

200

 

 

 

 

 

 

 

 

 

 

150

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

ug(t)

50

 

 

 

 

 

 

 

 

 

u(t)

 

 

 

 

 

 

 

 

 

 

 

0 2×10

− 3

− 3

 

− 3

 

− 3

0.01

0.012 0.014 0.016 0.018 0.02

 

4×10

6×10

8×10

 

 

 

− 50

 

 

 

 

 

 

 

 

 

 

− 100

 

 

 

 

 

 

 

 

 

 

− 150

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

Совпадение графических диаграмм этих функций указывает на верное

решение задачи.

 

 

 

 

 

 

 

 

 

103

ЗАДАЧА 26.

РАСЧЕТ ЦЕПИ НЕСИНУСОИДАЛЬНОГО ТОКА ГАРМОНИЧЕСКИМ МЕТОДОМ

Задана схема цепи и параметры отдельных элементов. Несинусоидальная функция ЭДС е(t) задана в виде таблицы координат точек в интервале полного периода T. Требуется определить действующие значения токов в ветвях схемы, напряжений на отдельных элементах, активные мощности источника и отдельных приемников энергии, проверить баланс мощностей.

 

 

 

 

 

 

 

 

i1

 

 

R1

 

L1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i2

 

 

 

 

 

 

 

i3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e(t)

 

 

 

 

L2

 

 

 

 

 

 

 

C3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2

 

 

 

 

 

 

R3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

:= 28×10− 6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

:= 21

R

:= 38

 

 

R

3

:= 26

 

L

:= .025

 

L

:= .092

 

C

3

 

 

 

1

2

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

 

 

 

 

 

tn

 

 

( 0 1 2

3

4 5 6

7

8 9 10

11

12 13

14

15 16

 

 

17 18 19 20 ) ×10− 3

 

 

 

en ( 52 122 167 138 73 28 14 15 26 38 29 −11 −61 −110 −147 −144 −82 −4 30 27 52)

Задача решается в 3 этапа

1.Гармонический анализ функции e(t)

1.1.Функция e(t) аппроксимируется кубическими сплайнами

tn := ( 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 )T×10− 3

en:=( 52 122 167 138 73 28 14 15 26 38 29 -11 -61 -110 -147 -144 -82 -4 30 27 52)T

cs := cspline (tn ,en) e(t) := interp (cs ,tn , en ,t)

104

1.2. Выполняется разложение функции e(t) в гармонический ряд Фурье

T := .02

ω := 2

π

j :=

 

 

-1

T

 

 

 

 

 

1 T

Eo := × e(t) dt = 9.978

T 0

 

 

 

 

 

 

T

 

 

 

 

 

T

 

 

 

 

 

 

2

 

 

 

 

2

 

Em

k

:=

 

×

e(t) ×sin (k ×ω ×t) dt + j ×

 

× e(t) ×cos(k

 

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

T

 

 

 

T

 

 

 

 

 

 

Emk

 

=

arg(Emk) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×deg

 

 

 

99.854

 

 

 

29.996

 

 

 

 

 

 

 

 

 

 

 

 

49.943

 

 

 

79.919

 

 

 

 

 

 

 

 

 

 

 

 

 

30.014

 

 

 

-79.525

 

 

 

 

 

 

 

 

 

 

 

 

 

19.898

 

 

 

-69.71

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9.742

 

 

 

-59.825

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k := 1 .. 5

×ω ×t) dt

5

(

 

 

 

×sin (k ×ω ×t + arg(Emk)))

eg(t) := Eo +

 

Emk

 

 

 

 

 

 

 

 

k = 1

 

 

 

 

 

2.Расчет токов и напряжений в схеме для отдельных гармоник

2.1.Расчет схемы для постоянной составляющей:

Ro := R

1

+ R

2

I1o :=

Eo

= 0.169

I2o := I1o

I3o := 0

 

 

 

 

 

 

 

 

Ro

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U1o := I1o×R1 = 3.552

U2o := I2o×R2 = 6.427

 

 

 

 

 

PEo := Eo×I1o = 1.688

P1o := U1o×I1o = 0.601

 

 

 

 

P2o := U2o×I2o = 1.087

P3o := 0

 

 

 

 

 

 

 

 

 

2.2. Расчет схемы для k-ой гармоники

 

 

 

 

Z1k := R1 + j ×k ×ω ×L1

 

Z2k := R2 + j ×k ×ω ×L2

Z3k :=

R3

- j ×

 

1

 

 

 

 

k ×ω ×C3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Zabk := Z2k ×

Z3k

k := Z1k + Zabk

 

 

 

Z2

+ Z3

 

 

 

 

 

 

 

 

 

k

 

k

 

 

 

 

 

 

105

I1mk

:=

Emk

 

 

I1mk

 

 

=

 

ψi1 k := arg(I1mk)

ψi1 k =

 

 

 

 

 

 

k

 

 

 

 

 

×deg

 

 

 

1.337

 

 

 

 

 

 

14.686

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.569

 

 

 

 

 

 

77.018

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.433

 

 

 

 

 

 

-79.035

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.334

 

 

 

 

 

 

-79.466

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.167

 

 

 

 

 

 

-81.56

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U1mk := I1mk ×Z1k

 

 

U1mk

 

 

=

ψu1 k := arg(U1mk)

 

ψu1 k =

×deg

 

 

 

 

 

 

 

 

 

29.969

 

 

 

35.192

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14.931

 

 

 

113.815

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13.659

 

 

 

-30.745

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12.616

 

 

 

-23.226

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7.453

 

 

 

-19.696

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U2mk := I1mk ×Zabk

 

U2mk

=

 

 

 

 

 

 

70.06

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

38.462

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

23.391

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14.47

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6.279

 

 

 

U2mk

 

 

 

 

 

 

 

I2mk

:=

I2mk

 

=

 

 

 

 

 

Z2k

 

 

 

 

 

 

1.467

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.556

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.247

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.119

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.042

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψu2 k := arg(U2mk)

ψu2 k =

×deg

 

27.776

 

 

 

 

67.415

 

 

 

 

 

-105.58

 

 

 

 

 

-108.928

 

 

 

 

 

-109.737

 

ψi2 k := arg(I2mk)

 

 

ψi2 k =

×deg

 

-9.481

 

 

 

 

10.735

 

 

 

 

 

-171.914

 

 

 

 

 

179.267

 

 

 

 

 

174.996

 

 

 

 

I3mk :=

U2mk

 

 

I3mk

 

=

ψi3 k := arg(I3mk)

ψi3 k =

 

 

 

 

 

 

Z3k

×deg

 

 

0.601

 

 

104.893

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.615

 

 

132.835

 

 

 

 

 

 

 

 

 

 

 

 

0.509

 

 

-50.034

 

 

 

 

 

 

 

 

 

 

 

 

0.376

 

 

-61.381

 

 

 

 

 

 

 

 

 

 

 

 

0.182

 

 

-68.568

 

 

 

 

 

 

 

 

 

 

 

106

3.Обработка результатов расчета

3.1.Мгновенные значения токов и напряжений, как функции времени

 

 

 

 

5

(

 

 

 

×sin (k ×ω ×t + ψi1 k))

 

 

 

i1 (t) := I1o +

I1mk

 

 

 

 

 

k = 1

 

 

 

 

 

 

 

 

 

 

 

 

5

(

 

 

 

×sin (k ×ω ×t + ψi2 k))

 

 

 

i2 (t) := I2o +

I2mk

 

 

 

 

 

k = 1

 

 

 

 

 

 

 

 

 

 

 

 

5

(

 

 

 

×sin (k ×ω ×t + ψi3 k))

 

 

 

i3 (t) := I3o +

I3mk

 

 

 

 

 

k = 1

 

 

 

 

 

 

 

 

 

 

 

 

5

 

(

 

 

 

×sin (k ×ω ×t + ψu1 k))

 

 

 

u1(t) := U1o +

U1mk

 

 

 

 

 

 

k = 1

 

 

 

 

 

 

 

 

 

 

 

5

 

(

 

 

 

×sin (k ×ω ×t + ψu2 k))

 

 

 

u2(t) := U2o +

U2mk

 

 

 

 

 

 

k = 1

 

 

 

 

 

 

 

 

3.2. Графические диаграммы функций напряжений

 

 

 

200

 

 

 

 

 

 

 

 

 

 

 

 

150

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

 

 

e(t)

50

 

 

 

 

 

 

 

 

 

 

 

u1(t)

 

− 3

− 3

− 3

 

 

− 3

 

 

 

 

 

0 2×10

 

 

0.01 0.012 0.014 0.016

0.018

0.02

u2(t)

4×10

6×10

8×10

 

− 50

 

 

 

 

 

 

 

 

 

 

 

 

− 100

 

 

 

 

 

 

 

 

 

 

 

 

− 150

 

 

 

 

 

 

 

 

 

 

 

 

− 200

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

107

 

3.3. Графические диаграммы функций токов

 

 

 

2.5

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

1.5

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

i1(t)

0.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i2(t)

 

 

− 3

− 3

− 3

− 3

 

 

 

i3(t) − 0.5 0

 

 

 

 

2×10

4×10

6×10

8×10

 

0.01 0.012 0.014 0.016

0.018

0.02

 

− 1

 

 

 

 

 

 

 

 

 

− 1.5

 

 

 

 

 

 

 

 

 

− 2

 

 

 

 

 

 

 

 

 

− 2.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

3.4. Действующие значения токов и напряжений:

 

 

 

2

 

5

 

(

 

 

 

)2 = 1.117

 

 

 

 

 

 

 

I1 :=

(I1o)2

+

1

×

k = 1

 

 

I1mk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

5

 

(

 

 

 

)2 = 1.139

 

 

 

 

 

 

 

I2 :=

(I2o)2

+

1

×

k = 1

 

 

I2mk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

5

 

(

 

 

 

)2 = 0.766

 

 

 

 

 

 

 

I3 :=

(I3o)2

+

1

×

k = 1

 

 

I3mk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

5

(

 

 

 

 

)2 = 83.833

 

 

 

 

 

 

 

E :=

(Eo)2 +

1

×

 

 

 

 

Emk

 

 

 

 

 

k = 1

108

 

 

1

5

(

U1 :=

(U1o)2 +

×

2

 

 

 

k = 1

 

 

 

1

5

(

U2 :=

(U2o)2 +

×

2

 

 

 

k = 1

 

U1mk )2 = 27.817

U2mk )2 = 60.276

3.5. Активная мощность источника и приемников энергии

 

 

5

PEk := Re(.5 ×Emk ×I1mk)

PE := PEo +

PEk = 90.78

 

 

k = 1

 

 

5

P1k := Re(.5 ×U1mk ×I1mk)

P1 := P1o +

P1k = 26.196

 

 

k = 1

 

 

5

P2k := Re(.5 ×U2mk ×I2mk)

P2 := P2o +

P2k = 49.338

 

 

k = 1

 

 

5

P3k := Re(.5 ×U2mk ×I3mk)

P3 := P3o + P3k = 15.246

k = 1

Pn := P1 + P2 + P3 = 90.78

Вывод: баланс мощностей выполняется.

109

ЗАДАЧА 27.

РАСЧЕТ ЦЕПИ НЕСИНУСОИДАЛЬНОГО ТОКА ДИФФЕРЕНЦИАЛЬНЫМ МЕТОДОМ

Задана схема цепи и параметры отдельных элементов. Несинусоидальная функция ЭДС е(t) задана в виде таблицы координат точек в интервале полного периода T. Требуется определить действующие значения токов в ветвях схемы, напряжений на отдельных элементах, активные мощности источника и отдельных приемников энергии, проверить баланс мощностей.

 

 

 

i1

R1

L1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i2

 

 

 

 

 

 

i3

 

 

 

e(t)

 

 

 

 

L2

 

 

 

 

 

 

C3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2

 

 

 

 

 

R3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T := .02

R1 := 21

 

R2 := 38

R3 := 26

 

L1 := .025

 

L2 := .092

C3 := 28×10− 6

Задача решается в 3 этапа 1. Аппроксимация функции в интервале нескольких периодов

Координаты точек заданной функции е(t) в интервале 1-го периода оформляются в виде столбцовых матриц, затем производится сшивание аналогичных матриц для последующих периодов, после чего создается кубический сплайн и производится аппроксимация.

tn1 := 10− 3×( 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 )T en1:=( 52 122 167 138 73 28 14 15 26 38 29 -11 -61 -110 -147 -144 -82 -4 30 27)T

tn2 := tn1 + T tn3 := tn1 + 2×T tn4 := tn1 + 3×T tn5 := tn1 + 4×T

tt := stack (tn1 ,tn2 ,tn3 ,tn4 ,tn5)

ee := stack (en1 ,en1 ,en1 ,en1 ,en1)

cs := cspline (tt ,ee)

e(t) := interp (cs ,tt ,ee ,t)

110

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]