Solving a State Space Model using Runge Kutta method, Mid Point or Heun Method to view system response

Question

Anas Rehman on 28 Oct 2022

1
Link

Direct link to this question

https://se.mathworks.com/matlabcentral/answers/1837878-solving-a-state-space-model-using-runge-kutta-method-mid-point-or-heun-method-to-view-system-respon

Commented: Benjamin Pommer on 28 Feb 2023

Hi,

I was designing a project where I want to develp a function [y,t] =step_response(Gs,method) that can solve a state space model for system responce. The function would accept Transfer function and the method to be used,converts the tf into state space representation

My state space equation looks like

dx(t) = Ax(t) + Bu(t)

y(t) = Cx(t) + Du(t)

{The result will be ploted similearly to what we get from using step(sys) where sys = ss(A,B,C,D) } I need to impliment this function using 4th order Runge Kutta, Mid Point or Heun Method.

Im taking

h = 0.1; %so the results becomes obvious

u = 1; %As Im trying to replicate the step() function

Using Stepinfo to determine the time vector

SS = stepinfo(sys);

Tf = SS.SettlingTime %Should i be using the stettling time?

To = 0; %Starting time

I dont want the complete solution Just the start, how should i make a function f() which i will use

I got this from online searching and for my case should u be constant if so then what is this I and how to get x0