«numerical-analysis» 태그된 질문

연속 문제의 대략적인 이산 솔루션을 계산하기위한 알고리즘의 구성 및 분석 표준적인 예는 차이 몫을 통한 미분의 근사입니다.

17
파이썬을위한 고품질 비선형 프로그래밍 솔버가 있습니까?
볼록하지 않은 전역 최적화 문제를 해결해야 할 몇 가지 문제가 있습니다. 현재 MATLAB의 Optimization Toolbox (특히 fmincon()algorithm = 사용 'sqp')를 사용하고 있습니다. 그러나 내 코드의 대부분은 Python이며 Python에서도 최적화를하고 싶습니다. 파이썬 바인딩과 경쟁 할 수있는 NLP 솔버가 fmincon()있습니까? 반드시 비선형 평등과 불평등 제약을 다룰 수있다 사용자가 Jacobian을 제공하도록 요구하지 않습니다. …

3
병렬 ODE 방법의 최신 상태는 무엇입니까?
현재 ODE 통합을위한 병렬 방법을 찾고 있습니다. 광범위한 접근 방식을 설명하는 새롭고 오래된 문헌이 많이 있지만, 일반적으로 주제를 설명하는 최근의 설문 조사 또는 개요 기사를 찾지 못했습니다. Burrage [1]의 저서가 있지만 거의 20 년이 되었기 때문에, 초현실적 알고리즘과 같은 더 현대적인 아이디어는 많이 다루지 않습니다. [1] K. Burrage, 정규 미분 …



3

1
숫자 코드를 어떻게 디버깅합니까?이 진동 오류의 원인은 무엇입니까?
조용히 많은 통찰력이 양식 경험을 얻을 수 있습니다. 누군가 전에 이와 비슷한 것을 보았는지 궁금했습니다. 플롯은 advection-diffusion 방정식의 초기 조건 (녹색), 반복 200에서의 용액 (파란색), 반복 400에서의 용액 (빨간색)을 보여줍니다. μ≈0.07μ≈0.07\mu\approx0.07C≈0.0015C≈0.0015C\approx 0.0015 배경. 이산화는 대류 및 확산 항의 중심 차이입니다. 나는 이것이 1 차 대류와 2 차 확산이라고 믿는다. 셀면의 …

1
유한 체적 방법을 사용할 때 경계 조건을 어떻게 적용해야합니까?
이전 질문에 이어이 비 균일 유한 체적 메쉬에 경계 조건을 적용하려고합니다. 도메인의 lhs ( Robin 유형 경계 조건을 적용하고 싶습니다 x=xL)x=xL)x=x_L). σL=(dux+au)∣∣∣x=xLσL=(dux+au)|x=xL \sigma_L = \left( d u_x + a u \right) \bigg|_{x=x_L} 여기서 σLσL\sigma_L 은 경계 값이며; a,da,da, d 는 경계, 이류 및 확산에 각각 정의 된 계수이고; 는 경계에서 …

3
과학 프로그래밍 콘테스트
저는 정기적으로 소위 "프로그래밍 콘테스트 (Programming Contests)"에서 경쟁합니다. 여기서 제한된 시간 동안 자신의 코드로 문제를 해결하고 문제 해결 기술을 사용하여 어려운 알고리즘 문제를 해결할 수 있습니다. 이러한 모양에 대한 참조 예를 보려면 Google Code Jam 또는 ACM-ICPC와 같은 경쟁을 검색하십시오. (프로그래밍 콘테스트가 무엇인지 아는 경우 다음 단락을 건너 뛸 수 …


3
행렬 제곱근 역의 효율적인 계산
통계에서 일반적인 문제는 대칭 양수 한정 행렬의 제곱근 역을 계산하는 것입니다. 이것을 계산하는 가장 효율적인 방법은 무엇입니까? 나는 우연히 어떤 문학 (나는 아직 읽지 않은), 일부 부수적 R 코드 여기에 내가 편의를 위해 여기에 재생됩니다, # function to compute the inverse square root of a matrix fnMatSqrtInverse = function(mA) { …

2
Newton-Raphson 반복을 사용하지 않고 비선형 PDE를 해결할 수 있습니까?
나는 일부 결과를 이해하려고 노력하고 있으며 비선형 문제를 해결하는 것에 대한 일반적인 의견에 감사드립니다. 피셔 방정식 (비선형 반응-확산 PDE), ut=duxx+βu(1−u)=F(u)ut=duxx+βu(1−u)=F(u) u_t = du_{xx} + \beta u (1 - u) = F(u) 신중한 형태로 u′j=Lu+βuj(1−uj)=F(u)uj′=Lu+βuj(1−uj)=F(u) u_j^{\prime} = \boldsymbol{L}\boldsymbol{u} + \beta u_j (1 - u_j) = F(\boldsymbol{u}) 여기서 은 미분 연산자이고 은 …

4
고정 된 수의 RHS 평가를위한 최적의 ODE 방법
실제로, 수치 해결 런타임 IVP X ( t 0 ) = X 0 종종 우측 (RHS)의 평가 기간을 지배한다 (F)를 . 그러므로 다른 모든 연산이 즉각적인 것으로 가정하자 (즉, 계산 비용이 없음). IVP를 해결하기위한 전반적인 실행이 제한되는 경우,이 평가 한 횟수 제한에 해당 F 일부 N ∈ N을 .x˙(t)=f(t,x(t)) for …

1
파이썬 함수에서 치명적인 취소를 피하는 방법은 무엇입니까?
함수를 숫자로 구현하는 데 문제가 있습니다. 큰 입력 값에서 결과는 매우 많은 수의 매우 작은 수라는 사실로 인해 어려움을 겪습니다. 치명적인 취소가 올바른 용어인지 확실하지 않으므로 수정하십시오. 무언가 잘못되었다는 증거 : 6의 더 큰 입력에 대해 진동과 0.0의 어설 션을 어떻게 피할 수 있습니까? 내 기능은 다음과 같습니다. import numpy …

2
유한 차분 법에 대한 폰 노이만 안정성 분석의 대안
나는 다음과 같이 주어진 결합 된 1 차원 다공성 탄성 방정식 (바이오 트 모델) 을 해결하기 위해 노력 하고 있습니다. ∂− ( λ + 2 μ ) ∂2유∂엑스2+ ∂피∂엑스= 0−(λ+2μ)∂2유∂엑스2+∂피∂엑스=0-(\lambda+ 2\mu) \frac{\partial^2 u}{\partial x^2} + \frac{\partial p}{\partial x} = 0 도메인Ω=(0,1)에서 경계 조건이있는 ∂ x 2 ]=q(x,t): ∂∂티[ γp + …

2
많은 변수의 수치 적분
하자 및 이 변수의 함수가 되십시오.x⃗ =(x1,x2,…,xn)∈[0,1]nx→=(x1,x2,…,xn)∈[0,1]n\vec{x} = (x_1, x_2, \dots, x_n) \in [0,1]^nf(x⃗ ):[0,1]n→Cf(x→):[0,1]n→Cf(\vec{x}): [0,1]^n \to \mathbb{C} 이 반복 적분에 대한 재귀 체계가 있습니까? ∫[0,1]n∏dxif(x⃗ )∫[0,1]n∏dxif(x→)\int_{[0,1]^n} \prod dx_i \;f(\vec{x}) 경우 및 I 브레이크 100 개의 세그먼트로, 우리가 까지 추가 포인트. 더 똑똑한 방법이 있어야합니다.n=10n=10n = 10[0,1][0,1][0,1]1020102010^{20} 실제로 통합하려는 기능 …

당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.