數(shù)學實驗matlab解常微分方程

《數(shù)學實驗matlab解常微分方程》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在學術(shù)論文-天天文庫。

1、GUIZHOUUNIVERSITY實驗課程名稱數(shù)學實驗實驗項目名稱MATLAB解常微分方程年級2010級專業(yè)信計(101)學生姓名成富學1007010167理學院實驗時間：20012年5月9日學生實驗室守則一、按教學安排準時到實驗室上實驗課，不得遲到、早退和曠課。二、進入實驗室必須遵守實驗室的各項規(guī)章制度，保持室內(nèi)安靜、整潔，不準在室內(nèi)打鬧、喧嘩、吸煙、吃食物、隨地吐痰、亂扔雜物，不準做與實驗N容無關(guān)的事，非實驗用品一律不準帶進實驗室。三、實驗前必須做好預(yù)習（或按要求寫好預(yù)習報告），未做預(yù)習者不準參加實驗。四、實驗必須服從教師的安排和指導，認真按規(guī)程操作，未經(jīng)教師允許不得擅

2、自動用儀器設(shè)備，特別是與本實驗無關(guān)的儀器設(shè)備和設(shè)施，如擅自動用或違反操作規(guī)程造成損壞，應(yīng)按規(guī)定賠償，嚴重者給予紀律處分。五、實驗中要節(jié)約水、電、氣及其它消耗材料。六、細心觀察、如實記錄實驗現(xiàn)象和結(jié)果，不得抄襲或隨意更改原始記錄和數(shù)據(jù)，不得擅離操作崗位和干擾他人實驗。七、使用易燃、易爆、腐蝕性、有毒有害物品或接觸帶電設(shè)備進行實驗，應(yīng)特別注意規(guī)范操作，注意防護；若發(fā)生意外，要保持冷靜，并及時向指導教師和管理人員報告,不得自行處理。儀器設(shè)備發(fā)生故障和損壞，應(yīng)立即停止實驗，并主動向指導教師報告，不得自行拆卸查看和拼裝。八、實驗完畢，應(yīng)清理好實驗儀器設(shè)備并放回原位，清掃好實驗現(xiàn)場，經(jīng)

3、指導教師檢杳認可并將實驗記錄交指導教師檢杳簽字后方可離去。九、無故不參加實驗者，應(yīng)寫出檢查，提出申請并繳納相應(yīng)的實驗費及材料消耗費,經(jīng)批準后，方可補做。十、自選實驗，應(yīng)事先預(yù)約，擬訂出實驗方案，經(jīng)實驗室主任同意后，在指導教師或?qū)嶒灱夹g(shù)人員的指導下進行。十一、實驗室內(nèi)一切物品未經(jīng)允許嚴禁帶出室外，確需帶出，必須經(jīng)過批準并辦理手續(xù)。姓名成富學號1007010167實驗組實驗時間2012-5-9指導教師劉軼中成繢實驗項目名稱MATLAB解微積分學生所在學院:理學院專業(yè):信計班級：101實驗?zāi)康募耙螅簩嶒災(zāi)康模?：學習MATLAB有關(guān)常微分方程計算指令；2:學會常微分方程(組)的

4、數(shù)値求解Euler法和剛性方程組問題:3：分析數(shù)值'汁算法求解數(shù)學問題的局限性：4：建模實驗研宂常微分方程。實驗要求：在MATLAB環(huán)境下熟練運用各知識作圖實驗(或算法〉原理：根據(jù)函數(shù)的表達式，建立M文件規(guī)則MATLAB中語句的功能法輸入相應(yīng)的句子。建立M文件。運用數(shù)學函數(shù)中基本函數(shù)作圖。、運用plot(x,y),ezplot(x,y)輸入相應(yīng)的變量。MATLAB將會根據(jù)輸入的變作出相應(yīng)圖案實驗硬件及軟件平臺：PCMATLAB實驗步驟：1打幵MATLAB,進入到窗口狀態(tài)。2根據(jù)要求輸入要求的題目。3運用MATLAB特殊的函數(shù)編程。4檢驗程序是否正確。5運用結(jié)果。6分析得到的

5、實驗數(shù)裾。實驗內(nèi)容(包括實驗具體內(nèi)容、算法分析、源代碼等等):1(1)：>>fun=inline(,x+y,,，x，，，y,);?[t，y]=ode45(fun,[0123],1)01231.00003.436611.778136.17111(2)。?fun=inline(,[-2*y(l)-3*y(2);2*y(l)+y(2)]*/t,y,);?[t，y]=ode45(fun，[010]，[-2.7;2.8]);plot(y(:，l)，y(:，2))?plot(y(:，l)，y(:，2))1(3)：?clear;close;fun=@(t,y)[y(2);0.01*y(

6、2)A2-2*y(l)+sin(t)];[t,y]=ode45(fun，[05]，[0;l]);plot(t，y(:，1))1(4)：?clear;close;fun=@(t,y)[y(2);2.5*y(2)+1.5*y(l)+45*exp(2*t)];[t,y]=ode45(fun,[02]，[2;1]);plot(t，y(:，1))1(5)：?clear;close;fun=@(t，y，mu)[y(2);-mu*(y⑴八2-1)*y(2)-y(l)];[t，y]=ode45(fun，[020]，[2;0]，[]，l);mu=l.plot(y(:，l)，y(:，2));h

7、oldon;[t,y]=ode45(fun,[020],[2;0],[],2);mu=2.plot(y(:，l)，y(:，2)，’r’);holdoff;mu=1mu=21(6)：?clear;close;sinit=bvpinit(l:0.5:3，[2;0])odefun=inline(’[y(2);-2/t*y(2)+(2*y(l)+10*cos(log⑴))/t/t]'，’t’，’y’);bcfun=inlineC[ya(l)-l;yb(l)-3],/ya,,,yb,);sol=bvp4c(odefun，bc