http://10.8.1.3/class/ 2026-05-09T22:31:42+00:00 http://10.8.1.3/class/ http://10.8.1.3/class/lib/exe/fetch.php?media=wiki:dokuwiki.svg text/html 2026-04-15T23:52:06+00:00 Anonymous (anonymous@undisclosed.example.com) http://10.8.1.3/class/doku.php?id=numericalipa:ans1&rev=1776297126&do=diff 【数値情報処理a】第1回の答え 課題 以下の問に答えなさい。 (1) RStudio の自分のユーザー名を16進数の整数とみなして10進数に変換しなさい。 例えば、ユーザー名が「c123456」だったら「c123456」を16進数の整数とみなす。$16^{7}$\begin{eqnarray} a_{7}&=&(\mathrm{c})_{16}=12\\ a_{6}&=&(\mathrm{1})_{16}=1\\ a_{5}&=&(\mathrm{2})_{16}=2\\ a_{4}&=&(\mathrm{3})_{16}=3\\ a_{3}&=&(\mathrm{4})_{16}=4\\ a_{2}&=&(\mathrm{5})_{16}=5\\ a_{1}&=&(\mathrm{6})_{16}=6 \end{eqnarray}\begin{eqnarray} &&a_{7}\times 16^{6}+a_{6}\times 16^{5}+a_{5}\times 16^{4}+a_{4}\times 16^{3}+a_{3}\times 16^{2}+a_{2}\times… text/html 2026-04-22T15:49:28+00:00 Anonymous (anonymous@undisclosed.example.com) http://10.8.1.3/class/doku.php?id=numericalipa:ans2&rev=1776872968&do=diff 【数値情報処理a】第2回の答え 課題 関数 $f(x)$ において $f(x)=0$ の解をニュートン法で全て求めなさい。 解の収束は $\delta=10^{-10}$ の精度とする。 それぞれの解と計算に用いた $x^{(0)}$ の初期値を答えること。 解答例 $f(x)$$f(x)=0$$f(x)$$f^{\prime}(x)$$\delta=10^{-10}$$f(x)=-x^{4}-x^{3}+10x^{2}+x-2$\[ f^{\prime}(x)=-4x^{3}-3x^{2}+20x+1 \]$f(x)=-x^{4}-2x^{3}+10x^{2}+x-2$\[ f^{\prime}(x)=-4x^{3}-6x^{2}+20x+1 \]$f(x)=-x^{4}-3x^{3}+10x^{2}+x-2$\[ f^{\prime}(x)=-4x^{3}-9x^{2}+20x+1 \]$f(x)=-x^{4}-4x^{3}+10x^{2}+x-2$\[ f^{\prime}(x)=-4x^{3}-12x^{2}+20x+1 \]$f(x)=-x^{4}-5x^{3}+1… text/html 2026-04-29T15:44:55+00:00 Anonymous (anonymous@undisclosed.example.com) http://10.8.1.3/class/doku.php?id=numericalipa:ans3&rev=1777477495&do=diff 【数値情報処理a】第3回の答え 課題 $y=f(x)$ と $x$ 軸ではさまれた領域の面積をシンプソンの公式で数値積分するプログラムを作って提出しなさい。 * 求める面積は下図の赤い部分である。 * $y$$f(x)=-x^{4}-x^{3}+10x^{2}+x-2$$f(x)=0$$f(x)<0$$[a,b]$$f(x)=0$$f(x)=-x^{4}-x^{3}+10x^{2}+x-2$$f(x)=-x^{4}-2x^{3}+10x^{2}+x-2$$f(x)=-x^{4}-3x^{3}+10x^{2}+x-2$$f(x)=-x^{4}-4x^{3}+10x^{2}+x-2$$f(x)=-x^{4}-5x^{3}+10x^{2}+x-2$$f(x)=-x^{4}-6x^{3}+10x^{2}+x-2$… text/html 2026-04-08T23:38:03+00:00 Anonymous (anonymous@undisclosed.example.com) http://10.8.1.3/class/doku.php?id=numericalipa:lecture1&rev=1775691483&do=diff 【数値情報処理a】第1回 数の表現と誤差 数値情報処理a について 数値計算とは、人間が解析的に解くのが難しかったり、解くのに時間がかかりすぎる問題をコンピュータ（計算機）に解かせることである。 コンピュータの発達に伴い、数値計算の手法も大きく発達してきた。 AI でも数値計算の手法が至るところで使われている。 本講義では R を使って数値計算の基本的な手法を学ぶ$10!=1\times 2\times 3\times \cdots \times 9\times 10$$n!$$0$$9$$0$$9$$\mathrm{a}$$\mathrm{f}$$51966$\[ 51966=(1100\ 1010\ 1111\ 1110)_{2}=(\mathrm{cafe})_{16} \]\(m\)$(a_{m}a_{m-1}\cdots a_{3}a_{2}a_{1})_{N}$\begin{eqnarray} &&\sum_{i=1}^{m}a_{i}\times N^{i-1}\\ &=&a_{m}\times N^{m-1}+a_{m-1}\times N^{m-2}+\cdots… text/html 2026-04-16T00:13:40+00:00 Anonymous (anonymous@undisclosed.example.com) http://10.8.1.3/class/doku.php?id=numericalipa:lecture2&rev=1776298420&do=diff 【数値情報処理a】第2回 代数方程式 前回の課題 【数値情報処理a】第1回の答え 数値計算で生じる誤差 * 丸め誤差 * 打ち切り誤差 * 桁落ち * 情報落ち ---------- 丸め誤差 浮動小数点数の仮数部の上限の桁数を超えたときに丸めることによって生じる誤差である。\[ S_{n}=\sum_{i=1}^{n}(0.5)^{i-1} \]$n\rightarrow\infty$\[ \lim_{n\rightarrow\infty}S_{n}=\lim_{n\rightarrow\infty}\frac{1-(0.5)^{n}}{1-0.5}=\frac{1}{0.5}=2 \]$S_{n}$$n=10$$2$$2-1.998047=0.001953$$n$$ax^{2}+bx+c=0$\[ x_{\pm}=\frac{-b\pm\sqrt{b^{2}-4ac}}{2a} \]$b$$a,c$$a=b=1, c=10^{-8}$$x_{+}$\[ x_{+}= -1.00000001000\cdots\times 10^{-8} \]$x_{+}… text/html 2026-05-01T03:54:40+00:00 Anonymous (anonymous@undisclosed.example.com) http://10.8.1.3/class/doku.php?id=numericalipa:lecture3&rev=1777607680&do=diff 【数値情報処理a】第3回 数値積分 前回の課題 【数値情報処理a】第2回の答え 積分 $x$ の区間 $[a,b]$ で関数 $y=f(x)$ と $x$ 軸ではさまれた領域の面積 $S$ は積分（定積分）で求められる。 \[ S=\int_{a}^{b}f(x)dx=[F(x)]_{a}^{b}=F(b)-F(a) \] $x$ を積分変数、 $f(x)$ を被積分関数、 $F(x)$ を $f(x)$ の原始関数$F(x)$$x$$f(x)$\[ \frac{dF(x)}{dx}=F^{\prime}(x)=f(x) \]$f(x)$$F(x)$$f(x)$$F(x)$$[a,b]$\[ f(x)=x^{n}\quad(n\neq -1) \]$f(x)$\[ F(x)=\frac{1}{n+1}x^{n+1} \]\[ S=\int_{a}^{b}x^{n}dx=\left[\frac{1}{n+1}x^{n+1}\right]_{a}^{b}=\frac{b^{n+1}-a^{n+1}}{n+1} \]$n=1$\[ f(x)=x \]\[ S=\int_{a}^{b}… text/html 2026-04-30T02:37:56+00:00 Anonymous (anonymous@undisclosed.example.com) http://10.8.1.3/class/doku.php?id=numericalipa:lecture4&rev=1777516676&do=diff 【数値情報処理a】第4回 乱数とモンテカルロ法 前回の課題 【数値情報処理a】第3回の答え 乱数 乱数はバラバラの数であり、理想的にはどの数が出るか予想ができず規則性も再現性も持たない。 例えばサイコロ1個を転がして出る目は 1,2,3,4,5,6 の数字のどれがでるか予想できない乱数である。$0<x<1$$<x<$$\mathrm{B}(m,p)$$\mathrm{N}(0,1)$$\mathrm{N}(m,s^{2})$$\Lambda(0,1)$$\Lambda(m,s^{2})$$\mathrm{t}(k)$$\chi^{2}(k)$$\mathrm{F}(k_{1},k_{2})$$0<x<1$$0<x<1$$x$$0<x<1$$1<x<7$$1<x<2$$1$$2\le x<3$$2$$3\le x<4$$3$$4\le x<5$$4$$5\le x<6$$5$$6\le x<7$$6$$2^{19937}−1$$\simeq4.3\times 10^{6001}$$\mathrm{N}(0,1)$$x=0$$x=0$$\mathrm{N}(0,1)$$\mathr… text/html 2026-04-05T15:32:21+00:00 Anonymous (anonymous@undisclosed.example.com) http://10.8.1.3/class/doku.php?id=numericalipa:top&rev=1775403141&do=diff 数値情報処理a * 【数値情報処理a】第1回 数の表現と誤差 * 【数値情報処理a】第2回 代数方程式 * 【数値情報処理a】第3回 数値積分 * 【数値情報処理a】第4回 乱数とモンテカルロ法 * lecture5 * lecture6 * 学内 RStudio サーバ1（学内限定） * 学内 RStudio サーバ2（学内限定） 質問がある場合は s4 かメールで連絡してください。 * s4 2026年度 数値情報処理a * e-mail:[email]…