book_bbyuan的博客

解决Ode45微分方程的实际问题

对meshgrid的数据挨个判断，如果 在xy平面的投影点 在区域内，那么保留，否则置为NaN。求K = 25456000*(1-4*(z/1000)^2)^(1/2);2、把范围之外的数据都设为NaN，mesh surf等直接就不显示了。K是 d1和d2的应变量，也间接是 x和z的二元函数。如果我们只想要交叠部分的mesh显示，该怎么办？我们有了 投影区域，还有meshgrid的每个点。先投影到二维，找到二维封闭曲线的边界。是一个有范围的空间曲面。1、需要知道交叠区域的范围。一定区域的三维数据，比如。

Matlab辨识工具箱里面的函数很多，这里用一个简单的例子来展示估算系统传递传递函数的使用方法。很遗憾，我们也不具体了解这个未知系统的 零极点个数（或者说多项式的阶数）系统未知，我们通过 功率谱估算 来求解系统响应。估算的系统响应，已经把两个峰值的特性 仿真出来了。互功率谱和自功率谱相除就得到了系统频响。title('传输tf幅度 小范围');title('传输tf相位 小范围');我们得到系统频响（幅度和相位）我们的系统输入输出数据如下。下面估算系统的传递函数。得到的 频响是一致的。...

最近要使用 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...

语音的常用处理就是分帧，分帧，fft变换，做一些处理，然后反变换回到时域。这里面常用到的就是enframe 和 overlapaddenframe分帧，我想大家都很熟悉了，这里不再介绍，overlapadd具体含义可以看这篇文章https://blog.csdn.net/u010592995/article/details/81001751上面这篇文章介绍的很清楚。一、问题的引入http://www.ee.ic.ac.uk/hp/staff/dmb/voicebox/mdoc/v_m

audacity的 VS2015工程网上已经有了，https://blog.csdn.net/weixin_34224941/article/details/93950373按照这个instructions就可以搭建一个比较稳定的vs2015（vs2017）的版本我这里主要是怀旧，用旧版本的vs2010 搭建一个轻量级的audacity工程供调试使用audacity的编译需要wxWid...

atan2 ，倒算 tan，输出范围 是 -pi 到 pi注意看 不同IDE下面的 y x 输入！！！matlabVC装13的微软就excel和其他不一样，爱玩个性化，不知道 ios上的这个函数是怎样的，有朋友知道告诉我下...

WebRTC VAD用了一些滤波器，这些滤波器的设置还是非常巧妙的，这篇文章分析他们所用滤波器之原理。下面分析的滤波器都是在文件vad_filterbank.c里面 在提取特征的时候，用的是滤波器组用来提取的log域上的能量，其中分band的频带是：（对于8khz采样的音频文件来说）3k-4k feature[5] ， 2k-3k  feature[4] ， 1k-2k  feat...

Aurora2 语音数据库 经常在学术上被用来作为评测数据库，很多语音识别的base数据库，HTK的训练模型数据库，VAD，基频检测都用Aurora2语音数据库。 Aurora2语音数据库的具体情况可以参考 http://aurora.hsnr.de/aurora-2.html,大致含义是 这个数据库是由 TIDigits 数据库 变换过来的，说白了就是 降采样之后...