0 
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;

(実行) 無効な添え字または範囲外の添え字を受け取り続けます。この問題を解決して正しい答えを得るにはどうすればよいですか? x[1:N-1]; の -1 を取り出してみました。と y[1:N-1] ですが、間違った答えが得られます。等間隔の間隔を想定する必要があるからですか？もしそうなら、どうすればいいですか？台形方程式: (x-x0)*(y+y0)/2 または (xi - xi-1) * (yi + yi-1) / 2.

4

1 に答える 1

0

下付き文字の問題は、xが列ベクトルではなく行ベクトルとして扱われているという事実から生じるため、nrow(x) = 1. SAS IML ベクトルは 0 ではなく 1 からインデックス付けされるため、添字エラーが発生します。これを修正するには、 を使用n = 1します。n-1n = 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;
于 2014-12-12T16:43:37.660 に答える