Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Categories

Welcome to the new platform of Programmer's Heaven! We apologize for the inconvenience caused, if you visited us from a broken link of the previous version. The main reason to move to a new platform is to provide more effective and collaborative experience to you all. Please feel free to experience the new platform and use its exciting features. Contact us for any issue that you need to get clarified. We are more than happy to help you.

Need help in matlab system identification for mimo system.

joe_mikaeljoe_mikael Posts: 2Member
given heater as the input and rel.humidity and temp. as the outputs.

To find this coefficients of this mimo system.This is my coding so far and i get confused .what should i find next?i also attached the .mat file .Please help me.May God bless you all. :)

[code]temp = A(:,1);rh = A(:,2); %temp % rh as the output
heater = A(:,3); %heater as the input
y =[temp rh];
u = heater;
trn_data_n = 100;
total_data_n = 218;
yn = dtrend(y);
un = dtrend(u);
ave_out = mean(y);
ave_in = mean(u);
T = 5; %minutes
hvac = iddata(yn,un,T);
hvac.InputName = {'Heater'};
hvac.OutputName = {'Temperature';'Rel.Humidity'};
plot(hvac(:,1,1)) % temp. vs. heater
pause
plot(hvac(:,2,1)) % rh vs. heater
ms = step(hvac); % step responses
step(ms)
impulse(ms,'sd',3) % impulse response with confidence regions corresponding to 3 standard deviations
%state space model
mp = pem(hvac(1:trn_data_n)); % Use only the first half of the data for estimation:
A=mp.A; B=mp.B; C=mp.C; D=mp.D;
[nem, den]=ss2tf(A,B,C,D,1); %convert ss to transfer function
[/code]
Sign In or Register to comment.