Duffing 진동의 세속항과 장시간 위상

1. 진폭이 맞아도 시계가 틀릴 수 있다

두 진동의 진폭과 에너지가 비슷해도 봉우리에 도달하는 시각은 다를 수 있습니다. 한 주기 동안의 아주 작은 주파수 차이가 여러 주기 쌓이면 큰 위상 오차가 됩니다. 고정한 시간에서 정확한 섭동 전개가 더 긴 관측 구간에서도 정확하리라는 보장은 없습니다.

매일 조금씩 어긋나는 두 시계를 떠올려 봅시다. 시곗바늘의 크기가 계속 커지지 않아도 시각은 크게 달라집니다. 근사식에서 시간에 따라 커지는 항도 실제 진폭 증가가 아니라 잘못 펼친 시계의 흔적일 수 있습니다.

이번 장에서는 단순 시간 전개, 두 단계 주파수 보정과 독립적으로 포함 구간을 계산한 비선형 주기해를 비교합니다. 작은 시간 간격으로 적분한 수치 궤적을 정확해라고 부르지 않습니다.

2. 방정식과 유효 결합을 먼저 정한다

무차원 초기값 문제는

$$x''+x+\lambda x^3=0,\qquad x(0)=A\ge0,\quad x'(0)=0,\quad\lambda\ge0$$

입니다. 질량과 선형 강성을 조정해 선형 각주파수가1이 되도록 했으며 시간도 그 단위로 잽니다. A>0이면 x=Aq로 두어

$$g=\lambda A^2,\qquad q''+q+gq^3=0,\quad q(0)=1,\quad q'(0)=0$$

를 얻습니다. 작은 λ만으로 약한 비선형성이라고 할 수 없습니다. λ=0.01일 때 A=1이면 g=0.01이지만 A=4이면 g=0.16입니다. 같은 결합계수에서도 진폭이 비선형성의 크기를 바꿉니다.

A=0이면 물리적 해는 x=0인 평형입니다. 선형화한 방정식의 주파수1과 실제 진동 궤적의 주기는 다른 개념입니다. 움직이지 않는 평형에 물리적 진동 주기를 부여하지 않습니다.

3. 정확한 운동이 할 수 있는 일을 확인한다

방정식에 x′를 곱하면 에너지가 보존됩니다.

$$E=\frac{x'^2}{2}+\frac{x^2}{2}+\frac{\lambda x^4}{4} =\frac{A^2}{2}+\frac{\lambda A^4}{4}.$$

λ≥0에서는 |x|가 커질수록 퍼텐셜이 증가하므로 |x|≤A, |x′|≤A√(1+g/2)입니다. 비영 전환점 진폭은 −A와 A 사이의 주기 궤적을 만듭니다. 근사를 보기 전부터 사용할 수 있는 범위 검사입니다.

따라서 근사식의 t에 비례하는 항을 실제 진폭이 무한히 증가하는 운동으로 해석할 수 없습니다. 이 보존계에서는 근사식의 적용 범위가 끝났다는 신호입니다. 외력·감쇠·다른 퍼텐셜을 추가하면 이 결론의 가정부터 다시 점검해야 합니다.

4. 단순한 1차 전개를 직접 구한다

원래 시간 좌표를 유지한 채 q=q₀+gq₁+O(g²)로 둡니다. 선도해는 q₀=cos t이고 다음 방정식은

$$q_1''+q_1=-\cos^3t=-\frac34\cos t-\frac14\cos3t, \qquad q_1(0)=q_1'(0)=0$$

입니다. 강제항에 같은 미분 연산자의 동차해인 cos t가 들어 있습니다. 이 공진 성분이 t sin t 형태의 항을 만듭니다. 초기조건까지 적용하면

$$q_{\mathrm N}(t)=\cos t+ g\left[\frac{\cos3t-\cos t}{32}-\frac38t\sin t\right]$$

입니다. 세 번째 조화파는 유계인 모양 보정이고, 시간에 비례하는 항이 세속항입니다. 세속항 자체와 그 포락선을 구분하세요. sin t가0이 되는 순간에는 보정값이 작아 보여도 긴 시간 구간 전체의 정확도가 복구된 것은 아닙니다.

5. 비균일 극한을 읽는 법

고정한 t에서 보정은 g와 함께0으로 갑니다. 하지만 관측 길이가1/g처럼 커지면 gt가 작지 않을 수 있습니다. 다음 상계는 즉시 얻을 수 있습니다.

$$|q_{\mathrm N}(t)-\cos t|\le g\left(\frac1{16}+\frac38|t|\right).$$

이는 표시한 보정과 선형 코사인의 차이에 대한 상계입니다. 정확한 Duffing 해와의 오차 상계로 이름을 바꾸면 안 됩니다. 물리적 변위 보정으로 읽으려면 A를 곱합니다.

더 간단한 식 cos((1+α)t)=cos t−αt sin t+…에서도 같은 일이 일어납니다. 이동한 위상을 펼치면 원래 함수가 유계여도 세속항이 생깁니다. 위상 이동을 코사인 안에 남기는 것은 근사의 구조를 재조직하는 방법입니다.

6. 시계를 펼치지 말고 주파수를 바꾼다

τ=ωt를 도입하고 τ에 대해 주기적인 해를 찾습니다. 물리 시간의 이차 미분은 ω²q_ττ입니다. ω=1+ω₁g+…와 q=cos τ+gY₁(τ)+…를 넣으면

$$Y_1''+Y_1= \left(2\omega_1-\frac34\right)\cos\tau-\frac14\cos3\tau$$

를 얻습니다. 주기적인 보정에는 주파수1의 공진 강제항이 없어야 하므로 ω₁=3/8로 정합니다. 남은 세 번째 조화파는 유계 보정으로 풀립니다.

이것이 Poincaré–Lindstedt 구성입니다. 그림을 보고 주파수를 맞추는 것이 아니라, 주기적인 보정 방정식이 풀릴 수 있도록 시계를 정합니다. 원래 시간에 대한 전개에서 피할 수 없었던 세속 성분이 주파수 이동으로 들어갑니다.

7. 1차 주파수 보정 파형

그 결과는

$$\omega_1^{\mathrm{app}}=1+\frac38g,\qquad q_{\mathrm{PL1}}(t)=\cos\tau+ \frac{g}{32}(\cos3\tau-\cos\tau),\quad \tau=\omega_1^{\mathrm{app}}t$$

입니다. 모든 τ에 같은 보정 시계를 적용합니다. 속도를 구할 때는 미분의 연쇄법칙으로 ω와3ω가 붙습니다. 주된 코사인만 보정하고 세 번째 조화파는 옛 시계로 두는 것은 다른 근사입니다.

파형이 유계이고 전환점 초기조건을 만족하는 것은 장점입니다. 그렇다고 주파수가 무한히 긴 시간까지 정확하다는 뜻은 아닙니다. 다음 단계에서는 남아 있는 작은 주파수 오차가 언제 드러나는지 봅니다.

8. 2차까지 계산을 이어 간다

다음 차수에서도 공진 코사인을 없애면

$$\omega_2^{\mathrm{app}}=1+\frac38g-\frac{21}{256}g^2, \qquad Y_2(\tau)=\frac{23\cos\tau-24\cos3\tau+\cos5\tau}{1024}$$

입니다. 선도항과1차항의 τ도 ω₂_app t로 바꾸고 g²Y₂를 더합니다. 23−24+1=0이므로 q(0)=1을 유지하며 모든 사인 미분이 처음에0이어서 q′(0)=0도 유지합니다.

2차 주파수는 약한 결합에서의 잘린 전개입니다. g=16이면 −14가 되지만 정확한 단일 우물 진동자의 주파수는 여전히 양수입니다. 계산기는 형식적인 식을 보존하면서 물리적 주파수 근사로는 부적합하다고 표시합니다. 실제 궤적이 없다는 결론으로 바꾸지 않습니다.

9. 복수 척도로 같은 현상을 설명한다

느린 시간 T=gt를 도입하고 전개 중에는 t와 T를 독립 변수로 취급합니다. q₀=a(T)cos(t+φ(T))로 두고 d/dt=D₀+gD₁을 적용하면 (D₀²+1)q₁=−2D₀D₁q₀−q₀³입니다.

공진 사인과 코사인 계수가 사라지는 조건은

$$\frac{da}{dT}=0,\qquad \frac{d\phi}{dT}=\frac38a^2$$

입니다. 정규화한 초기 진폭 a(0)=1에서는 선도 위상이 t+3gt/8이 되어1차 Lindstedt 시계와 일치합니다. 진폭은 일정하지만 위상이 천천히 이동합니다. 이 보존계에서는 에너지가 진폭을 고정하고 비선형성이 시계를 바꾸는 셈입니다.

이 계산은 가해 조건과 느린 변조를 보여 주는 예제입니다. 모든 동역학계와 모든 장시간에 적용되는 평균화 정리의 증명은 아닙니다. 감쇠·외력·추가 공진이 생기면 변조 방정식도 달라집니다.

10. 작은 잔차는 별도의 보증이다

유한 파형 Y(τ)에서 x=A Y(ωt)를 원래 방정식에 대입합니다. 정규화한 잔차는

$$R(\tau)=\omega^2Y''+Y+gY^3$$

입니다. 홀수 배각 코사인은 z=cos τ의 Chebyshev 다항식이므로,2차 파형의 잔차는 최대15차인 유한 다항식입니다. R(z)=∑rⱼzʲ이면 |z|≤1을 이용해 |R|≤∑|rⱼ|를 얻습니다. 물리적 방정식 잔차의 상계는 여기에 A를 곱한 값입니다.

1·2차 구성은 각각 g²와 g³ 미만의 항을 소거합니다. 이 상계는 방정식을 얼마나 잘 만족하는지 보증하지만, 위상 누적이나 전 시간 구간 궤적 오차를 단독으로 보증하지 않습니다. 작은 잔차에서 작은 장시간 오차로 가려면 안정성과 누적 효과에 관한 추가 논증이 필요합니다.

11. 위상 진단량을 먼저 정의한다

정확한 기본 각주파수를 ν라고 하면 관측 끝 H에서 비교하는 양은

$$\Delta\Phi_j(H)=(\omega_j^{\mathrm{app}}-\nu)H$$

입니다. 단순 전개의 시계는 ω_app=1로 둡니다. 이 값은 기본 주파수에서 얻은 펼친 위상 차이입니다. atan2(x,x′)로 읽은 각도도, 변위 오차도 아니며2π로 나눈 나머지를 취하지 않습니다. 비양수인 잘린 주파수는 형식적인 위상식일 뿐 적절한 물리 주파수 해석이 아닙니다.

g=0.01에서 기준 주파수는 약1.00374183617774233입니다. 1차식은1.00375,2차식은1.003741796875입니다. 다음은 반올림한 위상 차이 예시입니다.

H	단순 시계	1차 보정 시계	2차 보정 시계
100	−0.374184	0.000816382	−0.00000393027
10,000	−37.4184	0.0816382	−0.000393027
1,000,000	−3741.84	8.16382	−0.0393027

정확한 포함 구간은 계산 표를 읽으세요. 유계인 파형이라도 긴 시간에는 맞는 시계에서 멀어질 수 있음을 보여 줍니다. 같은 위상 차이에서도 관측 순간이 봉우리인지 영점을 지나는지에 따라 위치 오차는 다르게 나타날 수 있습니다.

12. 차수를 올리면 관련 시간 척도도 달라진다

정확한 약결합 주파수는1+3g/8−21g²/256+…로 시작합니다. 단순 시계의 주파수 오차는 g,1차 보정은 g²,2차 보정은 g³에서 시작합니다. 시간을 곱하면 각각1/g,1/g²,1/g³ 척도에서 사라지지 않는 위상 누적을 예상할 수 있습니다.

이는 다른 자료를 고정하고 g를 줄이는 점근적 척도에 관한 설명입니다. 계수는 실제 판단에서 중요합니다. 0.04라디안과4라디안은 둘 다 점근적 의미의 크기1일 수 있어도 실용적인 오차는 크게 다릅니다. “장시간 정확하다”는 문구보다 요청한 H에서 실제 포함 구간으로 구한 위상 차이가 더 많은 정보를 줍니다.

시간 척도를 말할 때 어떤 매개변수를 고정했는지도 써야 합니다. A를 바꾸면 g 자체가 달라지므로 λ에 대한 설명을 진폭 변화 실험에 그대로 옮기면 안 됩니다.

13. 정확한 주기 기준해를 만든다

이 초기값 계열은 다음 해석적 해를 갖습니다.

$$x(t)=A\operatorname{cn}\!\left(\sqrt{1+g}\,t\mid m\right), \qquad m=\frac{g}{2(1+g)}.$$

m은 타원함수의 parameter입니다. modulus k를 쓰는 자료에서는 m=k²입니다. 두 표기를 혼동하면 다른 값을 계산하게 됩니다. cn의 미분방정식에 위 m을 대입하면 Duffing 식이 되고, 초기값과 초기 미분이 A와0을 줍니다.

완전 타원적분 K(m)에서 주기 T=4K(m)/√(1+g)를 얻습니다. 기준해는 비선형 주기 파형입니다. 더 좋은 주파수를 넣은 코사인 하나로 바꿔 놓지 않습니다. 파형 보정과 시계 보정은 별개의 효과이기 때문입니다.

14. 산술·기하 평균으로 주파수를 포함한다

M=AGM(1,√(1−m))이면 K=π/(2M)이므로

$$\nu=\sqrt{1+g}\,M,\qquad T=\frac{2\pi}{\nu}$$

입니다. AGM은 양수 두 개를 산술평균과 기하평균으로 반복해서 교체합니다. 정확한 기하평균 수열은 극한 아래에서, 산술평균 수열은 위에서 접근합니다. 방향을 보존한 구간 연산으로 반올림 후에도 이 하한·상한을 유지합니다.

구현은18회 반복의 포함 구간을 사용합니다. 제곱근과 nome 계산에는512비트 구성 요소, 삼각함수에는768비트 작업 격자를 쓰며 최종 값은192비트 격자로 바깥쪽 반올림합니다. 이는 연산 설정이지 모든 결과에 정확한 상대192비트가 있다는 약속이 아닙니다.

15. 시간 간격을 누적하지 않고 파형을 포함한다

Nome Q=exp(−πK(1−m)/K(m))는 이 범위0<m<1/2에서1/16보다 작습니다. C=4M√(Q/m)로 두면 cn 파형을 수렴하는 급수로 쓸 수 있습니다.

$$\frac{x(t)}A=C\sum_{n=0}^{\infty} \frac{Q^n}{1+Q^{2n+1}}\cos((2n+1)\nu t).$$

N=64항만 남기면 위치 꼬리는 A C Qᴺ/(1−Q)로 제한됩니다. 속도를 구하면 (2n+1)ν가 추가되므로 그에 맞는 기하급수 꼬리를 따로 합합니다. 이 꼬리와 구간 산술 반올림을 함께 전파합니다.

따라서 H에 도달하려고 작은 시간 간격을 수없이 누적하는 기준해가 아닙니다. 그래도 주파수 불확실성에 따른 위상 구간은 시간과 함께 넓어질 수 있으므로 실제 반환 구간을 읽어야 합니다. 독립 검산은 같은 Fourier 계산을 되풀이하지 않고 고정밀 타원함수 cn/sn/dn을 직접 사용합니다.

16. 에너지가 맞아도 위상이 맞는 것은 아니다

같은 에너지 곡선 위에 있어도 그 곡선의 다른 위치에 있을 수 있습니다. 따라서 근사 에너지가 거의 보존되어도 시각이 크게 어긋날 수 있습니다. 반대로 큰 에너지 결함은 근사의 다른 성질도 깨졌음을 알려 주는 유용한 증거입니다.

관측점에서는 위치·속도, 주기 기준해에 대한 각각의 오차, 방정식 잔차, 에너지와 에너지 결함을 따로 표시합니다. 정확한 기준해는 초기 에너지를 보존합니다. 한 진단량을 다른 진단량으로 대체하거나 위상 차이에 곧바로 위치 오차라는 이름을 붙이지 않습니다.

학습된 적분기나 신경 ODE의 결과를 평가할 때도 같은 구분이 필요합니다. 에너지 그래프가 평평한지, 궤적 모양이 닫혀 있는지, 봉우리가 맞는 시각에 나타나는지는 서로 다른 검사입니다.

17. 긴 기록에는 국소적으로 해상한 그림이 필요하다

전체 기록은 [0,H]의 시간65점입니다. H가 매우 크면 모든 진동을 보여 줄 수 없습니다. 이 점들을 매끄러운 선으로 잇기만 하면 실제로 없는 느린 주기가 보일 수 있습니다.

두 번째65점 창은 관측점 근처에 두고 구간 밖으로 나가지 않게 자릅니다. 폭은 최대32/ceil(ν_upper)로 정해 몇 주기를 유용한 표본 수로 보게 합니다. 전체 격자의 표시는 주기당8점 기준을 보수적으로 검사합니다. 이는 표본화 규약이며 ODE 이산화 오차 정리가 아닙니다.

전체 시간의 위상 진단은 누적을, 국소 파형은 모양과 시각을 보는 데 사용하세요. 희소한 그림이나 닫혀 보이는 위상 초상만으로 연속 시간 전체의 궤적을 보증하지 않습니다. 표본 간 변화와 관측 창을 바꿨을 때의 변화도 구분해야 합니다.

18. 특수한 경우와 입력 범위를 존중한다

입력은 비음 유리수 λ,A≤16, 관측 길이 H≤999999999999와 관측 비율[0,1]입니다. 분자·분모는 각12자리 이내입니다. 유효 결합은4096까지 갈 수 있으며 강한 예제는 일부러 근사 실패를 보여 줍니다.

진폭0에는 주기가 없습니다. 진폭이 비영이고 λ=0이면 조화 진동이며 세 근사가 모두 정확해와 일치합니다. H=0이면 초기 상태와 근사 오차0을 대수적으로 처리합니다. 아주 작은 양의 결합은 긴 관측 구간에서도 정확한0으로 반올림하지 않습니다.

이 모형은 감쇠·주기 외력·음의 결합에서의 탈출·다중 우물·혼돈을 다루지 않습니다. 그런 현상은 다른 가정과 기준 문제가 필요합니다. 단순한 보존계 예제를 복잡한 강제 진동자의 완성된 설명으로 취급하지 않습니다.

19. 손계산 문제와 해설

1. 진폭의 역할. λ=0.01에서 A=1과4를 비교하세요. g는0.01과0.16입니다. λ가 같아도 뒤의 실험이 덜 약한 비선형성을 갖습니다.
2. 세속 강제항 찾기. cos³t=(3cos t+cos3t)/4로 전개하세요. cos t가 q₁″+q₁과 공진해 t sin t를 만듭니다. 세 번째 조화파만 있었다면 보정은 유계입니다.
3. 첫 시계 유도. 공진 계수2ω₁−3/4를0으로 놓으면 ω₁=3/8입니다. 그림에 주파수를 맞추는 대신 주기적 보정이 존재할 조건으로 결정합니다.
4. 초기조건 확인. 2차 파형의23,−24,1은 합이0이고 사인 미분은 처음에0입니다. 전환점 초기조건을 확인한 것이지 전 방정식을 정확히 푼 것은 아닙니다.
5. 위상 누적 계산. g=0.01, H=10,000에서1차 시계의 차이는 약0.08164라디안이고 H=1,000,000에서는 약8.164라디안입니다. 유계 진폭이 누적을 막지는 못합니다.
6. 음수 주파수 해석. g=16의2차 잘린 식은−14입니다. 정확한 궤적은 여전히 주기적입니다. 근사의 실패를 물리 진동자의 음수 각주파수나 궤적 부재로 해석하지 않습니다.
7. 잔차와 오차 구분. 작은 다항식 잔차만으로 전체 시간의 위치 오차를 제한할 수 있을까요? 안정성·누적 논증이 추가로 필요합니다. 잔차·위상 차이·실제 관측점 오차를 따로 비교하세요.
8. 그림 선택. H=1,000,000의65점은 모든 주기를 해상하지 못합니다. 국소 창으로 진동을 보고 펼친 위상으로 누적을 봅니다. aliasing된 점에서 느린 새 진동을 추론하지 마세요.

20. 연구 연결과 다음 범위

도구를 고르기 전에 네 섭동 문제를 비교하기

단원에서 같은 문자를 쓴다고 같은 물리적 매개변수인 것은 아닙니다. 평형 문제의 g=λh²/k³은 보정 전 변위에서, Duffing 문제의 g=λA²은 지정한 전환점 진폭에서 정의합니다. 사차 적분의 g는 gX⁴/4라는 자체 규약의 무차원 결합이고, 경계층의 δ=ε/|a|는 단위 구간에 대한 공간 폭입니다. 논증을 다른 문제로 옮길 때는 방정식·자료·극한 과정을 다시 적으세요.

문제와 복습 단원	작은 양과 극한	이 모형에서 얻는 보증	항을 하나 더 계산하는 것만으로 해결되지 않는 실패
볼록 평형	g=λh²/k³; k>0와 h 고정	근 포함 구간, 변위 오차 ≤ 잔차 크기/k	수렴 반경 밖의 Taylor 급수, k=0에서 기울기 하한 상실
사차 적분	남긴 차수 고정 후 g→0⁺	부호 있는 유한 나머지와 별도 구적 포함 구간	고정 g>0에서 항이 결국 성장함; 무한 합과 적분의 무근거 교환
경계층	비영 흐름과 자료를 고정하고 δ→0	합성해의 균일 값 오차 = 불일치 크기·exp(−1/δ)	축약식이 경계조건을 잃음; 늘린 좌표와 정합이 필요함
Duffing 운동	관측 시간 H를 명시하고 g→0	기준해 포함 구간, 잔차 상한과 기본 위상 누적	고정 시간 전개가 늘어나는 시간 구간에서 실패함; 시계를 재구성해야 함

네 계산기가 다음 결과를 제출했다고 생각해 보세요. 각각 무엇을 더 확인해야 할까요? 서로 다른 물리량과 단위의 숫자이므로 행 사이의 크기를 비교해 방법의 순위를 매기지 않습니다.

평형의 출력. k=h=1, λ=0.1에서 이차 근사 p=0.93의 잔차는0.0104357입니다. 기울기 하한 k=1이 있으므로 변위 오차 상한0.0104357로 바꿀 수 있습니다. 이 볼록 방정식의 가정이 필요하며 임의의 미분방정식으로 그대로 옮기지 않습니다.

적분의 출력. g=0.1에서 S₂=0.9578125의 첫 누락항 상한은0.0270703125입니다. 무한 급수는 발산해도 유한 나머지 보증은 유효합니다. 증명은 부호 있는 지수함수 나머지와 유한 가우스 모멘트에 의존하며, 평형 잔차를 k로 나누는 논증이 아닙니다.

경계 문제의 출력. ε=0.1, a=1, f=0, A=0, B=1에서 제안된 u=1은 내부 방정식 잔차가0이지만 경계 오차는1입니다. 반면 합성해1−exp(−10x)의 균일 값 오차는 exp(−10)으로 증명됩니다. 잔차 점수를 믿기 전에 어떤 경계조건을 검사했는지 명시해야 합니다.

장시간 운동의 출력. g=0.01, H=1,000,000에서 유계인 Lindstedt 1차 파형의 기본 위상 누적은 약8.16382라디안입니다. 유계성·정확한 초기조건·작은 잔차는 서로 다른 검사이며, 어느 하나만으로 긴 시간 뒤 봉우리의 시각이 맞는다고 결론낼 수 없습니다.

연구 비교에서는 먼저 관측량과 노름을 정하고, 매개변수와 시간·공간 범위, 증명의 가정과 기준해의 불확실성을 적으세요. 경계 결함·위상 누적·절단 나머지·안정성 상수를 보고서에 남깁니다. 이 네 문제의 진단은 입문 종합이며 WKB/Stokes, 일반 평균화 정리, 재합성, 연속체 장이론의 존재 증명을 제공하지 않습니다.

이 계열은 신경 ODE, 학습된 적분기, 축약 모형과 긴 rollout의 통제된 기준 문제입니다. 초기조건·불변량·국소 모양·미분·기본 주파수·누적 시각을 점검하세요. 짧은 rollout 성공이 결합을 줄이면서 늘어나는 관측 길이에서도 성공한다는 뜻은 아닙니다.

세속항과 시간 재조직의 배경은 Hunter Chapter5를 참고하세요. 비선형 기준은 NIST DLMF §22.19(ii), AGM 구성, cn Fourier 급수를 사용합니다. 자료의 modulus 표기와 이 장의 parameter 표기를 구분합니다.

이 예제까지 계획한 섭동 처리 유형을 다뤘더라도 구현·실제 reader 검증과 원주제 범위 감사를 마쳐야 입문 완료를 판정합니다. 일반 평균화 정리, WKB·Stokes, 비선형 공진과 장이론 재합성은 더 깊은 학습이 남습니다. 전체 교육과정에서는 이후 상관함수와 장이론·다중스케일의 후속 주제를 계속합니다.