求位移:遍历t,t在所有值上,杜哈梅积分的解: for i = 1:length(t) %{ 杜哈梅积分Duhamel Integral的方程式示意: yp(t) = (1/(m*wd)) * integral(F(tau) * exp(-zeta*wd*(t-tau)) * sin(wd*sqrt(1-zeta^2)*(t-tau)), tau=0 to t) 在此基础上设计此函数。此...
c_new = exp_decay * (sin_inc * s + cos_inc* c);s = s_new;c = c_new;u(i) = s / (omega_d * k);endfigure;plot(t,u, 'b-', 'LineWidth',2);xlabel('Time (s)');ylabel('Displacement (m)');title('Numerical Response Using Enhanced Duhamel Integral');gridon;...
求位移:遍历t,t在所有值上,杜哈梅积分的解: for i = 1:length(t) %{ 杜哈梅积分Duhamel Integral的方程式示意: yp(t) = (1/(m*wd)) * integral(F(tau) * exp(-zeta*wd*(t-tau)) * sin(wd*sqrt(1-zeta^2)*(t-tau)), tau=0 to t) 在此基础上设计此函数。此...