Fourier 모드·격자 분산·연속 근사

1. 복잡한 모양을 단순한 운동으로 나누기

현의 움직임이 복잡해 보이는 이유는 서로 다른 공간 패턴이 겹쳐 있기 때문일 수 있습니다. 정상 모드 기저를 쓰면 “각 점이 어떻게 움직이는가?”를 “각 독립 패턴이 얼마나 들어 있는가?”로 바꿉니다. 평형 근처의 선형 계에서는 각 모드의 좌표가 서로 독립된 스칼라 방정식을 따릅니다. 과학적 머신러닝의 여러 스펙트럼 표현도 이 관점에서 출발합니다.

스칼라장의 작용과 경계 변분 및 국소 보존 흐름과 Noether 변환을 먼저 복습하세요. 여기서는 같은 공간 이산화를 사용하고 시간은 연속으로 둡니다. 평형과 순간 기준장, 물리적 파수와 격자 인덱스, 고정 모드 수렴과 전체 스펙트럼의 정확성을 각각 구별하는 것이 목표입니다.

2. 변환을 고르기 전에 격자를 정한다

구간 길이를 $L$, 움직이는 값의 수를 $N$이라 하겠습니다. 지정한 무차원 모형에서 질량 행렬은 $M=hI$이고, 활성 점 또는 고정 외부 값을 연결하는 변을 사용해 에너지를 씁니다.

$$U(q)=\frac{c^2}{2h}\sum_{(i,j)}(q_j-q_i)^2 +h\sum_i\left(\frac a2q_i^2+\frac b4q_i^4\right),\qquad H=\frac h2\sum_i v_i^2+U(q).$$

경계	활성 점의 위치	간격	구간 밖의 처리
주기	$x_j=jh$	$h=L/N$	마지막 점을 0번 점과 연결
고정값 0	$x_j=(j+1)h$	$h=L/(N+1)$	두 외부 0은 자유도가 아님
무유속	$x_j=(j+1/2)h$	$h=L/N$	셀 중심과 내부 $N-1$개 변

$L,N$이 같아도 경계를 바꾸면 $h$가 달라질 수 있습니다. 노드 중심에서 정의한 Neumann 조건은 여기의 셀 중심 모형과 다릅니다. 그 경우에는 끝점의 가중치와 기저를 다시 유도해야 합니다. “Neumann이므로 코사인 변환”이라는 말만으로는 이산화가 정해지지 않습니다.

3. 선언한 기준장 주변에서 선형화하기

활성 점에서 $q_i=\bar q+\eta_i$로 놓고 기준 속도는 0으로 둡니다. 국소 퍼텐셜의 곡률은 다음과 같습니다.

$$\mu^2=a+3b\bar q^2.$$

$\mu^2$라는 표기를 썼지만 양수라고 가정하지 않습니다. $g=\nabla U(\bar q)$, $K=\nabla^2U(\bar q)$로 쓰면

$$h\ddot\eta=-g-K\eta+O(\eta^2)$$

입니다. 기준장이 유한 문제의 정적 평형이려면 $g=0$이어야 합니다. 그렇지 않으면 일정한 힘 $-g$를 남겨 두어야 합니다. 헤시안을 대각화하는 계산 자체는 가능하지만, 그 고유값이 정지한 해 주변의 실제 진동을 뜻하지는 않습니다.

주기·셀 중심 무유속 경계에서는 $a\bar q+b\bar q^3=0$이면 균일 기준장이 평형입니다. 고정값 0 경계는 양 끝 내부 점에 추가 가속도 $-c^2\bar q/h^2$를 줍니다. 따라서 퍼텐셜의 비영 최솟값을 모든 점에 넣었다고 해서 경계까지 포함한 평형이 되는 것은 아닙니다.

4. 질량 내적에 맞게 정규화한다

변의 에너지에서 얻은 무차원 양의 이차 차분 행렬을 $D$라 하겠습니다. 선형 연산자는

$$A=M^{-1}K=\frac{c^2}{h^2}D+\mu^2I, \qquad Ku_\alpha=\lambda_\alpha Mu_\alpha$$

이고, 기저는 $u_\alpha^TMu_\beta=\delta_{\alpha\beta}$로 정규화합니다. 유클리드 길이가 1인 고유벡터라면 $\sqrt h$로 한 번 더 나눠야 합니다. 이를 빠뜨리면 고유값은 맞아 보여도 모드 좌표·노름·에너지가 틀립니다.

여기서는 모든 점의 질량이 같아 $A$도 대칭입니다. 비균일 질량 문제에서 일반화 고유값 문제를 무조건 보통의 대칭 문제로 바꿔도 된다는 뜻은 아닙니다. 질량 변환이나 일반화 해법을 사용해야 합니다. 중복 고유값의 고유공간에서는 여러 직교 기저를 고를 수 있으며, 개별 계수는 달라도 그 공간 전체의 노름은 같습니다.

5. 주기 격자의 평면파로 분산 관계 유도하기

주기 차분식에 $u_j=e^{ikjh}$를 대입하면

$$2u_j-u_{j-1}-u_{j+1}=4\sin^2(kh/2)u_j$$

입니다. 따라서 격자 연산자의 고유값과 허용 파수는

$$\lambda_h(k)=\mu^2+\frac{4c^2}{h^2}\sin^2(kh/2), \qquad k_m=\frac{2\pi m}{L}$$

가 됩니다. 평형에서 선형 모드 좌표는 $\ddot Q_m+\lambda_hQ_m=0$을 만족합니다. 양의 $\lambda_h$에서만 실수 진동수 $\omega_h=\sqrt{\lambda_h}$를 사용합니다.

격자 표본은 정수 $r$에 대해 $k$와 $k+2\pi r/h$를 구별하지 못합니다. 이는 부동소수점 오차가 아니라 표본값이 정확히 같은 현상입니다. $k$와 $-k$도 같은 고유값을 갖습니다. 실수 표현에서는 두 방향을 표현하되 자유도를 중복 계산하지 않아야 합니다.

6. 실수 주기 모드의 수를 정확히 센다

다음은 질량 내적으로 정규화한 한 가지 실수 기저입니다.

$$u_0(j)=\frac1{\sqrt{Nh}},\quad u_{m,c}(j)=\sqrt{\frac2{Nh}}\cos\frac{2\pi mj}{N},\quad u_{m,s}(j)=\sqrt{\frac2{Nh}}\sin\frac{2\pi mj}{N}.$$

사인·코사인 쌍은 $1\le m\le\lfloor(N-1)/2\rfloor$에서 사용합니다. $N$이 짝수이면 Nyquist 모드 하나를 더합니다.

$$u_{N/2}(j)=\frac{(-1)^j}{\sqrt{Nh}}.$$

Nyquist 사인 벡터의 표본값은 모두 0이므로 별도의 자유도가 아닙니다. 상수 모드에도 독립된 사인 짝이 없습니다. 총 차원은 홀수에서 $1+2(N-1)/2=N$, 짝수에서 $1+2(N/2-1)+1=N$입니다.

같은 파수의 사인·코사인 모드는 고유값이 같지만, 주어진 상태의 계수는 일반적으로 다릅니다. 수치 고유값 해법의 결과를 정렬하기만 하면 여기서 선택한 사인·코사인 방향까지 일치할 것이라고 기대하지 마세요. 중복 고유공간에서의 기저 선택을 먼저 맞춰야 합니다.

7. 고정 벽과 셀 중심 무유속의 기저는 다르다

$N$개 내부 점을 가진 고정값 0 문제에서는

$$u_m(j)=\sqrt{\frac2L}\sin\frac{\pi m(j+1)}{N+1},\quad k_m=\frac{\pi m}{L},\quad 1\le m\le N$$

를 사용합니다. 셀 중심 무유속 문제에서는 상수 벡터 $1/\sqrt L$와 다음 코사인 벡터를 사용합니다.

$$u_m(j)=\sqrt{\frac2L}\cos\frac{\pi m(j+1/2)}N, \quad 1\le m\le N-1.$$

두 경우 모두 허용된 $k_m$를 같은 격자 분산식에 대입할 수 있습니다. 다만 차분 행렬의 끝점 대각 원소가 다릅니다. Dirichlet은 2이고 셀 중심 Neumann은 1입니다. 그 차이가 기저와 영 공간 모드의 유무를 결정합니다.

Gilbert Strang의 이산 코사인 변환 논문은 경계의 위치와 대칭 이차 차분 행렬의 고유벡터가 연결됨을 설명합니다. 여기의 식은 앞서 정한 변 에너지·질량·셀 배치에 맞춰 유도한 것입니다. 다른 배치를 같은 변환 이름으로 뭉뚱그리지 마세요.

8. 투영·재구성·에너지를 함께 검사한다

선택한 실수 기저에서 모드 좌표와 속도를 정의합니다.

$$Q_\alpha=h\sum_j u_\alpha(j)\eta_j, \qquad V_\alpha=h\sum_j u_\alpha(j)v_j, \qquad \eta_j=\sum_\alpha u_\alpha(j)Q_\alpha.$$

기저의 완전성으로 질량 가중 Parseval 항등식이 나옵니다.

$$\sum_\alpha Q_\alpha^2=h\sum_j\eta_j^2, \qquad \sum_\alpha V_\alpha^2=h\sum_jv_j^2.$$

고정한 헤시안에 대한 이차 에너지는

$$H_2=\frac12\sum_\alpha(V_\alpha^2+\lambda_\alpha Q_\alpha^2)$$

입니다. 재구성 오차, 두 Parseval 식, 에너지 합을 함께 검사하세요. 에너지 하나만 맞으면 부호나 기저 순서의 오류를 놓칠 수 있습니다. $b\ne0$이면 $H_2$는 일반적으로 전체 비선형 에너지가 아닙니다. 기준장이 비평형이면 에너지 전개의 일차항도 따로 고려해야 합니다.

9. 네 점의 상태를 직접 분해해 본다

주기 $N=L=4$, $h=c=1$, $a=b=\bar q=0$에서

$$q=(1,0,-1,0),\qquad v=(1,0,0,0)$$

을 선택하겠습니다. 상수·코사인·사인·Nyquist 기저의 고유값은 $(0,2,2,4)$이고, 그 순서로 모드 좌표를 구하면

$$Q=(0,\sqrt2,0,0),\qquad V=(1/2,1/\sqrt2,0,1/2)$$

입니다.

성분	운동 에너지 기여	퍼텐셜 에너지 기여
상수	$1/8$	$0$
코사인, $m=1$	$1/4$	$2$
사인, $m=1$	$0$	$0$
Nyquist	$1/8$	$0$

따라서 $K=1/2$, $U=2$, $H=5/2$입니다. 변위는 하나의 코사인 모드만 쓰지만, 한 점에 준 초기 속도는 세 모드의 속도에 나뉩니다. 변위 그림만으로 전체 상태가 정해지지 않는 이유입니다. 상수 모드의 속도는 다른 모드가 진동하는 동안에도 전체 평균의 이동을 만듭니다.

10. 양·영·음의 곡률은 서로 다른 운동이다

정지한 기준장의 동차 선형 방정식에서 다음 세 경우를 구별합니다.

• $\lambda>0$: $Q=A\cos(\sqrt\lambda t)+B\sin(\sqrt\lambda t)$.
• $\lambda=0$: $Q=A+Bt$. 유한 주기를 가진 진동이 아닙니다.
• $\lambda<0$: $Q=Ae^{\sqrt{-\lambda}t}+Be^{-\sqrt{-\lambda}t}$.

곡률이 0인 방향은 같은 운동 에너지를 유지하며 계속 이동할 수 있습니다. 에너지가 양의 준정부호일 때에는 에너지 보존만으로 위치의 유계를 보장할 수 없습니다. 이차 에너지가 부정부호이면 불안정한 좌표가 커지는 동안에도 에너지가 일정할 수 있습니다.

이는 일정한 선형 연산자에 관한 결론입니다. 비선형 문제의 영 곡률 방향에는 고차항이 중요할 수 있습니다. 기준장이 비평형이면 표시한 진동수·성장률은 순간 고정한 동차 연산자의 성질이며, 기준 힘을 제거하거나 실제 궤적의 성장을 증명하지 않습니다.

11. 경계는 음의 국소 곡률을 안정화할 수 있다

고정값 0 내부 점 세 개, $c=2$, $a=-1$, $b=\bar q=0$을 생각하겠습니다. $L=4$이면 $h=1$이고

$$A=\begin{pmatrix}7&-4&0\\-4&7&-4\\0&-4&7\end{pmatrix}, \qquad \lambda_{\min}=7-4\sqrt2>0$$

입니다. 기울기와 벽의 복원 효과가 음의 국소 곡률보다 강합니다. 같은 계수와 점의 수를 유지하되 길이를 $L=8$로 늘리면 $h=2$이고

$$\lambda_{\min}=1-\sqrt2<0$$

가 되어 불안정 모드가 생깁니다. 연속 문제의 최저 모드 조건은 $\mu^2+c^2(\pi/L)^2>0$이지만, 유한 격자의 고유값도 따로 검사해야 합니다. 국소 퍼텐셜의 판정, 유한 구간의 안정성, 연속 문제의 임계 길이는 같은 주장이 아닙니다.

12. 위상 속도와 기호적 군속도를 구별한다

주기 분산식의 양의 진동수 가지에서 $k>0$일 때

$$v_{\mathrm{ph}}=\frac{\omega_h(k)}k,\qquad v_{\mathrm g}=\frac{d\omega_h}{dk} =\frac{c^2\sin(kh)}{h\omega_h(k)}$$

입니다. 질량항이 없는 양의 가지 $0<kh\le\pi$에서는

$$v_{\mathrm{ph}}=c\frac{\sin(kh/2)}{kh/2},\qquad v_{\mathrm g}=c\cos(kh/2)$$

로 단순해집니다. Nyquist에서 위상 속도는 $2c/\pi$이지만 분산식의 군속도는 0입니다. 격자 위에서 빠르게 번갈아 나타나는 무늬가 포락선을 속도 $c$로 전달한다는 뜻은 아닙니다.

유한 구간의 사인·코사인 모드는 서 있는 공간 패턴입니다. 분산식의 미분을 파동 묶음의 전달로 해석하려면 충분히 큰 구간, 좁은 파수 분포, 적절한 가지 등의 조건이 필요합니다. 이 계산은 파동 묶음을 시간 적분하지 않습니다. $k=0$ 또는 비진동 모드의 위상 속도는 0으로 나누어 만들지 않고 정의하지 않은 값으로 둡니다.

13. 물리적 문제를 고정한 채 격자를 세분한다

$L,c,a,b,\bar q$와 물리적 모드 인덱스 $m$을 고정하면 대응하는 $k$도 고정됩니다. 주기·셀 중심 격자는 $N\to2N$, 내부 Dirichlet 격자는 $N\to2(N+1)-1$로 바꿔야 세 경우 모두 $h$가 절반이 됩니다.

$x=kh/2\in[0,\pi/2]$에서 다음 경계를 얻습니다.

$$0\le\lambda_{\mathrm{cont}}-\lambda_h =c^2k^2\left[1-\left(\frac{\sin x}{x}\right)^2\right] \le\frac{c^2k^4h^2}{12}, \qquad \lambda_{\mathrm{cont}}=\mu^2+c^2k^2.$$

이 구간의 $\sin x\le x$와 $\sin x\ge x-x^3/6$을 결합하면 위 경계가 나옵니다. 전개의 선도 고유값 오차는 $c^2k^4h^2/12$이므로 연속 세분화의 오차 비율은 4로 접근합니다. 어떤 해상도에서나 정확히 4인 것은 아닙니다.

고유값이 0에 가까우면 제곱근의 조건이 나빠져 진동수 오차를 별도로 봐야 합니다. 연산자의 고유값이 수렴했다는 사실만으로 시간 적분기나 전체 비선형 해의 수렴을 결론내리지 않습니다.

14. 고정 모드 수렴과 고주파 전체의 정확성은 다르다

매 격자에서 새 Nyquist 모드를 고르면 $k=\pi/h$도 계속 바뀝니다. 이때 공간 기호의 상대 오차는

$$1-\frac4{\pi^2}$$

로 남으며 0으로 가지 않습니다. 이는 앞 절의 고정 $k$ 실험이 아닙니다. 세밀한 격자가 새로운 주파수를 표현할 수 있게 되는 것과, 매 격자의 차단 주파수 근처까지 연속 모형처럼 정확해지는 것은 다릅니다.

세분화 계산에서는 처음 고른 Nyquist의 인덱스를 그대로 유지합니다. 따라서 $N$을 두 배로 늘린 뒤 그 인덱스는 더 이상 Nyquist가 아닙니다. 기준장의 연속 경계 적합성도 확인하세요. 활성 점을 비영 상수로 두고 외부는 0으로 고정하면 세분화할수록 경계에 급격한 점프가 남습니다. 그 순간 헤시안의 스펙트럼을 비교할 수는 있지만, 이는 연속 고정 벽을 만족하는 매끄러운 상수 기준장의 수열이 아닙니다.

15. 비선형 나머지를 실제로 계산한다

사차 퍼텐셜의 유한 에너지 전개는 다음과 같이 정확합니다.

$$H(q,v)=U(\bar q)+g^T\eta+H_2(\eta,v) +hb\sum_j\left(\bar q\eta_j^3+\frac{\eta_j^4}{4}\right).$$

주기 $N=L=4$, $a=-1$, $b=\bar q=c=1$에서 $q=(2,1,0,1)$, $v=(1,0,0,0)$을 선택하겠습니다. 기준장은 $\mu^2=2$의 평형이지만 에너지는

$$U(\bar q)=-1,\quad g^T\eta=0,\quad H_2=9/2, \quad R=1/2,\quad H=4$$

입니다. 모드 기저를 만들 수 있다고 해서 유한 섭동이 자동으로 작은 것은 아닙니다. $a=0,b=\bar q=1$로 바꾸면 같은 균일 기준장의 가속도는 모든 점에서 $-1$이 되어 비평형입니다. 고정값 0 경계에서는 이중 우물의 국소 최솟값도 양 끝의 힘을 받습니다. 각 모드를 독립된 물리 진동자로 해석하기 전에 비선형 나머지와 기준 힘을 함께 검사하세요.

16. 수치 검증이 증명한 범위를 읽는다

간격, 에너지 성분, 에너지 전개 잔차 등은 정확한 분수로 계산할 수 있습니다. 반면 대부분의 사인·코사인 계수, 고유값, 진동수는 무리수이므로 바깥쪽으로 반올림한 포함 구간으로 표시합니다. 0을 포함하는 좁은 구간과 정확히 0이라는 대수적 증거는 다릅니다. 고유값의 부호가 현재 정밀도로 결정되지 않으면 미결정으로 남깁니다.

계산에서는 질량 직교성, 모드 에너지, 두 Parseval 식, 재구성과 선택한 고유방정식을 대조합니다. 독립 참조는 변의 결합 행렬에서 에너지 행렬을 구성한 뒤 고정밀로 대각화하며 제품의 분산 공식을 호출하지 않습니다. 고정밀 수치 대조가 유향 산술의 포함 증명을 대신하는 것은 아닙니다.

이 검사들은 시간을 전진시키지 않습니다. 적분기를 사용하는 연구에는 시간 간격 의존성, 위상 오차, 안정성, 에너지 드리프트 검사가 추가로 필요합니다. 중복 고유값에서는 임의 해법의 고유벡터를 그대로 비교하지 말고 고유공간 또는 선언한 기저 규약을 비교하세요.

17. AI 연구에 연결할 때의 실험 설계

스펙트럼 신경 연산자는 계수 사이의 관계를 학습할 수 있지만, 그 계수는 구간·격자·경계·정규화에 의존합니다. 고정 벽 데이터에 주기 FFT 표현을 적용하면 경계에서 불연속을 만들 수 있습니다. 해상도가 바뀔 때 물리적 파수를 추적하지 않으면 주파수 배열의 같은 위치가 같은 대상을 뜻하지 않을 수도 있습니다.

$L$, 물리적 대역, 기준 조건을 고정한 실험에서 낮은 모드의 예측 오차와 차단 주파수 근처의 오차를 나눠 보세요. 재구성, 질량 노름, 에너지와 비선형 나머지도 기록합니다. 비선형 곱셈은 해상도 밖의 고조파를 만들고 그것이 남은 대역으로 alias될 수 있습니다. 차수와 이산화에 맞는 처리가 필요하며 하나의 de-aliasing 규칙을 설명 없이 모든 비선형식에 적용하지 않습니다.

Stamatopoulos의 격자 분산 연구는 공간 차분 방식에 따라 분산 관계가 달라짐을 보여 주는 참고 문헌입니다. 그 논문은 3차원 스칼라장 시뮬레이션을 다루며, 이 단원의 1차원 경계와 손계산은 별도로 유도했습니다. 한 격자의 수치 일치만으로 해상도 간 일반화나 새로운 물리 법칙을 주장하지 마세요.

18. 먼저 예측하고 계산으로 설명하기

1. 네 점의 변위가 하나의 모드이면 모드 속도도 하나일까요? 아닙니다. $Q=(0,\sqrt2,0,0)$이지만 $V=(1/2,1/\sqrt2,0,1/2)$입니다. 변위와 속도는 별도 데이터입니다.
2. 곡률이 0이면 좌표도 정지할까요? 동차 선형식은 $\ddot Q=0$이므로 $Q=A+Bt$입니다. $B\ne0$이면 이동합니다.
3. 국소 곡률이 음수이면 유한 문제도 불안정할까요? 세 점 고정 벽 예제는 $L=4$에서 $7-4\sqrt2>0$, $L=8$에서 $1-\sqrt2<0$입니다. 경계와 길이를 포함해 판정해야 합니다.
4. Nyquist의 위상 속도는 0일까요? 질량항 없는 양의 가지에서 위상 속도는 $2c/\pi$, 기호적 군속도는 0입니다.
5. 오차 비율이 4 근처이면 모든 모드의 정확성이 검증됐을까요? 아닙니다. 이차 오차 전개는 고정 $k$ 실험입니다. 매 격자에서 Nyquist를 새로 고르면 조건이 달라집니다.
6. 순간 헤시안이 양수이면 독립적인 비선형 진동자로 보아도 될까요? 기준장의 평형과 나머지의 충분한 작음을 먼저 확인해야 합니다. 비평형 기준장과 유한 진폭 이중 우물 예제는 서로 다른 조건에서 이 추론을 깨뜨립니다.

재현할 때는 경계의 배치, $L,N,h$, 기준 힘, 기저 정규화, 물리적 모드 인덱스, 고유값 부호, 정확한 에너지 전개를 함께 기록하세요. 유한 대수, 공간 극한, 시간 진화, 비선형 근사 중 실제로 무엇을 검사했는지 나누어 기술해야 합니다.