資源描述:
《數(shù)模競賽軟件選講ppt課件.ppt》由會員上傳分享,免費在線閱讀,更多相關內(nèi)容在教育資源-天天文庫。
1、數(shù)學軟件選講MathematicaMatlabSAS第一篇Mathematica基礎知識作為一門新的編程語言圖形處理(二維、三維及其參數(shù)方程的形式)極限、微分與積分求解方程(組)、微分方程(組)在線形代數(shù)方面的應用數(shù)值處理文件及其它高級操作第一章基礎知識一、Mathematica3.0界面及運行介紹二、基本數(shù)值運算1.整數(shù)運算:加、減、乘、除、冪、階乘2.數(shù)學常量:E、Pi、I、Degree、Infinity3.函數(shù)及數(shù)學函數(shù)4.浮點數(shù)及復數(shù)運算:N函數(shù)三、變量及表達式1.變量的定義及清除◆變量的特點(1)變量的默認作用域是全局的(2)全局變量不需事先定
2、義或聲明(3)盡量避免使用下劃線定義變量2.多項式及其操作(1)定義、替換符操作(2)常用操作:Expand、Factor、Together、PartSimplify、Collect、Coefficient、Exponent四、序列及其操作1.序列的定義2.序列的生成:Table函數(shù)3.序列的操作(1)添加刪除:Append、Prepend、Insert、Delete、DeleteCases(2)取元素:Part、Take、Drop、Select(3)檢測:Length、Count、Position五、表達式“頭”的概念:Head及Apply函數(shù)六、自定義函
3、數(shù)1.一元函數(shù)例:Clear[f,x]f[x_]:=x^2+4x-22.多元函數(shù)例:f[x_,y_]:=x^2+y^2-33.迭代函數(shù)例:f[n_]:=f[n-1]+f[n-2];f[0]=1;f[1]=1;1·條件語句◆邏輯判斷符==>=<=>
4、
5、&&◆/;運算符x=a/;test僅當test為True時才執(zhí)行賦值語句◆If語句語法:If[test,then,else]若test為True,則執(zhí)行then,若test為False,則執(zhí)行else.◆Which語句語法:Which[test1,value1,tes
6、t2,…]依次計算testi,給出對應第一個test為True的value◆Switch[expr,form1,value1,form2,…]比較expr與formi,給出與第一個form值匹配的value例1.定義如下的函數(shù):使用/;定義:f[x_]:=0/;x<=0f[x_]:=x/;x>0&&x<=2f[x_]:=x^2/;x>2②使用If定義:f[x_]:=If[x<=0,0,If[x>2,x^2,x]]③使用Which定義:f[x_]:=Which[x<=0,0,x>2,x^2,True,x]2·輸出語句Print3·循環(huán)語句◆Do語句語法:Do[
7、expr,{i,imin,imax,di}]計算expr,i=imin,…,imax,步長為di◆While語句語法:While[test,body]當test為True時,計算body◆For語句語法:For[start,test,incr,body]以start為起始值,重復計算body和incr,直到test為False時為止◆循環(huán)控制語句Break和ContinueBreak[]退出最里面的循環(huán)Continue[]轉(zhuǎn)入當前循環(huán)的下一步基本二維圖形①Plot[f,{x,xmin,xmax}],用于繪制形如y=f(x)的函數(shù)的圖形。當將多個圖形繪制在同一
8、坐標系上時,形如:Plot[{f1,…,fn},{x,xmin,xmax}]注意:有時需要使用Evaluate函數(shù)。第三章圖形處理例:在同一坐標系下繪出sinx,sin2x,sin3x,sin4x,sin5x的圖形。常用的選項:PlotStyle->Hue[a]設置線條顏色PlotRange->{a,b}控制顯示范圍DisplayFunction控制圖形顯示AspectRatio圖形的寬、高比AxesOrigin設置原點坐標程序:Clear[a,y,x]v=200;g=9.8;y[a_,x_]:=Tan[a]*x-g*x^2*Sec[a]^2/(2v^2)P
9、lot[Evaluate[Table[y[i,x],{i,Pi/12,5Pi/12,Pi/12}]],{x,0,4000}]例:有如下的拋物線簇:②ListPlot[List],用于繪制散點圖。注意,List的形式應為:例:在同一坐標系下繪制下列兩組散點圖p1={{0,0},{0,45},{5.3,89.6},{22.6,131.2}};p2={{0,0},{2.68,44.8},{12.57,88.28},{27,130.3}};程序:g1=ListPlot[p1,PlotJoined->True,DisplayFunction->Identity];g2
10、=ListPlot[p2,PlotJoined->T