Matlab入門教程很齊全

1、Matlab入門教程很齊全入門教程很齊全 Matlab入入門門教程很教程很齊齊全全 MatlabMatlab簡介簡介 數(shù)組和矩陣數(shù)組和矩陣 MatlabMatlab繪圖繪圖 MatlabMatlab在在微積分微積分中的應(yīng)用中的應(yīng)用 Matlab入門教程很齊全入門教程很齊全 三個代表性的計算機(jī)數(shù)學(xué)語言：三個代表性的計算機(jī)數(shù)學(xué)語言： MATLABMATLAB是是MATrix LABoratory MATrix LABoratory 的縮寫，由的縮寫，由MathWorksMathWorks公司推公司推 出。長于數(shù)值計算，編程方便。在各個領(lǐng)域都有領(lǐng)域?qū)＜揖帉懗?。長于數(shù)值計算，編程方便。在各個領(lǐng)域都有領(lǐng)

2、域?qū)＜揖帉?的工具箱，能高效、可靠地解決各種問題。的工具箱，能高效、可靠地解決各種問題。 MathematicaMathematica，Wolfram ResearchWolfram Research公司公司 MapleMaple，Waterloo MapleWaterloo Maple公司公司 有強(qiáng)大的解析運算和數(shù)學(xué)公式推導(dǎo)、定理證明能力，數(shù)值計有強(qiáng)大的解析運算和數(shù)學(xué)公式推導(dǎo)、定理證明能力，數(shù)值計 算能力比算能力比MatlabMatlab弱，更適合純數(shù)學(xué)求解。弱，更適合純數(shù)學(xué)求解。 Matlab入門教程很齊全入門教程很齊全 MATLAB的功能 MATLABMATLAB產(chǎn)品組是從支持概念設(shè)計、

3、算法開發(fā)、建模仿真，產(chǎn)品組是從支持概念設(shè)計、算法開發(fā)、建模仿真， 到實時實現(xiàn)的集成環(huán)境，可用來進(jìn)行：到實時實現(xiàn)的集成環(huán)境，可用來進(jìn)行： 數(shù)據(jù)分析數(shù)據(jù)分析 數(shù)值與符號計算數(shù)值與符號計算 工程與科學(xué)繪圖工程與科學(xué)繪圖 控制系統(tǒng)設(shè)計控制系統(tǒng)設(shè)計 數(shù)字圖像信號處理數(shù)字圖像信號處理 建模、仿真、原型開發(fā)建模、仿真、原型開發(fā) 財務(wù)工程、應(yīng)用開發(fā)、圖形用戶界面設(shè)計財務(wù)工程、應(yīng)用開發(fā)、圖形用戶界面設(shè)計 Matlab入門教程很齊全入門教程很齊全 MATLAB語言特點 編程效率高，允許用數(shù)學(xué)的語言來編寫程序編程效率高，允許用數(shù)學(xué)的語言來編寫程序 用戶使用方便，把程序的編輯、編譯、連接和執(zhí)行融為一體用戶使用方便，把

4、程序的編輯、編譯、連接和執(zhí)行融為一體 高效方便的矩陣和數(shù)組運算高效方便的矩陣和數(shù)組運算 語句簡單，內(nèi)涵豐富語句簡單，內(nèi)涵豐富 擴(kuò)充能力強(qiáng)，交互性，開放性擴(kuò)充能力強(qiáng)，交互性，開放性 方便的繪圖功能方便的繪圖功能 該軟件由該軟件由c c語言編寫，移植性好語言編寫，移植性好 Matlab入門教程很齊全入門教程很齊全 l MATLAB操作窗口操作窗口 雙擊桌面快捷鍵，啟動軟件。雙擊桌面快捷鍵，啟動軟件。 接受命令的窗口接受命令的窗口 Matlab入門教程很齊全入門教程很齊全 MATLAB的環(huán)境 菜單項；菜單項； 工具欄；工具欄； 【Command WindowCommand Window】命令窗口；命

5、令窗口； 【Launch PadLaunch Pad】分類幫助窗口；分類幫助窗口； 【W(wǎng)orkspaceWorkspace】工作區(qū)窗口；工作區(qū)窗口； 【Command HistoryCommand History】指令歷史記錄窗口；指令歷史記錄窗口； 【Current DirectoryCurrent Directory】當(dāng)前目錄選擇窗口；當(dāng)前目錄選擇窗口； Matlab入門教程很齊全入門教程很齊全 MatlabMatlab簡介簡介 數(shù)組和矩陣數(shù)組和矩陣 MatlabMatlab繪圖繪圖 MatlabMatlab在在微積分微積分中的應(yīng)用中的應(yīng)用 Matlab入門教程很齊全入門教程很齊全 1.1

6、.變量變量 MatlabMatlab不需要任何類型聲明和維數(shù)說明不需要任何類型聲明和維數(shù)說明, ,變量名的第一個字變量名的第一個字 符必須是字母。符必須是字母。 a=1a=1； num_students=25;num_students=25; 2.2.常用的常數(shù)常用的常數(shù) pipi：3.141592653.14159265 i,j: i,j:虛數(shù)單位；虛數(shù)單位；1i; 3-4j; 3e5i1i; 3-4j; 3e5i Inf Inf 無限值；無限值；NaNNaN空值空值 e e 以以1010為底的冪次。為底的冪次。1.602e-20, 6.532e121.602e-20, 6.532e12 M

7、atlab入門教程很齊全入門教程很齊全 3.3.運算符運算符 + +，- - ，* *，/ /， 左除左除: :23 = 1.500023 = 1.5000 冪冪: x=2; x3; x(-3): x=2; x3; x(-3) 復(fù)數(shù)共軛轉(zhuǎn)置復(fù)數(shù)共軛轉(zhuǎn)置 x=3+4i xx=3+4i x=3-4i=3-4i . .點運算點運算 (1)(1)當(dāng)當(dāng)x x是一個向量時，求是一個向量時，求 不能寫成不能寫成x2x2，而必須寫成，而必須寫成 x.2x.2 (2) (2)兩矩陣之間的點乘運算兩矩陣之間的點乘運算C=A.C=A.* *B B ijijij bac 2 i x Matlab入門教程很齊全入門教程

8、很齊全 . .常用函數(shù)常用函數(shù) sqrt(x)sqrt(x) 開平方開平方 abs(x) abs(x) 絕對值絕對值abs(3-4i)abs(3-4i) exp(x) eexp(x) ex x; log(x) ; log(x) 以以e e為底，為底，x x的對數(shù)的對數(shù) log(exp(2)log(exp(2) round(x) round(x)取整取整 syms x; syms x; 定義定義x x為符號變量為符號變量 . .幫助函數(shù)幫助函數(shù) help: help elfun; help specfun; help elmathelp: help elfun; help specfun; he

9、lp elmat Matlab入入門門教程很教程很齊齊全全 6.6.構(gòu)造數(shù)組構(gòu)造數(shù)組 (1)(1)直接構(gòu)造，用空格或逗號間隔數(shù)組元素直接構(gòu)造，用空格或逗號間隔數(shù)組元素 A = 2 3 5 1 A = 2 3 5 1 或或 A = sqrt(2),3e2,log(5),1+2i A = sqrt(2),3e2,log(5),1+2i (2) (2)用增量法構(gòu)造數(shù)組用增量法構(gòu)造數(shù)組 (first:last) (first:last) 或或 (first:step:last)(first:step:last) A = 10:15 A = 3:0.2:4 A = 9:-1:0 A = 10:15 A

10、= 3:0.2:4 A = 9:-1:0 (3) (3)用用linspacelinspace函數(shù)構(gòu)造數(shù)組函數(shù)構(gòu)造數(shù)組 x = linspace(first,last,num)x = linspace(first,last,num) x = linspace(0,10,5) x = linspace(0,10,5) Matlab入門教程很齊全入門教程很齊全 7.7.構(gòu)造矩陣構(gòu)造矩陣 (1)(1)簡單創(chuàng)建方法簡單創(chuàng)建方法 row = e1,e2,row = e1,e2,em; A = row1,em; A = row1；row2row2；rownrown A = 2 4 1;4 5 2;7 2 1

11、 A = 2 4 1;4 5 2;7 2 1 (2)(2)構(gòu)造特殊矩陣構(gòu)造特殊矩陣 ones ones 創(chuàng)建一個所有元素都為的元素創(chuàng)建一個所有元素都為的元素 zeros zeros 創(chuàng)建一個所有元素都為創(chuàng)建一個所有元素都為0 0的元素的元素 127 254 142 A Matlab入門教程很齊全入門教程很齊全 eye eye 創(chuàng)建對角元素為，其他元素為的元素創(chuàng)建對角元素為，其他元素為的元素 rand rand 創(chuàng)建一個矩陣或數(shù)組，其中的元素服從均勻分布創(chuàng)建一個矩陣或數(shù)組，其中的元素服從均勻分布 rand(5)rand(5)* *20;20; randnrandn創(chuàng)建一個矩陣或數(shù)組，其中的元素服

12、從正態(tài)分布創(chuàng)建一個矩陣或數(shù)組，其中的元素服從正態(tài)分布 diag diag 創(chuàng)建對角矩陣創(chuàng)建對角矩陣 C = 3 2 1; V=diag(C);C = 3 2 1; V=diag(C); (3)(3)聚合矩陣聚合矩陣 水平聚合水平聚合 C = A BC = A B 垂直聚合垂直聚合 C = A;BC = A;B n n diag 2 1 21 ),( Matlab入門教程很齊全入門教程很齊全 8.8.獲取矩陣元素獲取矩陣元素 A=2,3,3;4 9 4;6,3,0 A=2,3,3;4 9 4;6,3,0 取單個元素：取單個元素： 取多個元素：取多個元素： 獲取所有元素：獲取所有元素： 9.9.獲

13、取與矩陣有關(guān)信息獲取與矩陣有關(guān)信息 length length 返回最長維長度返回最長維長度 ndims ndims 返回維數(shù)返回維數(shù) numel numel 返回元素個數(shù)返回元素個數(shù) size size 返回每一維的長度返回每一維的長度 rows cols = size(A)rows cols = size(A) A(3,1) A(:,2) A(3,:) A(:) Matlab入門教程很齊全入門教程很齊全 數(shù)組和矩陣數(shù)組和矩陣 9.矩陣的基本運算矩陣的基本運算 例例 已知已知 422134 305 ,203 153211 AB a=4 -2 2;-3 0 5;1 5 3; b=1 3 4;-

14、2 0 -3;2 -1 1; a*b 12 10 24 7 -14 -7 -3 0 -8 ans =AB Matlab入門教程很齊全入門教程很齊全 ( )R A rank(a) ans = 3 數(shù)組和矩陣數(shù)組和矩陣 9.矩陣的基本運算矩陣的基本運算 例例 已知已知 422134 305 ,203 153211 AB Matlab入門教程很齊全入門教程很齊全 det(a) ans = -158 數(shù)組和矩陣數(shù)組和矩陣 9.矩陣的基本運算矩陣的基本運算 例例 已知已知 422134 305 ,203 153211 AB Matlab入門教程很齊全入門教程很齊全 求特征多項式 poly(A) A=sy

15、m(A); 將A轉(zhuǎn)換成符號矩陣 poly(A) 數(shù)組和矩陣數(shù)組和矩陣 9.矩陣的基本運算矩陣的基本運算 例例 已知已知 422134 305 ,203 153211 AB Matlab入門教程很齊全入門教程很齊全 1 AB A/B相當(dāng)于矩陣方程XB=A ans = 0 0 2.0000 -2.7143 -8.0000 -8.1429 2.4286 3.0000 2.2857 數(shù)組和矩陣數(shù)組和矩陣 9.矩陣的基本運算矩陣的基本運算 例例 已知已知 422134 305 ,203 153211 AB Matlab入門教程很齊全入門教程很齊全 1 A B AB ans = 0.4873 0.4114

16、 1.0000 0.3671 -0.4304 0 -0.1076 0.2468 0 數(shù)組和矩陣數(shù)組和矩陣 9.矩陣的基本運算矩陣的基本運算 例例 已知已知 422134 305 ,203 153211 AB Matlab入門教程很齊全入門教程很齊全 10.多項式求根多項式求根 例例 已知已知 p = 1 0 -2 -5; roots(p) ans = 2.0946 -1.0473 + 1.1359i -1.0473 - 1.1359i 52)( 3 xxxp 數(shù)組和矩陣數(shù)組和矩陣 Matlab入門教程很齊全入門教程很齊全 MatlabMatlab簡介簡介 數(shù)組和矩陣數(shù)組和矩陣 MatlabMa

17、tlab繪圖繪圖 MatlabMatlab在在微積分微積分中的應(yīng)用中的應(yīng)用 Matlab入門教程很齊全入門教程很齊全 1.1.二維圖形繪制二維圖形繪制 plot(t,y)plot(t,y) 例例1 1用用MatlabMatlab畫出的圖形。畫出的圖形。 x=-5:0.05:5; y=x.2; plot(x,y)plot(x,y,b+) 例例2 2繪制繪制y=sin(tan(x)-tan(sin(x)在在 區(qū)間內(nèi)的曲線。區(qū)間內(nèi)的曲線。 x=-pi:0.05:pi;y=sin(tan(x)-tan(sin(x); plot(x,y) 2 xy , Matlab入門教程很齊全入門教程很齊全 將多條曲

18、線畫在一個圖上：將多條曲線畫在一個圖上： plot(t1,y1,選項選項1,t2,y2,選項選項2,) plot(x,x.2,rO,x,x.3,b.) 2.三維圖形繪制三維圖形繪制 plot3(x,y,z) plot3(x1,y1,z1,選項選項1,x2,y2,z2,選項選項2,) 例例 試?yán)L制參數(shù)方程試?yán)L制參數(shù)方程 ， ， 的三維曲線。的三維曲線。 t=0:.1:2*pi; 注意點運算注意點運算 x=t.3.*sin(3*t).*exp(-t);y=t.3.*cos(3*t).*exp(-t); z=t.2; plot3(x,y,z),grid t etttx )3sin()( 3t ett

19、ty )3cos()( 3 2 tz Matlab入門教程很齊全入門教程很齊全 15.利用利用Mathematica作出數(shù)列作出數(shù)列 1 ) 1 1 ( n n n x的點圖，的點圖， 觀察當(dāng)觀察當(dāng) n時，時， n x的變化趨勢。并利用數(shù)值計算的命令計算當(dāng)?shù)淖兓厔?。并利用?shù)值計算的命令計算當(dāng)n 取很大的整數(shù)時，取很大的整數(shù)時， n x的取值。的取值。 n=1:10000; xn=(1+1./n).(n+1); plot(n,xn,.) Matlab入門教程很齊全入門教程很齊全 16.函數(shù)函數(shù)xxycos 在在內(nèi)是否有界？又問當(dāng)內(nèi)是否有界？又問當(dāng)),( 時，時，x這個函數(shù)是否為無窮大？為什么？

20、用這個函數(shù)是否為無窮大？為什么？用 Mathematica作圖并驗證你的結(jié)論。作圖并驗證你的結(jié)論。 x = -100:100; plot(x,x.*cos(x) Matlab入門教程很齊全入門教程很齊全 P168 20.利用利用Mathematica作出函數(shù)作出函數(shù) )45( 2 1 )( 2 x cxx xf 的圖形，分別取的圖形，分別取-1，0，1，2，3等等5個值，試比較作出的個值，試比較作出的5個個 圖，并從圖上觀察極值點、駐點，增加、減少區(qū)間，上凸、圖，并從圖上觀察極值點、駐點，增加、減少區(qū)間，上凸、 下凸區(qū)間以及漸近線。下凸區(qū)間以及漸近線。 x=-5:0.1:4 plot(x,1.

21、/(x.2+2*x-1) plot(x,1./(x.2+2*x) plot(x,1./(x.2+2*x+1) plot(x,1./(x.2+2*x+2) plot(x,1./(x.2+2*x+3) Matlab入門教程很齊全入門教程很齊全 17（1）在計算機(jī)屏幕上作出函數(shù)）在計算機(jī)屏幕上作出函數(shù) 1 . 0 )(xxf 和和 xxgln)( 的圖形，何時開始的圖形，何時開始 ?gf （2）再作出函數(shù)）再作出函數(shù))(/ )()(xfxgxh的圖形。選用的圖形。選用 適當(dāng)?shù)娘@示區(qū)域，展示適當(dāng)?shù)娘@示區(qū)域，展示 x 時，時， )(xh 的變化趨勢。的變化趨勢。 （3）確定正數(shù)）確定正數(shù)X，使當(dāng)，使當(dāng)X

22、x 時，時， ?1 . 0 )( )( xf xg Matlab入門教程很齊全入門教程很齊全 x=linspace(3.43063112146e15-1e6,3.43063112146e15 +1e6,100); plot(x,x.0.1,r+,x,log(x),b.) axis(3.43063112146e15-1e3 3.43063112146e15+1e3 35.77152063979 35.7715206398); x=1:0.05:5 plot(x,x.0.1,r.,x,log(x),b+) axis(XMIN XMAX YMIN YMAX) x=linspace(1e29,1e29

23、+1e16,1000); plot(x,log(x)./(x.0.1),r.) axis(1e29-1e20 1e29+1e20 0.08 0.12) Matlab入門教程很齊全入門教程很齊全 MatlabMatlab簡介簡介 數(shù)組和矩陣數(shù)組和矩陣 MatlabMatlab繪圖繪圖 MatlabMatlab在在微積分微積分中的應(yīng)用中的應(yīng)用 Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 1、求函數(shù)值、求函數(shù)值 例例1 在命令窗口中鍵入表達(dá)式在命令窗口中鍵入表達(dá)式 并求并求 時的函數(shù)值。時的函數(shù)值。 2 ln3, x y zxeyx 2,4xy x=2,y

24、=4 z=x2+exp(x+y)-y*log(x)-3 x = 2 y = 4 z = 401.6562 命令窗口顯示結(jié)果：命令窗口顯示結(jié)果： Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 2、求極限、求極限 極限問題：極限問題： 求單側(cè)極限：求單側(cè)極限： ),(limit 0 xxfunL ) ,(limit 0 rightorleftxxfunL )(lim 0 xfL xx )(lim 0 1 xfL xx )(lim 0 2 xfL xx Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 2、求極限、求極限

25、 例例2 求極限求極限 syms x; limit(sin(x)/x,x,0) ans = 1 定義符號變量定義符號變量 x x x sin lim 0 Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 2、求極限、求極限 例例3 求極限求極限 syms x; limit(exp(x3)-1)/(1-cos(sqrt(x-sin(x),x,0, right) ans = 12 定義符號變量定義符號變量 xx e x x sincos1 1 lim 3 0 Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 2、求極限、求

26、極限 例例4 求極限求極限 lim() n nnn syms n; limit(sqrt(n+sqrt(n)-sqrt(n),n,inf) ans = 1/2 定義符號變量定義符號變量 Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 3、求導(dǎo)數(shù)、求導(dǎo)數(shù) 調(diào)用格式：調(diào)用格式： ),(xfundiffy ),(nxfundiffy 求導(dǎo)數(shù)求導(dǎo)數(shù) 求求n階導(dǎo)數(shù)階導(dǎo)數(shù) 多元函數(shù)的偏導(dǎo)數(shù)：多元函數(shù)的偏導(dǎo)數(shù)： ),),(nymxfdiffdifff ),),(mxnyfdiffdifff 或或 Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中

27、的應(yīng)用中的應(yīng)用 3、求導(dǎo)數(shù)、求導(dǎo)數(shù) syms x diff(sin(x.3),x) ans = 3*cos(x3)*x2 定義定義X為符號變量為符號變量 求求 dy dx 習(xí)題習(xí)題P168 (1)sin( 3 xy Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 3、求導(dǎo)數(shù)、求導(dǎo)數(shù) syms x diff(atan(log(x),x) (2)arctan(lnxy ans = 1/x/(1+log(x)2) Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 3、求導(dǎo)數(shù)、求導(dǎo)數(shù) syms x diff(1+1/x)x,

28、x) ans = (1+1/x)x*(log(1+1/x)-1/x/(1+1/x) (3) x x y) 1 1 ( Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 3、求導(dǎo)數(shù)、求導(dǎo)數(shù) syms x (4)(2 2 xxfy Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 3、求導(dǎo)數(shù)、求導(dǎo)數(shù) 例例6 設(shè)設(shè) 10 10ln x yxx，求，求 y syms x y=10 x+x10+log(x) y = x10+10 x+log(x) diff(y) ans = 10*x9+10 x*log(10)+1/x 定義定義

29、X為符號變量為符號變量 求求 dy dx Matlab入門教程很齊全入門教程很齊全 例例7 設(shè)設(shè) ln(1),yx求求 2 1 2 x d y dx syms x; y=log(1+x); a=diff(y,x,2) a = -1/(1+x)2 x=1;eval(a) ans = -0.2500 求求 2 2 d y dx 求求 2 1 2 x d y dx 將符號表達(dá)式將符號表達(dá)式 轉(zhuǎn)換成數(shù)值表達(dá)式轉(zhuǎn)換成數(shù)值表達(dá)式 Matlab入門教程很齊全入門教程很齊全 例例6 設(shè)設(shè) 22 2 x zexyy，求，求 222 22 , zzzzz xyxyx y syms x y; z=exp(2*x)*

30、(x+y2+2*y); a=diff(z,x) b=diff(z,y) c=diff(z,x,2) d=diff(z,y,2) e=diff(a,y) z a x z b y 2 2 z c x 2 2 z d y 2 az e yx y Matlab入門教程很齊全入門教程很齊全 a =2*exp(2*x)*(x+y2+2*y)+exp(2*x) b =exp(2*x)*(2*y+2) c =4*exp(2*x)*(x+y2+2*y)+4*exp(2*x) d =2*exp(2*x) e =2*exp(2*x)*(2*y+2) 222 222 xx z aexyye x 2 22 x z be

31、y y 2 222 2 424 xx z cexyye x 2 2 2 2 x z de y 2 2 222 x z eey x y Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 4、求極值和零點、求極值和零點 例例7 已知已知 5432 ( )323f xxxxx ，求，求 （1）函數(shù)的零點；（）函數(shù)的零點；（2）函數(shù)在）函數(shù)在-1，2上的最小值上的最小值 fzero(3*x5-x4+2*x3+x2+3,0) ans = -0.8952 起始搜索點起始搜索點 函數(shù)函數(shù) 命令函數(shù)命令函數(shù) x,f=fminbnd(3*x5-x4+2*x3+x2+3,-1

32、,2) x = -1 f=-2 Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 4、求極值和零點、求極值和零點 ，求，求 例例7 已知已知 222 ( , , )2.5sinf x y zxyxy z 函數(shù)在點（函數(shù)在點（1，-1，0）附近的最小值）附近的最小值 X,FVAL= FMINSEARCH(x(1)2+2.5*sin(x(2)- x(3)*x(1)*x(2)2,1 -1 0) X = 0.0010 -1.5708 0.0008 FVAL =-2.5000 Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 5

33、、求積分、求積分 （）不定積分（）不定積分 ),int(xfunF （）定積分與無窮積分（）定積分與無窮積分 ),int(baxfI (a,b)為定積分的積分區(qū)間，求解無窮積分時為定積分的積分區(qū)間，求解無窮積分時 允許將允許將a，b設(shè)置成設(shè)置成-Inf或或Inf。 Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 5、求積分、求積分 例例8 求不定積分求不定積分 cos2 cos3xxdx int(cos(2*x)*cos(3*x) ans =1/2*sin(x)+1/10*sin(5*x) 例例9 求定積分求定積分 2 1 ln e xxdx Integ

34、rate：積分：積分 int(x2*log(x),x,1,exp(1) eval(ans) ans = 4.5746 Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 5、求積分、求積分 例例10 求二重積分求二重積分 2 2 1 ,2,12 2 y dxdyxy x syms x y; f=y2/x2; int(int(f,x,1/2,2),y,1,2) ans =7/2 Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 5、求積分、求積分 P251 18.(1) xdx 10 sin syms x int(sin(

35、x)10,x) ans = -1/10*sin(x)9*cos(x)-9/80*sin(x)7*cos(x)- 21/160*sin(x)5*cos(x)-21/128*sin(x)3*cos(x)- 63/256*cos(x)*sin(x)+63/256*x Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 5、求積分、求積分 P251 18.(2) dx xx 3 42 ) 1() 1( 1 int(1/(x+1)2*(x-1)4)(1/3),x) ans = -3/2*(x+1)/(x+1)2*(x-1)(1/3) Matlab入門教程很齊全入門教程

36、很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 5、求積分、求積分 P251 18.(3) 1 0 6 sindxx int(sin(x(1/6),x,0,1) ans = -606*cos(1)+390*sin(1) Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 5、求積分、求積分 P251 18.(4) 1 02 1 )1ln( dx x x int(log(1+x)/(1+x2),x,0,1) ans = 1/4*pi*log(2)-Catalan-1/2*i*dilog(1/2- 1/2*i)+1/2*i*dilog(1/2+1/2*i)

37、 vpa(ans) quadl(log(1+x)./(1+x.2),0,1) ans = 0.2722 變步長數(shù)值積分變步長數(shù)值積分 Matlab入門教程很齊全入門教程很齊全 例例11 求定積分求定積分 2 1 2 0 x edx int(exp(-x2/2),0,1) ans = 1/2*erf(1/2*2(1/2)*2(1/2)*pi(1/2) 2 0 2 ( ) x t erf xedt 22 2 0 2 2 t ansedt x=0:0.01:1; y=exp(-x.2/2); trapz(x,y) ans = 0.8556 y=exp(-x.2/2); quadl(y,0,1) an

38、s = 0.8556 變步長數(shù)值積分變步長數(shù)值積分 梯形法數(shù)值積分梯形法數(shù)值積分 Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 5、求積分、求積分 例例12 2 000 22 4 dzdydxxze zyx int(int(int(4*x*z*exp(-x2*y-z2),z,0,pi),y,0,pi),x,0,2) ans = (-Ei(1,4*pi)+Ei(1,4*pi)*exp(pi2)-2*log(2)+2*log(2) *exp(pi2)-eulergamma-log(pi)+eulergamma*exp(pi2) +log(pi)*exp(p

39、i2)/exp(pi2) vpa(ans,60) ans = 3.14127228346146476723483587 1 ),(dtteznEi nzt Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 6、解微分方程、解微分方程 ),( 21m fffdsolvey ) ,( 21 xfffdsolvey m 指明自變量指明自變量 i f既可以描述微分方程，又可以描述初始條件或邊界條件。既可以描述微分方程，又可以描述初始條件或邊界條件。 )( )4( ty:D4y 3)2( y :D2y(2)=3 Matlab入門教程很齊全入門教程很齊全 MATLAB

40、在在微積分微積分中的應(yīng)用中的應(yīng)用 6、解微分方程、解微分方程 例例13 計算初值問題：計算初值問題： 1)0(y xy dx dy dsolve(Dy=x+y,y(0)=1,x) ans =-x-1+2*exp(x) 一定要大寫一定要大寫 Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 6、解微分方程、解微分方程 例例14 求方程求方程 dsolve(Dy+2*y/x=sin(3*x)/(x2),x) ans = (-1/3*cos(3*x)+C1)/x2 2 3sin2 x x y xdx dy 的通解。的通解。 Matlab入門教程很齊全入門教程很齊

41、全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 6、解微分方程、解微分方程 例例15 設(shè)輸入信號為設(shè)輸入信號為5) 12cos()( 5 tetu t ，求，求 )(2)(4)(5 )(24)(50)(35)(10)( )3()4( tututu tytytytyty syms t u=exp(-5*t)*cos(2*t+1)+5; uu=5*diff(u,t,2)+4*diff(u,t)+2*u uu = 87*exp(-5*t)*cos(2*t+1)+92*exp(-5*t)*sin(2*t+1)+10 Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)

42、用 6、解微分方程、解微分方程 例例14 設(shè)輸入信號為設(shè)輸入信號為5) 12cos()( 5 tetu t ，求，求 )(2)(4)(5 )(24)(50)(35)(10)( )3()4( tututu tytytytyty syms t y; y=dsolve(D4y+10*D3y+35*D2y+50*Dy+24*y=87*exp(- 5*t)*cos(2*t+1)+92*exp(-5*t)*sin(2*t+1)+10) y = -547/520*exp(-5*t)*sin(2*t+1)-343/520*exp(- 5*t)*cos(2*t+1)+5/12+C1*exp(-2*t)+C2*e

43、xp(- 3*t)+C3*exp(-t)+C4*exp(-4*t) Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 6、解微分方程、解微分方程 若已知若已知0)0()0(, 2)0(, 3)0( )3( yyyy y=dsolve(D4y+10*D3y+35*D2y+50*Dy+24*y=,87*exp(- 5*t)*cos(2*t+1)+92*exp(- 5*t)*sin(2*t+1)+10,y(0)=3,Dy(0)=2,D2y(0)=0,D3y(0)=0) y = -547/520*exp(-5*t)*sin(2*t+1)-343/520*exp(-

44、 5*t)*cos(2*t+1)+5/12+(-51/13*sin(1)-445/26*cos(1)- 69/2)*exp(-2*t)+(5/8*sin(1)+179/8*cos(1)+73/3)*exp(- 3*t)+(97/60*sin(1)+133/30*cos(1)+19)*exp(- t)+(41/15*sin(1)-271/30*cos(1)-25/4)*exp(-4*t) Matlab入門教程很齊全入門教程很齊全 MATLAB在在微積分微積分中的應(yīng)用中的應(yīng)用 6、解微分方程、解微分方程 例例16 t t etytxty etytxtxtx 4)(3)(4)( )(2)()(2)( x,y=dsolve(D2x+2*Dx=x+2*y-exp(- t),Dy=4*x+3*y+4*exp(-t) x = -6*t*exp(-t)+C1*exp(-t)+C2*exp(1+6(1/2)*t)+C3*exp(-(- 1+6(1/2)*t) y = 6*t*exp(-t)-C1*exp(- t)+4*C2*exp(1+6(1/2)*t)+2*C2*exp