seinzumtode https://blog.hatena.ne.jp/seinzumtode/ how to code something https://seinzumtode.hatenadiary.jp/ clear; close all; zetas = [1,0.5,0.2,0.1,0.01]; w_wn = linspace(0,2,100); Mds=[]; figure(); hold on; for idx=1:length(zetas) zeta = zetas(idx); Md = 1./sqrt((1-w_wn.^2).^2+4*zeta^2.*w_wn.^2); plot(w_wn, Md,'DisplayName',sprintf('zeta=%.2g',zeta)); end legend(); xlabel('加振振動数比w_wn'); ylabel('変位の応答倍率M… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fseinzumtode.hatenadiary.jp%2Fentry%2F2023%2F05%2F01%2F141612" title="加振周波数ωと変位の応答倍率Md/位相角Φの関係 - how to code something" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://cdn-ak.f.st-hatena.com/images/fotolife/s/seinzumtode/20230501/20230501141523.png Hatena Blog https://hatena.blog 2023-05-01 14:16:12 rich https://seinzumtode.hatenadiary.jp/entry/2023/05/01/141612 1.0 100%