Matlab實驗大作業(yè).doc

數(shù)學實驗綜合實驗報告實驗名稱 Matlab基礎知識 學 院 計算機與通信工程學院 專業(yè)班級 通信F091 組 長 莊 原（40955008） 組 員 宋文軍（40955012） 吳克強（40955011）2011年 5月4一、 【實驗目的】1、 鞏固MATLAB的基本運算，學會綜合各種指令，合理的分配與組合，是問題變得容易解決。2、 學會把具體問題轉(zhuǎn)化為容易用MATLAB解決的數(shù)學模型，并且把數(shù)學模型轉(zhuǎn)換為MATLAB語言3學會綜合使用MATLAB并用它解決實際問題。二、 【實驗任務】問題：在一B=0.2T,方向垂直于軌道面向下的均勻磁場中,有一長L=1m質(zhì)量m=0.1kg的金屬桿,沿一傾角=30金屬滑道由靜止下滑,若滑道與自感系數(shù)L=0.5H的線圈相連,分析該金屬桿的運動速度及線圈內(nèi)電流隨時間的變化關系。L三、 【本實驗數(shù)學模型的構造】首先，明確本實驗的自變量和因變量。根據(jù)題目，自變量是線圈的下降速度，因變量也是要求的就是求一個變化關系，是一個變化關系函數(shù)。而顯然，這樣的變化函數(shù)是一個復合函數(shù)。那么對應到我們的程序中也應該是一個嵌套的模式。再根據(jù)物理知識,列出相關變量的物理方程,下降速度及線圈電流與時間關系的微分方程,建立起相應數(shù)學模型,用Matlab求解微分方程,得到最后的關系分析。四、 【實驗程序及解釋】%建立并定義子函數(shù)%感應線運動的微分方程:m(dv/dt)=mgsin-Bai;L(di/dt)=Bavfunction yp=indctn(t,y)global m a L g sita Byp=g*sin(sita)-B*a*y(2)/m B*a*y(1)/L%主函數(shù)clearglobal m a L g sita Bm=0.1;a=1;L=0.1;g=9.8;sita=pi/6;B=0.05tspan=0 100y0=0 0t,y=ode23(indctn,tspan,y0)subplot(2,1,1)plot(t,y(:,1),r)xlabel(時間(s);ylabel(下落速度(m/s)grid onsubplot(2,1,2)plot(t,y(:,2),b)xlabel(時間(s);ylabel(線圈電流)grid on五、 【實驗結(jié)果】六、 【實驗總結(jié)】 1、到目前，已經(jīng)大體上了解了Matlab的系統(tǒng)結(jié)構及其工作環(huán)境，對各個窗口的應用已有基本的認識。 2、通過設置數(shù)據(jù)類型可以得到不同數(shù)據(jù)精度，并且可以分析不同數(shù)據(jù)類型的相同與不同之處。用MATLAB畫各種數(shù)學圖形。有一個直觀的認識。