MATLAB/simulink

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

ウッダル数は \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^n + 1\end{align} で表される。 カレン数のうち素数のものをカレン素数という。 今回はMATLABでカレン素数を探す。 以下ソース n=10; count=1; fo...

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

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

Interface買ってみた。CD付きで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で分散を逐次計算しようとしてうまくいかなかった。 分散の逐次計算は \begin{align}\sigma_{n+1}^2=\dfrac{n(\sigma_n^2+\mu_n^2)+x_{n+1}^2}{n+1}-\mu_{n+...

平均は \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); ...

これの続き 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...

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

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

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

これのMATLAB版 n = 15; count = 0; for i = 1:n disp(); count = count + isprime(i); end disp(); function p = isprime(n) k = fa...

polynomialに生成した多項式を放り込んでsolveで解を求める。coefficientsには高い順に係数を入れればいい。今の例だと\(x^2+5x+6=0\)を解く。 coefficients = ; syms x; polynom...