Solving and Plotting Initial Value Response using EXPM() function

Nick Alden

2017 10 2

1 个回答

3 次查看（30 天）

0 个投票

How can I make matrix dimensions work when trying to solve and plot Initial Value Response to State Space model and using expm() function?

Code:

clear all
close all
clc
syms t
time = [0:0.01:10]
A = [0 1 ; -5 -2]
B = [0;1]
x_naut = [1;-1]
[eig_vector, eig_value] = eig(A)
v_inv = inv(eig_vector)
eAt = expm(A.*time)
x = expm(A).*x_naut
fplot(time,x)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Follow Question

回答（1 个）

Star Strider 2017-10-2

在 MATLAB Online 中打开

0 个投票

Try this:

time = [0:0.01:10];
A = [0 1 ; -5 -2];
B = [0;1];
x_naut = [1;-1];
for k1 = 1:numel(time)
    eAt = expm(A.*time(k1));
    x(:,k1) = [1 0; 0 1]*eAt*x_naut;
end
figure(1)
plot(time,x)
grid