Mathematic在數(shù)學(xué)中的應(yīng)用2

Mathematica在數(shù)學(xué)中的應(yīng)用

中國科學(xué)院研究生院

TheApplicationofMathematicainMath2Mathematica軟件的使用(1)用Mathematica軟件進(jìn)行各種簡單數(shù)學(xué)處理(2)用Mathematica軟件進(jìn)行作圖(3)用Mathematica軟件編寫程序本章教學(xué)目的Mathematica中的四種括弧Mathematica一般用到了四種括弧，每種括弧的功能都是不同的，下面總結(jié)一下：1、圓括號()：用于組合運(yùn)算，在計(jì)算中表示優(yōu)先計(jì)算關(guān)系2、方括弧[]：用在函數(shù)中，表示其中為函數(shù)的參數(shù)3、大括弧{}：用來表示集合4、雙括號[[]]：用來表示索引初等代數(shù)1.1有理式的展開(常用命令)命令說明Expand[poly]展開多項(xiàng)式polyFactor[poly]對多項(xiàng)式poly作因式分解FactorTerms[poly]提取數(shù)字公因子Exponent[poly,x]多項(xiàng)式poly中x的最高次數(shù)Cofficient[poly,expr]多項(xiàng)式poly中項(xiàng)expr的系數(shù)Poly[[n]]或Part[poly,n]多項(xiàng)式poly的第n項(xiàng)(1)f=Expand[(x+y+3)^2](2)Factor[f](3)Exponent[f,x](4)Coefficient[f,x](5)Factor[(x^3+2x+1)/(x^3+x^2+x+1)](6)Apart[%](*將表達(dá)式分解為最簡單分式*)1.2有理式運(yùn)算的常用命令命令說明ExpandNumerator[expr]展開分子ExpandDenominator[expr]展開分母Expand[expr]展開分子,每項(xiàng)除以分母ExpandAll[expr]分子,分母完全展開Together[expr]通分Apart[expr]分解為部分分式之和Cancel[expr]約分命令說明PolynomialQuotient[p,q,x]x的多項(xiàng)式p與q相除的商式PolynomialRemainder[p,q,x]x的多項(xiàng)式p與q相除的余式PolynomialGCD[p1,p2,….]多項(xiàng)式p1,p2,….的最大公因式PolynomialLCM[p1,p2,….]多項(xiàng)式p1,p2,….的最小公倍數(shù)1.3多項(xiàng)式的代數(shù)運(yùn)算(1)PolynomialQuotient[1+x^2,x+1,x](2)PolynomialGCD[x^2+2x+1,x^3+1,x^5+1]命令說明Solve[方程或方程組,{變量}]求方程(組)的精確解NSolve[方程或方程組,{變量}]求方程(組)的(全部)近似解FindRoot[方程,{變量,初值}]用Newton法求方程組的一個(gè)近似解1.4方程求解(1)Solve[a*x+b==0,x](2)Reduce[a*x+b==0,x](*推導(dǎo)方程組的解*)(3)FindRoot[Sin[x]==0,{x,3}]Mathematica對符號運(yùn)算的限制用戶在使用Mathematica對符號表達(dá)式進(jìn)行處理時(shí),也許會執(zhí)行Expand[(x+1)^(100^100)]的命令。這個(gè)命令執(zhí)行的結(jié)果，表達(dá)式將有100^100+1項(xiàng)。也許計(jì)算機(jī)可以很方便地執(zhí)行這個(gè)計(jì)算任務(wù)，但到計(jì)算的后期，每升高(x+1)的冪指數(shù)時(shí)，計(jì)算的增加量是非常大的。當(dāng)運(yùn)算超出計(jì)算機(jī)容量的限制時(shí)，Mathematica的計(jì)算只能停下來。即使計(jì)算下去，花費(fèi)時(shí)間也是相當(dāng)客觀(1):進(jìn)行包含10000位數(shù)的算術(shù)運(yùn)算(2):展開有1000項(xiàng)的多項(xiàng)式(3):用有幾千項(xiàng)的三個(gè)以上的變量分解多項(xiàng)式(4):使用一個(gè)遞歸計(jì)算1000次以上(5):找出100＊100矩陣的逆矩陣的計(jì)算(6):輸出結(jié)果超過10頁的計(jì)算用符號做標(biāo)記在Mathematica中，可以指定符號為不同類型的目標(biāo)作標(biāo)記。這個(gè)標(biāo)記就像數(shù)值的單位一樣。比如：微積分

Mathematica可以完成幾乎所有標(biāo)準(zhǔn)數(shù)學(xué)函數(shù)的積分,可以對有理多項(xiàng)式進(jìn)行積分,只要分母的冪指數(shù)不要太高.(1)D[Sin[n*x],{x,3}](2)Dt[Sin[n*x],x](3)Integrate[Log[x],x]有些積分可能無法用確定的式子表示,至少不能用標(biāo)準(zhǔn)的數(shù)學(xué)函數(shù)表示。比如像，這樣的積分根本不能用數(shù)學(xué)手冊中定義的標(biāo)準(zhǔn)函數(shù)的形式求解。Mathematica對于不可積分的數(shù)學(xué)表達(dá)式將不作任何處理。線形代數(shù)計(jì)算方法、擬和、最優(yōu)化計(jì)算Mathematica軟件的作圖給出一個(gè)一元函數(shù)及其作圖區(qū)間，用Plot語句可以立刻作出函數(shù)在相應(yīng)區(qū)間上的圖形(1)Plot[Sin[x],{x,0,2Pi}](2)Plot[Sin[x],{x,0,2Pi},

AspectRatio->Automatic,

PlotStyle->{GrayLevel[0.1], Dashing[{0.02,0.01}], Thickness[0.01]},

AxesLabel->{“x”,“y”}]二維函數(shù)作圖

修改坐標(biāo)的原點(diǎn)為(3,0),并設(shè)定圖形名稱修改x軸的刻度,y軸的刻度仍然不變在繪圖時(shí)，允許使用選項(xiàng)對繪制圖形的細(xì)節(jié)提出各種要求和設(shè)置。如果不設(shè)置任何選項(xiàng)，則Mathematica軟件作圖時(shí)選項(xiàng)取默認(rèn)值。選項(xiàng)默認(rèn)值意義AspectRatio1/GoldRatio=0.618圖形高度與寬度的比例。如果要圖形按實(shí)際情況顯示，設(shè)置的選項(xiàng)應(yīng)為AutomaticAxesAutomatic是否畫坐標(biāo)軸以及設(shè)置坐標(biāo)軸的中心位置.選項(xiàng)值為None時(shí),不設(shè)置坐標(biāo)軸AxesLabelNone設(shè)置坐標(biāo)軸上的標(biāo)記符號.用{“字符串1”,“字符串2”}的形式定義軸的橫坐標(biāo)和縱坐標(biāo)標(biāo)記PlotRangeAutomatic圖形中坐標(biāo)的范圍選項(xiàng)默認(rèn)值意義PlotPoints22采樣函數(shù)的點(diǎn)數(shù).對于函數(shù)值變化劇烈的表達(dá)式,應(yīng)設(shè)定較大的點(diǎn)數(shù)PlotStyleAutomatic設(shè)置曲線的樣式。取默認(rèn)值時(shí)畫出一條黑色、實(shí)的曲線。選項(xiàng)意義Graylevel[g]灰度比值，g取0到1之間的數(shù),0為黑色，1為白色RGBColor[r,g,b]紅、綠、藍(lán)三色的強(qiáng)度，0到1之間的數(shù)Thickness線的寬度值為tPlotStyle的可選項(xiàng)目二維參數(shù)作圖

使用Plot命令只能繪出一般的函數(shù)曲線，要繪制參數(shù)曲線，可以用ParametricPlot命令,其一般形式:

ParametricPlot[{x[t],y[t]},{t,tmin,tmax},選項(xiàng)](1)ParametricPlot[{Cos[t],Sin[t]},{t,0,2Pi}](2)ParametricPlot[{Cos[t],Sin[t]},{t,0,2Pi},

AspectRatio->Automatic]三維函數(shù)作圖

作出二元函數(shù)f(x,y)的立體圖形的命令是Plot3D,其格式為:Plot3D[f[x,y],{x,xmin,xmax},{y,ymin,ymax},選項(xiàng)](1)Plot3D[Sin[Sqrt[x^2+y^2]],{x,-5,5},{y,-5,5}]與Plot語句類似,Plot3D語句中也可以加入許多選項(xiàng)(2)Plot3D[Sin[Sqrt[x^2+y^2]],{x,-5,5},{y,-5,5}, Boxed->False,Axes->False,

PlotPoints->50,Mesh->False]三維函數(shù)作圖

圖形軸上加上標(biāo)記

學(xué)習(xí)過畫法幾何或工程制圖的都知道，制圖時(shí)通常用三視圖來表示一個(gè)物體的具體形狀特性。我們在生活中也知道從不同觀察點(diǎn)觀察物體，其效果是很不一樣的。在繪制立體圖形時(shí)，在系統(tǒng)默認(rèn)的情況下，觀察點(diǎn)在(1.3,-2.4,2)處。這個(gè)參考點(diǎn)選擇是具有一般性的，因此偶爾把圖形的不同部分重在一起也不會發(fā)生視覺混亂。改變視點(diǎn)

帶有陰影和網(wǎng)格的圖形對于理解曲面的形狀是很有好處的。在有些矢量圖形的輸出裝置中，你可能得不到陰影，但是當(dāng)有陰影時(shí)，輸出裝置可能要花很長時(shí)間來輸出它。通常情況下，Mathematica為了使圖形更加逼真而用明暗分布的形式給空間立體曲面著色。在這種情況下，Mathematica假定在圖形的右上方有三種光源照在物體上。但有時(shí)這種方法會造成混亂，此時(shí)你可用Lighting->False來采取根據(jù)高度在表面上涂以不同灰度的陰影的方法。選項(xiàng)默認(rèn)值意義AspectRatio1圖形高度與寬度的比率AxesTrue是否畫坐標(biāo)軸AxesLabelNone設(shè)置坐標(biāo)軸上的標(biāo)記符號BoxedTrue是否顯示外框MeshTrue是否在表面畫網(wǎng)格PlotRangeAutomatic圖形中坐標(biāo)的范圍PlotPoints15采樣函數(shù)的點(diǎn)數(shù)ShadingTrue表面是陰影還是空白Plot3D語句的各種常用的選項(xiàng)三維參數(shù)作圖

在Mathematica軟件中三維參數(shù)作圖有兩種形式,一種是空間曲線參數(shù)作圖,其命令為:ParametricPlot3D[{x[t],y[t],z[t]},{t,tmin,tmax},選項(xiàng)](1)ParametricPlot3D[{6Cos[t],6Sin[t],3*t}, {t,-8,8},AspectRatio->1]另外一種是空間曲面參數(shù)作圖,其命令為: ParametricPlot3D[{x[u,v],y[u,v],z[u,v]}, {u,umin,umax},{v,vmin,vmax},選項(xiàng)](2)ParametricPlot3D[{Cos[u]*Cos[v], Sin[u]*Cos[v],Sin[v]},{u,0,2Pi}, {v,-Pi/2,Pi/2},Axes->False,Boxed->False](2)ParametricPlot3D[{{Cos[u]*Cos[v],Sin[u]*Cos[v],

Sin[v]},{2Cos[u]*Cos[v],2Sin[u]*Cos[v], 2Sin[v]}},{u,0,Pi},{v,-Pi/2,Pi/2}, Boxed->False,Axes->False]數(shù)據(jù)作圖

在Mathematica軟件也可以根據(jù)數(shù)據(jù)作出圖形,其命令為:

ListPlot[數(shù)據(jù),選項(xiàng)](1)p=Table[{n,Prime[n]},{n,1,20}];L istPlot[p](2)ListPlot[p,PlotStyle->AbsolutePointSize[4]](3)ListPlot[p,PlotJoined->True]圖形的組合

上述的各種圖形命令中，ParametricPlot,ParametricPlot3D，Plot三個(gè)語句不僅可以畫出一個(gè)函數(shù)的圖形，而且可以同時(shí)畫出幾個(gè)函數(shù)的圖形。其一般形式為： 圖形命令[{函數(shù)1，函數(shù)2，…},變量范圍,選項(xiàng)](1)Plot[{Sin[x],x,x-x^3/6,x-x^3/6+x^5/120}, {x,-2Pi,2Pi}]

圖形元素作圖

如果要繪制一些最基本的圖形,可以用Graphics語句(三維圖形用Graphics3D)作出基本圖形元素,再用Show語句顯示圖形二維圖形元素集合意義Point[{x,y}]點(diǎn){x,y}Line[{x1,y1},{x2,y2},…]連接{x1,y1},{x2,y2},…的折線段Rectangle[{x1,y1},{x2,y2}]以{x1,y1},{x2,y2}為對角線兩頂點(diǎn)的填實(shí)矩陣Polygon[{x1,y1},{x2,y2},…]以{x1,y1},{x2,y2}為頂點(diǎn)的填實(shí)多邊形Circle[x,y,r]圓心在{x,y},半徑為r的圓Circle[{x,y},{rx,ry}]中心在{x,y},長短半軸分別為rx,ry的橢圓二維圖形元素集合意義Circle[{x,y},r,{t1,t2}]從弧度t1到弧度t2的圓弧Disk[{x,y},r]圓心在{x,y}，半徑為r的填實(shí)圓Text[“expr”,{x,y}]在{x,y}處的文本”expr”(1)執(zhí)行下列語句V1=Graphics[Circle[{0,0},{3.5,4}]];V2=Graphics[Line[{{-2,2.5},{-1,2.5}}]];V3=Graphics[Line[{{2,2.5},{1,2.5}}]];V4=Graphics[Circle[{-1.5,1.5},0.5]];V5=Graphics[Circle[{1.5,1.5},0.5]];V6=Graphics[Disk[{-1.65,1.5},0.15]];V7=Graphics[Disk[{1.35,1.5},0.15]];V8=Graphics[Polygon[{{-0.5,-1},{0.5,-1},{0,0}}]];V9=Graphics[Circle[{0,-2},{0.5,0.3}]];v10=Graphics[Text[“Hello”,{0,-5}]];Show[V1,V2,V3,V4,V5,V6,V7,V8,V9,V10,

AspectRatio->Automatic]三維圖形的元素的函數(shù)－略，感興趣的圖形查相關(guān)手冊編程

分枝結(jié)構(gòu)

在復(fù)雜的計(jì)算中經(jīng)常需要根據(jù)表達(dá)式的情況（它是否滿足一些條件）確定是否做某些處理，或在滿足不同的條件時(shí)做不同的處理。Mathematica軟件提供了一些描述條件分枝的結(jié)構(gòu)，它們在常用的程序里，用于控制程序的執(zhí)行過程。1：If語句

Mathematica軟件中If語句有三種形式

形式一：If[test,expr]當(dāng)test的值為True時(shí)，對expr求值，將它的值作為整個(gè)語句的值；當(dāng)test的值為False時(shí)，則給出空值Null

形式二：If[test,expr1,expr2]

當(dāng)test的值為True時(shí)，求expr1的值作為整個(gè)語句的值，當(dāng)test的值為False時(shí)，求expr2的值作為整個(gè)語句的值。

形式三：If[test,expr1,expr2,expr3]當(dāng)test的值為True時(shí)，求expr1的值作為整個(gè)語句的值，當(dāng)test的值為False時(shí)，求expr2的值作為整個(gè)語句的值;當(dāng)test求不出值為True與False時(shí),求expr3的值作為整個(gè)語句的值 (1)abs[x_]=If[x>=0,x,-x](2)f[x_]:=If[x>5,3,2,1](3)f[6](4)f[5](5)f[a]2：Which語句

Which[test1,expr1,test2,expr2,…]

該語句依次求出每一個(gè)條件的值，當(dāng)求出第一個(gè)值為True的條件時(shí)，求出對應(yīng)表達(dá)式的值作為整個(gè)語句的值。(例)g[x_]:=Which[x>=8,8,x>=6,6,x>=4,4,True,0]

用True作為Which語句的最后一個(gè)條件，可以處理“其他”情況。在此處即為，當(dāng)x<4時(shí)，g[x]取值為0

循環(huán)結(jié)構(gòu)

高級程序設(shè)計(jì)語言都提供了描述重復(fù)執(zhí)行的循環(huán)語句。在Mathematica軟件中也提供了一些類似的循環(huán)控制結(jié)構(gòu)。1、While[test,expr]在計(jì)算時(shí)，條件test先被求值。若求出值為True，則對表達(dá)式求值，然后再重復(fù)上述過程；一旦test的值不是True，整個(gè)循環(huán)結(jié)果計(jì)算結(jié)束。例如下面的程序可用來計(jì)算與100！k=1;s=0;p=1;While[k<=100,s=s+k;p=p*k;k++];Print[“s=”,s,“p=”,p]2、For[start,test,incr,body]在計(jì)算時(shí)，其初始表達(dá)式start首先求值，然后進(jìn)入循環(huán)，依次計(jì)算條件test，步長表達(dá)式incr與循環(huán)體body，一旦test的值不是True，整個(gè)循環(huán)結(jié)構(gòu)計(jì)算結(jié)束.s=0;p=1;For[k=1,k<=100,k++,s=s+k;p=p*k];Print[“s=”,s,“p=”,p]3、Do[expr,{i,imin,imax,di]在循環(huán)變量i依步長di從imin取到imax