1

\int_{-\infty}^a \int_{-infty}^b sum_{i,j}^K a_i*a_j*x^i*y^j*exp(- x^2 - y^2 + x*y)dx dy

a_i と a_j は定数です。積分は線形なので、和と積分を入れ替えることができますが、この場合は K^2 積分を評価する必要があり、時間がかかりすぎます。その場合、私は次のことを行います。

for i = 1:K
   for j = 1:K
       fun = @(x,y) x.^i.*y.^j.*exp(-2.*(x.^2 + y.^2 - 2.*x.*y))
       part(i,j) = alpha(i)*alpha(j)*integral2(fun,-inf,a,-inf,b)
   end
end

時間がかかりすぎるので積分を一つだけ評価したいのですが、sum_{i,j}^K a_i*a_j*x^i*y^j*exp(-x^2 - y^2 + x*y)、つまり、integral2 への供給方法です。| | どんな助けにも非常に感謝します。

4

1 に答える 1

0

どうですか:

I = [1:K]';
J = [1:K]';
funhandle = @(x,y)(sum(sum(...
repmat(alpha(I),1,K)...
.*repmat(alpha(J)',K,1)...
.*repmat(x.^I,1,K)...
.*repmat((y.^J)',K,1)))...
*exp(-x^2-y^2+x*y))

私はそれがあなたが正しく望む合計をするだろうと思います...

于 2013-03-22T10:10:37.617 に答える