- 博客(44)
- 资源 (17)
- 收藏
- 关注
RSS订阅原创 如何使用ipopt进行非线性约束求目标函数最小值（NLP非线性规划）内点法（Interior point method）
线性规划LP，二次规划（quadratic programming），如果要求更一般的非线性规划IPOT是个很好的选择，求解器很多，根据情况自己选择。eval_h: 计算目标函数和约束条件的二阶导数（即Hessian矩阵，二阶偏导数）。eval_jac_g: 计算约束条件的雅可比矩阵（两个约束条件的一阶偏导数）目标的黑森矩阵，注意走位，注意骚走位，注意下三角阵骚走位。n 是变量个数，m是约束条件个数，g是具体的约束函数。4、eval_jac_g 约束函数的jacobi矩阵。
2024-10-07 18:15:47
1159
原创 webrtc ns 降噪之粉红噪声参数推导
哇，我们分析了几乎所有的color noise，模型都是可以满足的，那么其实 webrtc种的pinknoise 参数估计应该叫color noise 估计，对吗？其实就是 S = C / f ^alpha , 其中f 和i 有个比率关系 f = i * fs/N , 常数比例，你可以认为，估算的时候全部塞到 C里面了。这两个参数是我们需要估算的， alpha 是零的时候退化成白噪声。恰好是 1/f 的，（狭义的粉红是正比 1/f 系数是1,）所以大多数是安装这个来的，这个更通用，（包含白噪声的情况）
2024-08-19 16:45:50
719
原创 传统降噪的原理
200帧以内的，应该用那个更新最慢的第三组lq，200帧以后的开始看谁先到200，用谁的？应该是 3 2 1 3 2 1 这样，第三组应该先到200，清零，第二组快到200了，再是第一组。噪声估计比较多的选择的，常用的有Martin的MS，还有Cohen的IMCRA，还有Gerkmann的SPP等方法。超过200，会强制清零，那么就是三个 counter，谁先达到200，就用谁的lquantile 的 offset。这样就是为了提高效率，按大 ，中，小，三个更新速率去接近magn，
2024-08-19 14:35:44
898
1
原创 利用GPT 将 matlab 内置 bwlookup 函数转C
这给地方还有一个加1，c代码中求和一样的code，没有加1，原因是这里做index的时候，c是从0开始，index本身比matlab就要少1。它说的很清楚，BW是二值图像矩阵，lut是一个512的查找表，对3*3邻居矩阵进行一些滤波（权重和）得到最后结果。这样计算出来发现还不对，仔细对比，中心数据都对上了，发现还是边缘不对。这个函数没有现成的m文件参考，内置已经打成库了，所以没有参考源代码。matlab中的原矩阵需要扩展，扩展的地方都要填0，然后只取中间。显然它给的矩阵滤波器顺序搞错了，它的顺序是。
2024-07-03 17:43:27
549
原创 如何使用Matlab进行三角剖分（自定义函数实现delaunayTriangulation 使用Bowyer-Watson 算法）
前一篇文章《三角剖分插值》，讲的是使用剖分后的结果进行不规则的散点插值，这篇文章主要是讲如何形成剖分。
2024-05-23 20:56:03
2078
3
原创 Matlab调C/C++简单模板例子
如果你是需要快速搭建一个matlab调c/c++环境，这篇文章可以参考。plhs[0] plhs[1]是输出，需要多少输出变量，那么就申请多少。所有的数字，数组，二位数字，到c这边都是一维数组，且按列优先访问。1、Matlab里面所有变量都是矩阵，包括单变量也是1*1的矩阵。2、Maltab矩阵按列优先访问，这个和fortran保持一致。有了c代码，想在matlab里面调用，可以参考我这个模板。
2024-04-09 15:13:21
337
原创 使用遗传算法解决微分方程的参数估算问题
我们常常使用GA算法得到某个目标函数的最小值，对于解决一般的方程，使用起来问题不大，但现在有个微分方程组的参数需要估算一下，今天给出个例子来演示一下如何使用ga来解决ode方程的参数估算问题。略。
2024-02-13 07:04:51
783
原创 Matlab三角剖分插值问题分析（二）
例如：随着人工智能的不断发展，机器学习这门技术也越来越重要，很多人都开启了学习机器学习，本文就介绍了机器学习的基础内容。提示：以下是本篇文章正文内容，下面案例可供参考例如：以上就是今天要讲的内容，本文仅仅简单介绍了pandas的使用，而pandas提供了大量能使我们快速便捷地处理数据的函数和方法。
2023-11-23 18:48:52
1953
原创 matlab二维曲面散点图插值方法
它将散点数据进行三角剖分，并使用插值方法对每个三角形进行插值计算，从而生成二维曲面。它支持多种插值方法，包括三次样条插值、最近邻插值、线性插值和自然邻近法插值。类创建一个插值函数对象，需要提供散点数据的坐标和值，然后可以使用该对象对新的点进行插值计算。函数需要提供散点数据的坐标和值，以及插值点的坐标，然后在插值点上生成插值结果。
2023-11-14 20:24:58
1659
原创 用ode45解一个带有积分的微分方程（integro-differential equations）
用ode45解一个带有积分的微分方程（integro-differential equations）
2023-01-11 18:43:57
1116
原创 Matlab三维数据区域显示问题
对meshgrid的数据挨个判断，如果 在xy平面的投影点 在区域内，那么保留，否则置为NaN。求K = 25456000*(1-4*(z/1000)^2)^(1/2);2、把范围之外的数据都设为NaN，mesh surf等直接就不显示了。K是 d1和d2的应变量，也间接是 x和z的二元函数。如果我们只想要交叠部分的mesh显示，该怎么办？我们有了 投影区域，还有meshgrid的每个点。先投影到二维，找到二维封闭曲线的边界。是一个有范围的空间曲面。1、需要知道交叠区域的范围。一定区域的三维数据，比如。
2022-08-23 11:35:39
1394
原创 使用matlab辨识工具来估算震动系统的传递函数
Matlab辨识工具箱里面的函数很多，这里用一个简单的例子来展示估算系统传递传递函数的使用方法。很遗憾，我们也不具体了解这个未知系统的 零极点个数（或者说多项式的阶数）系统未知，我们通过 功率谱估算 来求解系统响应。估算的系统响应，已经把两个峰值的特性 仿真出来了。互功率谱和自功率谱相除就得到了系统频响。title('传输tf幅度 小范围');title('传输tf相位 小范围');我们得到系统频响（幅度和相位）我们的系统输入输出数据如下。下面估算系统的传递函数。得到的 频响是一致的。...
2022-08-12 11:22:47
2665
原创 用好blas和lapack库，dtrsm举例
最近要使用 dtrsm 进行解方程，试用了一两个例子，但总是有误，归根结底还是对fortran的一些性质不熟悉，这里做一下记录1、 XTRSM(SIDE, UPLD, TRANSA, DIAG, M, N, ALPHA, A, LDA, B, LDB)主要用来解方程，而且花样很多op( A )*X = alpha*B, or X*op( A ) = alpha*B, op( A ) = A or op( A ) = A**T.我使用了别人x学习BLAS库 -- DTR...
2021-11-15 07:16:36
1215
原创 最基本的麦克风阵列延时累加和算法，DelayAndSum beamformer
最近从 phased.TimeDelayBeamformer 扒了它的代码，简单改造作为最基本的麦克风阵列算法学习refer to https://www.mathworks.com/help/phased/ref/phased.timedelaybeamformer-system-object.html一、延时累加和beamformer可以说是最基本的麦克风阵列算法，方法很简单，就是通过平移每个通道数据到合适的位置，然后累计调整系数即可。但问题是什么？一般采集到的麦克风数据..
2021-07-30 12:03:53
6750
3
原创 使用ffmpeg维修视频音量小的问题
不知道从哪里下来了一个视频，奇怪的是播放时音量非常小，使用增加音量的方法也不见得有效ffmpeg -i 01.mp4 -filter:a "volume=10dB" output.mp4仔细研究了一下，发现音频部分是 5.1声道的在PC端播放的时候应该进行了内部的downmix到 2声道stereo，但怎么downmix，很有可能和当前使用的媒体播放器有关系。想办法使用ffmpeg 修改 增加音量吧1、首先从视频剥离音频ffmpeg.exe -y -i "" -vn -sn ..
2021-06-27 18:46:20
949
原创 简单房间麦克风声压级测量仿真程序
之前一直想写一篇关于麦克风测量声压级的文章，一直没腾出时间，刚好前段时间有个项目用到了这方面的知识，趁项目结束，把这块的相关知识也总结出来。大体构思是在一个房间，由点声源发声，然后测量在房间任何一处位置的声压，给出仿真程序。这里面肯定用到 RIR（房间冲击响应），模拟声传输路径、声压级的计算这三大部分。1、房间冲激响应和T60房间冲击响应函数还是使用德国老牌的方法，镜面法https://www.audiolabs-erlangen.de/fau/professor/habets/softwa
2021-04-16 19:26:17
1634
5
原创 论wav文件格式和wav读写（c代码）
老调重谈的话题。基本上标准的44个字节的头，对照这个表格就可以了但也有较少的超过44个字节头的，比如这个：这篇文章说的很清楚https://sites.google.com/site/musicgapi/technical-documents/wav-file-format还有其他的chunkID多出来的，非44个自己的头，就是这些ID搞的鬼。我们搞信号处理的，不是很在乎其他的额外信息，那么直接对这些数据跳过去即可他们都遵循着 4个字节的 chunkID，4个字.
2020-11-26 15:23:55
771
原创 语音的合成和分解 enframe和overlapadd 用法
语音的常用处理就是分帧，分帧，fft变换，做一些处理，然后反变换回到时域。这里面常用到的就是enframe 和 overlapaddenframe分帧，我想大家都很熟悉了，这里不再介绍，overlapadd具体含义可以看这篇文章https://blog.youkuaiyun.com/u010592995/article/details/81001751上面这篇文章介绍的很清楚。一、问题的引入http://www.ee.ic.ac.uk/hp/staff/dmb/voicebox/mdoc/v_m
2020-11-25 17:33:01
5610
原创 平面和三维曲面图相交求交线和画图
一、三维散点图画图有的时候我们要画一些三维图，但图的方程没有，只有一些散点图，想要画这样的图mesh和surf是暂时不行的他们要求输入的是矩阵，但我们获得是 x y z这样的一些散点的结合，这时候画图只能用plot3比如我们这些数据和图：不好意思，具体描述什么物理量暂时不清楚，据说是一个亚/超高斯曲面现在有个需求，就是要从过 原点 和 过最大值（记为x0，y0，z0） 拉一个平面，对这个三维图进行切割，画出切割出出来的平面和曲面的交线，下面最大值，貌似在12 -90.5
2020-07-22 15:55:06
10970
17
原创 盛大哼唱检索前端提取算法分析
恒大开源头算法https://github.com/owensss/sdhumming两只老虎pitch提取, 求diff ，插值等五点均值化后处理 resample去除很小接近0的的，都给磨平除去毛刺降1/5采样率取log2去除slience地方vad卡限 6.35 到 8 算下来是 81.5719 到 256hz...
2020-04-11 14:27:47
679
2
原创 再议GMM和HMM的一些用法
推荐murphy的工具箱，这篇博客是根据大牛的工具箱使用写一些noteKevin Murphy 现在在Google brain工作，曾在MIT lab工作，他是 Machine Learning: a Probabilistic Perspective 的作者他写的关于机器学习的工具箱非常好，规范，注释也多，可供大家学习https://www.cs.ubc.ca/~murphyk/So...
2020-02-10 09:45:03
1983
3
原创 audacity VS2010 工程搭建
audacity的 VS2015工程网上已经有了，https://blog.youkuaiyun.com/weixin_34224941/article/details/93950373按照这个instructions就可以搭建一个比较稳定的vs2015（vs2017）的版本我这里主要是怀旧，用旧版本的vs2010 搭建一个轻量级的audacity工程供调试使用audacity的编译需要wxWid...
2019-11-07 16:13:19
520
原创 卡尔曼滤波之我见
最近简单看了一下kalman滤波，这里先写些笔记卡尔曼滤波说的是“两套系统“得到的计算量或者测量量，一个是根据系统状态预测出来的下一个状态的各种参量，另一个是实际测量值，那么我们怎么在这两套系统中找到更合适的值。说白了，它所谓的这两套系统都是服从高斯分布的，他们的乘也符合高斯分布，根据这两个高斯分布的方差的大小，调整到合适的方差选择F有的地方有人喜欢写成A。B是控制矩阵, x应该加...
2019-09-12 16:56:13
459
原创 Aurora2 语音数据库的建立
Aurora2 语音数据库 经常在学术上被用来作为评测数据库，很多语音识别的base数据库，HTK的训练模型数据库，VAD，基频检测都用Aurora2语音数据库。 Aurora2语音数据库的具体情况可以参考 http://aurora.hsnr.de/aurora-2.html,大致含义是 这个数据库是由 TIDigits 数据库 变换过来的，说白了就是 降采样之后...
2019-08-22 17:17:56
2352
2
原创 声卡内录工具之 Audio Record Wizard
这几天小弟在搞海康、大华的视频流分析，剥离其中的音频声道通道内容分析大厂有时候喜欢玩小装13，不用通用格式，用ffmpeg分析不出来音频流，后来想想，硬解决音频stream extract，还得用sdk，算了，找到内录工具直接录下来简单分析就好了Audio Record Wizard处理能录制mic in的还能录制speaker out的声音直接你在pc上放你要播放的声...
2019-07-15 10:56:24
3804
1
原创 再议Uniform FFT modulated filterbank
前面写过一篇文章，但是感觉对于FFT 的filterbak还是没有说清楚，这边就再次把自己的一些理解说一下1、这里面包含的几个概念uniform，说明就是基础滤波器就一个，其他滤波器是它的一个变体，也就是调制（modulated）出来的，polyphase分解，多相分解，主要目的是同一个多相组里面对应的“延迟“”是一样的，等会我后面会详细揭秘2、FFT干嘛？只为了在analysis...
2018-11-09 17:58:46
719
原创 Matlab坏了？国庆小修Matlab小毛病
奇了怪了，今天用Matlab，不能双击打开m文件了不用也没有什么啊，直接用Matlab editor 是可以打开的但是只能先打开Matlab中的editor再把m文件拖到editor里面去。。。。，用用就不爽了奇了怪了，马上Baidu，谷歌之，没有这个问题的答案，我去。。。 finfo文件应该是matlab系统函数自带的，马上在matlab系统目录中找找，的确有几个叫fin...
2018-10-06 17:28:56
400
原创 中秋无事，小破一个时间限制的软件
偶然从网上的一个某家语音合成的工具，效果还可以，只是软件试用期只到2013年：已经说明白了，2013年就不能用了对于一般人，改改系统时间，用用就得了，但是对于作为程序员的我来说是不可以接受的，一定要作一下，把它破解掉先用侦壳工具language.exe看看，软件没有加壳，大胆搞起来马上用 Ollydbg52pojie 进行反汇编：通过搜索找到关键词：双击...
2018-09-23 18:23:11
559
1
原创 fftw3/gsl/kissfft/OouraFFT库中傅里叶变换/反傅里叶变换函数和Matlab中的fft/ifft的对应关系
只分析一维度的一、fftw_plan_dft_1d正变换：fftw_complex *in = fftw_malloc ( sizeof ( fftw_complex ) * n );fftw_complex *out = fftw_malloc ( sizeof ( fftw_complex ) * n );plan_forward = fftw_plan_dft_1d ...
2018-09-04 17:26:26
7650
2
原创 WebRTC VAD 中所用滤波器之分析
WebRTC VAD用了一些滤波器，这些滤波器的设置还是非常巧妙的，这篇文章分析他们所用滤波器之原理。下面分析的滤波器都是在文件vad_filterbank.c里面 在提取特征的时候，用的是滤波器组用来提取的log域上的能量，其中分band的频带是：（对于8khz采样的音频文件来说）3k-4k feature[5] ， 2k-3k feature[4] ， 1k-2k feat...
2018-06-18 20:45:41
4271
11
原创 atan2 matlab VC excel 的不同输入
atan2 ，倒算 tan，输出范围 是 -pi 到 pi注意看 不同IDE下面的 y x 输入！！！matlabVC装13的微软就excel和其他不一样，爱玩个性化，不知道 ios上的这个函数是怎样的，有朋友知道告诉我下...
2018-06-04 10:50:31
455
原创 Filter Bank 浅谈一二
语音信号中子带处理技术用的较为频繁，涉及到子带技术，不得不说说Filterbank一、最简单的两通道的滤波器组如下图，它有最基本的一些单元，比如分析滤波器，合成滤波器，插值和抽取其中H0，H1可以理解为低通、高通滤波器。我对一段音频进行低高通滤波，然后降采样，升采样，最后合成。他们的频谱关系如下图，滤波器我就不说了，效果很明显，这里有个问题就是为什么降采样？不做行不行。实际上降采样（或者叫抽取）是...
2018-05-18 17:31:25
21825
5
原创 WebRTC VAD算法初探
这几天对webrtc中的VAD算法做了一些研究，有了一些心得VAD算法中核心的是使用了GMM进行分类，noise作为一类，speech作为一类，两类求后验概率，并且实时的更新GMM的参数但是这个GMM的具体情况是怎么样的，它的均值，方差，都是多少，使用特征的维度是多少？参数又是怎么更新的 一、GMM的 权重，均值，方差分析频率域上分成了6个子带（channel），对每个子带求能量，初始值有一个GM...
2018-03-23 10:37:40
7323
9
原创 高斯混合模型（GMM model）以及梯度下降法（gradient descent）更新参数
关于GMM模型的资料和 EM 参数估算的资料，网上已经有很多了，今天想谈的是GMM的协方差矩阵的分析、GMM的参数更新方法1、GMM协方差矩阵的物理含义涉及到每个元素，是这样求算：用中文来描述就是，注意后面的那个除以样本数-1，就是大括号外面的E （这叫无偏估计）上面公式也提到了，协方差本质上是就是很多向量之间的内积，内积是什么？ 举个例子说明，
2017-12-19 15:31:30
7359
5
原创 潮州方言的八个声调初探
提及潮汕，你脑子里面是不是潮州菜？要我能想到的就是周星驰的食神的台词：你家潮州人啊？从这里面可以看出潮州话的特殊。言归正传，博主从学校开始就做调式识别，一直想研究南方方言的调式。汉语是有调的语言，世界上用调表意的语言比较少，也比较难掌握（作为外语来说）我们知道的有：越南语有六种音调，普通话有四个调（不算轻声），粤语五个声调（阴平、阴上、阴去、阳平、 阳上），潮州话八
2017-08-10 09:41:19
23962
原创 Matlab 多线程解决的一些办法
Matlab 是不支持多线程的，貌似这是一个truth ，但是有些时候是非常需要Matlab能够多任务处理的。怎么使得Matlab能支持多线程呢，貌似要实现这个有点困难，不过总有一些办法，使得我们在不同的情况下，处于不同的考虑，能接近或者基本达到我们的目标。我大概总结了一下可能用的几种方法，下面给与分别介绍。
2017-07-13 16:03:29
36031
4
自定义函数实现delaunayTriangulation 使用Bowyer-Watson 算法
2024-05-23
Gunplot_win32_vs2010_编译版本
2019-02-18
Gunplot_win32_2010版本
2019-01-28
ActiveTcl8.4.11.2.201775-win32-ix86-threaded.exe
2018-10-28
libsndfile的vs2010工程
2018-04-22
NIST SPHERE shorten无损音频文件解码程序
2017-12-12
多相滤波器的实现
2017-08-05
FAAD2解码AAC的matlab代码
2017-06-29
FM stereo signal simulation
2017-03-02
带图形界面形式的gdb 在cygwin下使用
2014-09-22
带图形界面形式的gdb，使用在Cygwin下
2014-09-22
wma 解码vc6.0 环境
2014-09-11
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人