MM.1 · Coordinate operators, separated modes and stable updates

방사 함수의 라플라시안을 f″로 표기하면 항이 하나 빠집니다; 이 연산자는 기호에 속하는 것이 아니라 좌표에 속합니다.

선행 지식: 대수학, 미분 및 벡터 성분.

• 대칭성에 의해 선택된 좌표계에서 발산과 라플라시안을 사용하십시오.
• 단순한 분리 경계값 모드 문제를 풀이하십시오.
• 원본 수치 업데이트의 수렴, 단위 및 안정성을 확인하십시오.

연산자를 쓰기 전에 대칭성에 따라 좌표를 선택하십시오. 구면 좌표계의 방사 함수에 대해 ∇²f=f″+(2/r)f′이며, 방사장 A=A_r r̂의 발산은 (1/r²)d(r²A_r)/dr입니다. 확인: A=k r̂/r²의 발산은 (1/r²)dk/dr=0으로 원점에서 멀어질 때 일정하며, f=1/r은 ∇²f=f″+2f′/r=2/r³−2/r³=0로 원점에서 멀어질 때 일정함을 만족합니다. 카르테시안 형태는 카르테시안 종속성에만 적용됩니다.

변수 분리는 경계값 문제를 상미분 방정식 모드로 변환합니다. 양쪽 끝이 고정된 줄의 경우 u=X(x)T(t)로 씁니다. 파동 방정식은 X″/X=(1/c²)T″/T=−k²을 주며, 고정된 끝은 X=sin(nπx/L)을 강요하고, 허용되는 주파수는 ω_n=nπc/L, f_n=nc/(2L)입니다. 각 모드는 모든 모드 합에 앞서서 경계 조건을 만족해야 합니다.

An explicit Euler update for y′=−λy is y_{n+1}=(1−λh)y_n, bounded when |1−λh|≤1, i.e. 0<h≤2/λ for λ>0; asymptotic decay needs 0<h<2/λ. With λ=10 /s, h must not exceed 0.2 s; h=0.25 s multiplies by −1.5 each step and diverges with alternating sign. Stability is a property of the recurrence, not of the true solution.

수치 해를 신뢰하기 전에 세 가지 방법으로 확인하십시오: 모든 항의 차원(λh는 무차량이어야 함), 격자 세분화 시 수렴성(h를 반으로 줄이고 비교; 명시적 오일러의 오차는 1차임), 그리고 정확한 해가 알려진 한계 경우. y′=−4y, y0=8, h=0.2인 경우, 5단계로 계산하면 8×0.2⁵=0.00256이 되지만 정확한 값은 8e^(−4)≈0.1465입니다: 안정적이지만 부정확하므로 격자를 더 세밀하게 설정해야 합니다.

구면 방위 장은 발산을 가진다

$A_r=Cr$에 대해 이 식은 C가 아닌 3C를 산출한다. 고정된 끈의 끝에는 $X(0)=X(L)=0$이 필요하다, 따라서 양의 정수 n에 대한 $X_n=\sin(n\pi x/L)$. Euler for $\dot y=-\lambda y$ gives $y_{n+1}=(1-\lambda h)y_n$. 양수 h에 대해 유계 절대 안정성은 $h\leq2/\lambda$을 허용합니다; 영으로의 감쇠는 엄격한 부등식 $0<h<2/\lambda$을 요구합니다. 상한 경계에서 수치 해는 감쇠 없이 진동합니다. 안정성은 정확성을 보장하지 않습니다.

카르테시안 라플라시안을 방사 함수에 적용하거나, 경계 조건을 위반하는 모드를 합치거나, 단계의 세밀화 없이 겉보기로 유계적인 실행을 수렴한 것으로 간주하는 것은 잘못되었습니다. 실제 문제에 맞춰 연산자, 경계 조건 및 안정성 조건을 매칭하십시오.

변수 분리법: 미지 함수를 단일 변수 인자의 곱으로 표현하여 편미분방정식을 푸는 방법.

안정성 조건: 수치 반복 공식의 오차가 증폭되지 않도록 하는 단계 크기 제한.