利用matlab研究受空氣阻尼的拋體運(yùn)動應(yīng)用物理學(xué)專業(yè)課程設(shè)計(jì)(留存版)
【正文】 課程設(shè)計(jì)圖8張456 湖南工業(yè)大學(xué) 課程設(shè)計(jì)任務(wù)書2011—2012 學(xué)年第 1 學(xué)期 理學(xué)院 學(xué)院（系、部） 應(yīng)用物理學(xué) 專業(yè) 081 班級課程名稱： 計(jì)算物理 一、 設(shè)計(jì)題目：利用Matlab研究受空氣阻尼的拋體運(yùn)動二、 完成期限：自 2011 年 11 月 28 日至 2011 年 12 月 2 日共 1 周內(nèi)容及任務(wù)分析拋體運(yùn)動所受空氣阻力的三種情況，寫出質(zhì)點(diǎn)運(yùn)動的微分方程；用龍格庫塔方法求解質(zhì)點(diǎn)運(yùn)動的微分方程；畫出粒子的空間軌跡和速度隨時(shí)間變化的彗星圖，找出軌跡的最高點(diǎn)，到達(dá)最高點(diǎn)的時(shí)間和速度。例如，采用MATLAB編制自由振動的子程序時(shí)只需要調(diào)用兩個(gè)函數(shù)：求逆矩陣的inv函數(shù)及求特征值和特征向量的eig函數(shù)。一階常微分方程可以寫作：y39。px=[。 %J?8 amp。dx/dt39。通過這次課程設(shè)計(jì)，我掌握了測試；熟悉了Matlab運(yùn)行環(huán)境和如何編寫Matlab語言，了解了項(xiàng)目開發(fā)的基本思想和概念。果然是團(tuán)結(jié)就是力量，只有互相之間默契融洽的配合才能換來最終完美的結(jié)果。這次課程設(shè)計(jì)終于順利完成了，在設(shè)計(jì)中遇到了很多問題，最后在老師的指導(dǎo)下，終于游逆而解。 %yl7rfXsubplot(2,1,2)axis([0 10 0 4]) hold on xlabel(39。strdd{3}=39。程序global m d e m=1。由于此算法精度高，采取措施對誤差進(jìn)行抑制，所以其實(shí)現(xiàn)原理也較復(fù)雜。常用的數(shù)值方法有差分法、有限元法、加權(quán)殘值法、邊界元法登。在理論力學(xué)和普通物理學(xué)的范圍內(nèi)，一般說來，拋射體的運(yùn)動可看作是質(zhì)點(diǎn)的運(yùn)動，在假定空氣阻力可以忽略或阻力至于速度的一次放成正比的情況下，用解析法能夠很容易的求解運(yùn)動方程，但這種假定，只有在被拋射出的物體運(yùn)動速度很低時(shí)才能近似成立。本課題的中心是利用matlab來解決阻尼拋體運(yùn)動，當(dāng)然想到用建立模型的方法去求解。K2=f(Xn+h/2,Yn+(h/2)*K1)。]。 % ?pF6hold on %2Ixlabel(39。由Vx—t圖可知，水平方向速度不變。這對于我們的將來也有很大的幫助。就像中國提倡的艱苦奮斗一樣，我們都可以在實(shí)驗(yàn)結(jié)束之后
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1