MATLAB/simulink

MATLAB/simulink

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

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

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))...
MATLAB/simulink

MATLABでオイラーの定数を計算する

詳しくはここ MATLABの行列演算を使うと楽 N=10000; n=1:1:N; result=sum(1./n-log(1+1./n))
MATLAB/simulink

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/simulink

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/simulink

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

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

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

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

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

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

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

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

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

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

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

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

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
MATLAB/simulink

三角関数の近似式 2

十分小さい正の角度\(\theta\)について、\(\cos \theta\)は\(\tan \theta \)を用いて \begin{align}\cos \theta \approx 1 - \frac{\tan^2 \theta}{2...
MATLAB/simulink

三角関数の近似式

十分小さい正の角度\(\theta\)について、\(\tan \theta\)は\(\cos \theta \)を用いて \begin{align}\tan \theta \approx \sqrt{2(1- \cos\theta)} \e...
MATLAB/simulink

MATLABで分散を逐次計算しようとしてうまくいかなかった話

MATLABで分散を逐次計算しようとしてうまくいかなかった。 分散の逐次計算は \begin{align}\sigma_{n+1}^2=\dfrac{n(\sigma_n^2+\mu_n^2)+x_{n+1}^2}{n+1}-\mu_{n+...
MATLAB/simulink

MATLABで平均を逐次計算する

平均は \begin{align}\mu_{n+1} = \frac{1}{n+1} (n \mu_n + x_{n+1})\end{align} で逐次計算できる。 以下コード N=10; x=1:1:N; mu=zeros(1,N); ...
MATLAB/simulink

MATLABでベルヌーイ数を求める

これの続き N=6; B=zeros(1,N); B(1,1) = 1; disp(B(1,1)); for i = 2:1:N B(1,i) = getBernoulliNumber(i, B); disp(B(1,i)); end fu...
MATLAB/simulink

MATLABで素数の逆数和を計算する

素数の逆数和は \begin{align}P=\frac{1}{2} + \frac{1}{3}+ \cdots = \infty\end{align} となる。これを計算する。 以下コード N=100; P=zeros(1,N); for...
MATLAB/simulink

MATLABで素因数分解をする

これの続き num = input('数字を入力してください: '); fprintf('入力された数字 -> %d\n', num); PrimeFactorization(num); function PrimeFactorizatio...
MATLAB/simulink

MATLABのisprimeの速度を比較してみた

以下のソースコードで速度を比較。 n=15のとき 自作関数:0.0025 isprime:0.0016 n=150のとき 自作関数:0.0049 isprime:0.0064 なお自作のmyisprimeで計算できるのはn=170程度まで ...