プログラミング

MATLABでテイラー展開してグラフ化するプログラムを書いた。以下コード close all f = @(x) cos(x); a = 0; n = 15; x_range = ; =plotTaylorSeries(f, a, n, x_...

※本抽選は厳正に行われています。（+90kg固定） % ステップ1: 文字列入力 segments = cell(1, 6); segments{1} = '+50kg'; segments{2} = '+60kg'; segments{3...

ChatGPTにネルダーミード法を使った関数の最適解を求めてもらった あってるかは後日確認するつもり % 最小化する関数 func = @(x) (x(1) - 3)^2 + (x(2) - 2)^2; % 初期点 x0 = ; % 収束許...

マンデルブロ集合を書くだけ % パラメータ設定 maxIter = 5000; % 最大反復回数 xlim = ; % x範囲 ylim = ; % y範囲 resolution = 1000; % 解像度 % 複素数平面のメッシュグリッド...

matlabで振り子を動かしてみる 運動方程式などの細かい話は次回 % パラメータ設定 g = 9.81; % 重力加速度 (m/s^2) L = 1.0; % 振り子の長さ (m) theta0 = 0.1; % 初期角度 (ラジアン) ...

複素関数を使えば複素数を写像できる。 ディジタル制御では \begin{align}z=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))...

詳しくはここ 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...