MATLAB - Kaip atkurti signalo harmonikų koeficientas

[powersys]

Labas,

MATLAB kodas yra žemiau.

Prašome kreiptis PART-kode.Signalas, ty siga, yra apskaičiuota remiantis pateikta harmonikų koeficientai (ty H1, H3, ... H9).

Tada, iš dalies-B, Y = FFT (sigA.y) / L yra atliekamos.Naujus harmonikų koeficientai yra apskaičiuojami (koeficientus ženklas lemia kampo poliškumas (T)), ir išsaugoti kaip HC.

PART-C signalas iš naujo apskaičiuota remiantis harmonikų koeficientai (ty HC) apskaičiuojamas pagal PART-B.Signalas yra vadinamas sigB.

Prašome kreiptis į fig4.png žemiau.SigB yra nustatyta, kad inversija Siga, ty "out-of-faze su Siga.Ar galiu sužinoti, kodėl?Tai kaip man nustatyti koeficientą teisingas ženklas?

Kodas%%%%%%%%%%%%%%%

% KODAS

%%%%%%%%%%%%%%%aiškiai visiems; uždaryti visus; CLC;% #################################################

% PART -

% Sukurti signalas

% #################################################

% Koeficientai 1, 3, 5, 7, ir 9-ojo Harmonic

h1 = 1,00;

H3 =- 0,05;

H5 =- 0,04;

H7 =- 0,03;

H9 = 0,02;step = 1;j = 1;

for i = 0: žingsnis: 359

TR = I * pi/180;

sigA.x (j) = i;

sigA.y (j) = H1 * SIN (TR) H3 * SIN (3 * TR) H5 * sin (5 * TR) H7 * SIN (7 * TR) H9 * SIN (9 * TR);

j = j 1;

pabaigapav (1);

sklypo plotas (sigA.x, sigA.y);

% #################################################

% (END) Europos Sąjungos PART -

% #################################################% #################################################

% PART-B

% Compute FFT ir signalo sukurta Europos Sąjungos PART-.

% Priimti koeficientai 1,3,5,7 ir 9-ojo Harmonic kartą.

% #################################################

L = ilgis (sigA.y)

Y = FFT (sigA.y) / l;

Magy = ABS ;

angY = kampas (Y);paveikslas (2);

baras (2 * ABS (Y (2:20)));

% Baras (2 * abs (T (2: length / 2 1)));j = 1;

for i = 1: grindų (L / 2)

HC (j) = angY (i 1) / (abs (angY (i 1)) 1 E-100) * Magy (i 1) / Magy (2);

j = j 1;

pabaiga

% #################################################

% (END) Europos Sąjungos PART-B

% #################################################% #################################################

% Dalis-C

% Rekonstruoti signalas grindžiamas harmonikų koeficientus

% Gautas PART-B.

% #################################################

j = 1;

for i = 0: žingsnis: 359

TR = I * pi/180;

sigB.x (j) = i;

sigB.y (j) = HC (1) * SIN (tr ) ...

HC (3) * sin (3 * TR ) ...

HC (5) * sin (5 * TR ) ...

HC (7) * sin (7 * TR ) ...

HC (9) * SIN (9 * TR);

j = j 1;

pabaigapav (3);

sklypo plotas (sigB.x, sigB.y);

% #################################################

% (END) Europos Sąjungos PART-C

% #################################################% #################################################

% Dalis D

% Palyginkite signalus sukurta ne visą ir ne visą darbo B.

% #################################################

pav (4);

sklypo plotas (sigA.x, sigA.y, "r", sigB.x, sigB.y, 'b');

% #################################################

% (END) Europos Sąjungos PART-D

% #################################################

 

[MHanif]

Labas
Kodas

HC (j) = angY (i 1) / (abs (angY (i 1)) 1 E-100) * Magy (i 1) / Magy (2);