[問題] 如何畫出振動系統中時間-加速度的圖?已回收
tspan = [0: 0.001: 5];
y0 = [0; 0; 0; 0; 0; 0; 0; 0];
[t, y] = ode23 (@dfunc1, tspan, y0);
plot(t, y(:, 1));
xlabel('t');
ylabel('Zt(t)');
-----------------------------------------------
function f = dfunc1(t, y)
f = zeros(8, 1);
ct = 1360;
kt = 45005.3;
cc = 131.59;
kc = 2500;
cs = 980;
ks = 16000;
ky = 160000;
mt = 110.4/7;
mp = 276/7;
ms = 240;
mu = 36;
f(1) = y(2);
f(2) = (-ct*(y(2)-y(4))-kt*(y(1)-y(3)))/mt;
f(3) = y(4);
f(4) = (ct*(y(2)-y(4))+kt*(y(1)-y(3))-cc*(y(4)-y(6))-kc*(y(3)-y(5)))/mp;
f(5) = y(6);
f(6) = (cc*(y(4)-y(6))+kc*(y(3)-y(5))-cs*(y(6)-y(8))-ks*(y(5)-y(7)))/ms;
f(7) = y(8);
f(8) = (ks*(y(5)-y(7))+cs*(y(6)-y(8))-ky*(y(7)-0.02))/mu;
以上是我的程式碼
我現在只會畫出時間t-位移y(1) 和 時間t-速度y(2)的關係圖
可是時間-加速度我一直不知道要怎麼畫
麻煩大大幫我解答一下
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.118.199.92
※ 編輯: u8911059 來自: 140.118.199.92 (02/09 23:34)
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
