numerical-integration

2熱

1答えて

私は数値積分に関するプレゼンテーションの一部を行っています。トーク自体は数値積分（主に重要サンプリングと層別サンプリング）のより良い形になりますが、私は、モンテカルロ積分のサンプリングを一様分布からサンプリングしています。私がいることを発見しました： mean(sin(runif(1e8, 0, pi))) がかなり期待されている1より、0.636597の答えを与えています。この回答はサンプ

0熱

1答えて

Javaでの数値二変量積分

Javaで数値二変量積分（2つの変数を持つ関数の統合）を実行できますか？私はHow to solve an Integral in Java?を見てきましたが、単変量積分（1変数）についてのみ説明しています。

3熱

1答えて

ベクトル化Matlab - ループなしでガウス関数（コード）をベクトル化する方法

ガウスの直交ガードを使用して積分を計算する関数を作成しました。ループを使わずに最後の2行をベクトル化する必要があります。何か案は？ function Q = gauss5(f,a,b,n) % Divide the interval into n parts y = linspace(a,b,n); % Create function to use arrayfun

1熱

1答えて

浮動小数点除算（Python）の特定と修復

私は多項式に誤差を当てはめることによって台形則を使用した数値積分でエラーを見つけるためにスクリプトで使用されるコードセグメントを持っています。このコードセグメントは、浮動小数点除算ゼロエラーを投げているので、なぜそれを修正するかわかりません。誰かが私に答えを導くのに役立つでしょうか？ def trap(f,a,b,dx,exact): N = int(numpy.round(float

0熱

1答えて

フォワード - オイラー法

私は実験的なプログラムをやっています。一連の実験データから高度を求める必要があります。時間と加速度のデータのみを記録し、二重積分を解くためにForward Euler Methodを使用しました。私が書いたコードはこれです。 public void forward_method(double[] acceleration, double[] time){ double deltaT =

0熱

1答えて

Python vegasエラー

私は最近、Python 2.7を使用するWindows 10マシンにvegas Monte Carlo統合パッケージ（バージョン3.0.1）をインストールしました。すべてがうまくいくように見え、ベガスを輸入することができます。私はラスベガスのチュートリアルからこのコードをコピー：Jupyterノートからこれを実行する import vegas import math def f(x):

0熱

1答えて

変位への加速

私は実験的なプログラムをやっています。私はスマートフォンをサポートして車のダッシュボードに固定しました。旅行中、私はアプリケーションから慣性センサの値を読みました。正確には、各登録の加速度計データ時間（固定間隔で）を秒単位で読み取ります。ここで、垂直加速度から垂直変位に渡したいと思います。これを行うには、二重積分を行う必要があります。オイラー法を試しました。 v0=0.v0=0などの初期条件

1熱

1答えて

Pythonを使用した三重積分の数値計算

私はthis積分を深く計算する必要があります。私は数ヶ月間、Pythonのnumpyパッケージ、特にintegrate.tplquad関数を使用しています。上記のコードで from __future__ import division from math import * import numpy as np import scipy.special as special import s

1熱

1答えて

Matlabで数値的に統合されています。なぜ部門が増えると正確な出力が低下するのですか？

私は現在、分散ロードビームを与えられたMatlabの問題に取り組んでおり、数値的に（複合台形ルールを使用して）積分してせん断力と曲げモーメントを求めています。そこから、私は最大の瞬時値とこれが起こる位置を見つける。それから私は再び2倍に分けてせん断力と分散荷重を求めます。私は、数値積分と微分が正しく計算されていることを確認し、確認するために分析的に統合された関数を使用しています。私の問題は、計算