markdown
二階線形 PDE の分類md 424f35f
lecture/math/partial-differential-equations/classification-of-second-order-linear-pdes.lecture.n.md
Download PDF

二階線形にかいせんけい PDE の分類ぶんるい

date2026-07-16document_iddoc_81d376f88ff17ce0094b1c94c0e9b6e1description二変数の二階線形 PDE を主部の符号で点ごとに分類し、特性方向・標準形・典型的性質の適用条件を整理する。prerequisitesPDE とは何か / PDE の初期値問題と境界値問題 / 接平面・連鎖律・Jacobian / 二次形式と正定値行列type講義content_typelecturestatusactiverelateddata/lecture/math/partial-differential-equations/initial-and-boundary-value-problems.lecture.n.md / data/lecture/math/partial-differential-equations/method-of-characteristics.lecture.n.md / data/lecture/math/partial-differential-equations/heat-wave-and-laplace-equations.lecture.n.md / data/lecture/math/vector-calculus/laplacian-and-physical-applications.lecture.n.md / data/lecture/math/linear-algebra/quadratic-forms-and-positive-definite-matrices.lecture.n.md
mathpartial-differential-equationsclassificationlecture

1導入どうにゅう

この講義こうぎでは、二変数にへんすう二階線形にかいせんけい PDE を最高階微分さいこうかいびぶんからてんごとに分類ぶんるいし、特性方向とくせいほうこう標準形ひょうじゅんけい判定はんていする方法ほうほう説明せつめいする。楕円型だえんがた放物型ほうぶつがた双曲型そうきょくがたという名前なまえ典型的てんけいてき平衡へいこう拡散かくさん伝播でんぱ予想よそうする手掛てがかりになるが、かただけで存在そんざい一意性いちいせい安定性あんていせい保証ほしょうするものではない。

この判別量はんべつりょうによる三分類さんぶんるい独立変数どくりつへんすうが 2 場合ばあいあつかう。よりおおくの変数へんすうでは、主係数行列しゅけいすうぎょうれつ固有値こゆうち符号ふごうかぞえる一般化いっぱんか必要ひつようになる。

2一般形いっぱんけい主部しゅぶ

領域りょういき ΩR2 で、実数値じっすうち連続係数れんぞくけいすう方程式ほうていしき

Auxx+2Buxy+Cuyy+Dux+Euy+Fu=G

かんがえる。A,,F,G(x,y)既知関数きちかんすうであり、u やその微分びぶんには依存いぞんしない。このことが uかんする線形性せんけいせい意味いみする。分類ぶんるい使つか主部しゅぶPrincipal part

Auxx+2Buxy+Cuyy

であり、主係数行列しゅけいすうぎょうれつ

M(x,y)=(ABBC)

である。D,E,F,Gかい具体形ぐたいけい適切性てきせつせい影響えいきょうするが、かたM だけでまる。

3判別量はんべつりょう三分類さんぶんるい

てん (x,y)

δ=B2-AC=-detM

く。M0、すなわち A,B,C がすべて 0 ではないてんでは、つぎのように分類ぶんるいする。

判別条件はんべつじょうけんかたM符号構造ふごうこうぞう実特性方向じつとくせいほうこう
δ<0楕円型だえんがたElliptic正定値せいていちまたは負定値ふていちなし
δ=0 かつ M0放物型ほうぶつがたParabolicrank 1 の半正定値はんせいていちまたは半負定値はんふていち重複ちょうふくする 1 方向ほうこう
δ>0双曲型そうきょくがたHyperbolic不定値ふていちことなる 2 方向ほうこう

A=B=C=0 なら δ=0 でも二階主部にかいしゅぶがなく、放物型ほうぶつがたとはべない。係数けいすう変化へんかする場合ばあい、これは領域全体りょういきぜんたいではなく各点かくてんでの分類ぶんるいである。

4代表方程式だいひょうほうていしき計算けいさん

Laplace 方程式ほうていしき

uxx+uyy=0

では A=C=1,B=0 なので δ=-1<0 であり、楕円型だえんがたである。

熱方程式ねつほうていしき(x,t)じゅん

ut-κuxx=0,κ>0

くと、A=-κ,B=C=0δ=0M は rank 1 なので放物型ほうぶつがたである。ut一階項いっかいこうなのでかた判定はんていにははいらないが、時間発展じかんはってんきをめる重要じゅうようこうである。

波動方程式はどうほうていしき

utt-c2uxx=0,c>0

くと、A=-c2,B=0,C=1 なので δ=c2>0 であり、双曲型そうきょくがたである。

5座標変換ざひょうへんかんかたたもたれる理由りゆう

あたらしい座標関数ざひょうかんすう (ξ,η)=(ξ(x,y),η(x,y))C2 きゅうとし、その Jacobian 行列ぎょうれつJ とする。detJ0てんでは逆関数定理ぎゃくかんすうていりにより局所座標きょくしょざひょうとして使つかえ、新座標しんざひょうでの主係数行列しゅけいすうぎょうれつ

M~=JMJT

という合同変換ごうどうへんかんける。したがって

detM~=(detJ)2detM

となり、δ=-detM符号ふごうわらない。可変係数かへんけいすう方程式ほうていしき変換へんかんすると座標ざひょう二階微分にかいびぶんから低階項ていかいこうしょうじるが、主部しゅぶかた非特異ひとくい座標変換ざひょうへんかんたもたれる。特異とくい変換へんかんではこの結論けつろん使つかえない。

6特性曲線とくせいきょくせんとの関係かんけい

曲線きょくせん ϕ(x,y)=constant垂直すいちょく法線方向ほうせんほうこうは、勾配こうばい (ϕx,ϕy)あらわせる。この曲線きょくせん主部しゅぶたいして特性的とくせいてきCharacteristicである条件じょうけん

Aϕx2+2Bϕxϕy+Cϕy2=0

である。接方向せつほうこう(dx,dy)けば、おな条件じょうけん同次式どうじしき

A(dy)2-2Bdxdy+C(dx)2=0

になる。このかたち垂直すいちょく曲線きょくせんあつかえる。曲線きょくせんy=y(x)け、m=dy/dx とすれば法線ほうせん(-m,1)比例ひれいするので、

Am2-2Bm+C=0

る。A0 なら、これは有限ゆうげんかたむm二次方程式にじほうていしきであり、判別式はんべつしき4(B2-AC)=4δ である。A=0 では二次項にじこうえ、B0 なら一次式いちじしきB=0 なら定数式ていすうしきになるため、dx=0垂直方向すいちょくほうこうかたむmあらわせない場合ばあいがある。そのときは同次式どうじしきもどるか、x=x(y)表示ひょうじ使つかう。射影的しゃえいてき方向全体ほうこうぜんたいかぞえれば、うえ三分類さんぶんるい実特性方向じつとくせいほうこう本数ほんすう一致いっちする。

熱方程式ねつほうていしきでは m=dt/dx=0重根じゅうこんなので、t=constant二階主部にかいしゅぶたいする特性曲線とくせいきょくせんである。それでも t=0 に 1 初期値しょきちいて前向まえむ発展はってんかんがえられるのは、ut一階いっかいのこ放物型発展構造ほうぶつがたはってんこうぞうによる。前講義ぜんこうぎ一階準線形いっかいじゅんせんけい PDE における非特性条件ひとくせいじょうけんと、この二階主部にかいしゅぶ特性条件とくせいじょうけんおな規則きそくとしてあつかってはならない。

data/lecture/math/partial-differential-equations/method-of-characteristics.lecture.n.md

7標準形ひょうじゅんけいへの入口いりぐち

一点いってんでは、非特異ひとくい線形座標変換せんけいざひょうへんかんと、必要ひつようなら方程式全体ほうていしきぜんたいへの -1 ばいにより、主部しゅぶつぎかたち正規化せいきかできる。

かた主部しゅぶ標準形ひょうじゅんけい
楕円型だえんがたuξξ+uηη
放物型ほうぶつがたuξξ
双曲型そうきょくがたuξξ-uηη、または特性座標とくせいざひょうuξη

可変係数かへんけいすう方程式ほうていしき近傍全体きんぼうぜんたい標準形ひょうじゅんけいうつすには、係数けいすう正則性せいそくせいかた一定いっていであること、局所座標きょくしょざひょう存在そんざいべつ確認かくにんする。大域的たいいきてき標準化ひょうじゅんか自動的じどうてきにはしたがわない。

8かたから予想よそうできることとできないこと

かた問題設定もんだいせってい解法かいほうえら入口いりぐちである。ただしつぎ対応たいおうは、代表的だいひょうてき追加仮定ついかかていのもとで成立せいりつする。

かた典型像てんけいぞう必要ひつようになる追加条件ついかじょうけんれい
楕円型だえんがた境界値きょうかいち内部平衡ないぶへいこう、maximum principle一様楕円性いちようだえんせい低階係数ていかいけいすう符号ふごう領域りょういき境界条件きょうかいじょうけん
放物型ほうぶつがた前向まえむ時間発展じかんはってん平滑化へいかつか拡散係数かくさんけいすうせい下限かげん時間方向じかんほうこう係数けいすう初期境界条件しょききょうかいじょうけん
双曲型そうきょくがた有限速度ゆうげんそくど伝播でんぱと energyことなる実特性方向じつとくせいほうこうたもたれること、適切てきせつ時間面じかんめん係数けいすう正則性せいそくせい

領域りょういき一様楕円性いちようだえんせいUniform ellipticity仮定かていするとは、領域全体りょういきぜんたい固定こていした σ{1,-1}定数ていすう 0<λ0[PARSE ERROR: Undefined("Command(\"le\")")]Λ0< があり、すべての ζR2たいして

λ0|ζ|2[PARSE ERROR: Undefined("Command(\"le\")")]ζT(σM)ζ[PARSE ERROR: Undefined("Command(\"le\")")]Λ0|ζ|2

各点かくてん成立せいりつすることを意味いみする。これは固有値こゆうち符号ふごう領域全体りょういきぜんたいでそろえ、0 と無限大むげんだい両方りょうほうから一様いちようへだてる条件じょうけんである。このような一様性いちようせいは、かたごとの評価ひょうか領域全体りょういきぜんたい使つかうときに重要じゅうようである。

たとえば前向まえむ熱方程式ねつほうていしきうし熱方程式ねつほうていしき主部しゅぶおな放物型分類ほうぶつがたぶんるいつが、安定性あんていせいおなじではない。また楕円型だえんがたというだけで、任意にんい境界条件きょうかいじょうけん一意解いちいかいさだめるとはかぎらない。初期条件しょきじょうけん境界条件きょうかいじょうけん適切性てきせつせいは、かたとはべつ問題全体もんだいぜんたい確認かくにんする。

data/lecture/math/partial-differential-equations/initial-and-boundary-value-problems.lecture.n.md

9反例はんれいまたは限界げんかい

係数けいすう場所ばしょ変化へんかする方程式ほうていしきでは、領域全体りょういきぜんたいかた一定いっていとはかぎらない。たとえば

yuxx+uyy=0

では A=y,B=0,C=1 なので δ=-y である。したがって y>0楕円型だえんがたy<0双曲型そうきょくがたy=0 で rank 1 の放物型ほうぶつがたになる。y=0型変化線かたへんかせんであり、一様性いちようせいうしなわれるため、片側かたがわ使つかえる評価ひょうか解法かいほうをそのまま横断的おうだんてき適用てきようできない。

10よくあるあやま

  • 2Buxy という係数規約けいすうきやく確認かくにんせず、B2-ACべつB代入だいにゅうする。交差項こうさこうBuxy規約きやくでは、主係数行列しゅけいすうぎょうれつ非対角成分ひたいかくせいぶんB/2 になる。
  • δ=0 だけをて、A=B=C=0退化たいかしたしきまで放物型ほうぶつがたぶ。
  • かた確認かくにんしただけで、境界条件きょうかいじょうけん時間方向じかんほうこう係数けいすう正則性せいそくせい調しらべずに適切性てきせつせい結論けつろんする。
  • 一点いってんでの分類ぶんるい領域全体りょういきぜんたいひろげ、型変化かたへんか一様性いちようせい見落みおとす。

11関連かんれんリンク

data/lecture/math/partial-differential-equations/heat-wave-and-laplace-equations.lecture.n.md

Classification of Second-Order Linear PDEs

1Introduction

This lecture explains how to classify a second-order linear PDE in two variables pointwise from its highest-order derivatives and how to determine its characteristic directions and canonical form. The terms elliptic, parabolic, and hyperbolic suggest typical equilibrium, diffusion, and propagation behavior, but the type alone does not guarantee existence, uniqueness, or stability.

The discriminant classification presented here applies to two independent variables. In higher dimensions, it is generalized by counting the signs of the eigenvalues of the principal coefficient matrix.

2General form and principal part

On a domain ΩR2, consider an equation with real-valued continuous coefficients,

Auxx+2Buxy+Cuyy+Dux+Euy+Fu=G.

Here A,,F,G are prescribed functions of (x,y) and do not depend on u or its derivatives; this is the meaning of linearity in u. The {principal part} used for classification is

Auxx+2Buxy+Cuyy,

with principal coefficient matrix

M(x,y)=(ABBC).

The terms D,E,F,G affect the particular form and well-posedness of solutions, but the type is determined by M alone.

3Discriminant and the three types

At a point (x,y), set

δ=B2-AC=-detM.

At points where M0, equivalently where A,B,C are not all zero, the classification is as follows.

Discriminant conditionTypeSign structure of MReal characteristic directions
δ<0{Elliptic}positive or negative definitenone
δ=0 and M0{Parabolic}positive or negative semidefinite of rank 1one repeated direction
δ>0{Hyperbolic}indefinitetwo distinct directions

If A=B=C=0, then the second-order principal part is absent, so the equation is not called parabolic even though δ=0. When the coefficients vary, the classification is pointwise rather than automatically valid throughout the domain.

4Computations for representative equations

For Laplace's equation

uxx+uyy=0,

one has A=C=1 and B=0, so δ=-1<0 and the equation is elliptic.

Write the heat equation in the variable order (x,t) as

ut-κuxx=0,κ>0.

Then A=-κ and B=C=0, so δ=0 and M has rank 1; the equation is parabolic. The term ut is first order and therefore does not enter the classification, although it is essential in determining the direction of time evolution.

For the wave equation

utt-c2uxx=0,c>0,

one has A=-c2, B=0, and C=1, so δ=c2>0 and the equation is hyperbolic.

5Why the type is preserved under coordinate changes

Let the new coordinate functions (ξ,η)=(ξ(x,y),η(x,y)) be of class C2, and let J be their Jacobian matrix. At a point where detJ0, the inverse function theorem makes them local coordinates, and the principal coefficient matrix in the new coordinates is transformed by congruence:

M~=JMJT.

Consequently,

detM~=(detJ)2detM,

so the sign of δ=-detM is unchanged. Transforming a variable-coefficient equation also produces lower-order terms from second derivatives of the coordinate functions, but the type of the principal part is preserved by a nonsingular coordinate change. This conclusion is unavailable for a singular transformation.

6Relation to characteristic curves

The normal direction to a curve ϕ(x,y)=constant is the gradient (ϕx,ϕy). The curve is {characteristic} for the principal part when

Aϕx2+2Bϕxϕy+Cϕy2=0.

Writing the tangent direction as (dx,dy) gives the equivalent homogeneous condition

A(dy)2-2Bdxdy+C(dx)2=0.

This formulation also includes vertical curves. If the curve is written as y=y(x) and m=dy/dx, its normal is proportional to (-m,1), and hence

Am2-2Bm+C=0.

When A0, this is a quadratic equation for the finite slope m, with discriminant 4(B2-AC)=4δ. If A=0, the quadratic term disappears; the equation is linear when B0 and constant when B=0. A vertical direction dx=0 may then be invisible in the slope coordinate m. One must return to the homogeneous equation or use a representation x=x(y). Counting all projective directions yields the numbers of real characteristic directions in the classification table.

For the heat equation, m=dt/dx=0 is a repeated root, so t=constant is characteristic for the second-order principal part. Nevertheless, one initial value can be prescribed at t=0 for forward evolution because the first-order term ut remains and creates the parabolic evolution structure. The non-characteristic condition for a first-order quasilinear PDE in the preceding lecture must not be identified with this characteristic condition for a second-order principal part.

data/lecture/math/partial-differential-equations/method-of-characteristics.lecture.n.md

7Entry point to canonical forms

At a single point, a nonsingular linear coordinate transformation, together with multiplication of the entire equation by -1 when necessary, reduces the principal part to one of the following forms.

TypeCanonical form of the principal part
Ellipticuξξ+uηη
Parabolicuξξ
Hyperbolicuξξ-uηη, or uξη in characteristic coordinates

Reducing a variable-coefficient equation to canonical form throughout a neighborhood additionally requires sufficient coefficient regularity, constancy of type, and the existence of suitable local coordinates. A global reduction does not follow automatically.

8What the type does and does not predict

The type is an entry point for selecting a problem formulation and a method of solution. The following correspondences hold only under representative additional assumptions.

TypeTypical pictureExamples of additional conditions
Ellipticboundary values and interior equilibrium; maximum principlesuniform ellipticity, signs of lower-order coefficients, domain and boundary conditions
Parabolicforward time evolution and smoothinga positive lower bound for the diffusion coefficient, time direction, coefficients and initial-boundary data
Hyperbolicfinite-speed propagation and energy methodspreservation of distinct real characteristic directions, an appropriate time surface, coefficient regularity

{Uniform ellipticity} on a domain means that there are a fixed σ{1,-1} and constants 0<λ0[PARSE ERROR: Undefined("Command(\"le\")")]Λ0< such that, at every point and for every ζR2,

λ0|ζ|2[PARSE ERROR: Undefined("Command(\"le\")")]ζT(σM)ζ[PARSE ERROR: Undefined("Command(\"le\")")]Λ0|ζ|2.

This condition makes the eigenvalue signs consistent throughout the domain and keeps their magnitudes uniformly separated from both zero and infinity. Such uniformity is important when estimates associated with a type are applied over an entire domain.

For example, the forward and backward heat equations have the same parabolic principal-part classification, but their stability properties differ. Likewise, ellipticity alone does not imply that arbitrary boundary conditions determine a unique solution. Initial conditions, boundary conditions, and well-posedness must be assessed for the complete problem separately from its type.

data/lecture/math/partial-differential-equations/initial-and-boundary-value-problems.lecture.n.md

9Counterexample and limitation

For an equation whose coefficients vary with position, the type need not be constant throughout the domain. For example,

yuxx+uyy=0

has A=y, B=0, and C=1, so δ=-y. It is elliptic for y>0, hyperbolic for y<0, and parabolic of rank 1 for y=0. The line y=0 is a type-change line where uniformity is lost; estimates and methods valid on one side cannot automatically be applied across it.

10Common errors

  • Substituting the wrong coefficient into B2-AC without checking the convention 2Buxy. Under the convention Buxy, the off-diagonal entries of the principal coefficient matrix are B/2.
  • Calling a degenerate equation with A=B=C=0 parabolic merely because δ=0.
  • Concluding well-posedness from the type without examining boundary conditions, time direction, or coefficient regularity.
  • Extending a pointwise classification to an entire domain while ignoring type changes or loss of uniformity.

11Related links

data/lecture/math/partial-differential-equations/heat-wave-and-laplace-equations.lecture.n.md
raw .n.md をコピー
loc をコピー (filepath:line ~ line)
copy share link
copy encoded share link
path をコピー
copy share link
copy encoded share link
copy share link
copy encoded share link
タブを全て閉じる