結(jié)果：

改變相關(guān)變量后的變化情況：

請(qǐng)注意：該模擬結(jié)果僅僅具有參考意義而并無(wú)實(shí)際的物理意義

改變初始位置我們可以看到振幅發(fā)生了變化。

改變動(dòng)摩擦因數(shù)我們看到波動(dòng)向左側(cè)壓縮。

改變兩轉(zhuǎn)軸間的間距可以看到波動(dòng)向右側(cè)拉長(zhǎng)。

相關(guān)參數(shù)在代碼中給出。歡迎大家討論，如果這篇文章對(duì)你有幫助的話記得三連支持一下；

謝謝！

matlab代碼：

%%初步的動(dòng)力學(xué)基礎(chǔ)以及初步求解

clear all

global x l p d h g?

x=0.20;%兩軸的一半間距

l=0.4;%重物的長(zhǎng)度

p=1.5*10^3;%重物的密度

d=0.05;%重物的寬度

h=0.01;%重物的高度

g=9.8;%重力加速度

m=p*h*l*d;%重物的質(zhì)量

u=0.46;%動(dòng)摩擦因數(shù)

%——————————————————————————

a2=zeros(10000,1);

v2=zeros(10000,1);

x2=zeros(10000,1);

x2(1,1)=0.01;%重心偏移原點(diǎn)的位置

t=0.001;%時(shí)間步長(zhǎng)

G=m*g;%重力

for i=2:10000

G=m*g;%重力

N1=((x*G-G*x2(i-1,1))/2+x2(i-1,1)*G)/x;%受力分析的方程

N2=(x*G-G*x2(i-1,1))/(2*x);? ? ?%受力分析的方程

a2(i,1)=(N2*u-N1*u)/m;? ? ? ? ? ? %受力分析的方程

v2(i,1)=v2(i-1,1)+a2(i,1)*t;? %速度

x2(i,1)=x2(i-1,1)+v2(i,1)*t;? %位移

end

hold on

plot(x2,'Color','b')

%plot(v2,'Color','r')

%plot(a2,'Color','c')