數(shù)學實驗第一章

數(shù)學實驗第一章第1頁，課件共40頁，創(chuàng)作于2023年2月《數(shù)學實驗》第一章數(shù)學實驗與MATLAB向量創(chuàng)建與一元函數(shù)圖形矩陣創(chuàng)建與二元函數(shù)圖形數(shù)據(jù)顯示與字符串操作第2頁，課件共40頁，創(chuàng)作于2023年2月數(shù)學實驗——使用數(shù)學軟件快速獲取和處理有價值數(shù)據(jù)——觀察猜測驗證獲取信息實驗是獲取信息的一種活動以數(shù)據(jù)為載體！信息的表現(xiàn)形式：文字，聲音，圖片...MATLAB第3頁，課件共40頁，創(chuàng)作于2023年2月

數(shù)學實驗的學習內(nèi)容熟練掌握MATLAB的命令操作方式掌握MATLAB程序設(shè)計方法以MATLAB為操作平臺完成實驗作業(yè)第4頁，課件共40頁，創(chuàng)作于2023年2月思考問題，完成實驗報告,提升數(shù)學能力探月衛(wèi)星速度計算實驗第5頁，課件共40頁，創(chuàng)作于2023年2月一種解釋式語言.易學易用、代碼短效率高、具有強大的數(shù)值計算和繪圖功能、擴展性強.矩陣的數(shù)值運算、數(shù)值分析、數(shù)值模擬數(shù)據(jù)可視化、繪制2維/3維圖形可以與FORTRAN、C/C++做數(shù)據(jù)鏈接幾百個核心內(nèi)部函數(shù)幾十個工具箱(信號處理、自動控制、···

)MATLAB

(MATrixLABoratory)第6頁，課件共40頁，創(chuàng)作于2023年2月MATLAB桌面命令窗口工作空間瀏覽命令歷史窗口第7頁，課件共40頁，創(chuàng)作于2023年2月命令窗口、編輯窗口、圖形窗口第8頁，課件共40頁，創(chuàng)作于2023年2月命令窗口操作例1.1

A=magic(3)A=816357492例1.2

logoloadlogo;mesh(L)第9頁，課件共40頁，創(chuàng)作于2023年2月

例1.3

國際象棋發(fā)明人向印度國王求賜大麥,希望得到大麥數(shù)量由如下規(guī)則計算.在國際象棋棋盤的64個方格中,第一格放一粒麥粒,第二格放兩粒,第三格放四粒,……,以此類推.每格比前一格麥粒數(shù)多一倍,直到放滿64格為止.計算麥粒數(shù)表明這些大麥幾乎可以覆蓋地球表面.

級數(shù):S=1+2+22+23+···+263=方法二:n=0:1:63;S=sum(2.^n)S=1.8447e+019方法一:

N=2^64-1N=1.8447e+019第10頁，課件共40頁，創(chuàng)作于2023年2月程序窗口操作proportion=89.5963(覆蓋地球面百分比)錄入程序,調(diào)式成功后,將程序文件存盤、命名。在命令窗口中鍵入文件名并回車

chess第11頁，課件共40頁，創(chuàng)作于2023年2月例1.4

給定

=150,300,450,600,計算sin的值alpha=[15,30,45,60]*pi/180;sin(alpha)創(chuàng)建向量三種基本方法:

使用方括號、使用冒號、使用等分函數(shù)。ans=0.25880.50000.70710.8660方括號直接輸入法是創(chuàng)建向量(和矩陣)的常用方法，輸入時將向量元素用方括號“[]”括起來,元素之間用逗號（或空格）隔開.

alpha=(15:15:60)*pi/180;sin(alpha)第12頁，課件共40頁，創(chuàng)作于2023年2月

例1.5

用線性等分函數(shù)linspace()創(chuàng)建[0,2]上的等分點,繪出正六邊形和正十二邊形.冒號表達式創(chuàng)建向量使用格式

x=x0:step:xn當步長

step=1時可省略為

x=x0:xn當步長

step為負數(shù)時

x0應大于

xnalpha=linspace(0,2*pi,7)bata=linspace(0,2*pi,13);x1=cos(alpha);y1=sin(alpha);x2=cos(bata);y2=sin(bata);plot(x1,y1,x2,y2)————向量創(chuàng)建與一元函數(shù)圖形————第13頁，課件共40頁，創(chuàng)作于2023年2月sin正弦函數(shù)

asin反正弦函數(shù)

cos余弦函數(shù)

acos反余弦函數(shù)

tan正切函數(shù)

atan反正切函數(shù)

cot余切函數(shù)

acot反余切函數(shù)

sec正割函數(shù)

asec反正割函數(shù)

csc余割函數(shù)

acsc反余割函數(shù)

sinh雙曲正弦函數(shù)

asinh反雙曲正弦函數(shù)

cosh雙曲余弦函數(shù)

acosh反雙曲余弦函數(shù)

tanh雙曲正切函數(shù)

atanh反雙曲正切函數(shù)

sech雙曲正割函數(shù)

asech反雙曲正割函數(shù)

csch雙曲余割函數(shù)

acsch反雙曲余割函數(shù)

coth雙曲余切函數(shù)

acoth反雙曲余切函數(shù)

三角函數(shù)與雙曲函數(shù)第14頁，課件共40頁，創(chuàng)作于2023年2月abs(x)絕對值sqrt(x)開平方conj(z)共軛復數(shù)round(x)四舍五入floor(x)舍去正小數(shù)rat(x)

分數(shù)表示gcd(x,y)最大公因數(shù)exp(x)自然指數(shù)log(x) 自然對數(shù)Log10(x)以10為底對數(shù)angle(z)

復數(shù)z的相角real(z)

復數(shù)z的實部imag(z)

復數(shù)z的虛部fix(x)

舍去小數(shù)取整ceil(x)

加入正小數(shù)取整sign(x)

符號函數(shù)rem(x,y)

求x除以y的余數(shù)lcm(x,y)

最小公倍數(shù)pow2(x)

以2為底的指數(shù)log2(x)

以2為底的對數(shù)第15頁，課件共40頁，創(chuàng)作于2023年2月例1.6

用ezplot()命令繪衰減振蕩曲線函數(shù):

y=e-0.5x

sin5x

圖形.ezplot('exp(-0.5*x)*sin(5*x)',[0,10,-1,1])————向量創(chuàng)建與一元函數(shù)圖形————簡易繪圖命令第16頁，課件共40頁，創(chuàng)作于2023年2月在解決實際問題時,如果頻繁使用同一個數(shù)學表達式,則應該定義一個臨時函數(shù)以方便操作.

定義方法:

函數(shù)名=inline(‘表達式’)

例1.7

定義函數(shù)并分析函數(shù)性質(zhì)。fun=inline('x.*sin(1./x)')fplot(fun,[-0.15,0.15])N=1:5;x=2./(2*N+1)/pi;y=fun(x)y=-0.21220.1273-0.09090.0707-0.0579————向量創(chuàng)建與一元函數(shù)圖形————第17頁，課件共40頁，創(chuàng)作于2023年2月ezplot()

簡易繪圖方法，優(yōu)點：快速方便fplot()

函數(shù)繪圖方法,與簡易繪圖相似.要給定自變量變化范圍plot()

基本繪圖方法，利用一元函數(shù)自變量的一系列數(shù)據(jù)和對應函數(shù)值數(shù)據(jù)繪圖。具有很大靈活性例如plot(X,Y),plot(x1,y1,x2,y2)plot(X,Y,’r’),plot(x1,y1,’r’,x2,y2,’b’)————向量創(chuàng)建與一元函數(shù)圖形————MATLAB一元函數(shù)繪圖方法第18頁，課件共40頁，創(chuàng)作于2023年2月例1.8

用基本繪圖方法繪衰減振蕩函數(shù)

y=e–0.5xsin5x的圖形并用虛線表示振幅衰減情況。x=0:0.1:4*pi;y=exp(-0.5*x);y1=y.*sin(5*x);plot(x,y1,x,y,’--r’,x,-y,’--r’)第19頁，課件共40頁，創(chuàng)作于2023年2月例1.9

用直接方法創(chuàng)建3階希爾伯特矩陣formatratH=[1,1/2,1/3;1/2,1/3,1/4;1/3,1/4,1/5]bar3(H)H=

11/21/31/21/31/41/31/41/5第20頁，課件共40頁，創(chuàng)作于2023年2月注意事項矩陣元素必須在方括號[]之內(nèi);同一行相鄰元素間用逗號或空格分隔;

矩陣的行與行之間用分號分隔.————矩陣創(chuàng)建常用方法————1.直接輸入法;2.特殊矩陣函數(shù)法;3.

數(shù)據(jù)文件輸入B=9-3630-36192-18030-180180直接輸入法特殊矩陣函數(shù)法A=hilb(3)%用函數(shù)創(chuàng)建希爾伯特矩陣B=invhilb(3)%創(chuàng)建希爾伯特矩陣的逆陣A*B%驗證B為A的逆ans=100010001A=

11/21/31/21/31/41/31/41/5

第21頁，課件共40頁，創(chuàng)作于2023年2月——特殊矩陣函數(shù)表

——zeros(m,n) m×n階零矩陣

eye(m,n) m×n階單位矩陣ones(m,n) m×n階全1矩陣rand(m,n) m×n階隨機矩陣

randn(m,n)正態(tài)隨機數(shù)矩陣magic(n) n階魔方矩陣hilb(n) n階Hilbert矩陣

invhilb(n)逆Hilbert矩陣pascal(n) n階Pascal矩陣vander(C) 由向量C生成范德蒙矩陣第22頁，課件共40頁，創(chuàng)作于2023年2月

例1.10

創(chuàng)建4階幻方矩陣A,并驗證矩陣A各列元素之和、各行元素之和以及各對角元之和均為常數(shù)34。A=16231351110897612414151A=magic(4)sum(A)sum(A’)sum(diag(A))B=A(:,4:-1:1)sum(diag(B))ans=34343434%求列和%求行和%求A對角和%矩陣翻轉(zhuǎn)%求B對角和第23頁，課件共40頁，創(chuàng)作于2023年2月X=

-2-1012-2-1012-2-1012-2-1012-2-1012創(chuàng)建網(wǎng)格矩陣命令:[X,Y]=meshgrid(x,y)

例1.11

計算二元函數(shù)

z=xexp(–x2–y2)網(wǎng)格點值

[X,Y]=meshgrid(-2:2,-2:2)Y=

-2-2-2-2-2-1-1-1-1-1000001111122222二元函數(shù)圖形繪制方法：矩陣表示網(wǎng)格、計算網(wǎng)格(X,Y)上函數(shù)值、繪網(wǎng)面第24頁，課件共40頁，創(chuàng)作于2023年2月創(chuàng)建網(wǎng)格矩陣原理和方法x=1:6;y=1:8;%創(chuàng)建兩個向量[X,Y]=meshgrid(x,y)%將x和y分別擴充為8行6列X=123456123456123456123456123456123456123456123456Y=111111222222333333444444555555666666777777888888[X,Y]=meshgrid(1:6,1:8)%直接創(chuàng)建兩個矩陣X和Y第25頁，課件共40頁，創(chuàng)作于2023年2月第26頁，課件共40頁，創(chuàng)作于2023年2月例1.12

繪二元函數(shù)

z=xexp(–x2–y2)的圖形。

[x,y]=meshgrid(-2:0.2:2);z=x.*exp(-x.^2-y.^2);mesh(x,y,z)colormap([001])繪網(wǎng)面命令mesh()使用格式:mesh(x,y,z)

或

mesh(z)注記:x,y是維數(shù)相同的矩陣;表達式中“.*”和“.^”運算使得z是與x,y維數(shù)相同矩陣。功能相同的繪圖命令:surf,meshc,meshz,········第27頁，課件共40頁，創(chuàng)作于2023年2月例1.13

繪制一元函數(shù)y=sinx/x在[-8，8]上圖形。例1.14

繪二元函數(shù)圖形[x,y]=meshgrid(-8:0.5:8);r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;mesh(x,y,z)colormap([1,0,0])x=-8:8;y=sin(x)./x;Warning:Dividebyzero.plot(x,y)除零錯誤導致殘缺圖形分母加eps方法避免出錯第28頁，課件共40頁，創(chuàng)作于2023年2月[X,Y]=meshgrid(-2:0.2:2);Z=X.*exp(-X.^2-Y.^2);figure(1),contour(Z)%創(chuàng)建1號圖形窗口figure(2),contourf(Z,20)%創(chuàng)建2號圖形窗口例1.15

用contour()

命令繪二元函數(shù)z=xexp(–x2–y2)

等高線第29頁，課件共40頁，創(chuàng)作于2023年2月圓域上的復變函數(shù)圖形以復變函數(shù)的實部為二元函數(shù)繪圖例1.16

復變量滿足r=linspace(0,1,20);theta=linspace(-2*pi,2*pi,50);z=r'*exp(i*theta);u=r'.^(1/2)*exp(i*theta/2);x=real(z);y=imag(z);s=real(u);mesh(x,y,s)colormap([001])axisoffview(-74,0)第30頁，課件共40頁，創(chuàng)作于2023年2月r=[123];theta=linspace(0,2*pi,7);z=r’*exp(i*theta);x=real(z);y=imag(z);u=real(z.^2);mesh(x,y,u)第31頁，課件共40頁，創(chuàng)作于2023年2月復變函數(shù)圖形設(shè)有復平面上單位圓域內(nèi)變化的變量以u=z

的實部函數(shù)繪圖并輸出圖形文件r=linspace(0,1,20);theta=linspace(-pi,pi,25);z=r'*exp(i*theta);x=real(z);y=imag(z);mesh(x,y,x),holdoncolormap([001])mesh(x,y,-ones(size(x)))axisoff第32頁，課件共40頁，創(chuàng)作于2023年2月輸出圖形文件1.在圖形窗口用鼠標點擊窗口左上方的菜單欄“file”，選擇下拉菜單中的“Export”.(在高版本中選擇“SaveAs”)2.在對話框中選擇文件類型“bmp”格式，將圖形文件命名為riman.第33頁，課件共40頁，創(chuàng)作于2023年2月計算機浮點數(shù)表示法6.378137×103

6.378137e+0032.99792458×10-5

2.99792458e-005

例1.17

近似計算地球表面積,并以不同格式輸出數(shù)據(jù)

R=6378.137;S=4*pi*R^2formatlong,Sformatbank,Sformatshort短格式5.1121e+008formatlong

長格式5.112078933958109e+008formatbank銀行格式

511207893.40————數(shù)據(jù)顯示格式————第34頁，課件共40頁，創(chuàng)作于2023年2月MATLAB內(nèi)部常數(shù)常數(shù)返回值ans默認變量名，保存最近的結(jié)果eps浮點數(shù)相對精度(2.2204e-016)realmax最大浮點數(shù)(1.7977e+308)realmin最小浮點數(shù)(2.2251e-308)pi圓周率（3.1416）i,j虛數(shù)單位inf無限大NaN不合法的數(shù)值，非數(shù)值————數(shù)據(jù)顯示格式————第35頁，課件共40頁，創(chuàng)作于2023年2月

例1.18

立方倍積問題:求作一立方體,使其體積為已知立方體的二倍.如果已知立方體體積為V0=1,要作的立方體體積為V=2,則所求立方體高度為21/3在命令窗口中直接使用命令

h=vpa(2^(1/3),40)數(shù)據(jù)結(jié)果為

h=1.259921049894873190666544360283296555281MATLAB具有一種稱為VPA類型的數(shù)據(jù)顯示格式.VPA的全稱是可變精度浮點算法（Variableprecisionarithmetic）.顯示S的D位數(shù)方法如下R=vpa(S，D)第36頁，課件共40頁，創(chuàng)作于2023年2月字符串變量：將字符串用單引號括起來賦值給變量拼接字符串命令:strcat(T1,T2,T3)將整數(shù)轉(zhuǎn)換成字符命令:int2str()

例1.19

三個名人Euler，Elizabeth，Plato職業(yè)分別是mathematician，moviestar，philosopher