Tpx = 5 ((25 * т. ^ 2. + 68) * sin (t) + t * (25 * т. ^ 2. + 59) * cos (t)) / (25 * т. ^ 2. +34) ^ (3/2) Ошибка индекса в октаве - PullRequest
Новая
       72

Tpx = 5 ((25 * т. ^ 2. + 68) * sin (t) + t * (25 * т. ^ 2. + 59) * cos (t)) / (25 * т. ^ 2. +34) ^ (3/2) Ошибка индекса в октаве

0 голосов
/ 23 сентября 2019

Я пишу некоторые скрипты в Octave для определения таких вещей, как кривизна, тангенс единицы измерения и т. Д., Но я продолжаю получать эту ошибку.

error: index (67.2884): индексы должны быть либо целыми числамиОт 1 до (2 ^ 31) -1 или логическая ошибка: вызывается из Problem_2_Octave в строке 13 столбца 4 

clear;

t=pi/6; 

r=[5*t*cos(t)-3*t 5*t*sin(t)]
v=[5*cos(t)-5*t*sin(t) -3 5*sin(t)+5*t*cos(t)]
a=[-10*sin(t)-5*t*cos(t) 0 10*cos(t)-5*t*sin(t)]

magv = norm(v)

T=v/magv

Tpx=5((25*t.^2. + 68)*sin(t) + t*(25*t.^2. + 59)*cos(t))/(25*t.^2. + 34)^(3/2)
Tpy=(75*t)/sqrt((25*t^2)+34)
Tpz=((10*cos(t)-5*t*sin(t)) / (sqrt((25*t.^2) + 34))) - (25*t)*(5*sin(t)+5*t*sin(t))/((25*t.^2)+34)^(3/4)

Tp=[Tpx, Tpy, Tpz]
dottest=dot(T,Tp)

unitN = Tp/norm(Tp)
unitB = cross(T, unitN)

aT = ((25*t)/sqrt((25*t^2) + 34))
aN = magv*norm(Tp)
aA = aT*T + aN*unitN

k = (Tp/norm(v))
...