這一部分內容主要以學習截圖為主。方便自己回顧,給大家推薦乙個非常優秀的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...