スマホで数値計算（山崎純一）

% 2021-08-13 Ver 1.

% フーリエ変換するデーターを読込み、プロット
F0 = importdata('FFTTest-01.txt');
plot (F0)

% F0の1から300番目を取り出しプロット
F0s = F0(1:300,:);
plot (F0s)

% フーリエ変換し、プロット
t = 0.01; % サンプリング周期
l = 1024; % データー点数
F1 = fft(F0); % フーリエ変換
F2 = abs(F1/l);
F3 = F2(1:l/2+1);
ts = 1/t;
tt = ts*(0:(l/2))/l;
plot (tt,F3)

% 横軸を拡大
xlim([0,20])

% 2021-8-13 Ver 1.00

% 連立一次方程式 1
% x + 8y = 2
% 2x + 9y = -1
%
[1 8; 2 9] \ [2;-1]

% 連立一次方程式 2
% 2x + 8y = 2
% 2x + 3y = -1
%
[2 8; 2 3] \ [2;-1]

% 連立一次方程式 3
% 1x + 1y = 8
% 2x + 4y = 26
%
[1 1; 2 4] \ [8;26]

% 2021-8-13 Ver 1.00

% 式の三次元プロット
% x -5 5 y -5 5
fsurf(@(x,y) sin(x)-cos(y))

% 式の三次元プロット
% x -10 0 y -5 15
fsurf(@(x,y) sin(x)-cos(y),[-10 0 -5 15])

% 2021-8-13 Ver 1.00

% 関数を-5 5の区間でプロット
fplot(@(x) sin(x))

% 関数を0 5の区間でプロット
fplot(@(x) 3*sin(x)^2+cos(x),[0 5])

% 関数を２ポイントのラインでプロット
fplot(@(x) sin(x+pi/5),'Linewidth',2)

% 関数を２ポイントのラインでプロット
fplot(@(x) sin(x+pi/5),'Linewidth',2);
hold on
fplot(@(x) sin(x^2+pi/5),'--or');
fplot(@(x) sin(x^3+pi/3),'-.*c');
hold off

% 2021-8-13 Ver 1.00

% 円周率
fplot(@(x) sin(x))

% 円周率
pi

% 円周率
disp('digits(30)')
digits(30)
vpa(pi)

% 虚数
i
j
a = i
b = a^2

% ネイピア数
exp(1)

% オイラーの恒等式
% e^(π×i) = -1
exp(pi * i)