PROC IML трапециевидная линейка

proc iml;
start tr(x,y); * create function called tr;

N = nrow(x);
dx = x[2:N] - x[1:N-1];
ymean = (y[2:N] + y[1:N-1]) / 2;
return(dx` * ymean );
finish tr;

x = do(-2,5,0.01);
print "Integral of a over x is" ( tr(x,0#x+1) ); *Answer is 7;

Я продолжаю получать (исполнение) недопустимый индекс или индекс вне допустимого диапазона. Как решить эту задачу и получить правильный ответ? Я пробовал убрать -1 из x[1:N-1]; и y[1:N-1], но это дает мне неправильный ответ. Это потому, что мне нужно предполагать равные интервалы? Если да, то как бы я это сделал. Уравнение трапеции: (x-x0)*(y+y0)/2 или (xi - xi-1) * (yi + yi-1)/2.


sas calculus sas-iml
person user3547583    schedule 09.05.2014    source источник

Ответы (1)


arrow_upward
0
arrow_downward

Проблема индекса возникает из-за того, что x обрабатывается как вектор-строка, а не как вектор-столбец, поэтому nrow(x) = 1. Это означает, что вы определили n = 1, поэтому, когда индекс равен n-1, он фактически равен 0. Это вызывает ошибку нижнего индекса, поскольку векторы SAS IML индексируются с 1, а не с 0. Чтобы исправить это, используйте n = ncol(x).

Кроме того, вы можете сократить свой модуль tr, если хотите:

proc iml;
    start tr(x, y);
        i = 2:ncol(x);
        return( (x[i] - x[i-1])` * (y[i] + y[i-1]) / 2 );
    finish tr;

    x = do(-2, 5, 0.01);

    print 'The integral of a over x is' (tr(x, 0#x+1)) 'and should be 7.';
quit;
person Alex A.    schedule 12.12.2014