[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #65401] [octave forge] (control) A problem wit
From: |
Ahmed Salah Hammad Ahmed |
Subject: |
[Octave-bug-tracker] [bug #65401] [octave forge] (control) A problem with lsim command |
Date: |
Sat, 2 Mar 2024 14:01:07 -0500 (EST) |
Follow-up Comment #3, bug #65401 (group octave):
[comment #1 comment #1:]
> Can you please post your exact inputs so we can try to exactly recreate your
issue?
The first error message appears for this code (Sorry for the mess with how
matrix A looks..):
A = [0 1 0 0 0
0 0 0 0;
-4*w0^2*sx -Dx/Ix Dx/Ix 0 0
0 w0*(1-sx) 0 0;
0 Dx/Iwx -Dx/Iwx 0 0
0 0 0 0;
0 0 0 0 1
0 0 0 0;
0 0 0 -3*w0^2*sy -Dy/Iy
Dy/Iy 0 0 0;
0 0 0 0 Dy/Iwy
-Dy/Iwy 0 0 0;
0 0 0 0 0
0 0 1 0;
0 -w0*(1-sz) 0 0 0
0 -w0^2*sz -Dz/Iz Dz/Iz;
0 0 0 0 0
0 0 Dz/Iwz -Dz/Iwz;
];
B = [0 0 0; 1/Ix 0 0; 0 0 0; 0 0 0; 0 1/Iy 0; 0 0 0; 0 0 0; 0 0 1/Iz; 0 0 0];
C = eye(9);
D = zeros(9, 3);
x0 = [0 0 0 0 0 0 0 0 0].'; % initial conditions
sys = ss(A, B, C, D);
t = 0:.2:50;
u = 1 * ones(3,length(t)); % input
% u = [zeros(1, 100) 1000, 0*ones(1, numel(t)-100-1)] .* [1; 1; 1];
[Y, X] = lsim(sys, u, t, x0);
Changing the time vector into:
t = (0:.2:5).*[1;1;1];
produces the second error message:
time vector 't' must be real-valued or empty
_______________________________________________________
Reply to this item at:
<https://savannah.gnu.org/bugs/?65401>
_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/