14.09.2014 Views

Impact of fuel supply impedance and fuel staging on gas turbine ...

Impact of fuel supply impedance and fuel staging on gas turbine ...

Impact of fuel supply impedance and fuel staging on gas turbine ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

System Identificati<strong>on</strong><br />

The sec<strong>on</strong>d term <strong>on</strong> the right h<str<strong>on</strong>g>and</str<strong>on</strong>g> side can be estimated as the crosscorrelati<strong>on</strong><br />

vector <str<strong>on</strong>g>of</str<strong>on</strong>g> the input <str<strong>on</strong>g>and</str<strong>on</strong>g> output signals:<br />

1<br />

C x y =<br />

N− M+ 1<br />

N∑<br />

n=M<br />

x (i )<br />

n y n−k; k = 0,..., M; i = 1,2,3. (5.11)<br />

Compared to the basic CFD/SI method, which identifies SISO <str<strong>on</strong>g>and</str<strong>on</strong>g> MIMO<br />

models (see [34, 98, 102]), the auto-correlati<strong>on</strong> matrix c<strong>on</strong>tains now as well<br />

n<strong>on</strong>-diag<strong>on</strong>al entries as the input signals are correlated. This fact requires a<br />

discussi<strong>on</strong> about the choice <str<strong>on</strong>g>of</str<strong>on</strong>g> excitati<strong>on</strong> signals <str<strong>on</strong>g>and</str<strong>on</strong>g> the validati<strong>on</strong> methods<br />

to judge the quality <str<strong>on</strong>g>of</str<strong>on</strong>g> the results obtained. Equati<strong>on</strong> (5.9) can be rewritten as<br />

follows<br />

ˆθ= Γ −1<br />

xx C x y, (5.12)<br />

which is the optimal linear least square estimator for the unit impulse resp<strong>on</strong>se.<br />

It is also known as the Wiener-Hopf equati<strong>on</strong> or rather the Wiener-<br />

Hopf inversi<strong>on</strong> [40, 52, 73]. After the auto-correlati<strong>on</strong>s <str<strong>on</strong>g>and</str<strong>on</strong>g> cross-correlati<strong>on</strong>s<br />

are determined, the unknown parameters are solved using a method which is<br />

based <strong>on</strong> a least-square root algorithm for sparse linear equati<strong>on</strong>s [86].<br />

5.3 Excitati<strong>on</strong> signals<br />

The identificati<strong>on</strong> method presented above processes data from simulati<strong>on</strong>s<br />

with broadb<str<strong>on</strong>g>and</str<strong>on</strong>g> excitati<strong>on</strong>. Several types <str<strong>on</strong>g>of</str<strong>on</strong>g> broadb<str<strong>on</strong>g>and</str<strong>on</strong>g> excitati<strong>on</strong> signals are<br />

known, <str<strong>on</strong>g>and</str<strong>on</strong>g> it is worthwhile to investigate how the signal type influences the<br />

results <str<strong>on</strong>g>of</str<strong>on</strong>g> the identificati<strong>on</strong>.<br />

Insufficient excitati<strong>on</strong> strength will yield poor identificati<strong>on</strong> results, especially<br />

in the presence <str<strong>on</strong>g>of</str<strong>on</strong>g> noise. On the other h<str<strong>on</strong>g>and</str<strong>on</strong>g>, because the identificati<strong>on</strong> procedure<br />

is limited to linear systems, it has to be assured that maximum signal amplitudes<br />

do not generate n<strong>on</strong>-linear behavior. A maximum utilizati<strong>on</strong> <str<strong>on</strong>g>of</str<strong>on</strong>g> the<br />

amplitude limit over a wide range <str<strong>on</strong>g>of</str<strong>on</strong>g> frequencies is therefore beneficial. This<br />

can be judged in analyzing the so-called crest factor or peak-to-average ratio,<br />

proposed by Ljung [73]. The crest factor is a property <str<strong>on</strong>g>of</str<strong>on</strong>g> the waveform <str<strong>on</strong>g>and</str<strong>on</strong>g> is<br />

equal to the peak amplitude <str<strong>on</strong>g>of</str<strong>on</strong>g> the waveform divided by the mean square <str<strong>on</strong>g>of</str<strong>on</strong>g><br />

86

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!