Матричные процедуры. Процедура преобразования блочной матрицы в обычную. Процедура преобразования блочного вектора в обычный

ПРИЛОЖЕНИЕ 1риложение 2

Матричные процедуры

П.21.1.ПРОЦЕДУРА ПРЕОБРАЗОВАНИЕЯ БЛОЧНОЙ МАТРИЦЫ В ОБЫЧНУЮ

function [R]=CellArrayToMatrix(I)

[m n]=size(I);          Определяем размерность  блочной матрицы

[v c]=size(I{1,1});    Определяем размерность блока матрицы

R=zeros(m*v, n*c);         Обнуляем матрицу-результат

for i = 1:m                  Цикл по строкам блочной матрицы

for j = 1:n             Цикл по столбцам блочной матрицы

for k = 1:v           Цикл по строкам блоков

for l = 1:c         Цикл по столбцам блоков

if (I{i,j}==0);        Если элементы блока равны нулю, то

R(i*v-v+k,j*c-c+l)=0;       элементы матрицы R, адрес которых вычислен по данной формуле, будут равны нулю

else              иначе

R(i*v-v+k,j*c-c+l)=I{i,j}(k,l);            эти элементы будут равны элементам блоков матрицы блочной

end

end

end

end

end

П.21.2.ПРОЦЕДУРА ПРОБРАЗОВАНИЕЯ БЛОЧНОГО ВЕКТОРА В ОБЫЧНЫЙ

function [R]=CellArrayToVector(I)

r=size(I,1);              Определяем размерность  блочного вектора

s=size(I{1,1},1);      Определяем размерность блока вектора

R = zeros(r*s,1);             Обнуляем матрицу-результат

for i = 1:r                Цикл по строкам блочного вектора

for k = 1:s           Цикл по строкам блока

if (I{i}==0);  Если блок равен нулю, то

R(i*r-r+k)=0;  Элементы вектора R, адрес которых вычислен по данной формуле, будут равны нулю

else         иначе

R(i*r-r+k)=I{i}(k);    эти элементы будут равны элементам блоков блочного вектора

end

end

end

П.12.3. ПРОЦЕДУРА, ОСУЩЕСТВЛЯЮЩАЯ ШШАГ ЖОРДАНА ДЛЯ БЛОЧНОЙ 3х3 МАТРИЦЫ 3х3 УЗЛОВЫХ
 ПРОВОДИМОСТЕЙ С (3,3) – ЭЛЕМЕНТОМ (3,3) В КАЧЕСТВЕ ОПОРНОГО.

function y=Jord(Y)

y{3,3}=(Y{3,3})^-1;           Вычисляем обратную величину главного(опорного) элемента

for i=1:2             Перебор элементов главного столбца и главной строки

y{i,3}=Y{i,3}*y{3,3};        Прочие элементы главного столбца умножаем справа на обратную величину главного элемента

y{3,i}=-y{3,3}*Y{3,i};       Прочие элементы главной строки умножаются слева на обратную величину главного элемента и меняют свой знак на противоположный

for j=1:2                  Перебор остальных элементов

y{i,j}=Y{i,j}-Y{i,3}*y{3,3}*Y{3,j};        По данной формуле вычисляются прочие элементы, не принадлежащие ни к главному столбцу, ни к главной строке

end

end

П.21.4.ПРОЦЕДЕРВ ФОРМИРОВАНИЕ Я МАТРИЦЫ УЗЛОВЫХ ПРОВОДИМОСТЕЙ

function [Yyy]=SintezY(l,Z0,Y0,L,Ykz,num);

num – номер элемента топологического вектора, характеризующего тип КЗ

1 - Замыкание фазы А на землю

2 - Замыкание фазы В на землю

3 - Замыкание фазы С на землю

4 - Двухфазное замыкание АВ

5 - Двухфазное замыкание ВС

6 - Двухфазное замыкание АС

7 - Трёхфазное замыкание

8 - Трёхфазное замыкание на землю

9- Двухфазное замыкание АВ на землю

10 - Двухфазное замыкание ВС на землю

11 - Двухфазное замыкание АС на землю

FMY – топологический вектор, характеризующий тип КЗ

Y1k=(Z0*l)^-1;

Yk2=(Z0*(L-l))^-1;

Y11=Y0*l/2;                    формирование матрицы Yyy

Y22=Y0*(L-1)/2;

Yk=FMY{num}*Ykz;

Матрица узловых проводимостей

 

 

ПРОЦЕДУРА П.21.5. ФОРМИРОВАНИЕ ГЛОБАЛЬНЫХ КОНСТАНТ

globalFpsМатрица перехода от фазных составляющих к симметричным

globalFspМатрица перехода от симметричных составляющих к фазным

          Использован классический вариант с неунитарной матрицей, искажающий мощности. Коли так, то корня из трёх быть не должно!

 

globalf_nomНоминальная частота, Гц

globalw_nomЧастота, рад/сек

globalgradКоличество радиан в одном градусе

a=exp(2j*pi/3);

Fsp=Fps^-1;

f_nom=50;

w_nom=2*pi*f_nom;

grad=180/pi;