Question

我遇到一个问题，我知道这个问题非常小但它会破坏我的头脑，因为我无法看到错误的来源！我正在使用的代码（在Octave中）是：

      g = 9.814;
      Mend = 128;

      U = 15.0;

      wm = g/U;
      dw = 0.05;
      Hs = 0.16*(U*U)/g;

      sum = 0.0;
      for j = 1:Mend
      w(j) = dw + (2*j-1)*dw/2.0;
      Gw(j) = (5.0/16.0)*(wm^4/(w(j)^5))*(Hs^2)(exp(-1.25*(wm^4/(w(j)^4))));
      sum = sum + Gw(j);
      end

      plot(w,Gw);

我收到了这个错误：

 subscript indices must be either positive integers less than 2^31 or logicals

关于这一行：

Gw(j) = (5.0/16.0)*(wm^4/(w(j)^5))*(Hs^2)(exp(-1.25*(wm^4/(w(j)^4))));

我疯了，因为我看不出哪个下标索引认为它是零！

非常感谢任何帮助，并提前致谢！

Answer 1

修复！我刚刚意识到我在Hs和那个等式中的exp之间错过了一个*非常愚蠢的错误！

Octave中的指数问题

1 个答案: