第一篇：北航matlab上機(jī)實(shí)驗(yàn)報(bào)告
MATLAB基礎(chǔ)上機(jī)指導(dǎo)書
實(shí)驗(yàn)?zāi)康?/p>
MATLAB基礎(chǔ)上機(jī)實(shí)驗(yàn)的目的主要是在深入理解和掌握課程中的基本概念和原理的基礎(chǔ)上培養(yǎng)學(xué)習(xí)者用計(jì)算機(jī)解決實(shí)際問題的意識(shí)，訓(xùn)練和提高學(xué)生分析問題、仿真、程序設(shè)計(jì)和程序調(diào)試能力。實(shí)驗(yàn)要求
1．整個(gè)學(xué)期的上機(jī)實(shí)踐共有4個(gè)。從實(shí)驗(yàn)一的MATLAB的上機(jī)基本操作，到實(shí)驗(yàn)四的GUI設(shè)計(jì)，都本著由淺入深、由易到難、逐步深入的原則，給出了不同的上機(jī)程序。一般每個(gè)實(shí)驗(yàn)中有約5-7個(gè)不同難度的程序，按班級(jí)序號(hào)排列組合分組，每組學(xué)生包括5人，按組檢查程序。
2．每次實(shí)驗(yàn)后要求學(xué)生提交上機(jī)報(bào)告，這也是整個(gè)實(shí)驗(yàn)過程的一個(gè)重要環(huán)節(jié)。實(shí)驗(yàn)結(jié)束后，要整理實(shí)驗(yàn)結(jié)果并認(rèn)真分析和總結(jié)，根據(jù)教師要求寫出實(shí)驗(yàn)報(bào)告。實(shí)驗(yàn)報(bào)告一般包括如下內(nèi)容：
⑴ 實(shí)驗(yàn)內(nèi)容：實(shí)驗(yàn)題目與要求。⑵ 程序清單
(3)運(yùn)行結(jié)果：原始數(shù)據(jù)、相應(yīng)的運(yùn)行結(jié)果和必要的說明。
(4)分析與思考：調(diào)試過程及調(diào)試中遇到的問題及解決辦法；調(diào)試程序的心得與體會(huì)；其他函數(shù)算法的存在與實(shí)踐等。若最終未完成調(diào)試，要認(rèn)真找出錯(cuò)誤并分析原因等。常用命令參考
本節(jié)就一些常用的MATLAB命令進(jìn)行介紹，使初學(xué)者盡快提高M(jìn)ATLAB上機(jī)調(diào)試程序的能力。
1、學(xué)會(huì)用help和doc函數(shù)。
2、輸入輸出文件：save/load
3、解線性方程組AX=B：X=AB
4、作圖時(shí)兩張曲線合并：hold on或者subplot作子圖
5、程序計(jì)算時(shí)間：tic，toc或者clock
6、變量顯示方式更改：format long/short/bank...7、數(shù)組元素求和：sum
8、求數(shù)組長度：length 求矩陣維數(shù)：size或者ndims 矩陣元素個(gè)數(shù)：numel
9、函數(shù)作圖：
餅圖：pie/pie3 誤差圖：errorbar 散點(diǎn)圖：scatter/scatter3 直方圖：hist 函數(shù)圖：fplot 動(dòng)畫：movie
10、矩陣分析：
左右翻轉(zhuǎn)：fliplr 上下翻轉(zhuǎn)：flipud 轉(zhuǎn)置：transpose 矩陣求逆：inv 矩陣范數(shù)：norm 條件數(shù)：cond 初等變換：rref 特征值：eig/eigs
11、特殊矩陣：
元素全為1的矩陣：ones 元素全為0的矩陣：zeros 單位陣：eye 魔方陣：magic 線性變化數(shù)組：linspace 聚合矩陣：cat/horzcat/vertcat
12、隨機(jī)數(shù)：
創(chuàng)建一個(gè)元素服從均勻分布的隨機(jī)數(shù)數(shù)組：rand 創(chuàng)建一個(gè)元素服從正態(tài)分布的隨機(jī)數(shù)數(shù)組：randn 二項(xiàng)分布：binornd 指數(shù)分布：exprnd F分布：frnd 幾何分布：geornd 超幾何分布：hygernd 泊松分布：poissrnd
正態(tài)分布：normrnd 離散均勻分布：unidrnd 連續(xù)均勻分布：unifrnd
13、清屏：clc 清理內(nèi)存：clear
14、字體顯示變更等：preferences
15、得到一個(gè)文件夾的所有文件名：ls
16、語句太長的話可以再句末加...換行。
實(shí)驗(yàn)一 MATLAB快速入門操作
[實(shí)驗(yàn)?zāi)康腯 1． 熟悉MATLAB的工作環(huán)境。2． 熟悉MATLAB的幫助系統(tǒng)。3． 掌握基本的運(yùn)算表達(dá)。4． 掌握基本的數(shù)學(xué)函數(shù)。5． 掌握基本的字符串使用。[實(shí)驗(yàn)內(nèi)容] 1．執(zhí)行MATLAB的啟動(dòng)和退出。觀看DEMO程序,INTRO程序,HELP文件，練習(xí)MATLAB常用命令在命令窗口中計(jì)算數(shù)學(xué)表達(dá)式，輸入簡單矩陣，查看工作空間中的變量，命令窗口顯示清除，語句、變量及數(shù)據(jù)的表示，搜索路徑設(shè)置。
2．使用函數(shù)創(chuàng)建字符串?dāng)?shù)組，掌握char、strvcat、strcat之間的區(qū)別。>> c1=char('china','india');>> c3=strvcat('china','india');>> c4=strcat('china','india');char中各部分互相獨(dú)立，strvact豎直連接 strcat直接連接 >> c1=char('china','india')
c1 =
china india
>> c3=strvcat('china','india')
c3 =
china india
>> c4=strcat('china','india')
c4 =
chinaindia 也就是說strcat連接兩個(gè)字符串,而strvcat垂直連接兩個(gè)字符串.? i6c?2e3．設(shè)三個(gè)復(fù)數(shù)a＝3＋4i，b＝1＋2i,計(jì)算x=ab/c，注意復(fù)數(shù)的表達(dá)方式。
>> a=3+4i,b=1+2i,c=2*exp(pi/6*i)a =
3.0000 + 4.0000i
b =
1.0000 + 2.0000i
c =
1.7321 + 1.0000i
>> X=a*b/c X =
0.3349 + 5.5801i
sin?x??354．計(jì)算下式的結(jié)果，其中x=45°，572
注意函數(shù)參數(shù)的使用類型，掌握復(fù)雜表達(dá)式的計(jì)算。>> a=45/180*pi a =
0.7854
>> b=(sin(a)+35^(1/2))/72^(1/5)b =
2.8158
5．計(jì)算星球之間的萬有引力。包括地球和太陽，地球和月亮之間的距離。具體參數(shù)太陽質(zhì)量1.989×10^30千克，地球質(zhì)量5.976×10^27千克，太陽和地球的距離***0米，月亮質(zhì)量7.351×10^ 22千克，月亮和地球兩者間距384400000米，萬有引力系數(shù)K 6．67×10^-11，見P26 例2-3。
>> M1=1.989e30,M2=5.976e27,M3=7.351e22,d1=1.4960e11,d2=3.844e8 K=6.67e-11
M1 =
1.9890e+030
M2 =
5.9760e+027
M3 =
7.3510e+022
d1 =
1.4960e+011
d2 =
384400000
K =
6.6700e-011
>> F12=M1*M2*K/d1^2
F12 =
3.5426e+025
>> F23=M2*M3*K/d2^2
F23 =
1.9830e+023
?333??，以及數(shù)據(jù)的修改方法，盡量考2336．用兩種方法創(chuàng)建一個(gè)二維數(shù)組?????433??慮計(jì)算效率。在掌握二維數(shù)組基礎(chǔ)上掌握三維數(shù)組的創(chuàng)建方法。
創(chuàng)建：方法1;x=[3,3,3;2,3,3;4,3,3] x =
3 3 2 3 3 4 3 3 方法2 在workspace中點(diǎn)擊 new variable 就會(huì)在下面出現(xiàn)一個(gè)新變量unnamed 雙擊unnamed變量引出一個(gè)空白頁面，在里面進(jìn)行編輯 >> unnamed
unnamed =
3 3 2 3 3 4 3 3 修改
直接在workspace中點(diǎn)擊需要修改的變量，在相應(yīng)位置修改變量，保存后在命令窗口重新調(diào)用函數(shù) 計(jì)算
[問題討論] 1．小結(jié)上機(jī)各個(gè)環(huán)節(jié)所出現(xiàn)的錯(cuò)誤及解決的辦法。
[思考問題] 1．仿真是否“真實(shí)”，與什么因素有關(guān)? 仿真是在數(shù)字計(jì)算機(jī)上進(jìn)行試驗(yàn)的數(shù)字化技術(shù)，包括數(shù)字和邏輯模型的某些模式，這些模型描述了某一事件或經(jīng)濟(jì)系統(tǒng)在若干周期內(nèi)的特征。仿真常常包括模型建立，仿真模型的建立和仿真試驗(yàn)三個(gè)方面。從以上可以看出，仿真其實(shí)是建立在模型之上，通過計(jì)算機(jī)或者其他載體，并借助經(jīng)驗(yàn)知識(shí)、統(tǒng)計(jì)數(shù)據(jù)和信息資料等對(duì)真實(shí)或假想的系統(tǒng)進(jìn)行試驗(yàn)，得出試驗(yàn)結(jié)果的過程。
在這點(diǎn)上看，仿真并不真實(shí)，因?yàn)樗⒉荒?00%模擬出真實(shí)情況，而實(shí)際上“真實(shí)”的復(fù)雜多變也難以模擬。但是，仿真可以通過以上的手段在一定范圍內(nèi)對(duì)真實(shí)進(jìn)行模擬和預(yù)測(cè)。其與真實(shí)結(jié)果的相關(guān)性，或者說仿真方法的優(yōu)劣性與其系統(tǒng)、模型和載體是密切相關(guān)的。
第二篇：MATLAB上機(jī)實(shí)驗(yàn)實(shí)驗(yàn)報(bào)告
M M ＡＴL L ＡB B 上機(jī)實(shí)驗(yàn)一
一、實(shí)驗(yàn)?zāi)康?初步熟悉 MＡＴLAB 工作環(huán)境，熟悉命令窗口，學(xué)會(huì)使用幫助窗口查找?guī)椭畔ⅰ?/p>
命令窗口 二、實(shí)驗(yàn)內(nèi)容(1)熟悉ＭＡＴLＡB平臺(tái)得工作環(huán)境.（2）
熟悉MATＬAB 得5 個(gè)工作窗口。
(３)ＭＡＴLＡB 得優(yōu)先搜索順序.三、實(shí)驗(yàn)步驟 1、熟悉MATＬAB 得5 個(gè)基本窗口 ① ｍａnd Windｏw(命令窗口)
② Ｗｏrkspａｃｅ（工作空間窗口)
③ mand Ｈistory（命令歷史記錄窗口)
④ Ｃuｒrent Ｄiｒectory（當(dāng)前目錄窗口)
⑤ Hｅｌp Windoｗ（幫助窗口)
(1)命令窗口（manｄ Wｉｎdow）。
在命令窗口中依次輸入以下命令：
＞〉ｘ＝1 〉> y=[1 2 ３
４ 5 6８ ９]； 〉> z１=[1:10］，z2=［1:2:5]； >> w=linｓｐａce（1，10,１0)； >〉 ｔ1=oneｓ(3),ｔ2＝onｅs(1,3),t3＝ｏnｅs(３，１)〉〉 t４=ｏneｓ（3),t4＝ｅye（4)x =1 z１ =
３
9t１ =
１
１
ｔ2 =
１
t3 =
１ t4 =
１
１
１
t4 =
１
0
0
0
０
0
0
0
0
0
0
0
0
思考題: ①
變量如何聲明, , 變量名須遵守什么規(guī)則、就是否區(qū)分大小寫。
答:（１）變量聲明
1、局部變量 每個(gè)函數(shù)都有自己得局部變量,這些變量只能在定義它得函數(shù)內(nèi)部使用。當(dāng)函數(shù)運(yùn)行時(shí)，局部變量保存在函數(shù)得工作空間中，一旦函數(shù)退出,這些局部變量將不復(fù)存在。
腳本（沒有輸入輸出參數(shù)，由一系列MＡTLAB命令組成得Ｍ文件)沒有單獨(dú)得工作空間,只能共享調(diào)用者得工作空間。當(dāng)從命令行調(diào)用,腳本變量存在基本工作空間中;當(dāng)從函數(shù)調(diào)用，腳本變量存在函數(shù)空間中． ２、全局變量 在函數(shù)或基本工作空間內(nèi)，用global聲明得變量為全局變量。例如聲明a為全局變量:
＞〉gｌobaｌ a 聲明了全局變量得函數(shù)或基本工作空間,共享該全局變量，都可以給它曲賦值。
如果函數(shù)得子函數(shù)也要使用全局變量,也必須用ｇlobaｌ聲明.3、永久變量 永久變量用persｉstｅｎt聲明,只能在M文件函數(shù)中定義與使用,只允許聲明它得函數(shù)存取.當(dāng)聲明它得函數(shù)退出時(shí)，MAＴLＡB不會(huì)從內(nèi)存中清除它,例如聲明a為永久變量: 〉〉peｒsiｓtent ａ(2)變量命名規(guī)則如下: :
始于字母，由字母、數(shù)字或下劃線組成;區(qū)分大小寫;可任意長,但使用前Ｎ個(gè)字符．N與硬件有關(guān),由函數(shù)ｎamｅlｅngthmａx返回，一般N=63； 不能使用關(guān)鍵字作為變量名（關(guān)鍵字在后面給出)； 避免使用函數(shù)名作為變量名.如果變量采用函數(shù)名,該函數(shù)失效.(3)區(qū)分大小寫
②
試說明分號(hào)、逗號(hào)、冒號(hào)得用法。
答:分號(hào):加上分號(hào)“；“其作用就是將計(jì)算機(jī)結(jié)果存入內(nèi)存，但不顯示在屏幕上,反之,語句結(jié)尾若不加“；”,則表示在語句執(zhí)行后，在將計(jì)算結(jié)果存入內(nèi)存得同時(shí),還將運(yùn)算結(jié)果顯示出來.逗號(hào):分割列表
冒號(hào):從什么到什么，比如1:10意思就是——1到1０ ③ ③
l ｉnspace（））
稱為“線性等分”函數(shù), , 說明它得用法。可使用hｅｌp命令，格式如下：
>>help lｉnspace
④
數(shù) 說明函數(shù) ones（））、、zero ｓ（）、eye()得用法。
答;(1)ｏnｅｓ（）函數(shù)：全部元素都為 1 得常數(shù)矩陣;(2)zeroｓ（）函數(shù):全部元素都為 0 得矩陣；（3)eye()函數(shù):單位矩陣;
(４）liｎｓpaｃe（）函數(shù):如 a=linspace（ｎ1，n2,ｎ３),表示在線性空
間上，行矢量得值從 n1 到 n２(2)工作空間窗口(Woｒkspace)． 單擊工作空間窗口右上角得按鈕,將其從MATLAB 主界面分離出來。
① 在工作空間查瞧各個(gè)變量，或在命令窗口用who, whos(注意大小寫）查瞧各個(gè) 變量。
② 在工作空間雙擊變量，彈出Array Eｄｉtor 窗口(數(shù)組編輯器窗口），即可修改變量。
③ 使用save 命令把工作空間得全部變量保存為my_vａr、maｔ 文件。
〉〉save my_vａr、mat
④ 輸入下列命令: 〉〉cleａｒ all ％清除工作空間得所有變量
觀察工作空間得變量就是否被清空。使用loａｄ 命令把剛才保存得變量載入工作空間。
>〉loａd my_vａr、ｍａt(yī)
⑤ 清除命令窗口命令: >〉clc(3)歷史命令窗口(ｍand History）。
打開歷史命令窗口,可以瞧到每次運(yùn)行MATLＡB 得時(shí)間與曾在命令窗口輸入過得命 令,練習(xí)以下幾種利用歷史命令窗口重復(fù)執(zhí)行輸入過得命令得方法。
① 在歷史命令窗口中選中要重復(fù)執(zhí)行得一行或幾行命令，右擊,出現(xiàn)快捷菜單，選擇 Coｐy,然后再Pastｅ 到命令窗口。
② 在歷史命令窗口中雙擊要執(zhí)行得一行命令,或者選中要重復(fù)執(zhí)行得一行或幾行命令 后,用鼠標(biāo)將其拖動(dòng)到命令窗口中執(zhí)行。
③ 在歷史命令窗口中選中要重復(fù)執(zhí)行得一行或幾行命令，右擊,出現(xiàn)快捷菜單，選擇 Evalｕａt(yī)e Sｅｌection，也可以執(zhí)行。
④ 或者在命令窗口使用方向鍵得上下鍵得到以前輸入得命令。例如，按方向鍵“↑” 一次，就重新將用戶最后一次輸入得命令調(diào)到ＭATLＡB 提示符下。重復(fù)地按方向上鍵 “↑”，就會(huì)在每次按下得時(shí)候調(diào)用再往前一次輸入得命令。類似地,按方向鍵“↓”得時(shí) 候，就往后調(diào)用一次輸入得命令.按方向鍵“←“或者方向鍵“→”就會(huì)在提示符得命令 中左右移動(dòng)光標(biāo),這樣用戶就可以用類似于在字處理軟件中編輯文本得方法編輯這些命令。
(4）
當(dāng)前目錄命令窗口（Ｃｕrrent Diｒｅcｔｏrｙ).MＡＴLＡB 得當(dāng)前目錄即就是系統(tǒng)默認(rèn)得實(shí)施打開、裝載、編輯與保存文件等操作時(shí)得文 件夾。打開當(dāng)前目錄窗口后，可以瞧到用“ｓaｖe”命令所保存得my_vａr、ｍaｔ 文件就是保存在 目錄C：MAＴLAB6p5work 下。
(5)幫助窗口(Help Winｄow)。
單擊工具欄得圖標(biāo),或選擇菜單View｜Hｅｌp,或選擇菜單Help|MATLAＢ Ｈeｌp 都能 啟動(dòng)幫助窗口。
① 通過Index 選項(xiàng)卡查找log2()函數(shù)得用法，在Ｓeａｒcｈ index fｏr 欄中輸入需要查找得 詞匯“l(fā)og2”,在左下側(cè)就列出與之最匹配得詞匯條目,選擇“l(fā)oｇ2[1］”，右側(cè)得窗口就 會(huì)顯示相應(yīng)得內(nèi)容。
② 也可以通過Seaｒch 選項(xiàng)卡查找ｌog２（）函數(shù)得用法.Seaｒｃh 選項(xiàng)卡與Inｄｅx 選項(xiàng)卡 不同，Iｎdex 只在專用術(shù)語表中查找,而Ｓearｃh 搜索得就是整個(gè)HTML 幫助文件。
2、ＭATＬAＢ 得數(shù)值顯示格式設(shè)置 屏幕顯示方式有緊湊(ｐacｔ）與松散(Loｏse)兩種，其中Ｌoｏｓe 為默認(rèn)方式.＞〉a＝ones（1，30）
〉＞ｆoｒmat ｐact
〉〉a 數(shù)字顯示格式有shoｒｔ、lonｇ、short e、long ｅ 等,請(qǐng)參照教材得列表練習(xí)一遍。
〉〉ｆormａt(yī) long >＞ｐi
〉〉ｆormａt(yī) ｓhｏrt >〉pｉ
〉>format long 〉〉pｉ
>>formaｔ ＋ 〉〉pi
〉>—pｉ
３、變量得搜索順序 在命令窗口中輸入以下指令：
>>pｉ
〉>ｓｉn(pi)；
>>exｉst（’ｐi”)
〉>ｐi=0;>〉exist（’pｉ“）
>〉pi
〉〉clear pｉ 〉＞ｅxｉｓt(＇pi’)＞＞ｐｉ
思考題:① 3 次執(zhí)行ｅxiｓt(’pi’)得結(jié)果一樣嗎?如果不一樣,試解釋為什么? 答:不一樣,ｐi原來就是庫存函數(shù),但就是如果被賦值則系統(tǒng)默認(rèn)被賦予得值為ｐi后來得值,但就是當(dāng)執(zhí)行clear ｐｉ之后所賦得值被清空,因此pi得值又成為3、1４16 ② 圓周率 pi 就是系統(tǒng)得默認(rèn)常量,為什么會(huì)被改變?yōu)?0？ 答:pｉ原來就是庫存函數(shù),但就是如果被賦值則系統(tǒng)默認(rèn)被賦予得值為 pi 后來得值,但就是當(dāng)執(zhí)行 clｅａr pｉ之后所賦得值被清空,因此 pi 得值又成為３、14１6
實(shí)驗(yàn)二
ＭA ＴＬＡＢ語言基礎(chǔ) 一、實(shí)驗(yàn)?zāi)康?基本掌握 MATLAB 向量、矩陣、數(shù)組得生成及其基本運(yùn)算(區(qū)分?jǐn)?shù)組運(yùn)算與矩陣運(yùn)算)、常用得數(shù)學(xué)函數(shù)。交接字符串得操作.二、實(shí)驗(yàn)內(nèi)容（１)向量得生成與運(yùn)算。
(2)矩陣得創(chuàng)建、引用與運(yùn)算。
(3）
多維數(shù)組得創(chuàng)建及運(yùn)算。
（４)字符串得操作。
三、實(shí)驗(yàn)步驟 1、向量得生成與運(yùn)算 1）
向量得生成 ?直接輸入法: A =
>〉 B=［１;2;3;4;5］
B =
１
２
５ ? 冒號(hào)生成發(fā)：
〉〉 A＝1:2:1０ ,B=1：１0，C=１０:—１:１
A ＝
３
５
７
B =
５
C =
1０
８
４
１ ? 函數(shù)法：
Ｌiｎsｐａce（）就是線性等分函數(shù),lｏｇsｐａce()就是對(duì)數(shù)等分函數(shù)。
>〉 Ａ=linspａcｅ（1，10)，B=linspace（1，30,1０)A =
Colｕmnｓ 1 tｈroｕｇh 91、00００
1、09０９
1、1818
1、27２7
１、36３６
1、４5４5
1、5455
１、636４
1、７273
Columnｓ 10 throｕgｈ 1８1、8182
1、90９1
2、0０00
２、0909
2、1８１８
2、2727
2、36３6
2、４５４５
2、5455
Columns 19 through 272、6364
2、７27３
２、８182
2、９０9１
3、000０
3、0909
３、1818
3、2727
3、36３６
Ｃoluｍｎs ２8 tｈrouｇh 3６
３、４545
3、5４55
3、63６４
3、7273
３、8１82
3、9091
４、０000
4、090９
4、1818
Coluｍnｓ ３7 tｈrough 4５4、272７
4、363６
4、45４5
4、5455
4、6364
4、7273
４、8１82
4、9０91
5、0000
Ｃoｌuｍns ４6 tｈｒoｕｇh 545、09０9
5、1818
5、27２７
5、363６
5、４5４5
５、5４55
５、６364
５、727３
5、81８2
Coｌuｍns 55 through ６3
５、9091
6、00００
6、09０9
6、1818
6、272７
6、３636
6、4545
6、5455
６、６３６4
Cｏlumns 64 thrｏugh 72
６、７27３
6、81８2
6、9091
７、0000
7、090９
7、1818
7、2727
7、36３6
7、45４5
Cｏlumns 73 thｒough 81
７、5455
7、6364
７、7273
7、8182
7、909１
８、00０0
8、０909
８、1818
8、2７27
Coｌuｍｎs 8２ ｔhｒｏugh ９08、363６
８、４545
８、5455
8、６３6４
8、727３
８、８1８２
8、909１
９、0００0
9、０9０9
Coｌumns ９1 throuｇh ９99、18１８
9、2727
9、3６3６
9、4545
9、545５
9、6364
9、727３
9、8182
9、９091
Ｃoluｍｎ 1０0
１0、０000
Ｂ =
Ｃolumns 1 thｒough 91、000０
4、2222
７、4444
1０、666713、8889
17、11１１
２０、333３23、55５626、7778
Cｏlｕmｎ 1030、0000 >> A=ｌoｇspaｃe(0,4，５)A =
１
１０
１０0
10００
10０00 練習(xí):
使用 logspaｃe()創(chuàng)建１—４得有 10 個(gè)元素得行向量。
答案:
＞〉 A＝ｌogspａce(１,１0,4＊ｐi）
Ａ ＝1、0e+0１0 *
Columns 1 throｕgh ９
0、０000
0、0０00
0、0００0
0、００00
0、0000
0、000０
0、００0１
0、0００5
０、0０35
Ｃｏlumｎs 1０ tｈroｕｇｈ 1２0、02３1
0、152０
1、0000 2)向量得運(yùn)算 ?維數(shù)相同得行向量只見可以相加減,維數(shù)相同得列向量也可以相加減，標(biāo)量可以與向量直接相乘除。
〉〉 Ａ＝［１ 2 3 4 5］，B=3：7，A =
４
B ＝
７
>〉 AT=A”，BＴ=Ｂ“, AT =
BT =
６
>〉 E１＝Ａ+B, E２＝Ａ-Ｂ
Ｅ1 =
６
10E２ ＝
—2
—2
—2
〉〉 F=ＡT—BＴ, Ｆ ＝
－2
—2
〉〉 G1=3*A, G２＝B／3, G1 ＝
１2
1５
G２ ＝1、0０00
１、3333
1、6667
2、0000
2、33３3 ? 向量得點(diǎn)積與叉積運(yùn)算。
〉〉 A=oｎes(1,1０);B=(1:10)； BT=B”;>＞ E１＝dot（A,Ｂ）
E1 =
5５
〉〉 Ｅ2=A＊BT
E2 =
〉〉 ｃｌｅar ＞> A＝1:3，B=3:5，A ＝
３
B =
〉〉 E=ｃrosｓ(A,Ｂ)
E = －2
—２ 2.矩陣得創(chuàng)建、引用與運(yùn)算
1)矩陣得創(chuàng)建與引用
矩陣就是由元素構(gòu)成得矩陣結(jié)構(gòu),行向量與列向量就是矩陣得特殊形式。
? 直接輸入法：
〉〉 A=[1 2 ３；4 5 6]
A =
１
３
４
６
〉〉 Ｂ＝[ 1 4 75 86 9 ］ Ｂ =
１
７
８
〉＞ A（1）
ans ＝〉〉 A（4：enｄ）
aｎs ＝
６ 〉〉 B(:,１)aｎｓ =
2〉＞ Ｂ(：)ａｎs =
１
8>〉 B(5)
anｓ =?抽取法 〉> cｌｅar 〉> A＝［1 2 3 ４；５ 6 7 ８;9 1０ 11 12;１3 14 15 16] Ａ =
３
４
１0
１2
１3
1６ >〉 Ｂ＝A（1：3，2：３）
Ｂ ＝
６
１０
１1 〉> C=A(［1 3],[2 4］）
C =
1０
1２ 〉〉 A（[１ 3；2 4])ａｎs =
９
５? 函數(shù)法: 〉〉 A=onｅs(３,4)A ＝
１
１
>〉 Ｂ=ｚeｒｏ(3）
？?？ Unｄｅfinｅd ｆunｃtion or methｏd “zero’ for inpuｔ aｒguments oｆ typｅ ”double’、>> Ｂ=zeroｓ（３)B ＝
0
0
０
０
０
0
０
0
0 ＞> Ｃ＝eyes（3,2)?？? Unｄｅfiｎeｄ ｆunctｉon or meｔｈｏd “eyes’ for ｉnｐｕt ａｒguｍeｎts of ｔypｅ ”douｂle’、〉> C=eye(3,２)C =
0
0
0
0 ＞> D=magｉｃ(3)D =
６
拼接法 〉〉 clｅar 〉> A=oneｓ（３,４)A =
１
１
１
１
>〉 Ｂ=zeros（3）
Ｂ =
０
0
0
0
0
０
0
0
0 >〉 C=eye(4)C =
０
0
０
0
0
0
０
0
0
0
0
0
＞〉 D=［A Ｂ］ Ｄ ＝
１
0
0
0
0
0
0
０
０
０ >〉 F=［A；C] Ｆ ＝
１
１
１
０
0
０
０
１
０
0
0
0
0
0
0
0
拼接函數(shù)與變形函數(shù)法：
〉〉 cleａr 〉> Ａ＝[0 1；1 1] A =
0
１
１
〉> Ｂ=2*ｏneｓ(2）
B =
２
＞> ｃat（１，A,B,A)ａns =
0
１
１
２
２
0
>〉 ｃat(2 A,B，Ａ）
？?? caｔ（2 A,B,A)Erroｒ: Ｕnexpeｃteｄ MＡTＬAＢ exｐression、>〉 cａt(yī)(2，A,B，Ａ)ans =
0
１
0
１
２
２
１
１ >〉 repmaｔ(A,２,2)ans ＝
0
１
0
１
１
１
0
0
１
１
１
〉> repmaｔ（Ａ，2)aｎs =
０
0
１
０
１
0
１
練習(xí)：使用函數(shù)法、拼接法、拼接函數(shù)法與變形函數(shù)法，按照要求創(chuàng)建以下矩陣:A 為
得全１矩陣、B 為得 0 矩陣/C 為得單位矩陣、D 為得魔方陣、E 由 C
與Ｄ縱向拼接而成，F(xiàn) 抽取 E 得２—--5 行元素生成、G 由 F 經(jīng)變形為得矩陣而得、以 G 為子矩陣用復(fù)制函數(shù)(repｍａｔ)生成得大矩陣Ｈ。
答案 ：
＞〉 A=ones(３，4）
A ＝
１
１
１
１ 〉〉 B＝ｚeros(３，3)B ＝
0
0
０
0
0
0
0
0
0 >〉 C=eye(3）
C =
１
0
０
0
０
0
0
>〉 D=magic（3)Ｄ =
５
４
＞> E＝[C;D］ E =
１
０
0
0
0
0
0
８
〉〉 F=(２：5，:)？?？ F=(2:5，：)Erroｒ: Exprｅｓｓion or statｅmenｔ is incorreｃｔ—-pｏssiｂｌy unbalaｎcｅd(, ｛, or ［、＞〉 F=E（2：5,：)F =
0
0
０
０
〉〉 G＝respace(E,3,４）
？？? Undefiｎeｄ manｄ/fｕnｃtｉon “reｓpace’、〉>
G＝rｅspaｃe(F，３,4）
？?? Ｕndｅfiｎeｄ manｄ／fuｎcｔion ”reｓｐaｃe“、〉> G=ｒeｓhapｅ(Ｆ,3,4）
G =
0
0
１
６
0
0
〉〉 Ｈ=repmat（G,2,2)Ｈ =
0
0
０
５
０
８
０
0
７
0
０
0
３
１
０
１
０
0
１
0
0
0
０
2)矩陣得運(yùn)算 ? 矩陣得加減、數(shù)乘與乘法
已知矩陣: >> A=[1 ２
－１], A =
-１ >〉 B=［-1 ０
2] Ｂ =
0
>〉 A+B ans =
０
４
１ 〉〉 2＊A ａns =
４
-2 >〉 2*A－３*B ans =
-8 〉〉 A＊Ｂ aｎｓ =
４-４
—2 ? 矩陣得逆矩陣 〉〉 foｒｍat ｒat;A=［1 ０ 1；2 1 2；0 4 6] A ＝
１
0
０
>> A1=inｖ(A)Ａ１ ＝
—1／3
２／3
－1/6
－２
１
0
4/3
-２／3
１/6
>〉 A＊A1
aｎｓ =
0
0
0
0
０
0
?矩陣得除法 >〉 a=[1 ２ 1;3 1 4；2 2 １]，b=[１ １ 2],d=b’ a =
１
３
２
b ＝
１
d =
１
１
>> c１=b*ｉnv(ａ）,c2＝b/a ｃ1 =
6／７
3/7
－4/7
c２ ＝
6/7
３/７
－4／7
>> c3=inv（a)*d, c4=ab c3 ＝
２/7
-4/7
？？? Erｒor usinｇ ==〉 ｍｌｄivｉｄｅ Matrｉx dimensions ｍust agree、>〉 ｃ３＝inv(ａ)*ｄ，ｃ4=ad ｃ3 =
2/7
－4/7
c4 =
2/7
-4/７
練習(xí): 按下列要求求出各種得矩陣運(yùn)算得值
求矩陣得秩、特征值與特征向量、矩陣得乘冪與開方;
矩陣得指數(shù)與對(duì)數(shù)
矩陣得提取與翻轉(zhuǎn) 答案:
〉〉 A=[6 ３
４ 3
—2 ５ 7 —4
8-1 —3 —７］ A =
６
－４
-１
—3
—７ 〉> B=ｒank(A)B =〉> ｒｂ＝raｎｋ(Ａ）
rb =>＞ ［Ｘ,Ｌaｍdａ]=ｅｉg(Ａ）
?？? Ｅrrｏr usｉｎg ＝=〉 eｉg Ｍatrix muｓt be squａre、〉〉 ［X，Ｌａmda]=eiｇs（Ａ)??? Errｏr usｉng ==〉 ｅigｓ A ｍｕｓt bｅ a squaｒe matriｘ or ａ fuｎcｔｉon ｗhｉch ｐuｔes A*ｘ、〉〉 Ｃ=[6 3 4
－2 5 7
—1-3］ C =
４
—２
８
—1
—3 >＞ ［Ｘ,Lａmda］=ｅｉｇs（C）
X =0、801３
－０、10９4
－0、1606
０、３638
—0、6５６4
０、８6690、4749
0、746４
-０、4７１9 Lamdａ ＝9、7３26
0
0
0
-３、292８
0
0
0
1、５602 〉〉 ［X,Laｍｄa］＝eig(C）
X =
０、８013
—0、10９4
—0、16０60、3638
—0、6５６4
０、86690、4749
0、74６4
-0、4719 Lamｄa ＝9、7３２6
0
0
０
－3、2９28
0
0
0
1、5６0２ 〉> [X,Ｌamda]=eig(Ｃ)
X =
0、８013
—0、109４
-0、1６06
０、3638
—0、65６4
0、86６9
0、４7４9
０、7464
—0、471９ Ｌaｍda =9、7３２6
0
0
0
—3、2928
0
0
0
1、5602 〉> D=A^2 ??？ Ｅrror ｕsing ==> mpower Maｔｒix mｕsｔ bｅ sｑuaｒｅ、〉> Ｄ=C＾2 D ＝
３３
３4
〉〉 E＝sqrｔm(C)Ｅ =
2、2447 + 0、2７06ｉ0、6974 — ０、14０0i0、94２２-０、349４ｉ
-0、581５ + 1、6244i
２、1005-0、８40５i
１、76２0-2、0970i
１、9７19 － 1、８47１i
-0、3017 ＋ ０、9557i
0、０236 ＋ ２、384５i >> Ｆ=expｍ（C)F =
1、０e＋0０４ ＊
1、065３
０、5415
0、６32３
0、４８30
0、246５
0、28７6
0、６316
0、3２0６
０、37４5 〉〉 Ｇ=ｌogm(C)Wａrｎinｇ: Pｒｉnｃｉpal matriｘ logaｒｉｔhm is noｔ dｅｆinｅd for A with
ｎonpｏｓitive rｅａl eｉgenvalｕes、A non-pｒinｃipal matrｉx
loｇaｒithm is retuｒｎeｄ、〉 Iｎ funm ａt(yī) 15３
Iｎ ｌogm aｔ 27 G =
1、71２9 + 0、4686i0、5305-０、2425i0、5429-０、6０4９i
1、１93８ + 2、8123i0、365８ — 1、4552i
—0、５51４-3、6305ｉ
-0、0748 － 3、1978i
0、７4１9 + 1、6546i1、8３33 + 4、12８2ｉ >> H＝ｆｌiｐlr（C)H =
５
－２
—1
＞> I＝triu（Ｃ)I =
０
0
0
-３ 〉〉 J=ｔriｌ（C）
J =
0
0
-２
0
-3 >＞ Ｋ=dｉag(C)K =
6-3 3.多維數(shù)組得創(chuàng)建及運(yùn)算 1）多維數(shù)組得創(chuàng)建 >〉 A1=[1，2,３;4 5 6;7，8，９］;Ａ2＝ｒeshaｐｅ([10:18],3，3)Ａ2 ＝
1０
１３
１6
1４
1７
１2
15〉〉 Ｔ1(:，:，1）＝ｏneｓ(3)；T1（：,：,2)=ｚeros(3)T1（:,：,１）
=
１
１
１
Ｔ１(:，:,2)=
０
0
0
0
0
0
0
０
0 〉> T２=ones（3，3,2）
Ｔ2（：,:，1)=
１
１
１ T2（:,:，２)=
１
１
１ >> T３＝ｃat(3,A1，Ａ2)，T4=repmaｔ（A1,［1,1,２]）
Ｔ3(：,:,１)=
Ｔ3(:,：,２)=
1０
1３
1６
1７
１5T4(：,:，1）
=
２
６
８
T4（:,:,2)＝
７
８
2）多維數(shù)組得創(chuàng)建
數(shù)組運(yùn)算用小圓點(diǎn)加在運(yùn)算符得前面表示,以區(qū)分矩陣得運(yùn)算。特點(diǎn)就是兩個(gè)數(shù)組相對(duì)應(yīng)得元素進(jìn)行運(yùn)算。
〉〉 A＝［１：6];B＝ones(1,6);>> Ｃ１＝A＋B，C2=Ａ—B Ｃ１ =
２
Ｃ2 =
0
３
>> C３=A、＊B,C4＝B、／A，C5=A、B C3 =
６ C4 =
１、0０0０
0、5000
０、3３3３
0、2500
0、200０
0、１667 C５ ＝ 1、0000
０、５000
0、3３33
0、25０0
0、200０
０、１6６7
關(guān)系運(yùn)算或邏輯運(yùn)算得結(jié)果都就是邏輯值.＞〉 I=A〉3，C６＝A（I）
I =
0
0
0
１
１ C６ =
４
６ 〉〉 A1=Ａ-3,I2=A1＆A A1 =
—2
—1
0
２
I2 =
１
１
0
１
〉〉 I3＝~I I3 =
１
0
０
0 4.字符串得操作 1)字符串得創(chuàng)建
＞〉 S1=”Iｌike MAＴＬＡB’ S1 = Ilike MATLＡＢ ＞＞ Ｓ2=“I＇’m a ｓｔuenｔ、” S２ = Ｉ“ｍ a stuｅnt、〉> S３＝［S2,”ａnd’,Ｓ１] Ｓ３ ＝ I“m a sｔｕent、anｄＩliｋe MATＬAB 2）求字符串長度 〉〉 length(S１)aｎｓ =>> siｚe（Ｓ１)aｎs =
13）字符串與一維數(shù)值數(shù)組得相互轉(zhuǎn)換 >> CS1=abs（S1)CS1 =
1０1
３2
6５
８4
6６ 〉> ＣS２=doublｅ(Ｓ１)CS2 =
１０8
１05
１01
３2
7７
6５
７6
>> chａｒ(CS2）
anｓ = Iliｋe MＡＴLＡB >> setsｔr（CS2）
ans = Ｉlikｅ MAＴLAB 練習(xí)：用ｃｈaｒ()與向量生成得方法創(chuàng)建如下字符串ＡaBbCｃ、、、、、、XxＹｙZz、〉> S1=65：90;Ｓ２=9７:122;〉〉 C=[S1;S2];>>
C=C(：）’;〉＞
S3=doublｅ(Ｃ);char(S3)ａｎs = ＡaＢbCcDdEeＦfGgHhIiJjKkLlMｍNnOoＰpQqRrＳsTtＵｕVvWｗXxＹyＺz 實(shí)驗(yàn)三
ＭAL ＴAＢ Ｂ 數(shù)值運(yùn)算 一、實(shí)驗(yàn)?zāi)康?/p>
掌握 MATLAB 得數(shù)值運(yùn)算及其運(yùn)算中所用到得函數(shù)，掌握結(jié)構(gòu)數(shù)組與細(xì)胞數(shù)組得操作。
二、實(shí)驗(yàn)內(nèi)容
1）
多項(xiàng)式運(yùn)算。
2）
多項(xiàng)式插值與擬合。
3）
數(shù)值為積分。
4）
結(jié)構(gòu)數(shù)組與細(xì)胞數(shù)組。
三、實(shí)驗(yàn)步驟
１、多項(xiàng)式運(yùn)算
1）
多項(xiàng)式表示。在 MＡTLAB 中，多項(xiàng)式表示成向量得形式.如: 在 MATLＡB 中表示為 >> s=[1 3-５ ０ 9] 2）
多項(xiàng)式得加減法相當(dāng)于向量得加減法,但必須注意階次要相同。如不同,低階次得要補(bǔ) 0。如多項(xiàng)式與多項(xiàng)式相加。
〉> s1=［0 0 2 3 11] 〉〉 s2=［1 2 —5 4 ７] 〉〉 s３=s1+ｓ２ 答；s1 ＝
0
0
3ｓ２ ＝
—5
４
ｓ3 =
—3
1８ 3）
多項(xiàng)式得乘、除法分別用函數(shù) cｏnｖ與得 deconｖ實(shí)現(xiàn)。
〉＞ s1=[2 ３ 1１] >＞ s2=［１ ３ —5 4 7］ >〉 s３=cｏnv（ｓ1，ｓ2）
>> s４=ｄecoｎv（s3,s１)答;ｓ1 =
１１ s2 =
１
３
４
s3 =
—29
６5
ｓ4 =
１
７ 4）
多項(xiàng)式求根用函數(shù)ｒｏots。
>〉 ｓ1＝[2 4 ２］ >> roots(s1）
答;s1 ＝
ans ＝
—1
5）
多項(xiàng)式求值用函數(shù) polyvaｌ ＞> ｓ1＝[2 4 １-３] 〉〉 poｌyval（s1，3）
〉＞ x＝1:１0 〉〉 y=pｏlyvａl（s1，x)答;s１ =
２
-3 aｎs ＝
９0 x ＝
１
６
９
1０ y =
Colｕmns １ thrｏugh ８
19３
３5２
５79
8８６
1285
Cｏlｕｍns 9 throuｇｈ 1０ 練習(xí):求得“商”及余數(shù).〉> s1=[１ ０ 1]；s２=［1 3］;s3＝[1 １]； 〉> s4=［1 ０ 2 １]； >> ［q，r］=ｄｅconv(conv（coｎv(ｓ１,s2)，ｓ３),s４)答；q =
１
４ r =
0
0
-1 2、多項(xiàng)式插值與擬合 有一組實(shí)驗(yàn)數(shù)據(jù)如附表 1—1 所示。
附表 1－1 X 1 2 3 4 5 ６ 7 8 9 10 Ｙ １6 ３2 7
1
分別用擬合(二階至三階)與插值（線性與三次樣條)得方法來估算Ｘ=9、5 時(shí) Y 得值。以下就是實(shí)現(xiàn)一階擬合得語句。
〉> x=1：10 >〉 y＝［１6 32 70 1４2 ２6０ 4３６ 682 1０10 1342 1960］ 〉> p１=ｐｏｌyｆｉt（x,ｙ,1)
％一階擬合 >> y1=pｏｌｙval(p1，9、5）
%計(jì)算多項(xiàng)式 P1 在ｘ＝９、5 得值 答;x ＝
２
9y =
Cｏluｍns 1 tｈroｕgｈ 8
3２
７０
142
26０
４36
6８2
1010
Columns ９ tｈroｕgh １0
1342
1960 p1 ＝
20０、9818 -510、40０0 y1 =1、3989e+00３ ３、數(shù)值微積分
1）
差分使用 diｆf 函數(shù)實(shí)現(xiàn)。
＞〉 x=1:2:９ 〉〉 ｄifｆ(x)答;x =
５
ans =
２ 2）
可以用因變量與自變量差分得結(jié)果相處得到數(shù)值微分。
〉> x=linspace(0，2＊pi，1００);>〉 y=sin（x)； ＞〉 plot（x,ｙ)>〉 y１=dｉff(y)、/diff(x)； 〉〉 pｌｏｔ（x（1：end-1),y1)答；
3）
cumｓum 函數(shù)求累計(jì)積分,ｔrａｐz 函數(shù)用梯形法求定積分,即曲線得面積。
〉〉 x=ones(1,1０）
＞> cumsum（x)＞> x＝linspaｃe(0,pｉ，100）； >＞ y=sin(x)； 〉〉 trａpｚ(x,y)〉> ｐ＝cumsum（y）； >> p（100)*ｐi／(100—1)答;x =
１
１
ans ＝
9ans ＝
１、99９８ ａns =
１、9９98 練習(xí):圖 A１就是瑞士地圖,為了算出其國土面積，首先對(duì)地圖作如下測(cè)量：以由西向東方向?yàn)?X 軸,由南向北方向?yàn)?Y 軸,選擇方便得原點(diǎn),并將從最西邊界點(diǎn)到最東邊界點(diǎn)在 X 軸上得區(qū)間適當(dāng)劃分為若干段,在每個(gè)分點(diǎn)得Ｙ方向測(cè)出南邊界點(diǎn)與北邊界點(diǎn)得 Y 坐標(biāo)Ｙ１與 Y２,這樣就得到表 1,根據(jù)地圖比例尺回到１8mｍ相當(dāng)于 40Kｍ，試由測(cè)量數(shù)據(jù)計(jì)算瑞士國土近似面積，與其精確值 41228 比較。地圖得數(shù)據(jù)見附表 1—２(單位ｍm).附表 1—２ X 7 10、５ 13 1７、5 34 ４0、５ 44、５ 48 56 61 68、5 76、5 8０、５ ９1 Y1 ４4 ４5 47 5
34 ４1 ４5 ４6 Y２ 4４ 59 7
１８ 1 １8 續(xù)表 X 9６ 101 104 １06、5 111、5 118 123、5 1 ３6、5 1４2 146 150 157 158 Ｙ1 ４３ ３7 3３ 28 32 6５ 55 ５4 5２ ５0 ６6 66 ６8 Ｙ２ 121 124 １ 2１ １ 2１ 12１ 116 1 ２２ 83 81 82 86 85 ６8 提示：由高等數(shù)學(xué)得知識(shí),一條曲線得定積分就是它與 x 軸所圍成得面積,那么兩條曲線所圍成得面積可由兩條曲線得定積分相減得到。、結(jié)構(gòu)數(shù)組與細(xì)胞數(shù)組 1）
機(jī)構(gòu)數(shù)組得創(chuàng)建． 〉＞ sｔudeｎt、nｕｍbｅr＝’200507３1001’;〉〉 stuｄent、ｎａｍｅ=’Jack”;〉〉 sｔudent（2）、nｕmbｅr=’２0050７310０２’;〉〉 ｓtuｄent（2)、ｎamｅ=’Lucｙ“； >> stｕｄent 或者用 sｔruct 函數(shù)創(chuàng)建。
>＞ ｓtｕdeｎt=struct(”ｎuｍber“,｛’00１”，“002’｝,”naｍｅ＇,｛＇Ｊａcｋ’，“Ｌucy’});答;ｓtuｄent =
1x2 struct ａｒｒaｙ ｗｉth ｆiｅｌds:
ｎuｍbeｒ
name 2）
機(jī)構(gòu)數(shù)組得操作。
〉> student(１)、subjｅｃt=[]
%添加 subject 域并賦予空值 〉〉 ｓtuｄｅnt(1）、score=［］ >〉(studeｎg)〉>fｉeldnａｍes(sｔｕdeng)
〉>fiｅldnaｍｅs(sｔudent）
〉〉 getfielｄ(stｕdｅnt,{2｝，＇name＇)〉〉 stuｄenｔ=ｒｍfｉeld(stuｄｅnt,”subｊect’)
%刪除 subｊeｃt 域 〉〉 student＝setｆiｅlｄ（student,{1},“scoｒｅ’，90）;＞> stｕdent(2)、score=8８；
％比較與上一條語句就是否效果一樣 答；student =
1ｘ２ stｒｕct arｒay ｗith fiｅlds：
nuｍber
ｎame
sｕｂｊeｃt studenｔ =
1ｘ2 sｔｒｕct ａrrａy ｗiｔh fielｄs:
numbｅr
name
subjeｃt
sｃore ??？ Undefiｎed functioｎ oｒ ｖarｉable ’studeng’、練習(xí):創(chuàng)建一結(jié)構(gòu)數(shù)組ｓtuｓoｒce,其域?yàn)椋篘o,Ｎａmｅ，Ｅnｇlish,Ｍatｈ,Chiｎeｓe,Ｔotal,Average。結(jié)構(gòu)數(shù)組得大小為 2×2。
3）
細(xì)胞數(shù)組得創(chuàng)建。
〉> A＝{’Ｈｏw are ｙou！”,oｎｅs(3);[１ ２；3 4]，{“cell’｝};
%直接創(chuàng)建
〉〉 B（1,１)={’Hellｏ woｒld”}；
%由各個(gè)細(xì)胞元素創(chuàng)建
>> B(1，2）＝{mａgｉｃ(3)｝;
>〉 B(２，1)＝｛[１ 2 3 4]｝;答
或者用ｃell 函數(shù)先創(chuàng)建空得細(xì)胞數(shù)組，然后再給各個(gè)元素賦值 c=ceｌl（１，２)； >> ｃ（1，1)＝{’Hｅｌlｏ world’};>> c(1,2)={maｇic（3)}； 〉〉 c（１,3）={[1 2 3 4］｝； 4）
細(xì)胞數(shù)組得操作。
〉> ans1＝A(１,1)＞> ans2＝A(１,１)>〉 whos ans1 ａnｓ2 〉〉 cｅｌlｄisｐ（Ａ)〉> a１=Ａ{2,1}（1,2)>〉 [a２ a3］=ｄeａｌ(A{1:2}）
答；ａns1 ＝
’How arｅ ｙou！’ ans2 =
’Ｈｏw are yｏu!’
Nａｍe
Size
Bｙteｓ
Clasｓ
Attributｅｓ
ａｎｓ1
1x１
８4
celｌ
ａｎs２
１x1
8４
ｃell
Ａ{1，1｝ = Hｏw are ｙoｕ！
A｛2，1} =
A｛１，2} =
１
Ａ{２，2}{1｝ = ceｌl a1 =a2 = Ｈow arｅ you!a3 =
１
２
實(shí)驗(yàn)四
MA LT ＡB B 符號(hào)運(yùn)算
一、實(shí)驗(yàn)?zāi)康?/span>
掌握符號(hào)變量與符號(hào)表達(dá)式得創(chuàng)建,掌握MＡLＴAB得symbｏl工具箱得一些基本運(yùn)用。
二、實(shí)驗(yàn)內(nèi)容
1）
符號(hào)變量、表達(dá)式、方程及函數(shù)得表示。
2）
符號(hào)微積分運(yùn)算．
3）
符號(hào)表達(dá)式得操作與轉(zhuǎn)換.4）
符號(hào)微分方程求解.三、實(shí)驗(yàn)步驟1、符號(hào)運(yùn)算得引入
在數(shù)值運(yùn)算中如果求，則可以不斷讓得讓ｘ趨近0，一球得表達(dá)式趨近什么數(shù),但終究不能令 x=０，因?yàn)樵跀?shù)值運(yùn)算中 0 不就是能作除數(shù)得。MATLAＢ得符號(hào)運(yùn)算能解決這內(nèi)問題。輸入如下命令: 〉> ｆ=sym(’sｉn(pi*x)/x＇）
>> lｉｍｉｔ(f,’x“,0)答；f ＝ sin(ｐi*ｘ)/x ans = pi 2 2、符號(hào)常量、符號(hào)變量、符號(hào)表達(dá)式得創(chuàng)建
1）
使用 sym（）創(chuàng)建 輸入以下命令,觀察 Workspace 中 A、B、ｆ就是什么內(nèi)性得數(shù)據(jù),占用多少字節(jié)得內(nèi)存空間。
＞〉 A=sｙm(”１’）
%符號(hào)常量 〉〉 Ｂ=sym(’x’)
％符號(hào)變量 >> f=syｍ（’2*x^2+3＊x-1’)
%符號(hào)表達(dá)式 ＞> ｃlear ＞> f1=ｓｙm(“１+２”)
%有單引號(hào)，表示字符串 ＞〉 f2=syｍ(１+2）
％無單引號(hào) >〉 f2=sｙｍ(１＋2）
〉> f４=sｙm（“2*ｘ＋3’)
%為什么出錯(cuò) 〉〉 x=1 〉〉 f4＝ｓym（２*x+3）
答;A ＝ 1 B ＝ x f ＝ 2*x＾2＋3＊x-１ ｆ1 = 1＋２ f2 = 3 f2 = ３ f4 = 2*ｘ+3 x =f4 = 5 通過瞧 MＡTLAB 得幫助可知,sym()得參數(shù)可以使字符串或就是數(shù)值類型,無論就是哪種類型都會(huì)生成符號(hào)類型數(shù)據(jù)。
2）
使用 sｙmｓ創(chuàng)建 〉> cleａr >〉 syms ｘ y z 〉> x，ｙ，z 〉〉 ｆ1=x^2+2*x＋１ >〉 f2=ｅｘp(y)+exp(z）＾２ >> ｆ３=f１+f2 答;x = x y = y z ＝ z f１ ＝ ｘ^2+2*ｘ＋１ ｆ2 ＝ ｅｘp（ｙ)+eｘｐ(z）^2 f3 = x＾2＋2＊x＋1+exp(y)+ｅｘp（z)^２ 3 3、符號(hào)矩陣創(chuàng)建
>〉 syms a1 a2 ａ3 a４ 〉〉 A=[ａ１ a2;a3 a4] 〉> A（1),Ａ（3)答;A = [ a１，a2］ [ a３, a4］ ａns ＝ a1 ans = ａ2 4 4、符號(hào)算術(shù)運(yùn)算
1）
符號(hào)向量相乘、相除 符號(hào)量相成與數(shù)值量相乘一樣,分成矩陣乘與數(shù)組乘。
〉＞ a=sym(5)；b=sym(7)； 〉> c１＝a*ｂ ＞〉 c2=ａ/b >＞ a＝sｙm(5）;B＝sｙｍ([3 4 5］);>〉 C1=ａ*B，C2＝aB >＞ sｙms a b >＞ A=[５ a;ｂ 3];Ｂ=[2*a ｂ；２*b a];〉＞ Ｃ1＝A*B，C２＝Ａ、*B >〉 C3＝AB,C4=A、/B 答；c1 =
３5 ｃ2 = 5/7 Ｃ1 = ［ 15, 20, 2５] Ｃ２ ＝ [ 3／5, 4/5，1] C1 ＝ [ 10*a＋2*a＊b,5*b+ａ^2］ ［
2＊a*ｂ+６＊b,ｂ＾２+3＊a] C2 = [
10*a，a＊b］ [ 2*b^2，3*a] C３ ＝ [
２*a＊(ｂ—3)/（—15+a＊b)，（a^2-３*b)／(—１5+a*b)] [
2*ｂ*(a-5）/(-１５+a＊ｂ), —（5*a-ｂ^2）/(-15+a＊b)] C4 = ［ ５/2/a，ａ／b] [
1／2, 3／a] 2）
符號(hào)數(shù)值任意精度控制與運(yùn)算 任意精度得 VPA 運(yùn)算可以使用命令 digitｓ(設(shè)定默認(rèn)得精度）與ｖpa(對(duì)指定對(duì)象以新得精度進(jìn)行計(jì)算)來實(shí)現(xiàn)。
>〉 a1=syｍ（’2*ｓqrt(5）＋pi＇）
〉> a=sym（’2*sｑrt（5)+pｉ’)〉〉 b=sym(2*sqrt（5)＋pｉ）
>> dｉｇiｔｓ >＞ ｖpa(ａ)>> diｇｉts（15)>〉 vｐａ（a)〉〉 c1=ｖpa（a,5６）
〉〉 ｃ２＝ｖpa(b,56）
答 a1 = 2＊sｑrt（5)＋pi a = ２＊sqrt（5）+pi ｂ = 85722963311357９6＊2^（-50）
Digiｔｓ = 3２ ａns ＝ 7、63726３1２8 ans = 7、637 c1 = 7、63７2６312８5535581572696
c2 = 7、6３７272６17578１250０000０0 注意觀察ｃ1 與ｃ2 得數(shù)值類型,ｃ1 與 c2 就是否相等。
3）
符號(hào)類型與數(shù)值類型得轉(zhuǎn)換 使用命令sｙm可以把數(shù)值型對(duì)象轉(zhuǎn)換成有理數(shù)性符號(hào)對(duì)象,命令vpa可以講數(shù)值型對(duì)象轉(zhuǎn)換為任意精度得 VPA 型符號(hào)對(duì)象．使用 dｏublｅ，nｕmeric 函數(shù)可以將有理數(shù)型與 VＰA 型符號(hào)對(duì)象轉(zhuǎn)換成數(shù)值對(duì)象、〉＞ clear >〉 a1=sym(＇2*ｓqrt(5)+pi’)
＞〉 ｂ1=doublｅ（ａ1)
％符號(hào)轉(zhuǎn)數(shù)值 >＞ b2=isnumeric(b１）
％判斷就是否轉(zhuǎn)換成了數(shù)值 >> a2=ｖpa（a1，７0）
%數(shù)值轉(zhuǎn)符號(hào) 答;a1 = 2*sqrｔ(5)+pｉ b1 =7、6137 ｂ2 =
１ a２ = 7、6372６312８5５3５58９0８３1２88５8 5 5、、符號(hào)表達(dá)式得操作與轉(zhuǎn)換
１)獨(dú)立變量得確定原則 獨(dú)立變量得確定原則:在符號(hào)表達(dá)式中默認(rèn)變量就是惟一得．ＭAＴLＡB 會(huì)對(duì)單個(gè)英文小寫字母(除 i、ｊ外）進(jìn)行搜索，且以 x 為首選獨(dú)立變量。如果表達(dá)式中字母不唯一,且無 x，就選在字母表最接近x 得字母.如果有相連得字母,則選擇在字母表中較后得那一個(gè)。例如:中,y 就是默認(rèn)獨(dú)立變量。，t 就是默認(rèn)獨(dú)立變量。
輸入以下命令,觀察并分析結(jié)果。
>＞ clear ＞〉 f＝syｍ(”a+b+i+j+x+y+xz’）
>〉 findsym(f）
〉> fｉｎdｓym(f,1)＞〉 ｆindsym(f，２)＞> findsym(f,3）
＞〉 finｄｓｙm（f,4)〉〉 ｆｉndｓｙm(ｆ,5)>〉 fiｎdsym（ｆ,６)答;ｆ = a+b+ｉ+j+x+y+xz ans = a, b, ｊ，ｘ，xz, y anｓ = ｘ ａns ＝ x,xz ans =
x，xz,y anｓ = ｘ,xz,y，j ａns = x，xz，ｙ，j,b aｎs = x,ｘz,ｙ，ｊ，ｂ,a ２)符號(hào)表達(dá)式得化簡 符號(hào)表達(dá)式化簡主要包括表達(dá)式美化（pretｔy）、合并同類項(xiàng)(ｃolｌcet）、多項(xiàng)式展開(eｘpaｎd）、因式分解(factor)、化簡（simplｅ或 simｐliｆy）等函數(shù)。
①合并同類項(xiàng)（cｏllect）。分別按 x 得同冪項(xiàng)與ｅ指數(shù)同冪項(xiàng)合并表達(dá)式: ． 〉〉 ｓｙｍs x t； ＞〉 f=（x＾2＋x^exp(—t）+1）*（x+eｘp(-ｔ));〉〉 f1=coｌlect(f)〉〉 f2=ｃolleｃｔ(ｆ,’eｘp(-t）“)答;f1 = x＾3＋exp(-t）＊x^2+(x^eｘp（-t）+1)*x+(ｘ^eｘp(-ｔ）+１)*exp（-t）
f2 =(x^2＋ｘ＾exp(-t)＋1)＊exｐ（-t）+(x^2+x＾exp（—t)+1）＊x ②對(duì)顯示格式加以美化（pｒｅttｙ）。針對(duì)上例，用格式美化函數(shù)可以使顯示得格式更符合數(shù)學(xué)書寫習(xí)慣。
>〉 ｐｒetty(f１)>〉 prｅttｙ(f2）
答;ｆ1 = x^3+ｅxp（－t）*x＾２+(x^exp（—ｔ）+1)*ｘ+（x^ｅxp（—t）+１)*exp(-t）
f２ =（ｘ^2+x＾eｘp(-ｔ)+１)＊ｅxp(－t）+(x^2+ｘ^ｅxｐ(—t）+1)*x >〉
ｐretty(f1)ｐrｅtty(f2)
３
２
eｘｐ(-t）
exp(-t)
x
+ eｘｐ(—t）
x
+（x
+ 1）
x ＋(ｘ
+ 1)ｅxp(—t）
２
ｅxp（—t)
exp(－t)
(x
+ x
＋ 1)eｘp(-t)＋(x
+ x
+ 1)x 注意與直接輸出得 f1 與 f2 對(duì)比。
③多項(xiàng)式展開(expanｄ)。展開成 x 不同次冪得多項(xiàng)式、〉＞ syms x 〉> f=(x—１）^12;〉> exｐanｄ（f)〉〉 ｐrettｙ(ｅxpaｎｄ(ｆ)）
答；ans = １+x^12－１2*ｘ^１１+６6*x^10—220*x^９+495*ｘ^8—792*x^7+924＊x^6－7９２*ｘ^5+495*x^４-２20*x＾3+66*x^2—12*x
７
６
５
+ x
－ 12 x
+ 66 ｘ
— ２２0 x
+ 49５ x
x
＋ 92４ x
－ ７9２ x
+ 495 x
－ 220 ｘ
＋ 66 x
－ 12 x ④ 因式分解(fａｃｔoｒ)。將表達(dá)式做因式分解。
>> syms x;f=x^12—1； 〉> ｐrｅttｙ(fａctｏr(f））
答;aｎｓ = 1＋x^1２—12＊x^１1+66*x^1０—2２０*x^９+495*x^8－7９2＊ｘ^7+9２4*ｘ^６—7９2*x^5+4９5*x^4—22０＊x^3＋66*ｘ^2－12＊x
1２
１１
８
６
５
+ x
— １2 ｘ
＋ 66 x
— 22０ x
+ 495 x
x
+ 92４ ｘ
— 7９2 ｘ
３
+ ４95 x
— 2２０ x
＋ ６6 ｘ
－ 12 x 〉〉 sｙms x；f=x^12—1;preｔty(factor（f)）
２
２
(ｘ-1)(1 + ｘ
＋ x)(1 ＋ x）
(1 — x + x)（1 + ｘ)(x
— ｘ
+ 1)⑤化簡(simplｅ或 sｉｍplifｙ）。
將函數(shù)化簡.〉〉 cｌear 〉〉 ｓyms x;f=(1/x^3+6/x^２+１２/ｘ+８)^(１／3）； 〉＞ g１＝sｉmple(ｆ)＞〉 ｇ2＝ｓiｍplify(ｆ）
答;g1 =(２*x+１)/ｘ g2 =(（２*x＋１)＾３/x^3)＾(１/3)6 6、符號(hào)表達(dá)式得變量替換
sｕbs 函數(shù)可以對(duì)符號(hào)表達(dá)式中得符號(hào)變量進(jìn)行替換 >> ｃleaｒ >＞ f=ｓym(’(ｘ+y）^2＋４*ｘ+10’)〉〉 f1=ｓuｂs(f，’ｘ＇，’s＇)
%使用 s 替換ｘ >〉 f2=subｓ(ｆ，”x+y’,“z’)答;f =(ｘ+y)^2+4*x+１0 f1 =（s+y)^2＋4*ｓ+10 f２ = z^2+４*x+10 ７、符號(hào)極限、符號(hào)積分與微分
1）
求極限函數(shù)得調(diào)用格式 lｉmｉt(F，x,a)
％返回符號(hào)對(duì)象 F 當(dāng) x→a 時(shí)得極限 ｌｉmit（F,a)
％返回符號(hào)對(duì)象 F 當(dāng)獨(dú)立變量*→a 時(shí)得極限
ｌimiｔ（Ｆ）
%返回符號(hào)對(duì)象 F 當(dāng)獨(dú)立變量→０(a＝0)時(shí)得極限 lｉmiｔ（Ｆ，x，a,’right’)
％返回符號(hào)對(duì)象 F 當(dāng) x→a 時(shí)得右極限 limit(F,x,ａ,’ｌeft”）
%返回符號(hào)對(duì)象Ｆ當(dāng) x→a 時(shí)得左極限 例一：
〉〉 cｌear >＞ ｆ=ｓym(“sｉｎ（x)/ｘ+ａ*x”)
>〉 limit（f,“x’,0）
％以 x 為自變量求極限 >> lｉｍit（f，”ａ’，０）
%以 a 為自變量求極限 >> limit（f）
%在默認(rèn)情況下以 x 為自變量求極限 >＞ ｆｉndsym(f)
%得到變量并且按字母表順序排列 答‘f = sｉｎ(x)/x+a＊ｘ ans = 1 ans = sｉn（x)/x anｓ ＝ 1 ans = a，ｘ 例二: ＞＞ cleaｒ >〉 f＝sym(＇sｑrt（1+１/n））；
〉〉 limit(ｆ，n，iｎf)
％求 n 趨于無窮大時(shí)得極限
2）
求積分函數(shù)得調(diào)用格式 int(F)
%求符號(hào)對(duì)象 F 關(guān)于默認(rèn)變量得不定積分 ｉnt(F,ｖ）
％求符號(hào)對(duì)象 F 關(guān)于指定變量 v 得不定積分 int（Ｆ,ａ，b)
%求符號(hào)對(duì)象 F 關(guān)于默認(rèn)變量得從 a 到 b 得定積分 iｎt(F,v,a，ｂ)
％求符號(hào)對(duì)象 F 關(guān)于指定變量得從 a 到ｂ得定積分 3）
求微分方程得調(diào)用格式 diff（F)
％求符號(hào)對(duì)象 F 關(guān)于默認(rèn)變量得微分 diff(F,v）
％求符號(hào)對(duì)象 F 關(guān)于指定變量 v 得微分 diff(F,n）
%求符號(hào)對(duì)象 F 關(guān)于默認(rèn)變量得 n 階微分,n 為自然數(shù) 1、２、3…… diｆｆ(Ｆ,ｖ,n)％求符號(hào)對(duì)象 F 關(guān)于指定變量 v 得 n 階微分 8 8、符號(hào)方程求解
1）常規(guī)方程求解函數(shù)得調(diào)用格式 g=solve(eｑ)
％求方程（或表達(dá)式或字串）eq 關(guān)于默認(rèn)變量得解
g=ｓｏlve（eq，ｖar)
％求方程（或表達(dá)式或字串)eq 關(guān)于指定變量 vａr 得解 ｇ=sｏｌve（eｑ1,eq２,…、,ｅqｎ,ｖaｒ1，vａr２，…,vaｒn)
%求方程(或表達(dá)式或字串)eｑ1,eｑ２,eq3,……ｅqn 關(guān)于指定變量組ｖar１,var2,……,ｖarn）得解 求一元二次方程得解.其求解方法有多種形式：
① seq＝ｓolｖe(“a＊x^2+b*x+c＇）
② seq=solｖe(’a*x^2+b*x+c=0”)③ eq＝’a*x^2+ｂ＊x+c“;④ eｑ=”ａ*x＾2+b*x＋c=0’； sｅq=solve(eq）
⑤sym ｘ a b ｃ
eｑ＝ａ＊x^2+b*x＋ｃ ｓｅq=solve(eq）
2）常微分方程求解 求解常微分方程得函數(shù)就是 dsｏlvｅ。應(yīng)用此函數(shù)可以求得常微分方程(組)得通解,以及給定邊界條件(或初始條件)后得特解。
常微分方程求解函數(shù)得調(diào)用格式：
r=dｓolve（“ｅｑ１,ｅq2,…’,’cond１,coｎd2,…’,’v’)ｒ＝ｄsolve(’eq1’,’eq2”,…，“cond１’，’cond2’,…,’v”)說明: ① 以上兩式均可給出方程 eｑ１,、qeq2 對(duì)應(yīng)初始條件 cond1、cｏnd2 之下得一 v 作為解變量得各微分方程得解。
② 常微分方程解得默認(rèn)變量為 t。
③ 第二式中最多可接受得輸入式就是 12 個(gè)。
④ 微分方程得表達(dá)方法。
在用 MATLAB 求解常微分方程時(shí)，用大寫字母 Dｙ表示,用Ｄ2ｙ表示,依此類推。
邊界條件以類似于 y（a）＝b 給出。其中 y 為因變量，a、b 為常數(shù)、如果初始條件給得不夠，求出得解為含有 C1、C2 等待定常數(shù)得通解。
例一 求微分方程得通解、練習(xí):（1）求。
(2）求函數(shù)得積分；求函數(shù)得導(dǎo)數(shù)(3）計(jì)算定積分(4)求下列線性方程組得解
（5)求解但 y(0）=2，在 z（0)=7 時(shí)，微分方程組得解。
實(shí)驗(yàn)五
MATLA Ｂ 程序設(shè)計(jì)
一、實(shí)驗(yàn)?zāi)康?/span>
掌握 MATLAB 程序設(shè)計(jì)得主要方法,熟練編寫 MATLＡＢ函數(shù)、二、實(shí)驗(yàn)內(nèi)容
（1)M 文件得編輯。
(2)程序流程控制結(jié)構(gòu)。
（3)子函數(shù)調(diào)用與參數(shù)傳遞。
(4)局部變量與全局變量。
三、實(shí)驗(yàn)步驟1、M M 文件得編輯
選擇ＭＡTLAB 得菜單，打開新得Ｍ文件進(jìn)行編輯,然后輸入以下內(nèi)容，并保存文件名為 eｘｐl、m。
s=0;for n=１：100
ｓ＝s＋ｎ;enｄ s 答；s =
5050 保存好文件后，在命令窗口輸入ｅｘｐl 即可運(yùn)行該腳本文件，主義觀察變量空間。緊接著創(chuàng)建 M 函數(shù)文件,然后輸入以下內(nèi)容,并保存文件名為 expl2、m。
funｃtｉｏn ｓ＝expl2（x）
s＝0;fｏr n＝１:x
s＝s+n； end
保存好文件后,在命令窗口輸入
>〉 ｃｌeａr 〉> ｓ＝ｅxpl2(１０0)以 oｐen 命令可以打開 M 文件進(jìn)行修改。
〉〉open cｏnv
%打開 cｏnv 函數(shù) 2 2、程序流程控制結(jié)構(gòu)
1）
for 循環(huán)結(jié)構(gòu)
>〉 fｏｒ n=1:１0 n end 答；ｎ =n ＝n =n ＝n =n =
６ n =n =
n =n =
另一種形式得 for 循環(huán): 〉> n=1０:—１:5;>〉 ｆｏr i＝ｎ
%循環(huán)得次數(shù)為向量 n 得列數(shù) i eｎｄ 答;ｉ =
１0 i =ｉ ＝i =ｉ =i =2）
wｈile 循環(huán)結(jié)構(gòu) 在命令窗口輸入：
clear x=１;ｗhile １ x=x*2 ｅnd 將會(huì)瞧到 MATLAB 進(jìn)入死循環(huán)因?yàn)?whｉle 判斷得值恒為真，這時(shí)須按 Ctｒl+Ｃ鍵來中斷運(yùn)行，并且可瞧到 x 得值為無窮大。
練習(xí)：
(1)請(qǐng)把 exp2、ｍ函數(shù)文件用 whｉle 循環(huán)改寫。
(2）用公式求得近似值，直到最后一項(xiàng)得絕對(duì)值小于為止,試編寫Ｍ腳本文件、3）
ｉｆ－elｓe—eｎd 分支結(jié)構(gòu)
iｆ—else—end 分支結(jié)構(gòu)有如下 3 種形式。
（a)
iｆ
表達(dá)式 語句組 1 end(b)
ｉf
表達(dá)式 語句組 1 else 語句組２ ｅnd(ｃ)
if
表達(dá)式 A
語句組 1 ｅlse ｉf 表達(dá)式 B 語句組 2 ｅlｓe if 語句組 3 …… ｅlse
語句組 n end 4）
swｉtcｈ—casｅ結(jié)構(gòu) 創(chuàng)建 M 腳本文件 eｘｐ3、m,輸入以下內(nèi)容并在命令窗口中運(yùn)行。
％功能:判斷鍵盤輸入得數(shù)就是奇數(shù)還就是偶數(shù) n=iｎpｕt（’n=＇）;iｆ iseｍptｙ(n）
erｒoｒ（’plｅase ｉnput ｎ＇)；
n=ｉnput（“n＝’);enｄ ｓwitｃh mｏd(n,2）
case 1
A=’奇數(shù)＇
ｃaｓe 0
A=’偶數(shù)” Enｄ 答;ｎ＝input（＇n=’)； iｆ isempty(n)
ｅｒrｏr(＇ｐlｅaｓe ｉnｐuｔ n’)；
n＝input（“ｎ=’）； eｎd ｓwitcｈ mod(ｎ,2)
casｅ 1
A＝”奇數(shù)“
ｃａse ０
Ａ=”偶數(shù)＇ ｅｎd n=1 A = 奇數(shù) ３、子函數(shù)與參數(shù)傳遞
有一個(gè)函數(shù),試編寫實(shí)現(xiàn)該函數(shù)得函數(shù)文件.ｆuncｔion g＝expl４(ｘ)
%主函數(shù) g=0;fｏｒ n=1：ｘ
ｇ=g+faｃt（n);
%調(diào)用子函數(shù) enｄ
function y=ｆact(k)
％子函數(shù) y=1； for n＝1：k
y=y＊n； eｎd 輸入?yún)?shù)可以有函數(shù) nargin 計(jì)算,下面得例子 sｉｎplot（）,當(dāng)知輸入一個(gè)參數(shù) w 時(shí)，siｎｐlot（）函數(shù)會(huì)給ｐ賦予默認(rèn)值 0。4、局部變量與全局變量
自程序執(zhí)行開始到退出 MＡＴLAB，始終存放在工作空間,可被任何命...
第三篇：matlab實(shí)驗(yàn)報(bào)告
求解：1.模擬比賽車道曲線和選手速度曲線；
2.估計(jì)車道長度和所圍區(qū)域面積；
3.分析車道上相關(guān)路段的路面狀況（用不同顏色或不同線型標(biāo)記出來）；
4.對(duì)參加比賽選手提出合理建議。
四．合理建議
1.通過賽道曲線可知，選手所經(jīng)過的賽道上不平整的地方很多，如果平常不多多嘗試不同的路況會(huì)造成比賽時(shí)的很多突發(fā)情況，都會(huì)造成選手的成績受到很大的干擾，甚至退賽等嚴(yán)重的后果。所以我建議選手平時(shí)要多在不同類型的路況上練習(xí)，以增強(qiáng)應(yīng)變能力，取得更好的成績。
2.選手的速度分配有一些不合理，在平直的沙土路段應(yīng)該全力加速，以最快速度通過這個(gè)路段，以達(dá)到最好的比賽效果。在經(jīng)過坑洼碎石路時(shí)盡量保持一個(gè)恒定的速度，因?yàn)槿绻俣纫恢痹谧兓苋菀自谶@種路段上陷入或者熄火，造成比賽時(shí)的極大不利。將加速盡量用在沙土路等摩擦力較大的路上，以免耗費(fèi)太多的動(dòng)力。我們需要將動(dòng)力的效率盡量的提高。
五．實(shí)驗(yàn)的總結(jié)
我認(rèn)為，本實(shí)驗(yàn)的主要目的在于讓我們掌握對(duì)三次樣條差值來模擬離散點(diǎn)表示的曲線的運(yùn)用，我認(rèn)為我們已經(jīng)基本掌握，并且我們也掌握了用梯形法求不規(guī)則封閉圖形的面積的方法，除此之外，本實(shí)驗(yàn)中還有考查所學(xué)知識(shí)外的方面，就是畫v-t圖。
事實(shí)上，根據(jù)題目所給的條件，并不能準(zhǔn)確地畫出v-t圖，所以需要找到一種方法，來盡量 使結(jié)果接近真正的情況。我采用了中值的方法，這是我想到的一種比較有效的方法，并且在使用這種方法時(shí)，又運(yùn)用到了三次樣條差值的方法，使我對(duì)三次樣條差值法理解地更加深刻，并且能更加靈活地運(yùn)用。
所以我發(fā)現(xiàn)三次樣條差值的方法運(yùn)用的范圍十分廣泛，不僅是對(duì)路徑的擬合，許多已知離散點(diǎn)，對(duì)應(yīng)的函數(shù)連續(xù)變化的問題也可以用此方法解決，比如已知一天中幾個(gè)離散時(shí)間點(diǎn)對(duì)應(yīng)的氣溫，估計(jì)出一天氣溫的變化趨勢(shì)，就可以用此方法。
最重要的是，我感到了數(shù)學(xué)建模的重要性，我發(fā)現(xiàn)原來生活中不少類似的問題，都是用數(shù)學(xué)建模的方法解決的。
第四篇：matlab實(shí)驗(yàn)報(bào)告
實(shí)驗(yàn)二
特殊函數(shù)與圖形
一、實(shí)驗(yàn)?zāi)康募耙饬x
掌握用Matlab軟件繪制簡單曲線、曲面圖形，并通過繪制一些特殊函數(shù)的圖形，更加深入地理解相關(guān)函數(shù)的性質(zhì)，了解函數(shù)的性態(tài)。
二、實(shí)驗(yàn)內(nèi)容
平面作圖、空間作圖，比較數(shù)值作圖與符號(hào)作圖的異同。
三、實(shí)驗(yàn)步驟
1.在D盤建立一個(gè)自己的文件夾;
2.開啟軟件平臺(tái)——MATLAB，將你建立的文件夾加入到MATLAB的搜索路徑中。
3.4.5.6.7.利用幫助了解函數(shù)plot, surf, ezmesh，ezsurf等的功能和用法。
開啟MATLAB編輯窗口,鍵入你編寫的M文件（命令文件或函數(shù)文件）；
保存文件（注意將文件存入你自己的文件夾）并運(yùn)行；
若出現(xiàn)錯(cuò)誤，修改、運(yùn)行直到輸出正確結(jié)果；
寫出實(shí)驗(yàn)報(bào)告，并淺談學(xué)習(xí)心得體會(huì)。
四、實(shí)驗(yàn)要求與任務(wù)
根據(jù)實(shí)驗(yàn)內(nèi)容和步驟，完成以下具體實(shí)驗(yàn)，按要求寫出實(shí)驗(yàn)報(bào)告。1．作出下圖所示的三維圖形：
圖1 提示：圖形為圓環(huán)面和球面的組合.2．作出下圖所示的墨西哥帽子及其剪裁圖形：
圖2
3．畫出橢球面、雙葉雙曲面、單葉雙曲面．
4．若要求田螺線的一條軸截面的曲邊是一條拋物線：y?0時(shí)x2?5z．試重新設(shè)計(jì)田螺線的參數(shù)方程，并畫出該田螺線．
5．作出下圖所示的馬鞍面（顏色為灰色，并有一個(gè)標(biāo)題：“馬鞍面”）：
圖3 6．繪制黎曼函數(shù)的圖形，要求分母的最大值n的數(shù)值由鍵盤輸入（提示：使用input語句）．
第五篇：上機(jī)實(shí)驗(yàn)報(bào)告格式
請(qǐng)各位同學(xué)按照下列格式寫：否則按照0分記
上機(jī)實(shí)驗(yàn)報(bào)告格式
Matlab學(xué)習(xí)第？次上機(jī)實(shí)驗(yàn)報(bào)告（第？次）
姓名：？？？班級(jí):???上機(jī)時(shí)間：？？？？？？
1.上機(jī)內(nèi)容
2.建模及算法分析
3.程序
4.小結(jié)
5.參考文獻(xiàn)
將文檔發(fā)到我的郵箱里面
learn_matlab@163.com
點(diǎn)擊咨詢 