跳到主要内容

神戸大学 システム情報学研究科 2025年度 第一期 数学 3

Author

祭音Myyura

Description

独立変数 xx,従属変数 yy に関する微分方程式

y+a(x)y+b(x)y=0\begin{align} y'' + a(x)y' + b(x)y = 0 \tag{*} \end{align}

の二つの解 y=y1(x),y2(x)y=y_1(x),y_2(x) に対し,ロンスキアン W(x)W(x) を次式で定義する:

W(x)=y1(x)y2(x)y1(x)y2(x)W(x)= \begin{vmatrix} y_1(x) & y_2(x) \\ y_1'(x) & y_2'(x) \end{vmatrix}

(1) W(x)W(x) が微分方程式

W=a(x)W W' = -a(x)W

を満たすことを示せ。

(2) x0x \geq 0 とする。a(x),b(x)a(x), b(x) および微分方程式 ()(*) の解の一つ y1(x)y_1(x)

a(x)=3x+43x+1,b(x)=6x+53x+1,y1(x)=ex a(x)=-\frac{3x+4}{3x+1},\qquad b(x)=-\frac{6x+5}{3x+1},\qquad y_1(x)=e^{-x}

で与える。また y2(0)=0, y2(0)=1y_2(0)=0,\ y_2'(0)=1 とする。W(x)W(x)y2(x)y_2(x) を求めよ。

题目描述

对关于自变量 xx、因变量 yy 的微分方程

y+a(x)y+b(x)y=0(*)y''+a(x)y'+b(x)y=0 \tag{*}

的两个解 y1(x),y2(x)y_1(x),y_2(x),定义 Wronskian

W(x)=y1(x)y2(x)y1(x)y2(x).W(x)= \begin{vmatrix} y_1(x)&y_2(x)\\ y_1'(x)&y_2'(x) \end{vmatrix}.
  1. 证明 W(x)W(x) 满足

    W=a(x)W.W'=-a(x)W.
  2. x0x\geq0

    a(x)=3x+43x+1,b(x)=6x+53x+1,a(x)=-\frac{3x+4}{3x+1},\qquad b(x)=-\frac{6x+5}{3x+1},

    且方程 ()(*) 的一个解为

    y1(x)=ex.y_1(x)=e^{-x}.

    另一解满足

    y2(0)=0,y2(0)=1.y_2(0)=0,\qquad y_2'(0)=1.

    W(x)W(x)y2(x)y_2(x)

考点

  • Wronskian 与 Abel 恒等式:对行列式求导并用原微分方程消去二阶导数。
  • 一阶线性方程:由初始 Wronskian 解 W=aWW'=-aW
  • 降阶法:已知一个非零解时,利用 W=y1y2y1y2W=y_1y_2'-y_1'y_2 建立并求解 y2y_2 的一阶方程。

Kai

(1)

W(x)=y1y2y1y2W(x)=y_1y_2'-y_1'y_2

なので

W(x)=y1y2+y1y2y1y2y1y2=y1y2y1y2.W'(x) = y_1'y_2' + y_1y_2'' - y_1''y_2 - y_1'y_2' = y_1y_2'' - y_1''y_2.

微分方程式 ()(*) から

yj=a(x)yjb(x)yj(j=1,2)y_j'' = -a(x)y_j' - b(x)y_j \quad (j=1,2)

だから

W=y1(ay2by2)(ay1by1)y2=a(y1y2y1y2)=a(x)W.\begin{aligned} W' &= y_1(-a y_2' - b y_2) - (-a y_1' - b y_1)y_2 \\ &= -a(y_1y_2' - y_1'y_2) \\ &= -a(x)W. \end{aligned}

したがって W=a(x)WW'=-a(x)W が示された。

(2)

まず (W(x))

(1) より

W=aW,a(x)=3x+43x+1.W'=-aW,\qquad a(x)=-\frac{3x+4}{3x+1}.

したがって

WW=a(x)=3x+43x+1.\frac{W'}{W}=-a(x)=\frac{3x+4}{3x+1}.

積分して

WWdx=3x+43x+1dxlnW=x+ln(3x+1)+C.\int \frac{W'}{W}dx = \int \frac{3x+4}{3x+1}dx \quad\Rightarrow\quad \ln W = x+\ln(3x+1)+C.

よって

W(x)=Cex(3x+1).W(x)=C e^{x}(3x+1).

初期条件から定数 CC を決める。 x=0x=0 のとき

W(0)=1011=1.W(0)= \begin{vmatrix} 1 & 0\\ -1& 1 \end{vmatrix} =1.

一方

W(0)=Ce0(30+1)=C.W(0)=C e^{0}(3\cdot 0+1)=C.

ゆえに C=1C=1

W(x)=ex(3x+1).\boxed{W(x)=e^{x}(3x+1)}.

次に y2(x)y_2(x)

ロンスキアンの定義

W=y1y2y1y2W = y_1y_2' - y_1'y_2

y1=ex, y1=ex, W=ex(3x+1)y_1=e^{-x},\ y_1'=-e^{-x},\ W=e^{x}(3x+1) を代入すると

exy2(ex)y2=ex(y2+y2)=ex(3x+1).e^{-x}y_2' - (-e^{-x})y_2 = e^{-x}(y_2'+y_2) = e^{x}(3x+1).

両辺に exe^{x} を掛けて

y2+y2=e2x(3x+1).y_2' + y_2 = e^{2x}(3x+1).

これは一次線形微分方程式。解いていく:

y2+y2=e2x(3x+1).y_2' + y_2 = e^{2x}(3x+1).

左辺の積分因子は exe^{x} なので

(exy2)=exe2x(3x+1)=e3x(3x+1).(e^{x}y_2)' = e^{x}e^{2x}(3x+1)=e^{3x}(3x+1).

積分すると

exy2=e3x(3x+1)dx+C=xe3x+C,e^{x}y_2 = \int e^{3x}(3x+1)dx + C = x e^{3x} + C,

よって

y2(x)=xe2x+Cex.y_2(x) = x e^{2x} + C e^{-x}.

初期条件 y2(0)=0y_2(0)=0 から

0=0e0+Ce0C=0,0 = 0\cdot e^{0} + C e^{0} \Rightarrow C=0,

したがって

y2(x)=xe2x.\boxed{y_2(x)=x e^{2x}}.

(このとき y2(x)=e2x(2x+1)y_2'(x)=e^{2x}(2x+1) なので y2(0)=1y_2'(0)=1 も満たしている。)