LTI衝激響應卷積實驗 附加音效

2021-10-04 02:23:42 字數 1833 閱讀 9934

這一部分內容主要以學習截圖為主。方便自己回顧,給大家推薦乙個非常優秀的b站uper,所有內容均出自他的課程。**鏈結分享在此dr_can,ph.d. in dynamics and control 機械人工程師。

本實驗是關於筆記第三節現代控制理論-工程數學基礎(3)-lti衝激響應卷積。

利用家裡條件,在洗手間模擬衝激響應環境。嘗試使用木棍,鏟子及鐵具相互碰撞,對比實驗。

具體步驟如下:

1、選取一段原聲**

2、用手機錄下衝激響應

3、用處理軟體去除部分雜音

4、在matlab中跑以下**

**如下:

close all

clear

clc;[y1

,fs1]

=audioread

('zhuimeng.***');

%獲取音訊 1 檔案,並返回取樣資料到向量 y中, fs 表示取樣頻率[y2

,fs2]

=audioread

('tc.***');

%獲取音訊 1 檔案,並返回取樣資料到向量 y中, fs 表示取樣頻率

ft1 =y1(

:,1)

;siglength1 =

length

(ft1)

;%獲取聲音長度

t1=(

0:siglength1-1)

/fs1;

%求出音訊 1 對應的時間座標

subplot(3

,1,1

);plot

(t1,ft1)

;%畫出音訊 1 的波形圖

title

(' 追夢人的波形圖 ');

xlabel

('時間(s)');

ylabel

(' 振幅 ');

grid;

ft2 =y2(

:,1)

;siglength2 =

length

(ft2)

;%獲取聲音長度

t2=(

0:siglength2-1)

/fs2;

%求出音訊 2 對應的時間座標

subplot(3

,1,2

);plot

(t2, ft2)

;%畫出音訊 2 的波形圖

title

(' 衝擊響應的波形圖 ');

xlabel

('時間(s)');

ylabel

(' 振幅 ');

grid;

ft=conv

(ft1,ft2)

;%將兩段音訊進行卷積

siglength =

length

(ft)

;%獲取聲音長度t=(

0:siglength-1)

/fs1;

subplot(3

,1,3

);plot

(t,ft)

;%畫出卷積後的波形圖

title

('卷積後的波形圖 ');

xlabel

('時間(s)');

ylabel

(' 振幅 ');

圖上和**只對一種情況進行分析。

筆者對比三種物體相互碰撞產生的衝激響應,及卷積情況,排除誤差外,卷積後的結果基本上相同。

致敬dr_can!畢。

衝激串函式的傅利葉變換還是衝激串函式!?

我們都知道,衝激函式的傅利葉變換是乙個單頻的三角函式,所以從書上看到 衝激串函式的傅利葉變換還是衝激串函式 這個結論的時候非常震驚。從直覺上看,衝激串函式的傅利葉變換應該是一串三角函式,知乎上也有乙個發出了相同的疑問,但是以 不構成疑問 的方式被關閉了。那麼到底是怎麼回事呢?我們先來觀察一下這四個式...

全自動LTI部署OS

全自動lti部署os 零 通過adk製作winpe 需包含有imagex.exe工具,用來捕獲映像 一 使用winpe中的imagex捕獲映像 install.wim 二 使用mdt製作啟動映像 boot.wim 三 通過wds批量部署 一 捕獲wim映像 install.wim 1.準備模板機 安...

VHDL脈衝雷射測距

vhdl 這是乙個雷射測距機中接收模組和顯示模組之間的小模組 求解這段 的含義和作用 an highlighted block library ieee use ieee std logic 1164 all entity counter isport st clk,tr rec in std lo...