哥們，我心血來潮為妳努力了壹下午，求加分！
程序代碼如下:
β= 0.1；
w0 = 1；
delta _ t = 25/1000；
t(1)= 0；
因為i=1:1000
t(I+1)= delta _ t * I；
結束
閾值= 0.9；
x(1)= 100；%為了方便解決問題，最後改壹下尺度。
temp _ v 1 = 0；
n = 1；%x的角落
因為i=1:1000
對於j =-100:0.001:100
temp _ x = j；
temp _ v2 =(temp _ x-x(n))/delta _ t；
temp _ a =(temp _ v2-temp _ v 1)/delta _ t；
如果abs(temp_a+2*beta*temp_v2+w0^2*temp_x)<；閾值
n = n+1；
x(n)= temp _ x；
temp _ v 1 = temp _ v2；
打破；
結束
結束
結束
對於i=1:1001%返回到真實比例。
x(I)= x(I)/100；
結束
plot(t，x)；
xlabel(' Time ')；
ylabel('X距離')；
題目('彈簧振子阻尼振動曲線')；
彈簧振子的阻尼振動曲線可以通過直接運算得到。