李江龙的博客

冲击函数只有在时候才有非零值，其余时间值均为零，对于时刻来讲，系统并没有任何输入，因此系统方程的受迫响应部分为零，只有自由响应，也就是说系统的齐次解就是系统的冲击响应完全解。然而求系统方程的齐次解，需要知道系统方程的边界条件（初始状态），因此通过方程以及冲激函数确定方程的初始状态是第一步，将齐次解的形式解带入方程求齐次解的系数是第二步，求出后的结果就是冲激响应。第一部分】冲激函数匹配法求初...

【第一步】列出微分方程的特征方程【第二步】求取特征方程的特征根【第三步】根据特征根写出微分方程的通解【第四步】写出微分方程的特解【第五步】写出微分方程的完全响应表达式【第六步】将完全解分别代入微分方程以及微分方程的初始条件，得到有关于完全响应系数的线性方程组【第七步】求解该线性方程组，代入完全响应，得到微分方程的完全响应PS：微分方程的通解部分与激励信号...

记泰勒公式舍去高阶项在导数为0时候取得极值，令，以上式子可以变形为继续整理可以得到通过观察上式，我们可以发现，当给出一个初始点a的时候，我们可以通过泰勒公式用a点的一阶导数与二阶导数还有a的值近似的表达出使的点x。但是我们的这个近似是通过舍去泰勒公式中的高阶项实现的，因此我们要保证高阶项的舍弃对x的求值影响较小，也就是说初始点a的取值不能离x点太远，否则会出现较大...

用最速下降法求函数的极小值，迭代两次，计算各迭代点的函数值、梯度及其模，设初始点。以下代码在求解的时候并没有严格按照题目中要求的迭代两次进行计算，而是设置了最大容许误差e，当迭代过程中g的范数小于e时候停止迭代，最终迭代了五次，如果按照题目中仅迭代两次的话误差就会稍微大一些。clcclear allclose all% 定义允许误差e = 0.1e = 0.1;Q = [2 ...

拉格朗日乘子法的目的就是将有等式约束极值问题转化为无约束非线性优化问题，通常这种无约束极值问题我们称之为经典极值问题。这类问题通常有参数拟合（包括阵列中的MVDR最优系数加权系数）、选址问题等。可行域内极值点称为全局最优解，邻域内极值点称为局部最优解。一阶导数 = 梯度；二阶导数 = 海塞矩阵（hesse）一阶导数是列向量，二阶导数是一阶倒数转置再求导，每一个函数求导后还是按...

1、生产组织与计划问题某电视机厂生产I、II、III型号的电视机，按照市场需求，三种型号的电视机每天产量最少为200、250、100，该厂每天可以利用的总工时为1000个时间单位（小时），可利用的材料每天有2000个单位，每个产品工时消耗为2、1.2、1，原材料消耗为1、1.5、4，每台电视产生的利润分别为10、14、12，问每种型号的产量如何调配能够使工厂的利润获得最大化？2、解答2...

某工厂生产A\B\C三种型号的产品，三种型号产品收益分别为2、1、3元，其生产受下列条件限制：人工<=1,材料<=3，每种产品的人工消耗分别为1/3、1/3、1/3；每种产品的用料消耗分别为1/3、4/3、7/3，试决定该生产的最优方案。按照单纯形法，计算如下：clcclear allclose allsyms x1 x2 x3 y1 y2 y3 y4 y5; %...

IIR滤波器还是FIR滤波器?1】当需要线性相位（h(n)关于M/2对称或反对称，M是滤波器长度）时候，考虑FIR滤波器2】不需要线性相位时候，优先选择IIR滤波器；IIR滤波器有更短低的旁瓣非因果滤波器系统不具有物理可实现性1】∫ln(|H(w)|)dw<∞，那么系统就是因果的（帕雷-维纳定理）2】在时域内，如果h（n） = 0，当n<0时候，系统就是因果的以...

如果对于自变量x1、x2以及参数λ，有则认为f是凸函数，进一步，如果则认为f是严格凸函数。R向量空间中，如果集合 S 中任两点的连线上的点都在 S 内，则称集合 S 为凸集。在凸集范围内，求凸函数在约束条件下的最小值成为凸规划,即：s.t 如果C∈R^n向量空间，对于任意x∈C，有ax∈C，则称C为锥集（类似于空间绝对可...

本文使用协方差矩阵特征值分解的方法求解声音数目，对比两种主特征值数目的判定方法，代码如下clcclearclose allframeLength = 160;audioDataArray = zeros(frameLength,5);audioDataArray2 = zeros(frameLength,5);signalSourceInfo = zeros(4,1);sign...

从今天起我们开始学习算法开发的单元测试框架使用方法，我们以之前写的求解一个范围内素数个数的算法为例，介绍单元测试框架的基本使用方法。先将之前的代码写成function的形式，如下：function [divNumberVectorLength] = primeNumberSolute(calculationRange) numEnter = calculationRange;if...

本文介绍如何使用matlab自带的高阶累积量函数求取一个随机过程的高阶累积量，运行demo之前确保matlab中已经安装了HOSA的工具箱。安装过程中经常会碰到的两个问题：1、info.xml出问题，在安装的目录下找到这个文件，将里面的<area>改成<type>，将</area>改成</type>2、找不到定义的函数cum4est.m，...

本文展示如何绘制一个高斯分布过程的第一特征函数&amp;第二特征函数clcclearclose allu = 3;a = 0.5;Nr = 100;Nc = 1000;w = 1:Nr;x = random('norm',u,a,1,Nc);Q = zeros(Nr,Nc);Qabs = zeros(Nr,Nc);for iloop1 = 1:Nc Q(:...

会了第二种卷积手算方法，可以递推圆周卷积方法，原理是一样的的，只是延时变成了循环。同理的根据相关性的定义，将其中一个序列反向，然后按照本文这种方法进行求解，就是手算序列相关性了。对如下x、y两个序列进行卷积计算，根据定义可以知道卷积结果长度为 5首先，看第一种计算方法，根据定义式先将y序列翻转，然后逐位相乘第二种方法相对来说比较简单，转到z域后，卷积变为多项式相乘...

离散信号重采样（采样率转换）通常有两种实现办法1】一种是通过DA转换成模拟信号然后再通过AD重新对模拟信号进行采样2】一种是通过插值的办法在对输入信号进行插值或者抽取的方式，完成纯数字域的采样率转换本文主要记录第二种方法1.通过插值公式实现对原始信号的重构观察以上公式，gt函数中sin(pi*t/T)的周期为1/2*Fx=1/(2T)，即gt频率是原信号采样频率的1/2...

有期望信号d（n），纯净信号x（n），以及噪声信号g（n）;有滤波器h（m）,以及滤波器输出信号y（n），滤波器输出纯净信号x的估计值y，并且使y与x误差的期望为0. （相关性）根据最小均方差准则设计滤波器,滤波器设计问题转为函数求驻点的问题，在驻点处能够取得极值，如果能证明该式有连续二阶偏导，且得到的hesse矩阵正定...