numerical-integration

-1열

1답변

, 나는 확실 해요 하나의 변수 t를 처리해야했기 때문에 두 개의 변수 함수에 적응하는 데 잘못된 것이 있다고 가정합니다. 감사합니다! %%Calculation of dH/dt for mode m=1 for the entire sphere, NH and SH clear all %%Radius of photosphere r = 6.957*(10^5);

1열

1답변

MATLAB - sinusoidal forcing 함수를 ode45에 전달

저는 Matlab을 처음 사용하고 있으며 기본 사항을 파악하는 데 정말로 어려움을 겪고 있습니다. 저는 댐핑과 구동력을 가진 질량/스프링 시스템의 위치와 속도를 해결하는 함수 myspring을 가지고 있습니다. ode45을 실행하기 전에 명령 창에서 스프링 강성 (k), 감쇠 계수 (c) 및 질량 (m)에 대한 값을 지정할 수 있습니다. 내가 할 수없는 것

-3열

1답변

특정 기능의 0을 찾을 수있는 프로그램을 작성할 수 있습니까?

Gauss-legendre 통합에서 우리는 legendre 함수의 제로를 찾아야하지만 그렇게 할 수있는 코드를 작성하는 방법을 찾을 수 없습니까? 이 함수가 0 인 "xi"목록이 있다는 것을 알고 있지만이 "xi"를 소유하고있는 프로그램을 작성할 수 있습니까?

1열

1답변

더 많은 인수와 통합 기능이 있습니다.

간단한 문제가 있지만 어디서나 솔루션을 찾을 수 없습니다. 함수를 통합해야합니다 (예 : Simpson의 규칙 서브 루틴 사용).하지만 하나 이상의 인수에 함수를 전달해야합니다. 하나는 나중에 통합하려는 변수이고 다른 하나는 함수에서 오는 값입니다. 내가 함수 안에서 수행 할 수없는 다른 계산. Simpson 서브 루틴은 f (x) 만 받아 들여 f (x,

2열

1답변

k3_1의 출력은 -3.1445e + 24로 제한됩니다.

저는 RK4를 사용하여 ODE 시스템을 해결하고 있습니다. 저는 k3_1이 -3.1445e + 24로 제한되어 있기 때문에 직선 플롯을 생성합니다. 나는 그것이 왜 모자란 지 이해하지 못한다. function RK4system_MNModel() parsec = 3.08*10^18; r_1 = 8.5*1000.0*parsec; % in cm z_1

1열

2답변

커브 아래의 영역 - 1D 배열 (Java)

저는 파이썬과 같은 대부분의 언어에서 간단한 질문을합니다. 고정 소수점의 1D 배열에서 정수 (곡선 영역)를 얻으려고합니다. 자바는 분명히 많은 수치 적 통합 라이브러리를 가지고 있는데, 이들 모두는 함수 (f {double (x)})를 입력으로 요구하는 것 같습니다. 그러나 [1,4,10,11]와 같은 배열 (double [])을 수용 할 수있는 것을 찾

2열

1답변

Rcpp를 사용하여 C++에서 수치 적분을 수행하는 방법 R

Rcpp은 R에서 C++을 호출하여 수치 적분을 수행하는 방법에 대해 궁금해했습니다. 현재 설치에는 시간이 오래 걸리고 오류가 발생하기 쉽습니다. 기본 R 수치 통합 패키지보다 더 나은 것이 필요하다고 생각합니다. R 내의 수치 적 통합을 통해 C++에서 이러한 문제를 해결할 수 있을까요? funk <- function(x,b) { 10^b * exp(-x

2열

1답변

R의 적분

R의 적분 기능에 문제가 있습니다. 보봇 (Boe-Bot)을 플로팅하려하고 있지만 제대로 작동하지 않는 것 같습니다. t <- seq(0, 0.04, 0.0001) vi <- function(x) {5 * sin(2 * pi * 50 * x)} vo <- function(x) {integrate(vi, lower=0, upper=x)$value} t

-2열

1답변

원주의 Matlab 통합

간단한 질문 일 수 있습니다. Matlab에서 수치 적 통합을 이해하려고합니다. 나는 그것이 2*pi 것으로 예상하는 반면, 단위 원의 경우, fun = @(x) sqrt(cosd(x).^2+sind(x).^2); q = integral(fun,0,360); 대답은 360.0000입니다. 각도는 각도로 유지됩니다. 나는 혼란 스럽다. 내가 도대체

0열

1답변

Scipy를 사용하여 구적 구에 대한 다 변수 가우스 함수 정의

scipy를 사용하여 구적법에 대한 다 변수 가우스 pdf를 정의하는 데 문제가 있습니다. 평균 벡터와 공분산 행렬을 입력으로 받아 가우스 함수를 반환하는 함수를 작성합니다. 내가 가우스를 정의하고 가우스에 색인을 시도 make_mvn_pdf를 사용하는 경우 def make_mvn_pdf(mu, sigma): def f(x): retur