- 博客(32)
- 收藏
- 关注
RSS订阅原创 音乐信号处理基础
设置文件格式例如，将Windows系统音频文件转移到Mac上，需要把格式从wav转化成aiff。通过命令psf_format format;format=psf_getFormatExt("soundtrack.wav")getFormatExt就是提取字符串并分类的程序。执行这个函数就可以了。它返回psf_format类型文件，可以直接将原始二进制音频文件分配到合适的PSF_PRO...
2020-06-09 21:53:25
3593
1
原创 随机过程基础(8)---随机过程的线性变换
判断是否存在固定音高如果音高值几乎是恒定的，那么音高信息就能够找到。对于清音的片段，可以发现获得的音高因素很小，并且变化很大。那么可以说轻音片段的音高信息是缺失的。如果能找到重复的音高，这个片段就是浊音。预相位能量比可以求归一化的预相位能量比区分浊音和清音片段。对于浊音，相邻样本之间的变化值很小，而对于清音部分这个分子很大。求所谓预相位能量的程序如下：...
2020-05-10 00:29:06
2137
原创 数据结构 图论DFS、BFS、最短路径树、最小生成树 python实现
广度优先遍历(BFS)图的广度优先遍历可定义如下：首先访问出发点v，并将其标记为已访问过；然后从v出发依次访问v的每个未被访问的邻接点w。当v的所有邻接点全部访问完后，再对v的每个邻接点w继续进行广度优先遍历，直至图中所有和源点v有路径相通的顶点均已被访问为止。若此时图中仍有未访问的顶点，则另选一个尚未访问的顶点为新的源点重复上述过程，直至图中所有的顶点均已被访问为止。如下代码只考虑无向连通图...
2020-04-28 23:56:32
1538
原创 随机过程基础(7)
目录互功率谱非平稳过程的功率谱广义功率谱密度时变功率谱例子matlab互相关函数估计matlab功率谱估计例子：数字图像直方图随机过程应用：判断说话中的清音浊音互功率谱类似从自相关函数推出互相关函数，也可以从单个随机过程的功率谱推出互功率谱。两个随机过程X(t)和Y(t)，它们的互功率谱定义是GXY(ω)=E[limT→∞12TXT(ω)YT∗(ω)]G_{XY}(\omega)=E[\l...
2020-04-28 01:06:50
1514
原创 随机过程基础(6)--应用随机过程分析音乐（语音）信号（1）、随机序列功率谱（PSD）
随 机 序 列 功 率 谱对于平稳随机序列X(n)，如果相关函数满足∑m=−∞+∞∣RX(m)∣<∞\sum_{m=-\infty}^{+\infty}|R_X(m)|<\inftym=−∞∑+∞∣RX(m)∣<∞那么它的功率谱定义为自相关函数RX(m)R_X(m)RX(m)的傅里叶变换：GX(ejω)=∑m=−∞+∞∣RX(m)∣e−jmωG_X(e^{j\o...
2020-04-27 03:03:01
2983
原创 高级数据结构---B树、红黑树 python实现
B树一棵 2t (t>=2)阶（此处阶数表示每个节点最大的孩子数量）B树是一棵平衡的 2t 路搜索树。它或者是空树，或者是满足下列性质的树：1、根节点至少有两个子女；2、每个非根节点所包含的关键字个数j满足：t-1<=j<=2t-1；3、每个节点都包含了目前节点内key数量+1个孩子指针，叶子节点除外；4、节点孩子树中的key与当前节点中key的值存在大小关系；5、所...
2020-04-25 21:55:17
771
原创 电磁场与电磁波(2)
第一章部分习题：1：试采用直角坐标中∇⋅A=∂Ax∂x+∂Ay∂y+∂Az∂z\nabla\cdot\boldsymbol A=\frac{\partial A_x}{\partial x}+\frac{\partial A_y}{\partial y}+\frac{\partial A_z}{\partial z}∇⋅A=∂x∂Ax+∂y∂Ay+∂z∂Az相似的方法推出圆柱坐标下的...
2020-04-24 21:33:29
1617
原创 python--音乐播放小程序
1.网易云公开可下载内容爬取器（声明：未触犯软件使用规则及条款法律）伪装成Edge浏览器进入界面，爬取列表中第一个音频地址并存入相应文件夹中。这里有一个最简单的爬虫程序（注意，要先在网易云音乐网页（html语言）中获得第一个链接的地址），和一个最简单的tkinter GUI编程；from tkinter import *from selenium import webdriverglobal...
2020-04-24 02:23:44
614
原创 通信电子电路(5)
上章一些题儿：1：这题不要求定量计算参数，所以简单，把已有的模块合起来就可以了。只是所谓组合式偏置电路和零偏置电路没有提到。零偏置本来是PN结中的概念。在本征半导体基础上，n型半导体中添加5价（磷等）杂质（pentavalent impurities）而p型中添加3价（硼、铝等）杂质（trivalent impurities）。杂质的在P/N型半导体中添加是不均匀的，因此载流子（charg...
2020-04-22 22:35:47
3267
原创 通信电子电路(4)----高频功率放大器(2)
RC Ec Eb Ubm变化的影响ucemin>Uces，欠压u_{cemin}>U_{ces}，欠压ucemin>Uces，欠压ucemin=Uces，临界u_{cemin}=U_{ces}，临界ucemin=Uces，临界ucemin<Uces，过压u_{cemin}<U_{ces}，过压ucemin<Uces，过压ucemin=Ec−U...
2020-04-21 22:33:11
3907
原创 通信电子电路(3)---高频功率放大器
基本原理电路小信号调谐放大器的输入信号很小，在微伏和毫伏数量级，晶体管工作于线性区域；它的功率很小，但通过阻抗匹配，可以获得很大的功率增益（30~40dB）。小信号放大器一般工作在甲类状态，效率较低。而调谐功率放大器的输入信号要大很多，几百毫伏到几伏，晶体管工作延伸到非线性区域——截止和饱和区。这种放大器输出功率大，效率较高，一般工作在丙类状态，主要技术指标是输出功率、效率和谐波抑制度。如下，...
2020-04-20 22:20:49
4987
原创 通信电子电路(2)---使用multisim分析电路小技巧、晶体管等效电路
RLC电路瞬态分析在瞬态分析中设置参数如下，输出参数为信号源两端和R2两端的电压输出结果如下。示波器：示波器经常使用。右侧缩略栏中有4通道、2通道示波器，按照实际实验中的方法接就可以。失真分析器失真分析器在音频中常用，设计硬件效果器和功放时常用它来测试。其中，THD是total harmonic distortion，指(快速傅里叶变换后)信号不大于某特定阶数H（默认取1...
2020-04-18 22:51:21
5356
原创 随机过程基础(5)---各态历经性(2)、联合分布、连续信号功率谱密度（PSD）
各态历经性可以证明，随机过程X(t)具有均值遍历性的充要条件是limY→∞1T∫02T(1−τ2T)[RX(τ)−mX2]dτ=0\lim_{Y\rightarrow \infty}\frac{1}{T}\int_0^{2T}(1-\frac{\tau}{2T})[R_X(\tau)-m_X^2]d\tau =0Y→∞limT1∫02T(1−2Tτ)[RX(τ)−mX2]dτ=0具...
2020-04-17 21:48:38
4322
1
原创 随机过程基础(4)---各态历经性、典型随机过程matlab仿真
各态历经性部分源自《随机信号分析与处理》，罗鹏飞等著，清华大学出版社；公式旁标注*；时频转换部分证明源自《信号与系统分析》，吴京等著，国防科技大学出版社，证明处有标注。部分参考资料https://en.wikipedia.org、http://www.maths.qmul.ac.uk/~bb/TS_Chapter4_3&4.pdf各态历经部分代码来自https://www.math...
2020-04-16 22:57:26
5391
原创 C++贪吃蛇 最简单实现
目的就是实现贪吃????程序。按照面向过程的思想来想的话，搭建边框—>开辟蛇的数组并画出蛇—>随机生成第一个食物—>保持读取更改蛇头方向指令，同时每隔一段时间让蛇按照某方向前进（通过蛇所在的坐标更新实现），同时更新显示食物位置，同时检查：蛇头是否碰到边界、碰到自己（这样游戏结束）；是否碰到食物（这样就更新蛇长—>加分—>判断加等级、酌情加快蛇前进速度。）按照面向对象(OO...
2020-04-15 22:54:28
2792
原创 随机过程基础(3)---随机过程的统计描述、平稳随机过程
随机过程的概率分布一维连续概率分布：Fx(x,t)=P{X(t)≤x}F_x(x,t)=P\{X(t)\le x\}Fx(x,t)=P{X(t)≤x}对x偏导，就是X(t)的一维概率密度。显然，这两个二元函数将随机变量的分布\密度铺到整个时间轴上，一个时间对应一个变量分布情况。因此只考虑某个t的结论全部和随机变量结论相同，例如∫−∞∞f(x,t)dx=1\int_{-\infty}^{\inf...
2020-04-14 22:51:38
1222
原创 机器学习入门(1)---以手写数字识别为例
(部分代码源自pytorch mnist数据集文档)机器学习的核心组件包括：层(相当于转移函数)、模型(层网络)、损失函数(衡量某模型参数选择好坏的函数)、优化器(使损失函数最小)。简单的手写数字识别使用mnist数据集，其中70000张用于手写数字图片，60000张用于训练，10000张用于检验。(test data和training data应该分开，并且不能根据training data...
2020-04-13 21:27:53
1076
原创 通信电子电路multisim仿真(1)---单调谐放大器
单调谐放大器以上是经过我调试的单调谐放大器最基本的结构。由于没用变压器。这个电路当信号源和负载变化时无法阻抗匹配)。运行如下。这个结果是可以接受的。以下是加入变压器进行阻抗匹配后得到的标准的电路图。经我调试，得到如下结果：(由于multisim中抽头变压器的QL等参数不好调，这里的信号源姑且根据前几次的运行结果调整到LC回路的谐振频率。由于未知multisim抽头变压器线圈的品质因数...
2020-04-11 23:36:24
6540
原创 数值计算方法 线性方程组的数值解法(4)---向量和矩阵范数(norm) 高斯-赛德尔(Gauss-Seidel)迭代、共轭梯度(Conjugate Gradient)迭代
向量范数设x∈Rn\boldsymbol x\in \boldsymbol R^nx∈Rn则范数||x||满足：∣∣x∣∣≥0||\boldsymbol x|| \ge 0∣∣x∣∣≥0当且仅当x=0,||x||=0。任意实数λ∣∣λx∣∣=∣λ∣∣∣x∣∣||\lambda\boldsymbol x||=|\lambda|||\boldsymbol x||∣∣λx∣∣=∣λ∣∣∣x∣∣且对任意...
2020-04-10 22:26:30
1464
原创 随机过程基础(2)---多维随机变量常用性质、随机过程的引入
二维随机变量分布函数和概率密度二维随机变量的联合概率密度通过如下方式定义:fXY(x,y)f_{XY}(x,y)fXY(x,y)且∫−∞∞∫−∞∞fXY(x,y)dxdy=1\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}f_{XY}(x,y)dxdy=1∫−∞∞∫−∞∞fXY(x,y)dxdy=1可以想象，二维随机变量的联合概率密度占据了三维...
2020-04-09 22:16:26
1881
原创 数值计算方法--线性方程组的数值解法(3) 追赶法(Thomas)，选择主元(Pivoting)，求逆(Inversion)
追赶法追赶法适用于三对角矩阵，形如(XX000XXX000XX0000XXX000XX)\begin{pmatrix}X&X&0&0&0\\X &X& X& 0& 0\\0&X&X&0&0\\0&0&X&X&X\\0&0&0&X&X\en...
2020-04-08 21:54:09
3937
原创 随机过程基础(1)---一维随机变量标准分布函数、常见性质MATLAB实现
基础概念大部分在高中数学、概率论和离散数学等中已经涉及。本文生成正态分布变量部分代码摘自《Random Signal Processing》，Shaila Dinkar Apte著，Taylor Francis Group出版。正态分布老办法实现：sigma=1;mean=0;x=(-10:0.1:10);a=(x-mean).*(x-mean);y=(1/(sigma*sqrt(...
2020-04-07 22:00:20
3468
原创 矢量场散度(divergence)和旋度(curl)及Helmholtz定理MATLAB计算
旋度的定义 已给向量F⃗=P i⃗+Q j⃗+R k⃗\vec F = P\,\vec i + Q\,\vec j + R\,\vec k F=Pi+Qj+Rk 则该向量的旋度表示为rot F⃗=(Ry−Qz)i⃗+(Pz...
2020-04-06 22:38:02
13440
原创 数值计算方法(Numerical Methods)MATLAB实现(2)--choleski分解
1. 乔列斯基分解 乔列斯基分解当且仅当矩阵为正定对称阵，运算量大约是LU分解的一半(n3/6)，另有n次取平方根计算(n为方阵维度)。 对于3*3矩阵A = LLT。根据对称性，表示为 &...
2020-04-05 17:10:32
1006
1
原创 数值计算方法(Numerical Methods)MATLAB实现(1)---Gauss消元法、Doolittle分解
（1）—高斯消元法 1.1 消去阶段 假设系数矩阵前k行已被转化为上三角矩阵形式。当前枢轴方程(作为被减量的方程) 是第K行的方程，其下所有方程都待转换为下三角形式。 假设现要消去第i行的方程，也即系数Aik要被消去。将这行减去A...
2020-04-04 18:05:05
2573
原创 人工智能在音乐领域的应用概述(1)---听歌识曲
人工智能将成为很长一段时间热门话题。本文将概括其在音乐领域的垂直应用。AI+机器听觉(Content-Based MIR，基于内容的检索)从根据语言数据(直接用文字标定每一首歌的风格，歌手等信息)到直接根据音乐本身进行推断的转变。1.1 Music Fingerprinting（音乐指纹）H(X)理解为哈希函数，将音频数组X映射到哈希值Y，Y有更小的数据量，同时保留歌曲之间差异，具有辨...
2020-03-12 23:06:57
6745
3
原创 C++简易计算器_理解面向对象
为更好理解面向对象的过程，举一例，如输入表达式3+2*5，解析与求值过程如下： 一个字符一个字符的读，根据字符的种类去决定组合成哪种单词（token），这个过程叫词法分析，由词法分析器...
2020-03-10 15:52:17
970
原创 python AVL树(详细注释)
AVL树本质是二叉搜索树，但它又具有以下特点：只能是一棵空树或它的左右两个子树的高度差的绝对值不超过1，并且左右两个子树也分别满足同样条件。在AVL树中任何结点的两个子树的高度最大差别为一，所以它也被称为平衡二叉树。#AVL树（二叉平衡树）#节点类。AVL树相对一般二叉搜索树，节点增加树高属性，便于...
2020-03-10 15:32:05
1650
2
原创 python括号匹配(详细注释)
实际编程中，我们经常嵌套使用括号(“{}”、“[]” 、 “()”)，如果括号太多，可能会出现括号不匹配的情况，比如“(as))”、“{(bcd})”等。现用栈判断一串字符中的括号是否匹配。此代码不检查小括号是否一定被大级别括号包着，只检查左括号与右括号是否对应。#括号匹配,用栈实现#只用判空、压栈...
2020-03-10 11:38:49
5999
1
原创 python按照层序遍历建立二叉树(详细注释)
按照层序遍历建树是树的基本操作之一。通常用创建队列，逐层加入元素检查某层是否满，直到发现空位置的方法解决。本代码增加一种输入’#'字符的情况，表明该节点存储的数据为空。#构造二叉树并实现访问操作，输入#代表树中相应节点的数据是none#本题的二叉树是用层序构建的#节点类class Node(): ...
2020-03-10 10:34:00
4024
原创 python散列表 链地址法解决冲突 除留余数法散列 插入删除查找操作（详细注释）
散列是数据结构中较为重要的内容。两种基本方法之一的链地址法虽稍繁琐，但思路较简单，过程清晰。散列表（Hash table，也叫哈希表），是根据关键码值(Key value)而直接进行访问的数据结构。也就是说，它通过把关键码值映射到表中一个位置来访问记录，以加快查找速度。这个映射函数叫做散列函数，记录的数组叫做散列表。散列函数：除留余数法（取关键字被某个不大于散列表表长m的数p除后所得的余数为散...
2020-03-05 14:59:59
2509
3
原创 python 01背包问题（详细注释）
#首先定义01背包函数#输入参数p为一个列表，顺次表示物品价值的列表；输入参数W，顺次表示物品的体积；输入参数v是整型，表示背包容量的大小def knapsack(p, w, v): #n来存储物件的总数目 n = len(p) #建立动态规划的过程划分矩阵。lists是应该放进去的物品的集合，arr为初始的动态规划矩阵。动态规划矩阵有v+1列，n+1行。v代表的是背包...
2020-03-04 18:15:41
1888
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人