歐陽法創編陽法其中建模局部包括對天花板的建模、對連接直線的建模、對彈簧的建模、對振動小球的建模。模擬彈簧振動的過程主要包括：①更新小球和彈簧的位置。②更新水平線的位置和長度。③實時地畫出彈簧高度與時間的關系曲線。 1歐陽法創編021.03.09紹...................................1作.................................2畫................................4程序功能分析.....................................4模............................6天花板的建模...................................7直線的建模......................................8彈簧的建模......................................9球的建歐陽法創編歐陽法創編彈簧振動的實時動彈簧振子簡諧振動的物理模結19結..........................................23考文獻............................................24致謝..............................................歐陽法創編歐陽法創編附錄-彈簧阻尼振動程序...................................26程序功能分析%畫直線歐陽法創編歐陽法創編%彈簧tl=r*sin(C);%球%句柄[黃線]td=[];yd=[];歐陽法創編歐陽法創編Y=-(y-2)*Dy+7;Yf=Y(end)+tl;ididTydydYend)];程序功能介紹：模擬彈簧振子簡諧運動的過程，同時畫出相應圓球質心隨時上面的程序為模擬彈簧振子的簡諧振動，假設要模擬彈簧阻尼振動過程的程歐陽法創編歐陽法創編])；%畫頂板%畫直線歐陽法創編歐陽法創編%彈簧tl=r*sin(C);%球天花板的建模歐陽法創編歐陽法創編歐陽法創編歐陽法創編直線的建模彈簧的建模歐陽法創編歐陽法創編的曲線。如取決于是向量的行還是列的尺寸與矩陣相匹配。球白勺建模歐陽法創編歐陽法創編tl=r*sin(C);l歐陽法創編021.03.09笫一個頂點連接來封通過上幾步的建模我們可以得到如下所示的為下一步實時動畫做好準備我們需要重新設plot([0,15],[3.3,3.3],black);歐陽法創編歐陽法創編軸上刻度記號的位置以及每隔多少值會到下一個刻度記號置。對平彳新位置軸的設置為:plot([0,l5],[3.3,3.3],’black),plot為畫直線的常用函數彈簧振動的實時動畫彈簧振子簡諧振動的物理模型彈簧振子簡諧振動的物理模型彈簧振子簡諧振動的物理模型子的簡諧振動歐陽法創編歐陽法創編現在我們來定量地分析上圖所示的彈簧振子設彈簧的勁度系數為k,小球的質量為m,x彈性回復力Fx可M是兩個積分常數，他們的物理意義和確定方法將在后面討論。將式(1.4)對時間求一階和二階導數，得到小歐陽法創編歐陽法創編球運動速度和加速度sin()vxAl:o(1.4)的解即可寫成遇險函數的形式，也可以寫成正弦函x對平衡未知的位移)按余變化。因此，只在線性彈性回復力作用下的彈簧振作諧振動物體的速度和加速度也是按余弦或正弦函數規律隨時間變化的。見下的加速度大小總是與其位移大小成正比，二歐陽法創編歐陽法創編者符號相反，這一結論被視為諧振動利用利用實現彈簧的實時振動實現彈簧的實時振動實現彈簧的實時振動由節可知彈簧振動是的變化主要表現在彈簧和小球化，而且小球的位移是按正弦規律變化的。據此我們先設計一個算法。現用流程Matlab是通過以下程序實現彈簧振動的實時動畫歐陽法創編歐陽法創編%畫頂板%畫直線%彈簧C二():.l:2*pi;i?二0.35;tl=r*sin(C);Fl=fill(l3+r*cos(C),2+tl,V);%球歐陽法創編歐陽法創編plot([0,15],[3.3,3.3],’black)%句柄潢線]歐陽法創編歐陽法創編Yf=Y(end)+tl;該程序的建模局部已經在前面介紹過了，現%句柄潢線]歐陽法創編歐陽法創編位置為〔2,9〕,標注為'理想中的彈通過上面的可得到如下列圖所示的界面。模擬彈簧振子簡諧振動的過程包括以下內容：①更新小球和彈簧的位置。歐陽法創編歐陽法創編②更新水平線的位置和長度。③實時地畫出彈簧高度與時間的關系曲線。Yf=Y(end)+tl;歐陽法創編歐陽法創編Dy=(3/2-l/2*sin(pi*T))*l/2;計算T時刻彈簧對平衡位ld=[ld,l];yd=[yd,Y(cncl)];更新運動曲線的數函數是返回一個的隨機矩陣，這里用來改變小球的d黃線的數據。這里句柄黃線的主要作用為將振動彈簧的位移以運動曲線的形式變現出來，以便我們跟好的了歐陽法創編歐陽法創編調試結果調試結果調試結果0123理理理理理理理理理理理理012歐陽法創編歐陽法創編理理理理理理理理歐陽法創編歐陽法創編0123理理理理理理理理理理理理00123理理理理理理理理理理理理2021.03.09歐陽法創編2021.03.09歐陽法創編0123理理理理理理理理理理理理設計總結設計總結計總結更多的了解，對一些基2021.03.09歐陽法創編021.03.09參考文獻參考文獻參考文獻[1]徐東艷.MATLAB函數庫查詢辭典.中國鐵道出[2]蘇金明.MATLAB實用教程.電子工業出版[3]李麗.MATLAB工程計算及應用.人民郵電出版[4]蔡旭輝.MATLAB根底與應用教程.人民郵電出版[5]周建興.MATLAB從入門到精通.人民郵電出版[6]張錚.MATLAB程序設計與實例應用.中國鐵道出2021.03.09歐陽法創編歐陽法創編致謝致謝在本次設計的整個過程中，得到了敖老師和同學的大力幫助，借此時機向他感謝在百忙之中對我的論文進行評審并提-彈簧阻尼振動彈簧阻尼振動彈簧阻尼振動彈簧阻尼振動程序程序程序下面來模擬彈簧振子在阻尼力的作用下進歐陽法創編歐陽法創編行阻尼運動的過程，同時畫出相應圓球質心隨時間變化的曲線。相應的歐陽法創編歐陽法創編MATLAB程序如下：生成彈簧的橫坐標數據值C=0:.l:2*pi;%生成圓球的角度數據y值歐陽法創編歐陽法創編r,);%畫出運動曲線;td=[];%記錄時間的變量T=();%設置初始時的時間值T<12;%利用循環處理阻尼運動過程的模擬彈簧對平衡-位置的位移Y二(y-2)*Dy+7;%計算彈簧的縱坐標數值d歐陽法創編歐陽法創編更新跟蹤線的數據I<d=find(diff(sign(diff(yd)))==-2)+l;%計算極大值的X=td(Kd);