跳到主要内容

東京工業大学 工学院 システム制御系 2020年度 筆答専門試験(2019年8月実施)問題4

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

問1​

関数 x(t)x(t) に関する微分方程式の初期値問題

dx(t)dt=ax(t)+f(t),x(0)=x0,t≥0\frac{dx(t)}{dt}=ax(t)+f(t),\qquad x(0)=x_0,\qquad t\ge0

を考える。ただし a,x0a,x_0 は実定数、f(t)f(t) は tt の連続な関数である。

(1) f(t)=0f(t)=0 のとき、x(t)x(t) を a,x0a,x_0 を用いて表せ。

(2) f(t)=sin⁡tf(t)=\sin t、x0=0x_0=0 のとき、x(t)x(t) を求めよ。

(3) T>0T>0、zz を実定数として、

f(t)=ea(T−t)(∫0Te2aτ dτ)−1(z−eaTx0)f(t)=e^{a(T-t)}\left(\int_0^Te^{2a\tau}\,d\tau\right)^{-1}(z-e^{aT}x_0)

のとき、x(T)x(T) を求めよ。

問2​

関数 x(t)x(t) に関する微分方程式の初期値問題

d2x(t)dt2+3dx(t)dt+2x(t)=c,x(0)=x0,x′(0)=v0,t≥0\frac{d^2x(t)}{dt^2}+3\frac{dx(t)}{dt}+2x(t)=c, \qquad x(0)=x_0,\qquad x'(0)=v_0,\qquad t\ge0

を考える。ただし x0,v0,cx_0,v_0,c は実定数である。

(1) x0=3,v0=0,c=0x_0=3,v_0=0,c=0 のとき、x(t)x(t) を求めよ。

(2) c=4c=4 のとき、lim⁡t→∞x(t)\lim_{t\to\infty}x(t) を求めよ。

题目描述​

问1:考虑 x′=ax+f(t)x'=ax+f(t)、x(0)=x0x(0)=x_0、t≥0t\ge0,其中 a,x0a,x_0 为实常数,ff 连续。(1) f=0f=0 时求 x(t)x(t)。(2) f=sin⁡t,x0=0f=\sin t,x_0=0 时求解。(3) T>0,zT>0,z 为实常数,且 f(t)=ea(T−t)(∫0Te2aτ dτ)−1(z−eaTx0)f(t)=e^{a(T-t)}(\int_0^Te^{2a\tau}\,d\tau)^{-1}(z-e^{aT}x_0),求 x(T)x(T)。

问2:考虑 x′′+3x′+2x=cx''+3x'+2x=c、x(0)=x0x(0)=x_0、x′(0)=v0x'(0)=v_0、t≥0t\ge0。(1) x0=3,v0=0,c=0x_0=3,v_0=0,c=0 时求解。(2) c=4c=4 时求 t→∞t\to\infty 的极限。

Kai​

問1​

積分因子 e−ate^{-at} を用いると、初期値問題の解は

x(t)=eatx0+∫0tea(t−τ)f(τ) dτx(t)=e^{at}x_0+\int_0^te^{a(t-\tau)}f(\tau)\,d\tau

である。

(1) f=0f=0 を代入して、

x(t)=x0eat.\boxed{x(t)=x_0e^{at}}.

(2) f(t)=sin⁡t,x0=0f(t)=\sin t,x_0=0 より、

x(t)=eat∫0te−aτsin⁡τ dτ=eat−cos⁡t−asin⁡ta2+1.x(t)=e^{at}\int_0^te^{-a\tau}\sin\tau\,d\tau =\boxed{\frac{e^{at}-\cos t-a\sin t}{a^2+1}}.

(3) I=∫0Te2aτ dτI=\int_0^Te^{2a\tau}\,d\tau とおくと、T>0T>0 より I>0I>0 である。したがって、

x(T)=eaTx0+z−eaTx0I∫0Te2a(T−τ) dτ=eaTx0+z−eaTx0I∫0Te2as ds=z.\begin{aligned} x(T) &=e^{aT}x_0+\frac{z-e^{aT}x_0}{I} \int_0^Te^{2a(T-\tau)}\,d\tau\\ &=e^{aT}x_0+\frac{z-e^{aT}x_0}{I} \int_0^Te^{2as}\,ds =\boxed{z}. \end{aligned}

ここで s=T−τs=T-\tau と置換した。この計算は a=0a=0 の場合にも成立する。

問2​

(1) 特性方程式は λ2+3λ+2=(λ+1)(λ+2)=0\lambda^2+3\lambda+2=(\lambda+1)(\lambda+2)=0 であるから、x=Ae−t+Be−2tx=Ae^{-t}+Be^{-2t}。初期条件より

A+B=3,−A−2B=0A+B=3,\qquad -A-2B=0

なので A=6,B=−3A=6,B=-3 を得る。したがって、

x(t)=6e−t−3e−2t.\boxed{x(t)=6e^{-t}-3e^{-2t}}.

(2) 定数特解は x=2x=2。一般解と初期条件から、

x(t)=2+(2x0+v0−4)e−t+(2−x0−v0)e−2t.x(t)=2+(2x_0+v_0-4)e^{-t}+(2-x_0-v_0)e^{-2t}.

両指数項は t→∞t\to\infty で 00 となるので、初期値によらず

lim⁡t→∞x(t)=2.\boxed{\lim_{t\to\infty}x(t)=2}.