matlab - State space system gives different bode plot then transfer function matrix -

-

I have a state space system with matrix A, B, C and D.

I have either the sys2 = c * inv (z * I - A) * B + D , sys1 = ss (A, B, C, D) code >

However, when I draw the board plot of both systems, they are different, when they are the same.

What's going wrong here? Does anyone have a clue? I know that bodeplot generated by sys1 is correct.

The system can be downloaded here: 

  clear all; close all; CLC; Ts = 0.01; Z = TF ('Z', TS); % Discrete system A = [0 0 0; 0 0 1; 0.41 -1.21 1.8]; B = [0; 0; 0.01]; C = [7 -73 170]; D = 1; Set as state% sys1 = ss (A, B, C, D, SS); % Compute Transfer function sys2 = C * inv (z * eye (3) - a) * B + D; Calculate% actual transfer function [num, den] = ss2tf (a, b, c, d); Sys3 = tf (num, den, ts); Show% bode bode (sys1, 'b', sys2, 'r -', sys3, 'g--'); Edit: I made a small mistake, the transfer function matrix  sys2 = C * inv (z * I - A) * B + D  , Instead of  sys2 = C * inv (z * I - A) * b - d , which I did before. The problem still holds. Edit 2: I have indicated that when I calculate all the way, it is correct. 
 
  sys z; Collected ((z * eye (3) - a), z)

its notion Is that sys2 = C * inv (z * I-A) * B + D is incorrect. Your state-space system (A, B, C, D) is equal to sys2 = C * inv (S * I-A) * B + D . If you want to express it in the context of z , then you have to reverse the relationship z = exp (s * T) . sys1 is the correct representation of your state-space system What do I suggest for sys2 : 

  sys1 = ss (mjlsCE.A , MjlsCE.B, mjlsCE.C, mjlsCE.D, ts); Sys1_c = d2c (sys1); S = TF ('s'); Sys2_c = sys1_c.C * inv (s * eye (length (sys1_c.A)) - sys1_c.A) * sys1_c.b + sys1_c.D; Sys2_d = c2d (sys2_c, ts);

You should give the right result.


Comments

Post a Comment

Popular posts from this blog

mysql - How to enter php data into a html multiple select box -

-
I am working on a page that will create a team. I just want to create a selection box which will list all the users of all the users. After that, the user can select more than one friend (i.e. 5) and can send the request to these 5 people to join the team but I am having problems I did this: Select $ sql = "Select i. *, Join * Connect with M * m.id = i.clique_id where adder_id =: id"; $ Stmt = $ db- & gt; Query ($ sql, array ('id' = & gt; $ _SESSION ['id'])); While ($ record = $ stmt-> fetch ()) {$ surname = $ record ['surname']; } and this: & lt; Select multiple = "true" name = "category" well id = "category" name = "category" category = "" & gt; & Lt; Option value = "& lt ;? php echo $ surname ;? & gt;" & Gt; & Lt ;? Php ek $ nickname; ? & Gt; & Lt; / Options & gt; & Lt; / Select & gt; But only a friend's nickname...
Read more

java - Can't add JTree to JPanel of a JInternalFrame -

-
Image
There are two panels in my JInternalframe. I want to add TopPanel to Jtree to nominate another. But I can not add Jaitari to the top of the panel. Please help me in this piece of code: DefaultMutableTreeNode root = new DefaultMutableTreeNode ("deck"); DefaultMutableTreeNode Item Clubs = New DefaultMutableTreeNode ("Club"); AddAllCard (itemClubs); Root.add (itemClubs); DefaultMutableTreeNode item = New DefaultMutableTreeNode ("Diamonds"); AddAllCard (itemDiamonds); Root.add (itemDiamonds); DefaultMutableTreeNode Itemspads = New DefaultMutableTreeNode ("Spades"); AddAllCard (itemSpades); Root.add (itemSpades); DefaultMutableTreeNode item increases = new DefaultMutableTreeNode ("heart"); AddAllCard (itemHearts); Root.add (itemHearts); Default Trimodel Tree Model = New DefaultTreeModel (Route); Tree = new jetty (tree modal); ScrollPane = new JScrollPane (tree); // scrollPane.setViewportView (tree); . ScrollPane.getViewport () add (tr...
Read more

c++ - Cassandra datastax cpp driver - avoiding unnecessary copies -

-
I use the following method to create my insert statement for Cassandra (using datastax cpp driver) : ssStmt1 & lt; & Lt; "Insert some table (...) value (....)" string sStatement1 = ssStmt1.str (); Castestring cass1 = cass_string_init (sStatement1.c_str ()); CassStatement * pStmt1 = cass_statement_new (cass1, 0); I have also tried the following: CassString cass1 = cass_string_init ("Enter some table (id, date, c) value (? ,?,?); "); CassStatement * pStmt1 = cass_statement_new (cass1, 3); Cass_statement_bind_uuid (pStmt1, 0, uuidKey); Cass_statement_bind_int64 (PSTMT1, 1, timestamp); Cass_statement_bind_string (pStmt1, 2, sSomething); The first pattern works with mixed keys, the second is not. The key has a timestamp and a UUID. Timestamp is one of two errors in the field (based on the used bond statement). If bound in this way: cass_statement_bind_int64_by_name (pStmt1, "date", timestamp); I get the following error: Inva...
Read more