パチンコの確立計算機

パチンコの確立計算機なるものがあるみたい どうやって計算してるんだろう・・・

伝達関数とパルス伝達関数の相互変換について

\(s\)平面から\(z\)平面への変換式は \begin{align}\label{S-T transform}z=e^{sT}\end{align} で与えられる．\(z\)平面上の点および\(s\)平面上の点を \begin{alig

計算するのが大変な積分に用いる置換積分で何が起きるのか

計算するのが大変な積分に用いる置換積分で何が起きるのか 下の積分の例で見る \begin{align}\displaystyle \int x(2-x)^4 dx\end{align} \(t=2-x\)とおくと \begin{align}

MATLABを使って複素数を写像してみる

複素関数を使えば複素数を写像できる。 ディジタル制御では \begin{align}s=e^{sT}\end{align} を使うので\(T=1\)として写像してみる 例えば下のプログラムの例では虚軸が円に写される。 x=0; y=-5:0

MATLABで台形近似

台形近似で積分を計算してみる 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型非線形モデルの非線形ブロックの取り扱い

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でオイラーの定数を計算する

詳しくはここ 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制御とは

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枚の中から指定の五枚を引く確率

40枚の中から指定の五枚を引く確率は \begin{align}\frac{1}{{}_{40} C_{5}={1}{658008}\end{align} となる

MATLABでウッダル数を求める

ウッダル数は \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

MATLABで作ったカレン素数を探すプログラムを改造してみる

カレン数は \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

MATLABでカレン数とカレン素数をさがす

カレン数は \begin{align}n \times 2^i + 1\end{align} で表される。 カレン数のうち素数のものをカレン素数という。 今回はMATLABでカレン素数を探す。 以下ソース n=10; count=1; fo

RNNとLSTM

参考 RNNは入出力を等しく学習→長期的な依存性の学習が苦手

NumbersとLSTM

NUMBERSには横滑り現象なるものがあるらしくLSTMで学習して当てる試みがほそぼそとあるらしい Qiitaだとこれとか Github 機械学習に興味あるのでやってみようと思う

ポケモン図鑑を和訳してみる ニャオハ

和訳してみるGrassCatPokemonくさねこポケモンItsfluffyfurissimilarincompositiontoplants.ふわふわした毛皮は植物と成分が似ている。fluffyふわふわしたfur毛皮ThisPokemon

電力と同期化力の計算

電力と同期化力の関係は\begin{align}P_e=\frac{\partialP_e}{\partial\delta}=\frac{E_sE_r}{X}\cos\delta\end{align}

風力発電の出力

風車の回転断面積を\(A\)、風速を\(V\)、空気の密度を\(\rho\)、ロータの係数を\(C\)とすると風力発電の出力\(P\)は\begin{align}P=\frac{1}{2}C\rhoV^3A\end{align}

matlabのtanとtandの違いを簡単に見てみる

matlabのtanとtandの違いを簡単に見てみる。まずは89から90どの範囲で重ねて比較。ほとんど同じ値が得られている。2つの方法の誤差。通常誤差は0であるが微妙に生じている。90度に近くなるにつれて大きくなるようだ。以下コードtd=l

ローレンツ力とは

電場\(\boldsymbolE\)、磁場\(\boldsymbolB\)中を移動する電荷\(\boldsymbolq\)の荷電粒子に加わる力は荷電粒子の位置を\(\boldsymbolx\)、速度を\(\boldsymbolv\)とすれば

matlabでn以下の最大の素数を求める

100以下であればprevprime(100)を使えばいい

レギュレータ問題における操作量の表現

PID制御則が\begin{align}u(t)=K_{p}e(t)+K_{i}\inte(\tau)d\tau+K_{d}\frac{de(t)}{dt}\end{align}で与えられているとき、この制御則に対するレギュレータ問題を考え

PID制御とレギュレータ問題

PID制御は\begin{align}u(t)=K_{p}e(t)+K_{i}\inte(t)d\tau+K_{d}\frac{de(t)}{dt}\end{align}のような問題を言い、\(requiv0\)をレギュレータ問題という。

Lyapunov方程式を証明する

システムの安定性を調べるにはLyapunov方程式\begin{align}PA+A^{T}P=-Q\end{align}を調べればいい。\(P\)は\(A\)の固有値の実部が負であれば\begin{align}P=\int_0^\inft

急にMATLABのロゴ画像が欲しくなったとき

コマンドラインにlogoと打つとが出力される。

初等数学の基本定理について

ある正の整数\(a\)は\(a\)以下の素数\(p_1 p_2 \cdots p_{n-1} p_n\)の積で表すことができる。つまり\begin{align}a=p_1 p_2 \cdots p_{n-1} p_n\end{align}が

相加相乗平均

相加相乗平均とは\begin{align}\frac{a+b}{2} \geq \sqrt{ab}\end{align}のことを言う。

同期発電機の励磁方式

同期発電機には界磁巻線に励磁電流を与えるための励磁装置が必要である。励磁方式は次の3つに大別できる。・直流励磁機方式・交流励磁機方式・静止系励磁方式

書籍版電験王 電験2種二次試験 過去問徹底解説 令和5年度版を買ってみた

電験二種の二次試験対策に書籍版電験王 電験2種二次試験 過去問徹底解説 令和5年度版を買ってみた。良かった点・解説が分かりやすい・図が見やすい悪かった点・本が大きい使いやすい参考書だと思います。

collectを使って係数をまとめる

matlabでシンボリック演算を行った際の係数をまとめるときはcollectを使う

n個の数字からm個を取り出すときの組み合わせの数

\begin{align}{}_{n}C_{m}==\frac{n(n-1)(n-2) \cdots (n-m+2)(n-m+1)}{m(m-1)(m-2) \cdots 1}=\frac{n!}{m!(n-m)!}\end{align}

電流と電流密度

導体に電場\(E\)が与えられた時の電流密度は\begin{align}\boldsymbol{J}=\sigma \boldsymbol{E} \mathrn{}\end{align}電流との関係は\begin{align}I=\bold

MATLABで遊ぶ例が載ってるInterface買ってみた

Interface買ってみた。CD付きでMATLABでできる例があるらしい。今後試してみる。

電験二種二次試験の参考書

これを買ってみた

ESP32用の通信基板を作った

CANとI2Cで簡単に通信できるシールドを作った。コミケで頒布するかも

電束の定義

1の電荷から1の仮想の線が出るとする。この線を電束という。

It means a lotの意味

It means a lot は「大きな意味がある」だが、「とても大切にする」感情に使える。例It means a lot to me.(とても大切にする、とても嬉しい)It means a lot to you.(とても大切にしている)な

LaTeXでいい感じのPDF内リンクを作る

くわしくはプリアンプルに以下を追加すれば動く。\usepackage{hyperref}\usepackage{pxjahyper}\hypersetup{setpagesize=false, bookmarksnumbered=true,

ESP32のGPIOについて

GPIO34以降は入力のみ

MATLABで等比数列を実装する

MATLABで等比数列は次のように実装すればいい以下コードa=2;r=4;m=10;list=ones(1,m).*a;for n=2:1:m list(1,n)=a*r^(n-1);end

TeXWorksを使うときのbiberの設定

これを追加して、pLaTeX→biber→pLaTeXでうまくいく

メモ：複数のマイコンがバス上に存在するときのCAN通信

CANは衝突を起こす通信方式なので衝突回避が必要。マスターから送信指示を出すのが簡単

階差数列の定義

ある数列\(\{a_n\}\)について、隣合う2つの項の差\begin{align}b_n=a_{n+1}-a_n\end{align}で作られる数列\(\{b_n\}\)を\(\{a_n\}\)の階差数列という。

メモ：CAN割り込み中のシリアル通信

CANの割り込み処理中のシリアル通信をやりすぎると落ちる