why do I get --- ??? Subscript indices must either be real positive integers or logicals.----- on my code

Question

Mariela Flores 2017-10-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/364180-why-do-i-get-subscript-indices-must-either-be-real-positive-integers-or-logicals

编辑： KSSV 2017-10-31

this is my code:

    % code
  e
clc;
m=1;
c=0;
k=10;
wn=sqrt(k/m);
xi=c/(2*sqrt(k*m));
t0=0;
tf=5;
muestras=1000;
for i=1:muestras+1
    t(i)=t0+(i-1)*(tf-t0)/muestras;
    if t(i)<=tf
    x(i)=1/m((t(i)*cos(t(i)))*(((exp(-wn*t(i)))*(sin(t(i))-wn*cos(t(i)))+wn)/wn^2+1)+t(i)*(sin(t(i)))*(((-exp(-wn*t(i)))*(wn*(sin(t(i)))+(cos(t(i)))/wn^2+1)+(1/(wn^2+1))))-(1/(wn^2+1)^3)*((-exp(-wn*t(i))*-(wn^5)*(t(i)^2)+2*(wn^4)*t(i)^1+2*(wn^3)*(t(i)^2+1)+wn*(t(i)^2-6)-(2*wn^3-6*wn)*(cos(t(i))))));
end 
end
plot(t,x,'r')
xlabel('t')
ylabel('x(t)')
grid
hold on

the display shows the next error: ----??? Subscript indices must either be real positive integers or logicals.

Error in ==> EXACTOprueba at 15 x(i)=1/m((t(i)*cos(t(i)))*(((exp(-wn*t(i)))*(sin(t(i))-wn*cos(t(i)))+wn)/wn^2+1)+t(i)*(sin(t(i)))*(((-exp(-wn*t(i)))*(wn*(sin(t(i)))+(cos(t(i)))/wn^2+1)+(1/(wn^2+1))))-(1/(wn^2+1)^3)*((-exp(-wn*t >>

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

请先登录，再进行评论。

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

Answer 1

KSSV 2017-10-31

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/364180-why-do-i-get-subscript-indices-must-either-be-real-positive-integers-or-logicals#answer_288602

编辑：KSSV 2017-10-31

在 MATLAB Online 中打开

     % code
    e
  clc;
  m=1;
  c=0;
  k=10;
  wn=sqrt(k/m);
  xi=c/(2*sqrt(k*m));
  t0=0;
  tf=5;
  muestras=1000;
 t = zeros(muestras+1,1) ;
  for i=1:muestras+1
      t(i)=t0+(i-1)*(tf-t0)/muestras;
      if t(i)<=tf
      x(i)=1/m*((t(i)*cos(t(i)))*(((exp(-wn*t(i)))*(sin(t(i))-wn*cos(t(i)))+wn)/wn^2+1)+t(i)*(sin(t(i)))*(((-exp(-wn*t(i)))*(wn*(sin(t(i)))+(cos(t(i)))/wn^2+1)+(1/(wn^2+1))))-(1/(wn^2+1)^3)*((-exp(-wn*t(i))*-(wn^5)*(t(i)^2)+2*(wn^4)*t(i)^1+2*(wn^3)*(t(i)^2+1)+wn*(t(i)^2-6)-(2*wn^3-6*wn)*(cos(t(i))))));
  end 
  end
  plot(t,x,'r')
  xlabel('t')
  ylabel('x(t)')
  grid
  hold on