octave-bug-tracker
[Top][All Lists]
Advanced

[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/




reply via email to

[Prev in Thread] Current Thread [Next in Thread]