領會跟做第一個MATLAB程序

1、李文東中國海洋大學物理學2015.03.16第第2 2章章 領會、跟做第一個領會、跟做第一個MATLABMATLAB程序程序1 課程概述2 領會、跟做第一個MATLAB程序3 三種代碼執(zhí)行方式 （第6章）4 矩陣分析 （第3章）5 數值計算 （第4章）6 程序設計 （第6章）7 圖形用戶界面（GUI） （第8章）8 符號代數 （第2章）9 Simulink仿真 （第7章）10 數字圖像處理圖形顯示文件I/O數據文件矩陣工作空間查看繪圖處理查看繪圖處理后數據數據保存圖形保存5文件I/O符號計算GUI28矩陣變化數值分析程序設計工具箱6347數據類型圖形顯示文件I/O數據類型數據文件矩陣工作空間查

2、看繪圖處理查看繪圖處理后數據數據保存圖形保存5文件I/O符號計算GUI28矩陣變化數值分析程序設計工具箱63472.1 2.1 數據數據讀取與讀取與查看查看文件I/O數據文件矩陣工作空間查看繪圖處理查看繪圖處理后數據數據保存圖形保存圖形顯示5文件I/O符號計算GUI28矩陣變化數值分析程序設計工具箱63472.2 2.2 數據處理與分析數據處理與分析數據類型實際問題：一個振動測出（時間，位置）坐標，想知道加速度隨時間變化規(guī)律？數據：tx10.99219820.96891230.93050840.87758350.81096360.73168970.64099780.54030290.43117

3、7990.9817441000.997798文件I/O數據文件矩陣工作空間文件存放位置文件I/O數據文件矩陣工作空間數據文件看到了，數據怎么讀進來？文件I/O數據文件矩陣工作空間help xlsread參考書參考書：清空清空Command Windowdata=xlsread(data_SH.xls);沒有反應？沒有反應？理解：data=xlsread(data_SH.xls)文件I/O數據類型數據文件矩陣工作空間繪圖處理查看工作空間雙擊可編輯數據！以查看為主，以查看為主，編輯慎重！編輯慎重！y1=data(1,1)y2=data(1,1);分號表示不顯示結果。分號表示不顯示結果。y3=dat

4、a(2,3)y4=data(3,2)y5=data(2,1:2)y6=data(1:3,2)y7=data(1:2:6,1)矩陣處理查看繪圖查看查看help plot關閉矩陣處理查看繪圖txt=data(:,1)x=data(:,2)plot(t,x)help：處理查看繪圖圖形保存0102030405060708090100-1-0.8-0.6-0.4-0.200.20.40.60.81t/sx/cm振 動 位 置 隨 時 間 變 化 關 系 圖嘗試保存成bmp、jpg，對比效果圖形保存查看繪圖處理矩陣圖形顯示文件I/O數據類型數據文件矩陣工作空間查看繪圖處理查看繪圖處理后數據數據保存圖形保存

5、5文件I/O矩陣變化數值分析程序設計工具箱63472.1 2.1 數據數據讀取與讀取與查看查看2.2 2.2 數據處理與分析數據處理與分析圖形保存文件I/O數據類型數據文件矩陣工作空間查看繪圖處理查看繪圖處理后數據數據保存圖形保存矩陣理論數值計算程序設計63工具箱4圖形顯示5文件I/O7任務：一個振動測出（時間，位置）坐標，想知道加速度隨時間變化規(guī)律？矩陣工作空間處理處理后數據tx10.99219820.96891230.93050840.87758350.81096360.73168970.64099780.54030290.431177990.9817441000.997798ta1234

6、5678999100矩陣變化數值分析程序設計工具箱6347計算N階微分1009810098plot(t(1:98),a)100980102030405060708090100-0.02-0.015-0.01-0.00500.0050.010.0150.02t/(s)a/(cm/s2)振 動 加 速 度 隨 時 間 變 化 關 系 圖保存emf圖片處理查看繪圖處理后數據數據保存圖形保存文件I/O查看保存：t-a 坐標保存到excel表格中，t一列，a一列查看保存：t-a 坐標保存到excel表格中，t一列，a一列t-a 坐標合并為一個矩陣data_ta矩陣data_ta存儲到excel表格中為何

7、報錯為何報錯保存：t-a 坐標保存到excel表格中，t一列，a一列t-a 坐標合并為一個矩陣data_tat-a 坐標合并為一個矩陣data_ta第第1列列第第2列列參考書參考書沒有反應？矩陣data_ta存儲到excel表格中data文件I/O數據類型數據文件矩陣工作空間查看繪圖處理查看繪圖處理后數據數據保存圖形保存矩陣變化數值分析程序設計工具箱圖形顯示5文件I/Odata=xlsread(data_SH.xls);dataplot(t,x)t=data(:,1);x=data(:,2);0102030405060708090100-1-0.8-0.6-0.4-0.200.20.40.60

8、.81t/sx/cm振 動 位 置 隨 時 間 變 化 關 系 圖a=diff(x,2);0102030405060708090100-0.02-0.015-0.01-0.00500.0050.010.0150.02t/(s)a/(cm/s2)振 動 加 速 度 隨 時 間 變 化 關 系 圖t, adata_ta(:,1)=t(1:98);data_ta(:,2)=a;xlswrite(data_ta.xls,data_ta)plot(t(1:98),a)6347*讀寫讀寫excel:-xlsread xlswrite*數據查看：數據查看：-data(1,1) data(3,2) data(1:2,1) data(1:2:6,2)*微分：微分：-diff*繪圖：繪圖：-plot xlabel ylabel title*幫助：幫助：-help作業(yè)：一個振動測出（時間，位置）坐標，速度隨時間變化規(guī)律又如何？數據：tx10.99219820.96891230.93050840.87758350.81096360.73168970.64099780.54030290.431177990.9817441000.997798已知1