Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- A=[0,1;-3,-7];
- B=[0;1];
- C=[1,0;0,1];
- D=[0;0];
- h=0.1;
- WPC=0;
- syms s;
- Ad=expm(h*A);
- Bd=double(int(expm(s*A)*B,s,0,h));
- Cd=C;
- open('modellab3');
- sim('modellab3');
- figure;
- plot(tout,x);
- %Jan Zasadny
Add Comment
Please, Sign In to add comment
