計算するのが大変な積分に用いる置換積分で何が起きるのか 下の積分の例で見る \begin{align}\displaystyle \int x(2-x)^4 dx\end{align} \(t=2-x\)とおくと \begin{align}
ルジャンドル予想とはある任意の自然数\(n\)について \begin{align}n^2 \leq p_1 \leq \cdots \leq p_n\leq(n+1)^2\end{align} となるような素数\(p_1 […]
マクスウェルの方程式は \begin{align}\mathrm{div} \boldsymbol{B} (t,\boldsymbol{x}) &= 0 \[1.5ex]\mathrm{rot} \boldsy […]
これの続き。 オイラーの公式 \begin{align}e^{i \theta} = cos \theta + i \sin \theta\end{align} から加法定理を計算することができる。いま\(\theta_ […]
\( y=f(x)^{g(x)} \)などの形をとる関数の微分を行う場合対数微分法を使うと簡単になることがある。 この形の対数微分を考える。両辺の対数をとり \begin{align}\log{y(x)}&=g( […]
炭素原子から陽子の重さを計算する。1mol中に存在する炭素原子数をアボガドロ定数個とし、炭素12の陽子と中性子数は12、電子の重さを無視し、陽子と中性子が同じ重さであるとすれば \begin{align}m_{p} \a […]
遠藤がカイジに吹っ掛けた10分3割複利がどのくらいやばいか計算してみる。 繰り返し回数を\(n\)、元金を\(a\)とすればこの計算は \begin{align}y = 1.3^n a \end{align} 指数のグラ […]
これの続き。 刻み量が少なく正規分布に見えないので刻み量を増やす。ヒストグラムを生成する部分 のbinsを変えればいい 増やすとこうなる。なぜか段々になっている点が気になるが正規分布に近づいた。
サイコロのある面が出る確率はどの目でも一様であると考えれば \begin{align}P(X)=\frac{1}{6}\end{align} となる。いまサイコロを \(N\) 回振り、その平均を求めることを考える。 例 […]
転置行列は \begin{align}A=\begin{pmatrix} a & b \ c & d\end{pmatrix}\end{align} の時 \begin{align} ^{t} A =\ […]
三層のニューラルネットワークをPythonで実装する。 \(X\)を入力 、 \(W\)を重み 、 \(B\)をバイアスとすれば各層の計算は行列を使って \begin{align}A=XW+B\end{align} と計 […]
numpyを使って行列を定義するには とすればいい 出来てるか確認するには 大きさを確認するには とするといい。
シグモイド関数の一つをPythonを使って描画する 関数は \begin{align}y=\frac{1}{1+e^{-x}}\end{align}
あるシステムを表す線形非同次微分方程式 \begin{align}\dot{x} (t) = Ax(t) +Bu(t) \hspace{5mm} x(t_0)=x_0\end{align} の解を求める。 \begin{ […]
AmazonにRaspberryPi用のカメラがあったので買いました。カメラはこれ↓ 組み立てはプラねじでアクリルを止めるだけ。組み立て手順の説明書はないが何とかなる。組み立てると画像のようになった。 付属品に長さの異な […]
これの続き。RC直列回路の回路方程式は \begin{align}E=Ri(t)+\frac{1}{C} \int i(t) dt\end{align} これを解けば \begin{align}q(t)&=CE( […]
電流の関係式 \begin{align}i=\frac{dq}{dt}= C \frac{dv}{dt} \end{align} より出力電圧は \begin{align}\frac{dv}{dt} = \frac{E- […]
テント写像は \begin{align}x_{n+1}=\left 1 – 2 x_{n} \right \end{align}
これの続き。zpk2tfを使うと指定した極、零点、ゲインを持つ伝達関数を簡単に設計できるようになる。極、零点、ゲインは次のように指定する。 結果 ソース全体
Pythonで古典制御と現代制御の双方の視点からばねマスダンパ系を解析する
これの続き。Pythonで同じ解析をした。
「ブログリーダー」を活用して、しろねこさんをフォローしませんか?
計算するのが大変な積分に用いる置換積分で何が起きるのか 下の積分の例で見る \begin{align}\displaystyle \int x(2-x)^4 dx\end{align} \(t=2-x\)とおくと \begin{align}
複素関数を使えば複素数を写像できる。 ディジタル制御では \begin{align}s=e^{sT}\end{align} を使うので\(T=1\)として写像してみる 例えば下のプログラムの例では虚軸が円に写される。 x=0; y=-5:0
台形近似で積分を計算してみる Nが刻み数 minが下限、maxが上限 funcが被積分関数 N=100; min=0; max=1; t=linspace(min,max,N); dt=t(2)-t(1); S=zeros(size(t))
博士とったのでブログも再開します!
Hammerstein型非線形モデルの非線形ブロックによるゲインを\(\alpha\)とすると \begin{align}B(q^{-1}) &=\alpha b_{1} q^{-1} +\alpha b_{2} q^{-2} +
\(f(x)=1\)とする。この関数を\(a\)から\(b\)まで複数回積分すると \begin{align}\int_a^b 1 dx=a-b\end{align} \begin{align}\int_a^b \int_a^b 1 dx
\begin{align}\int_{-a}^{a} f(x) dx = 2 \int_{0}^{a} f(x) dx\end{align}
奇関数の定積分には \begin{align}\int_{-a}^af(x)dx=0\end{align} が成り立つ
詳しくはここ MATLABの行列演算を使うと楽 N=10000; n=1:1:N; result=sum(1./n-log(1+1./n))
オイラーの定数とは \begin{align}\gamma=\lim_{n \to \infty} (1 + \frac{1}{2} + \cdots + \frac{1}{n} - \log n)\end{align} の極限値のことであ
クロネッカー積には次の関係が成り立つ。 \begin{align}x_1 \otimes (y_1+y_2)&=x_1 \otimes y_1 + x_1 \otimes y_2 \\(x_1 + x_2 ) \otimes y_1
\(X,\mathcal{O}\)を位相空間とする。 \begin{align}{}^{\forall} x_1,x_2 \in X (x_1 \neq x_2) \hspace{2mm} {}^{\exists} \mathcal{O}_
PID制御とは比例・積分・微分の3つを組み合わせて行う制御方式である。 PID制御は次のように与えられる。 \begin{align}u(t)=K_P e(t) + K_{I} \int_0^{t} e(\tau) d\tau + K_D
正五角形の1辺の長さを1とすると正五角形の対角線の長さ\(a\)は余弦定理より \begin{align}a^2&=1^2 + 1^2 - 2 \times 1 \times 1 \times \cos 108\\&= 2
40枚の中から指定の五枚を引く確率は \begin{align}\frac{1}{{}_{40} C_{5}={1}{658008}\end{align} となる
ウッダル数は \begin{align}n \times 2^n -1\end{align} の形をしている数である。 MATLABでは次のように計算できる。 n=10; count=1; p=2; for i=1:1:n K(i)=i*p
カレン数は \begin{align}n \times 2^n + 1\end{align} であるが \begin{align}n \times p^n + 1\end{align} を考える。 n=10; count=1; p=3; f
カレン数は \begin{align}n \times 2^i + 1\end{align} で表される。 カレン数のうち素数のものをカレン素数という。 今回はMATLABでカレン素数を探す。 以下ソース n=10; count=1; fo
参考 RNNは入出力を等しく学習→長期的な依存性の学習が苦手
NUMBERSには横滑り現象なるものがあるらしくLSTMで学習して当てる試みがほそぼそとあるらしい Qiitaだとこれとか Github 機械学習に興味あるのでやってみようと思う
ベクトル関数がスカラー関数のときの微分は\begin{align}\frac{d \boldsymbol{K}}{dt}=0\end{align}となる。
ベクトル\(\boldsymbol{A}(t),\boldsymbol{B}(t)\)について\begin{align}\frac{d \boldsymbol{A}(t)}{dt}=\boldsymbol{B}(t)\end{align}の
次の式を連続の式という。\begin{align}\frac{\partial \rho}{dt} + \mathrm{div} (\rho \boldsymbol{v})=0\end{align}
ベクトルの微分は各成分ごとに微分したものと等しい。即ち\begin{align} \frac{d \boldsymbol{A}(t)}{dt}=\frac{dA_x(t)}{dt} \boldsymbol{i}+\frac{dA_y(t)}
静電場\(\boldsymbol{E}(\boldsymbol{r})\)について、\begin{align}\boldsymbol{E}(\boldsymbol{r})=- \nabla \phi (\boldsymbol{r})\end
ニュートンの運動方程式\begin{align}m \frac{d^2x(t)}{dt^2} =F\end{align}および自由落下を行っている物体に掛かる力\begin{align}F=-mg\end{align}より\begin{al
ベクトル関数の微分\(A(t)\)の微分係数は\begin{align}\frac{dA(t)}{dt}=\lim_{\delta t \to 0} \frac{A(t + \delta t)-A(t)}{\delta t}\end{ali
畳み込み積分のラプラス変換は次のようになる。\begin{align} \mathcal{L}&=\int_0^{\infty}e^{-st}\int_0^tf(u)g(t-u)dudt \\&=\int_0^{\infty
定義に従い計算すれば良い。\begin{align}\mathcal{L} & =\int_0^\infty e^{-st} (a f(t) + b g(t)) dt \\& =a \lim_{p \to \infty}
ある実数\(t\)によってベクトル\(A\)が定まる時、これをベクトル関数といい\(A(t)\)と書く。\(A(t)\)の変数が\(A_x,A_y,A_z\)であれば\begin{align}A(t)=A_x(t) \boldsymbol{
合同数の定義\begin{align} \begin{cases}X^2+Y^2=Z^2\\\displaystyle \frac{XY}{2}=n\end{cases}\end{align}楕円曲線の関係を求める。合同数の定義を平方完成す
合同数の定義\begin{align} \begin{cases}X^2+Y^2=Z^2\\\displaystyle \frac{XY}{2}=n\end{cases}\end{align}楕円曲線の関係を求める。合同数の定義を平方完成す
回路に電気エネルギーを供給する素子を電源という。外部にどんな負荷を接続しても一定の電流を出力する電源を電圧源という。理想電圧源の内部抵抗は零である。電圧源に接続された抵抗を小さくすることを考える。オームの法則より、\begin{align}
区間\((0,\infty]\)で定義された関数\(f(t)\)について次の無限積分\begin{align}\lim_{T \to \infty} \int^{T}_{0} e^{-st} f dt = \int_0^\infty e^{
合同数の定義\begin{align} \begin{cases}X^2+Y^2=Z^2\\\frac{XY}{2}=n\end{cases}\end{align}楕円曲線の関係を求める。合同数の定義を平方完成すれば\begin{align
\(n\)が合同数であるとは\begin{align}\begin{cases}x^2+y^2=z^2\\\frac{xy}{2}=n\end{cases}\end{align}となる有理数\(x,y,z\)が存在することである。
3辺の辺の長さが有理数であるような直角三角形の面積を合同数という。例:底辺を\(\frac{3}{2}\)、高さを\(\frac{20}{3}\)とすると斜辺は\begin{align}c&=\sqrt{\left ( \frac{
もし~だったらどうなるか
\(5 \times 5\)が\(25\)であることから\(5\)
集合の演算において、次の分配率が成り立つ。\begin{align}A \cup (B \cap C)\end{align}証明\begin{align}x \in A \cup (B \cap C) & \Leftrightarr