院試hub

東大 学際情報学府 学際情報学専攻 先端表現情報学コース 専門科目 2016年度 院試 過去問 解答例・解説(全5問)

全5問の分野は線形代数1問、微分積分・解析1問、力学1問です。固有値・固有ベクトル、フーリエ変換、固有振動数と振動系などが出ています。2015年度と共通するテーマはフーリエ変換です。解説・解答・最終答まで、このページですべて公開しています。問題本文は含みません。

最終更新:

第1問 — ロジスティックモデル

ロジスティック解の検算

t=0t=0 を代入すると x(0)=x0x(0)=x_0、t→∞t\to\infty で x(t)→Kx(t)\to K となる。問題の条件 0<x<K0<x<K と、増加して環境収容力 KK に近づく挙動が一致している。

平衡点の扱い

連立系には複数の平衡点がある。減衰振動の議論は、正の x,yx,y をもつ共存平衡点まわりで行うのが自然である。そのため、まず rb>cdrb>cd により y∗>0y^*>0 が必要であることを明示する。

線形安定性の条件

2次元線形系では、跡が負、行列式が正なら安定である。さらに振動を伴うには固有値が複素数、つまり判別式 τ2−4Δ\tau^2-4\Delta が負でなければならない。この3点を分けて書くと、条件の意味が明確になる。

解答

  1. 方程式は dxdt=rx−rKx2=rx(1−xK) \frac{dx}{dt}=rx-\frac{r}{K}x^2 =rx\left(1-\frac{x}{K}\right) である。変数分離すると dxx(1−x/K)=r dt \frac{dx}{x(1-x/K)}=r\,dt であり、初期値 x(0)=x0x(0)=x_0 を用いると x(t)=K1+(K−x0x0)e−rt=Kx0ertK+x0(ert−1) x(t)=\frac{K}{1+\left(\frac{K-x_0}{x_0}\right)e^{-rt}} =\frac{Kx_0e^{rt}}{K+x_0(e^{rt}-1)} となる。

  2. 連立方程式は dxdt=x(r−cx−ay),dydt=y(bx−d) \frac{dx}{dt}=x(r-cx-ay),\qquad \frac{dy}{dt}=y(bx-d) である。平衡点は (0,0),(rc,0) (0,0),\qquad \left(\frac{r}{c},0\right) および、正の共存平衡点が存在する場合 x∗=db,y∗=r−cd/ba=rb−cdab x^*=\frac{d}{b},\qquad y^*=\frac{r-cd/b}{a} =\frac{rb-cd}{ab} である。共存平衡点が第1象限にある条件は rb>cdrb>cd である。

  3. Jacobi 行列は J(x,y)=(r−2cx−ay−axbybx−d). J(x,y)= \begin{pmatrix} r-2cx-ay&-ax\\ by&bx-d \end{pmatrix}. 共存平衡点 (x∗,y∗)(x^*,y^*) では J∗=(−cd/b−ad/b(rb−cd)/a0). J^*= \begin{pmatrix} -cd/b&-ad/b\\ (rb-cd)/a&0 \end{pmatrix}. 跡と行列式は τ=tr⁡J∗=−cdb<0,Δ=det⁡J∗=d(rb−cd)b. \tau=\operatorname{tr}J^*=-\frac{cd}{b}<0,\qquad \Delta=\det J^*=\frac{d(rb-cd)}{b}. 減衰振動しながら収束するには、固有値が負の実部をもつ複素共役であればよい。したがって Δ>0,τ2−4Δ<0 \Delta>0,\qquad \tau^2-4\Delta<0 が必要十分である。すなわち rb>cd,(cdb)2<4d(rb−cd)b rb>cd,\qquad \left(\frac{cd}{b}\right)^2 <\frac{4d(rb-cd)}{b} である。

最終答

x(t)=K/{1+((K−x0)/x0)e−rt}x(t)=K/\{1+((K-x_0)/x_0)e^{-rt}\}。共存平衡点は (d/b,(rb−cd)/(ab))(d/b,(rb-cd)/(ab))。減衰振動収束条件は rb>cdrb>cd かつ (cd/b)2<4d(rb−cd)/b(cd/b)^2<4d(rb-cd)/b。

第2問 — 動的計画法

貪欲法ではない

大きい数から順に選ぶ貪欲法は一般には失敗する。例えば大きな値を選ぶことで、その両隣の合計がさらに大きい場合を失うことがある。この問題では、各位置で「選ぶ・選ばない」の2通りを最適部分構造として残す動的計画法が確実である。

解答

  1. 列 3,5,7,6,2,43,5,7,6,2,4 に対して、隣接しないように選ぶ。最大値は 5+6+4=15 5+6+4=15 であり、選ぶ部分集合は {5,6,4}\{5,6,4\} である。

  2. 先頭から ii 個まで見たときの最大合計値を SiS_i とする。空集合を許すなら S0=0,S1=max⁡(0,a1) S_0=0,\qquad S_1=\max(0,a_1) とおく。i≥2i\ge2 では、aia_i を選ばない場合は Si−1S_{i-1}、選ぶ場合は ai−1a_{i-1} を選べないので Si−2+aiS_{i-2}+a_i である。したがって Si=max⁡(Si−1,Si−2+ai)(i=2,…,n) S_i=\max(S_{i-1},S_{i-2}+a_i) \qquad (i=2,\ldots,n) を順に計算すれば、求める最大合計値は SnS_n である。

    擬似コードで書くと次のようになる。 s0←0,s1←max⁡(0,a1)for i=2 to n:s2←max⁡(s1, s0+ai)s0←s1,s1←s2return s1 \begin{array}{l} s0\leftarrow 0,\quad s1\leftarrow \max(0,a_1)\\ \text{for }i=2\text{ to }n:\\ \quad s2\leftarrow \max(s1,\ s0+a_i)\\ \quad s0\leftarrow s1,\quad s1\leftarrow s2\\ \text{return }s1 \end{array}

  3. 各 ii で定数回の比較と加算だけを行うため、時間計算量は O(n) O(n) である。上の擬似コードのように直前2個の値だけを保持すれば、追加メモリは O(1) O(1) で済む。選んだ要素自体も復元するなら、選択の有無を記録するため O(n)O(n) メモリを使う。

最終答

最大部分集合は {5,6,4}\{5,6,4\}、合計 1515。一般には Si=max⁡(Si−1,Si−2+ai)S_i=\max(S_{i-1},S_{i-2}+a_i) を用いる。時間計算量は O(n)O(n)。

第3問 — 無線LANとフーリエ変換

正規化の確認

本答案では逆変換側に 1/(2π)1/(2\pi) を置く約束を使った。この約束では、時間領域の積は周波数領域で 1/(2π)1/(2\pi) 倍の畳み込みになる。余弦の変換に π\pi が出るため、最終的に 1/21/2 のシフト和になる。

切片と零点

矩形パルスの X(0)X(0) は面積に等しいので 2t02t_0 である。三角波の Z(0)Z(0) も面積に等しいので t0t_0 である。原点の値を極限または面積で検算できる。

窓関数としての比較

矩形窓は時間幅の端で不連続なのでスペクトルの裾が大きく、漏れが強い。三角窓は端で連続的に0へ落ちるため漏れを抑えるが、主ローブが広がる。このトレードオフが窓関数の比較で最も重要である。

解答

  1. 2.4GHz帯は回折しやすく、壁などを通って遠くまで届きやすい。また古い機器を含め対応機器が多い。一方、5GHz帯は利用できるチャネルが多く、電子レンジやBluetoothなどとの干渉が比較的小さいため、近距離では高速・安定な通信を得やすい。

    1. ここでは角周波数 ω\omega に対して X(ω)=∫−∞∞x(t)e−jωt dt X(\omega)=\int_{-\infty}^{\infty}x(t)e^{-j\omega t}\,dt をフーリエ変換とする。対応する逆変換は x(t)=12π∫−∞∞X(ω)ejωt dω x(t)=\frac{1}{2\pi}\int_{-\infty}^{\infty} X(\omega)e^{j\omega t}\,d\omega である。

    2. 矩形パルスについて X(ω)=∫−t0t0e−jωt dt=2sin⁡(ωt0)ω=2t0sin⁡(ωt0)ωt0. X(\omega)=\int_{-t_0}^{t_0}e^{-j\omega t}\,dt =\frac{2\sin(\omega t_0)}{\omega} =2t_0\frac{\sin(\omega t_0)}{\omega t_0}. ω=0\omega=0 での値は極限により X(0)=2t0 X(0)=2t_0 であり、零点は ω=±nπt0(n=1,2,…) \omega=\pm \frac{n\pi}{t_0}\qquad (n=1,2,\ldots) である。概形は偶関数の sinc 形で、原点から振動しながら 1/∣ω∣1/|\omega| 程度で減衰する。

    3. y(t)=x(t)cos⁡ω0t y(t)=x(t)\cos\omega_0t である。時間領域の積は周波数領域の畳み込みに対応する。さらに F[cos⁡ω0t]=π{δ(ω−ω0)+δ(ω+ω0)} \mathcal F[\cos\omega_0t] =\pi\{\delta(\omega-\omega_0)+\delta(\omega+\omega_0)\} なので Y(ω)=12{X(ω−ω0)+X(ω+ω0)}. Y(\omega) =\frac12\{X(\omega-\omega_0)+X(\omega+\omega_0)\}. つまり、矩形パルスのスペクトルが ±ω0\pm\omega_0 だけ平行移動し、半分ずつ重なったものになる。

    4. 三角波は偶関数なので Z(ω)=2∫0t0(1−tt0)cos⁡ωt dt=2(1−cos⁡ωt0)t0ω2=t0(sin⁡(ωt0/2)ωt0/2)2. \begin{aligned} Z(\omega) &=2\int_0^{t_0}\left(1-\frac{t}{t_0}\right)\cos\omega t\,dt\\ &=\frac{2(1-\cos\omega t_0)}{t_0\omega^2}\\ &=t_0\left(\frac{\sin(\omega t_0/2)}{\omega t_0/2}\right)^2 . \end{aligned} ω=0\omega=0 では Z(0)=t0Z(0)=t_0 で、零点は ω=±2nπt0(n=1,2,…) \omega=\pm\frac{2n\pi}{t_0}\qquad(n=1,2,\ldots) である。矩形窓に比べて主ローブは広いが、スペクトルの裾は 1/ω21/\omega^2 程度で速く減衰し、サイドローブが小さい。したがって周波数分解能は下がるが、漏れは抑えやすい。

最終答

X(ω)=2sin⁡(ωt0)/ωX(\omega)=2\sin(\omega t_0)/\omega、Y(ω)={X(ω−ω0)+X(ω+ω0)}/2Y(\omega)=\{X(\omega-\omega_0)+X(\omega+\omega_0)\}/2、Z(ω)=t0{sin⁡(ωt0/2)/(ωt0/2)}2Z(\omega)=t_0\{\sin(\omega t_0/2)/(\omega t_0/2)\}^2。

第4問 — 材料力学

軸力の符号

外力全体は左向き 50+20=70 kN50+20=70\,\mathrm{kN}、右向き 30+40=70 kN30+40=70\,\mathrm{kN} でつり合っている。区間ごとの軸力は、切断面の片側だけを取り出してつり合いを書くと間違えにくい。

単位換算

N/mm2\mathrm{N/mm^2} は MPa\mathrm{MPa} と同じである。したがって 50 000 N/75 mm2=666.7 MPa50\,000\,\mathrm{N}/75\,\mathrm{mm^2}=666.7\,\mathrm{MPa} と直接読める。伸び計算では E=70 000 N/mm2E=70\,000\,\mathrm{N/mm^2} を使えば、長さの単位は mm\mathrm{mm} のままそろう。

体積変化の近似

微小変形では体積変化率を3方向ひずみの和で近似できる。通常の金属のポアソン比は 0<ν<0.50<\nu<0.5 なので、軸方向に伸びる量が横方向2方向の縮みを上回り、体積が増える。

解答

  1. 重ね合わせの原理とは、材料が線形弾性範囲にあり、変形が微小で幾何学的非線形性を無視できるとき、複数の荷重による応力・ひずみ・変位が各荷重単独の応答の和として表せるという原理である。弾性限度を超える塑性変形、大変形、接触条件の変化がある場合には一般に成立しない。

  2. 引張を正とする。各区間の軸力は、切断面の左側に作用する外力のつり合いから求める。 断面軸力垂直応力B (A–C)50 kN50 000/75=666.7 MPaD (C–E)20 kN20 000/75=266.7 MPaF (E–G)40 kN40 000/75=533.3 MPa \begin{array}{c|c|c} \text{断面} & \text{軸力} & \text{垂直応力}\\ \hline B\ (A\text{--}C) & 50\,\mathrm{kN} & 50\,000/75=666.7\,\mathrm{MPa}\\ D\ (C\text{--}E) & 20\,\mathrm{kN} & 20\,000/75=266.7\,\mathrm{MPa}\\ F\ (E\text{--}G) & 40\,\mathrm{kN} & 40\,000/75=533.3\,\mathrm{MPa} \end{array} いずれも引張応力である。

  3. 長さは AC=500 mmAC=500\,\mathrm{mm}、EG=300 mmEG=300\,\mathrm{mm}、全長 1000 mm1000\,\mathrm{mm} なので CE=200 mm CE=200\,\mathrm{mm} である。断面積 A=75 mm2A=75\,\mathrm{mm^2}、ヤング率 E=70 GPa=70 000 N/mm2 E=70\,\mathrm{GPa}=70\,000\,\mathrm{N/mm^2} を用いると、伸びは δ=∑NiLiEA=50 000⋅500+20 000⋅200+40 000⋅30070 000⋅75=7.81 mm. \begin{aligned} \delta &=\sum \frac{N_iL_i}{EA}\\ &=\frac{50\,000\cdot500+20\,000\cdot200+40\,000\cdot300} {70\,000\cdot75}\\ &=7.81\,\mathrm{mm}. \end{aligned}

  4. 一軸引張で軸方向の垂直ひずみを ε\varepsilon とすると、横ひずみはおよそ −νε -\nu\varepsilon である。微小変形では体積ひずみは3方向の垂直ひずみの和で εv≃ε−2νε=(1−2ν)ε \varepsilon_v\simeq \varepsilon-2\nu\varepsilon =(1-2\nu)\varepsilon となる。通常の金属ではポアソン比 ν<0.5\nu<0.5 なので、εv>0\varepsilon_v>0、すなわち体積は増加する。

  5. 鋼の応力ひずみ線図は、初期に直線的な弾性域を持ち、その後、降伏、加工硬化、最大応力、くびれ、破断へ進む。

    弾性限度は、荷重を除けば元に戻る範囲の上限である。降伏点は、応力がほぼ増えないまま塑性ひずみが進む点または領域である。引張強さは公称応力の最大値で、最大荷重を元の断面積で割った値に対応する。

最終答

断面応力は B:666.7 MPaB:666.7\,\mathrm{MPa}、D:266.7 MPaD:266.7\,\mathrm{MPa}、F:533.3 MPaF:533.3\,\mathrm{MPa} の引張。棒の伸びは約 7.81 mm7.81\,\mathrm{mm}。

第5問 — エレベータ待ち行列

待ち時間の定義

本問では、1階で停止している時間は待ち時間に含めない。したがって「出発時刻」ではなく「1階で扉が開く時刻」までを待ち時間として計算する。これにより、(3) では 10m10m ではなく 10m−210m-2、(4) では 11m−311m-3 が基準になる。

連続流としての扱い

到着が一様であるため、到着率 λ=2.2\lambda=2.2 人/分の連続流として計算している。60分後の待ち人数は、累計到着人数から、その時点までに出発したエレベータの輸送人数を引けばよい。

容量比の直感

需要容量比が1を超えると、待ち行列は解消せず時間とともに増える。(4) は定員が増えてもサイクルが長くなれば容量が下がる例であり、エレベータや交通施設では「1回の輸送量」だけでなく「周期」を同時に見る必要がある。

解答

  1. 1サイクルは、上昇3分、展望階停止2分、下降3分、1階停止2分の合計10分である。測定開始直後に出発するので、次に1階で扉が開くのは8分後である。待ち時間には1階停止中の2分を含めないため、到着時刻をサイクル内の t∈[0,10)t\in[0,10) とすると W(t)={8−t(0≤t<8),0(8≤t<10) W(t)= \begin{cases} 8-t & (0\le t<8),\\ 0 & (8\le t<10) \end{cases} である。一様到着より平均待ち時間は 110∫08(8−t) dt=3210=3.2 分. \frac1{10}\int_0^8(8-t)\,dt =\frac{32}{10} =3.2\,\text{分}.

  2. 到着率は λ=13260=2.2 人/分 \lambda=\frac{132}{60}=2.2\,\text{人/分} である。容量は10分ごとに20人なので μ=2010=2.0 人/分=120 人/時. \mu=\frac{20}{10}=2.0\,\text{人/分} =120\,\text{人/時}. したがって需要容量比は 132120=1.1 \frac{132}{120}=1.1 である。1サイクルあたり到着は 2.2×10=222.2\times10=22 人、輸送は20人なので、待ち人数は1サイクルごとに2人増える。60分後は6サイクル後であり、 6×2=12 人 6\times2=12\,\text{人} が1階で待っている。

  3. mm 回目の出発時刻は 10m10m 分、1階で扉が開く時刻は 10m−210m-2 分である。FIFOで乗るので、mm 回目のエレベータに乗る最後の客は、連続流として見ると累計 20m20m 人目であり、到着時刻は 20m2.2 \frac{20m}{2.2} 分である。全員の待ち時間が10分以上となるには、この最後の客でも (10m−2)−20m2.2≥10 (10m-2)-\frac{20m}{2.2}\ge 10 を満たせばよい。これより 0.909…m≥12 0.909\ldots m\ge 12 なので最小の整数は m=14m=14 である。したがって最初のエレベータは 10m=140 分後 10m=140\,\text{分後} に出発する。

    1. この場合、1サイクルは上昇3分、展望階停止2分、下降3分、1階停止3分で 11 分 11\,\text{分} である。容量は11分ごとに21人なので、1時間あたり容量は 2111×60=114.545… 人/時. \frac{21}{11}\times60=114.545\ldots\,\text{人/時}. したがって需要容量比は 132114.545…=24.221≃1.152. \frac{132}{114.545\ldots} =\frac{24.2}{21} \simeq 1.152.

    2. 60分までに出発するのは 11,22,33,44,5511,22,33,44,55 分の5回であり、輸送人数は 5×21=105 人 5\times21=105\,\text{人} である。一方、60分間の到着人数は 132 人 132\,\text{人} なので、待ち人数は 132−105=27 人 132-105=27\,\text{人} である。

    3. mm 回目の出発時刻は 11m11m 分、1階で扉が開く時刻は 11m−311m-3 分である。mm 回目の最後の乗客は累計 21m21m 人目であり、到着時刻は 21m/2.221m/2.2 分である。全員が10分以上待つ条件は (11m−3)−21m2.2≥10 (11m-3)-\frac{21m}{2.2}\ge10 である。すなわち 1.4545…m≥13 1.4545\ldots m\ge13 となるので、最小の整数は m=9m=9 である。出発時刻は 11m=99 分後 11m=99\,\text{分後} である。

    4. 1回あたりの定員は20人から21人に増えているが、1階停止時間が1分延びたため、サイクル長が10分から11分に増える。容量は 120120 人/時から約 114.5114.5 人/時へ下がり、需要容量比は 1.11.1 から約 1.1521.152 へ悪化する。少し定員を増やしても、停止時間が延びると総輸送力は低下し、待ち行列の増加が速くなる。

最終答

(1) 平均待ち時間 3.23.2 分。(2) 需要容量比 1.11.1、60分後の待ち人数 1212 人。(3) 140140 分後。(4) 需要容量比約 1.1521.152、60分後 2727 人、全員10分以上待つ最初の出発は 9999 分後。

年度別の解答例(11年度)

  • 2020年度(全7問)

    基礎数学 / 画像検索 / 重力走行と振動

  • 2019年度(全7問)

    基礎数学 / オペアンプ回路 / 流体の無次元数

  • 2018年度(全5問)

    二次曲線 / CとRPN / 色空間とDCT

  • 2017年度(全5問)

    回帰と正則化 / Cと正規表現 / 二色性反射

  • 2016年度(このページ・全5問)最終答まで公開

    ロジスティックモデル / 動的計画法 / 無線LANとフーリエ変換

  • 2015年度(全5問)最終答まで公開

    ガウス関数のフーリエ変換 / 乱数生成アルゴリズム / レイトレーシング

  • 2014年度(全5問)最終答まで公開

    積分とアステロイド / 水量計測と状態探索 / 画像処理

  • 2013年度(全5問)最終答まで公開

    QR分解と最小二乗 / 集合データ構造 / 信号処理

  • 2012年度(全5問)最終答まで公開

    線形最小二乗 / ソートの計算量 / 線形システムと相互情報量

  • 2011年度(全7問)最終答まで公開

    固有値と漸化式 / 命題論理 / 信号解析

  • 2010年度(全11問)最終答まで公開

    線形微分方程式 / 漸化式と高速計算 / z 変換とデジタルフィルタ

2016年度の出題内容

全5問

この5問の分野は線形代数1問・微分積分・解析1問・力学1問です。

大問分野主題解説の小見出し最終答
第1問—ロジスティックモデルロジスティック解の検算 / 平衡点の扱いあり
第2問—動的計画法貪欲法ではないあり
第3問微分積分・解析無線LANとフーリエ変換正規化の確認 / 切片と零点あり
第4問力学材料力学軸力の符号 / 単位換算あり
第5問線形代数エレベータ待ち行列待ち時間の定義 / 連続流としての扱いあり

同じテーマを出した他大学・他年度

6テーマ

前年度(2015年度)との違い

大問数
2015年度 5問 → 2016年度 5問
両年度に出たテーマ
フーリエ変換
2015年度の解答例

東京大学の公式過去問

公開年度範囲
2010〜2020(先端表現情報学コースとして公開されている専門科目)
公式の解答例
非公開(問題冊子のみ)
英語の扱い
TOEFL/TOEIC 等の扱いは募集要項で年度別に確認。公式過去問ページは専門科目PDFを公開。
院試hubの最終確認
2026-07-05
東京大学の公式過去問PDF