qd-suriken.com

「基礎から学べる機械力学」（伊藤）を参考に clear; close all; clc; e = 1; wn =1 ; t = linspace(0,10,30); x = e*wn/2.*t.*sin(wn*t); y = -e*wn/2.*t.*cos(wn*t)+e/2*sin(wn*t); plot(t,x); hold on; plot(t,y); legend('x(t)','y(t)') xlabel('time …

clear; close all; clc; E = 206e9; d = 0.1; I = pi/64*d^4; m = 30; l = 2; a = l/2; b = l/2; k = 6*E*I*l/(a*b*(l^2-a^2-b^2)); nc = 60/(2*pi)*sqrt(k/m); fprintf('危険速度: %.0f [rpm]\n',nc); 実行結果 危険速度: 4294 [rpm]

等温壁 Illingworth_Stewartson_transformation.m clear; close all; clc; eta = linspace(0, 20, 50); fprimeprime_tmps = linspace(0.45,0.48,1000); errors = []; for idx=1:length(fprimeprime_tmps) fprimeprime_tmp = fprimeprime_tmps(idx); [eta2, x…