estimation of markov parameter by recursive least square method

16 次查看（过去 30 天）
Nehul Thakkar
Nehul Thakkar 2011-2-15
回答： k k d 2019-4-30
I have to estimate Markov parameter recursively. Markov Parameters are defined by h0 = D, hi = C*Ae(i-1)*B , A,B,C,D are state-space matrices.
I have input signal u = idinput(8191,'prbs'), Gs = 2/s2+3s+2 Genertae Output Y = lsim(G,u,t)
Now How can I estimate Markov Parameter????????

请先登录，再回答此问题。

回答（2 个）

Rajiv Singh
Rajiv Singh 2011-3-12
Would something like this work? z = iddata(Y,u, t(2)-t(1)); m = n4sid(z,2); h = impulse(m)
k k d
k k d 2019-4-30
% G(s)=(s+4)/(s^3+6s^2+11s+6) Example from 'Modeling and identification of
% dynamic sysytem' by Sinha and Kuszta [Page-142]
clc
clear all
s=tf('s');
num=[10 50 200 400]; % Enter numerater coefficients from highest power of s
den=[1 33 292 660 400]; % Enter denomerater coefficients from highest power of s
G=tf(num,den)
[num,den]=tfdata(G,'v');
numf=fliplr(num) ; % For Time-moment
denf=fliplr(den) ; % For Time-moment
N=7 ; % N is number of parameters required
for i=1:N
numf=[numf 0]; % For Time-moment
TimeMoments=deconv(numf,denf); % For Time-moment
num=[num 0]; % For Markov Parameters
MarkovParameter=deconv(num,den); % For Markov Parameters
end
TimeMoments
MarkovParameter

类别

在 Help Center 和 File Exchange 中查找有关 Markov Chain Models 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by