SCILAB-"Index invalide"

[invite3e1e58b5]

voila comme dit dans le titre, voici l'erreur qu'il me marque et malheureusement je ne trouve pas l'erreur :/

t=linspace(-1,1,500);
function [H]=Heaviside(t)
n=prod(size(t));
for i=1 : n
if t(i)>=0 then
H(i)=1;
else
H(i)=0;
end
end
endfunction
H=Heaviside(t);

[ok,h,m]=getvalue("choix du pas et du coefficient d amortissement",["h";"m"],list("vec",1,"vec",1),["0.05";"0.5"]);
t=[0:h:50]';
N=prod(size(t));
A=[0,1;-1,-1];
B=[0;1];
I=eye(2,2);
f=cos(t+%pi/3)*H(t);
Xexpl=zeros(2,N);Xtrap=zeros(2,N);Ximpl=zeros(2,N);
for n=2:N
    Xexpl(:,n)=(I+h*A)*Xexpl(:,n-1)+h*B*f(n);
    Ximpl(:,n)=inv(I-h*A)*(Ximpl(:,n-1)+h*B*f(n));
    Xtrap(:,n)=inv(I-h*A/2)*((I+h*A/2)*Xtrap(:,n-1)+h*B/2*(f(n)+f(n-1)));
end


Xexact=f-exp(-m*t).*cos(sqrt(1-m*m)*t)-m*exp(-m*t).*sin(sqrt(1-m*m)*t)/sqrt(1-m*m);
plot2d(t,[Xexpl',Ximpl',Xtrap',Xexact,H],leg="expl@impl@trap@exact");

l'erreur se situe sur la fonction f
merci d'avance =)
-----

[invite3e1e58b5]

pardon, deja Bonjour , et voici un peu plus de précisions

Je dois réaliser un programme afin de résoudre l'équation différentielle suivante : y"(t)+y'(t)+y(t)=cos(t+pi/3).H(t) (avec H(t)=heaviside)
après avoir résolue cette équation a l'aide des transformées de Laplace, j’obtiens le programme ci-dessous :

lorsque je compile celui-ci , voici ce que m'affiche scilab :

f=cos(t+%pi/3)*H(t);
                    !--error 21 
Index invalide.

merci d'avance =)

[Tryss]

H est un vecteur, donc il faudrait que t soit entier (compris entre 1 et la taille de H).

Il faudrait donc écrire f=cos(t+%pi/3)*Heaviside(t);

Mais dans Scilab, * est la multiplication matricielle : il ne va pas apprécier de multiplier deux vecteurs ensemble comme ça (dimensions incompatibles). Il faut donc utiliser la multiplication éléments par éléments .* La bonne syntaxe est alors :

f = cos(t+%pi/3) .* Heaviside(t);

[invite3e1e58b5]

alors déjà merci de votre aide =)

alors après avoir essayé ce que vous m'avez dit il me marque toujours la même erreur.
de plus lorsque je change le H(t) en Heaviside(t) ,celui ci me dit "Opération non définie pour les opérandes données."
voila merci

[A voir en vidéo sur Futura]

[Tryss]

Ceci compile sans erreur :

t=linspace(-1,1,500);
function [H]=Heaviside(t)
n=prod(size(t));
for i=1 : n
if t(i)>=0 then
H(i)=1;
else
H(i)=0;
end
end
endfunction
H=Heaviside(t);

[ok,h,m]=getvalue("choix du pas et du coefficient d amortissement",["h";"m"],list("vec",1,"vec",1),["0.05";"0.5"]);
t=[0:h:50]';
N=prod(size(t));
A=[0,1;-1,-1];
B=[0;1];
I=eye(2,2);
f=cos(t+%pi/3).*Heaviside(t);
Xexpl=zeros(2,N);Xtrap=zeros(2,N);Ximpl=zeros(2,N);
for n=2:N
    Xexpl(:,n)=(I+h*A)*Xexpl(:,n-1)+h*B*f(n);
    Ximpl(:,n)=inv(I-h*A)*(Ximpl(:,n-1)+h*B*f(n));
    Xtrap(:,n)=inv(I-h*A/2)*((I+h*A/2)*Xtrap(:,n-1)+h*B/2*(f(n)+f(n-1)));
end

Xexact=f-exp(-m*t).*cos(sqrt(1-m*m)*t)-m*exp(-m*t).*sin(sqrt(1-m*m)*t)/sqrt(1-m*m);

plot2d(t,Xexpl',leg="expl@impl@trap@exact");
plot2d(t,Ximpl',leg="expl@impl@trap@exact");
plot2d(t,Xtrap',leg="expl@impl@trap@exact");
plot2d(t,Xexact,leg="expl@impl@trap@exact");
plot2d(t,H,leg="expl@impl@trap@exact");

[invite3e1e58b5]

ah ben ... que dire merci beaucoup
tu as donc changé le heaviside et les plot
merci encore =)

ps : ne pas clôturer le sujet svp car je risque encore peut-être d'avoir besoins de votre gentille aide ^^

[invite3e1e58b5]

bon je vous remercie j'ai réussi a tout faire
sujet clos