Howdy, Stranger!

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

Sign In with Facebook Sign In with Google Sign In with OpenID


We have migrated to a new platform! Please note that you will need to reset your password to log in (your credentials are still in-tact though). Please contact if you have questions.
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
plot(hvac(:,2,1)) % rh vs. heater
ms = step(hvac); % step responses
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
Sign In or Register to comment.