制御工学 自律システムの定義
微分方程式系が\(t\）を含まないとき、すなわち、ある微分方程式系が \begin{align}\dot{x} = f(x(t)) \hspace{5mm} f(x_e) = 0\end{align} のとき、このシステムを自律システムとい...
2023-03
制御工学 制御工学 安定性の定義
入力\(u(t)\)と状態ベクトル\(x(t)\)を用いて記述される次のシステムがあるとする。 \begin{align}y=f(x(t),u(t))\end{align} このシステムの内部安定性を調べるために\(u(t)\)を時間関数に...
MATLAB/simulink MATLABでベルヌーイ数を求める
これの続き。 N=10; B=zeros(1,N); B(1,1)=1; for i = 1:1:N B(1,i+1)=getBernoulliNumber(i, B); end B function y = getBernoulliNum...
C/C++/C# C++の正規表現を使ってみる
regexを使えばできる。 以下コード #include <iostream> #include <string> #include <regex> int main() { std::string str = "12345あいうえお"; ...
C/C++/C# C++でベルヌーイ数を求める
ベルヌーイ数を求めるには漸化式を解けばいい。 漸化式は \begin{align}B_0&=1\\B_n&=-\frac{1}{n+1} \sum_{k=0}^{n-1} {}_{n+1} C_{k} B_{k}\end{align} とな...
代数 冪乗の和公式
冪乗の和公式は次式で与えられる。 \begin{align}\sum_{i=1}^{n} i^k = \sum_{j=0}^{k} i^k \begin{pmatrix}k \\ j\end{pmatrix}B_j \frac{n^{k+1...
物理 理想気体の状態方程式
理想気体の状態方程式は圧力を\(P\)、体積を\(V\)、物質量を\(n\)、モル気体定数を\(R\)、熱力学温度を\(T\)とすると \begin{align}PV=nRT\end{align} で与えられる。
C/C++/C# 多倍長ライブラリのメモ
多倍長ライブラリmpirを含むデータをtupleにまとめvectorに格納するとうまくいかなくなる。
数学 小数の二乗とルートの外し方
小数\(0.5\)の二乗は\(0.25\)となるが、今回は\(\sqrt{0.25}\)を考える。 これは \begin{align}\sqrt{0.25}&=\sqrt{\frac{25}{100}}\\&=\sqrt{\frac{5^2...
代数 二次方程式の複素解と共有点の場所
二次方程式 \begin{align}y=ax^2+bx+c\end{align} について複素解になるのは \begin{align}b^2-4ac<0\end{align} のときである。このときの共有点の場所を調べる。 \(x=p+q...
物理 圧力の定義
単位面積\(S\)当たりに働く力\(F\)を圧力といい、圧力\(P\)は \begin{align}P=\frac{F}{S}\end{align} となる。単位はパスカル
数学 三角関数のマクロリーン展開を使ったバーゼル問題の解法
三角関数の級数展開を使ってバーゼル問題の値を求める。 \(\sin x\)と\(\frac{\sin x}{x}\)をマクロリーン展開する。 \begin{align}\sin x &=x - \frac{x^3}{3!} + \cdots...
数学 ディリクレ核とは
次の関数で得られる数列 \begin{align}D_n(x)=1+2\sum_{k=1}^{n} \cos (kx) = \frac{\sin \left( \frac{2n+1}{2}\right )x}{\sin \frac{x}{2...