步驟1：啟動Matlab

第一步是啟動并運行matlab，以便我們可以開始使用它。初次啟動matlab時，它應類似于以下屏幕截圖。第一步是為matlab分配一個工作目錄。這是程序?qū)闹刑崛∷形募奈恢茫彩菓４嫠衜atlab工作的位置。

我建議在一個您會記住它的位置創(chuàng)建一個新文件夾，并為它命名，以便您識別。創(chuàng)建新文件夾后，單擊屏幕右上方的“ 。..”，如第二張圖片中突出顯示的那樣。如第三張圖片所示，這將彈出一個瀏覽框。找到您在計算機上創(chuàng)建的新文件夾，然后選擇它。對于此示例，該文件名為“ 370”，位于桌面上。

步驟2：創(chuàng)建M文件

現(xiàn)在我們需要做的是創(chuàng)建一個新的M文件。 M文件的功能與直接在matlab中鍵入代碼完全一樣，但是您可以保存和修改代碼，然后重復運行。將代碼直接輸入到matlab中時，您分別鍵入每一行代碼。在M文件中，您可以編寫整個代碼，然后立即運行它。

要打開新的M文件，請單擊文件。將光標放在“新建”上，然后單擊“空白M文件”，如第一張圖片所示。打開的內(nèi)容應該看起來像第二張圖片。

由于此代碼可以重復運行，因此最好在每次運行之前關(guān)閉所有內(nèi)容并清除所有變量。這是通過兩行代碼完成的：

關(guān)閉所有

清除所有

如第三張圖所示，確保清除并關(guān)閉了所有內(nèi)容。

步驟3：創(chuàng)建時間向量

我們要做的第一件事要做的是在matlab中創(chuàng)建函數(shù)圖。第一步是創(chuàng)建獨立變量。在這種情況下，我們稱其為“ t”。我們將用于創(chuàng)建此變量的方法是創(chuàng)建向量。向量基本上是一系列數(shù)字。例如，1，2，3，4將是一個短向量。創(chuàng)建此向量的代碼是：

t = 0.1：0.01：10;

第一個數(shù)字0.1表示起點。第二個數(shù)字0.01是步長。第三個數(shù)字10表示終點。因此，此向量一直對應于0.1，0.11，0.12 。..一直到10。

要查看創(chuàng)建向量是否起作用，請單擊第二張圖中突出顯示的綠色運行按鈕。這將運行程序。要查看我們的矢量，請轉(zhuǎn)到主Matlab窗口。單擊桌面，然后將鼠標懸停在桌面布局上，然后單擊第三張圖片中概述的默認。現(xiàn)在您的屏幕應該看起來像第四張圖片。在右側(cè)，您將看到我們新創(chuàng)建的變量t。雙擊它，就像在第五張圖片中一樣，您將看到創(chuàng)建的一系列數(shù)字。

步驟4：運行和繪制函數(shù)圖

現(xiàn)在，我們將繪制在matlab中創(chuàng)建的函數(shù)的圖形。第一步是創(chuàng)建函數(shù)。這就像寫出所需的數(shù)學函數(shù)一樣簡單。第一幅圖中顯示了一個示例。用于此功能的代碼為：

y = sin（t）+ 4 * cos（5. * t）。^ 2;

在余弦中乘法之前的時間，并且在余弦的平方告訴Matlab之前，僅對時間向量的有價物執(zhí)行這些功能，而不是將時間向量視為矩陣并嘗試對其進行矩陣函數(shù)。

下一步是創(chuàng)建圖本身。這是使用第二幅圖中所示的代碼完成的。 plot命令中變量的順序非常重要，因此請確保按照下面的設(shè)置來設(shè)置代碼。

圖

h = axes（‘fontsize’，14）;

plot（t，y，‘linewidth，2）

xlabel（’Time（s）‘） ylabel（’Y Value‘）

Title（’Y Value vs Time‘）

網(wǎng)格上

最后，只需再次單擊綠色的運行箭頭，圖形便會像第三張圖片一樣彈出。

步驟5：從Excel中提取數(shù)據(jù)

我們現(xiàn)在將創(chuàng)建與之前相同的圖形，但是要從中導入函數(shù)數(shù)據(jù)一個Excel電子表格。第一張圖片是將要使用的excel電子表格的屏幕截圖。它與在上一步中在matlab中創(chuàng)建的數(shù)據(jù)點完全相同，只是在excel中創(chuàng)建的。

首先，我們可以從上一步中刪除創(chuàng)建時間向量的代碼和函數(shù)的代碼。現(xiàn)在，您的代碼應該看起來像第二張圖片。

插入代碼，如第三張圖片的頂部紅色框所示。這是讀取excel文件的代碼。 “ A”是指將包括電子表格中所有數(shù)字的矩陣，而“ B”包括電子表格中所有的文本。 t和y變量從第一個和第二個列中提取，如代碼所示。

［A，B］ = xlsread（’excelexample.xlsx‘）;

t = A（：，1 ）;

y = A（：，2）;

還可以修改圖形代碼，如第三張圖片的下部紅色框所示。實際上，這將從電子表格中提取圖表標題和軸標簽，并將它們放在圖形上。

xlabel（B（2））

ylabel（B（3））

Title（B（ 1））

要做的最后一件事是再次運行該程序，您將看到彈出與最終圖片相同的圖形。

步驟6：創(chuàng)建頻譜圖

在這一步中，我們將使用matlab通過讀取wav聲音文件來創(chuàng)建頻譜圖。頻譜圖有時稱為“ 2.5D圖形”，因為它使用二維圖形，并附加了顏色以顯示振幅。顏色提供了比簡單的2D圖形更多的細節(jié)，但沒有提供3D圖形的細節(jié)，因此提供了術(shù)語“ 2.5D”。matlab的specgram函數(shù)從wav文件中獲取一組數(shù)據(jù)點并執(zhí)行傅立葉對這些點進行變換以確定信號中存在的頻率。對于這個可指導的人員，了解傅里葉變換的工作原理，知道頻譜圖將顯示存在的頻率以及它們相對于時間的強度并不重要。該函數(shù)在X軸上繪制時間，在Y軸上繪制頻率。每個頻率的強度通過顏色顯示。

在這種情況下，wav文件是敲擊一塊金屬的聲音記錄，然后將金屬的振動記錄為聲音。使用該頻譜圖，我們可以輕松確定金屬片的共振頻率，因為該共振頻率將隨時間而持久。

要執(zhí)行此任務，請首先讓matlab使用以下代碼讀取wav文件：

［x，fs］ = wavread（’flex4.wav‘）;

在這種情況下，flex4.wav是我們wav文件的標題，變量x是數(shù)據(jù)點在文件中，fs表示采樣頻率。要執(zhí)行頻譜圖，只需鍵入以下代碼：

specgram ［x（：。1），256，fs］;

256對應于分析數(shù)據(jù)時執(zhí)行FFT的頻率。 Matlab基本上是將聲音文件切成塊，并對每個塊進行FFT。256告訴它每個塊應該有多大。細節(jié)并不重要，對于大多數(shù)應用程序來說，256是一個安全值。現(xiàn)在，如果您運行代碼，您將看到一個彈出的圖，如第二張圖所示。由此很容易看出，諧振頻率對應于該圖右下角的紅色峰值。這是相對于時間而言最長的峰。
責任編輯：wv