原创 線性規劃問題的靈敏度分析

某工廠生產A\B\C三種型號的產品,三種型號產品收益分別爲2、1、3元,其生產受下列條件限制:人工<=1,材料<=3,每種產品的人工消耗分別爲1/3、1/3、1/3;每種產品的用料消耗分別爲1/3、4/3、7/3,試決定該生產的最優方案。

原创 單純形法求解線性規劃問題

1、生產組織與計劃問題 某電視機廠生產I、II、III型號的電視機,按照市場需求,三種型號的電視機每天產量最少爲200、250、100,該廠每天可以利用的總工時爲1000個時間單位(小時),可利用的材料每天有2000個單位,每個產品工時消

原创 利用SVM實現噪聲分類

通過支持向量機實現對噪聲的分類,提升語音增強以及VAD的效果。 當樣例是以輸入/輸出對的形式給出時,學習方法稱爲有監督學習。 當樣例僅以輸入的形式給出,而沒有輸出結果時候,成爲無監督學習,包括密度估計、分佈類型學習、聚類等領域的學習。 另

原创 端到端音頻信號處理流程

本篇博客主要總結端到端語音信號處理流程,對於單獨每一個節點的概述將會在後續的時間裏逐步的更新。 常規端到端信號處理流程  

原创 數字濾波器的設計問題

IIR濾波器還是FIR濾波器? 1】當需要線性相位(h(n)關於M/2對稱或反對稱,M是濾波器長度)時候,考慮FIR濾波器 2】不需要線性相位時候,優先選擇IIR濾波器;IIR濾波器有更短低的旁瓣 非因果濾波器系統不具有物理可實現性 1】

原创 幾個基本的matlab指令

化最簡行階梯矩陣:rref >> A = magic(10) A = 92 99 1 8 15 67 74 51 58 40 98 80 7

原创 通過放置零極點設計濾波器

用取Z變換在單位圓上界面就是傅里葉變換的思想,通過放置零極點位置設計濾波器。傳輸函數0Hz位置的放大倍數通過zp2tf的第三個參數k設置,對應增益20*log10(k)。 %% 本程序DSP 243頁圖5.4.2中通過放置零極點設計低通

原创 寫在轉行一週年

轉行的前三個月花了些時間瞭解了陣列的技術棧,最後選了一條嵌入式路線實現。這兩個月中主要學習了單片機,做了一些基本的控制功能,刷完了幾本C語言的教程與語音處理的算法教程。因爲沒什麼c編程與語音算法基礎,這一段起步的起的非常困難,學習效率相對

原创 凸規劃問題與二階錐規劃

如果對於自變量x1、x2以及參數λ,有 則認爲f是凸函數,進一步,如果 則認爲f是嚴格凸函數。 R向量空間中,如果集合 S 中任兩點的連線上的點都在 S 內,則稱集合 S 爲凸集。   在凸集範圍內,求凸函數在約束條件下的最小值成爲凸

原创 貝塞爾函數法實現束控加權

使用貝塞爾函數法對各個頻率點的權值進行計算調整,實現恆定波束的形成,算法的效果不是很好,應該是算法的應用場景上有限制。本質上這就是一種時域內通過濾波器組實現恆定波束的方法。 驗證代碼如下 %% 本程序延時如何使用bessel法實現恆定波

原创 空域的匹配濾波

對於一個固定的陣型的陣列來說,理想遠場情況下,聲源處發出聲音到MIC收到聲音信號之間存在延遲,如果陣元位置位於相同的位置,則陣元與陣元之間的延時差就等於零,此時如果將聲源按照固定的距離爲半徑繞陣元位置一週,並掃描各個方向上信號的累加值,可

原创 使用matlab進行算法開發時的單元測試寫法

從今天起我們開始學習算法開發的單元測試框架使用方法,我們以之前寫的求解一個範圍內素數個數的算法爲例,介紹單元測試框架的基本使用方法。 先將之前的代碼寫成function的形式,如下: function [divNumberVectorL

原创 第一特徵函數&第二特徵函數

本文展示如何繪製一個高斯分佈過程的第一特徵函數&第二特徵函數 clc clear close all u = 3; a = 0.5; Nr = 100; Nc = 1000; w = 1:Nr; x = random('norm',u

原创 聲源數目測定

本文使用協方差矩陣特徵值分解的方法求解聲音數目,對比兩種主特徵值數目的判定方法,代碼如下 clc clear close all frameLength = 160; audioDataArray = zeros(frameLength

原创 高斯分佈估計子的性能與克拉默勞下界的討論

本文演示估計子的定義以及評估方法 % 假設一個N長度的高斯隨機過程x,按照平均數u、方差a分佈,然後根據隨機過程估計他的均 % 值u1,並對評估的有效性進行檢驗 clc clear close all %% 生成一個隨機變量,並將