1、 Matlab程序設(shè)計(jì)仿真訓(xùn)練實(shí)驗(yàn)報(bào)告設(shè)計(jì)題目：物體碰撞運(yùn)動(dòng)建模專業(yè)；通信工程班級(jí)：三班姓名：張樂(lè)學(xué)號(hào)：201100805524報(bào)告時(shí)間：2012.06指導(dǎo)老師：蔡益紅實(shí)驗(yàn)?zāi)康模毫私鈓atlab的基本特點(diǎn)與功能，基本掌握matlab的功能解決數(shù)學(xué)物理的相關(guān)問(wèn)題。實(shí)驗(yàn)環(huán)境MATLAB7.9I、題目：質(zhì)量為m的小球以速度珂正面碰撞質(zhì)量為M的靜止小球，假設(shè)碰撞是完全彈性的，即沒(méi)有能量損失，求碰撞后兩球的速度，及它們與兩球質(zhì)量比K=M/m的關(guān)系并對(duì)結(jié)果進(jìn)行分析。II、程序簡(jiǎn)介根據(jù)小球碰撞過(guò)程中的機(jī)械能守恒和動(dòng)量守恒，寫出兩個(gè)方程，然后解出兩個(gè)小球各自的速度表達(dá)式，再把已知參量代入到速度表達(dá)式即可求的小

2、球完全彈性碰撞后的速度；其中，跟據(jù)兩個(gè)小球的質(zhì)量關(guān)系K=M/m,可以分析兩個(gè)小球的質(zhì)量跟碰撞后的速度u和V的方向和大小關(guān)系。III、程序依據(jù)(機(jī)械能守恒)0.5*m*VA2-0.5*m*uA2-0.5*M*vA2=0(動(dòng)量守恒)m*V-m*u-M*v=O求解方程的：u=V*(m-M)/(m+M)；v=2*V*m/(m+M)；代入K=M/m得：u=V*(l-k)/(l+k);v=2*V/(l+k);IV、程序代碼symsuv;x,y=solve(0.5*m*VA2-0.5*m*uA2-0.5*M*vA2=0,m*V-m*u-M*v=0);x=vpa(x,4);y=vpa(y,4);clear v

3、=0 symsuv;x,y=solveC0.5*m*VA2-0.5*m*uA2-0.5*M*vA2=0,m*V-m*u-M*v=0);x=vpa(x,4)x=V-(1.0*(M*V-1.0*V*m)/(M+m)y=vpa(y,4)y=o(2.0*V*m)/(M+m)當(dāng)以K=M/m計(jì)算時(shí)：symsuv;u,v=solve(V-u-K*v=0,VA2-uA2-K*vA2=0);u=vpa(u,4)u=V(V-1.0*K*V)/(K+1.0)v=vpa(v,4)(2.0*V)/(K+1.0)VI、函數(shù)圖象x=0.1:0.1:10;y=2./(l+x)-l;z=2-2./(l+x);plot(x,y,x,z)V、數(shù)據(jù)分析當(dāng)K=1時(shí)，兩小球的質(zhì)量相等，由公式的V1=0,V2=V0，即碰撞后質(zhì)量為m的小球靜止，質(zhì)量為M的小球以速度大小V0,方向不變往前運(yùn)動(dòng)；當(dāng)K1時(shí)，碰后由速度表達(dá)式的VI為負(fù)，即速度反向，V2沿速度V0的方向運(yùn)動(dòng)；當(dāng)Kvl時(shí)，碰后由速度表達(dá)式得，兩個(gè)小球均沿V0的方向運(yùn)動(dòng)，且V1的速度小于V2的速度。實(shí)驗(yàn)總結(jié)：Matlab的功能很強(qiáng)大，對(duì)物理、數(shù)學(xué)的幫助很大。大大簡(jiǎn)便了物理