诚朴勇毅

原创 AD9854+STM32正弦波信号发生器

https://pan.baidu.com/s/1oBLHT23-HahRMiAKpZWaZQAD9854使用20M晶振，MCU为STM32F103ZET6，该程序可以在百度网盘免费下载。

原创 matlabR2013B许可证到期激活问题

我安装的Matlab R2013b在使用了不到两个月之后，一直提示要激活。本来想重新安装一个，后来找到了解决办法。在Matlab根目录的license文件夹里，下载以下链接的license.lic文件，替换原来的license.lic文件，重启Matlab即可。R2013a亲测有效。链接: https://pan.baidu.com/s/1mixMy84 密码: 3g59

原创 顺序存储结构的队列（C语言）

队列是一种特殊的线性表，它的特殊性体现在，它只能够从一段进，从另一端出，遵循先入先出的原则。这种独特的规则，可以运用在程序设计中的很多地方，非常的巧妙。下面，我用C语言来实现了循环队列。/*******************************************************************///描述：队列的顺序存储结构，该队列为循环队列，队列的创建，插入，

原创 matlab绘图背景为透明

【代码】matlab绘图背景为透明。

原创 【无标题】

【代码】【无标题】

原创 采用Otsu（最大类间方差法）大津算法进行信号检测

【代码】采用Otsu（最大类间方差法）大津算法进行信号检测。

原创 峰值检测方法

findpeaks高斯平滑预处理是一种常用的图像处理和信号处理技术，用于抑制噪声。在边缘检测的上下文中，我们通常使用高斯平滑来滤除高频噪声，同时尽量保留真实的边缘特征。

原创 多尺度自适应差分盒滤波

理论支撑：利用干扰切片的脉冲特性（陡峭边沿）与幅度波动的低频特性差异。self.scales = scales # 多尺度检测窗口。动态阈值：基于局部窗口（如50采样点）的极差自适应计算阈值。方案1：多尺度自适应差分盒滤波器（MS-ADBF）多尺度检测：小尺度捕捉锐利边沿，大尺度抑制幅度波动。

原创 脉内频率捷变LFM信号

可以看到，频率捷变之后旁瓣升高了。

原创 带约束的高斯牛顿法求解多音信号分离问题

将待估计参数定义为向量形式：$$$$

原创 小波变换检测信号突变位置

【代码】小波变换检测信号突变位置。

原创 word如何插入高清晰的matlab绘图

在matlab中画好的图另存为emf格式，保存到本地，然后在word中选择插图图片，注意不要复制粘贴。

原创 Latex公式排版，大括号cases环境使用

就是采用cases命令，它可以使用大括号包裹两行，且将条件对其，很方便，特此记录，供自己使用，也分享给给人参考。

原创 Latex参考文献中作者姓名显示为横线解决办法

从网上寻找原因，原来是因为有一篇参看文献的作者和这一篇一模一样导致的。

原创 恒虚警检测自适应门限计算方法

最近看到一篇文献，他给出的自适应门限计算方法效果还不错，特此记录。

原创 线性调频脉冲雷达信号的匹配滤波输出公式推导

原创 时域匹配滤波、频域匹配滤波以及自相关

上图可以看到经过脉压之后LFM信号被压缩成了一个尖峰。LFM信号的实部和匹配滤波器的实部是一样的。今天闲来无事，想看一下雷达信号匹配滤波的过程中匹配滤波器的时频特性。上面是频谱图，可以看大3者的带宽是一样的。假设，雷达发射的线性调频脉冲信号为。表示共轭，匹配滤波的频谱为。

原创 匹配滤波器的实现-时域卷积法和频域法

【代码】匹配滤波器的实现-时域卷积法和频域法。

原创 不同窗口宽度对应的短时傅里叶变换的结果

根据测不准原理，可知，时间窗越宽，频率分辨率越高，但是时间分辨率越低；时间窗越窄，频率分辨率越低，时间分辨率越高。下面给出一个例子，展示不同窗口宽度下的时频分布，窗函数统一采用hanning窗。图3 128点的时窗下时频分布。图4 256点的时窗下时频分布。图1 32点的时窗下时频分布。图2 64点的时窗下时频分布。

原创 短时傅里叶变换加窗和不加窗的区别

今天编了好长时间的代码，发现短时傅里叶变换后的时频图有点怪怪的，最后才知道原来是没有加窗的原因，特此记录一下，下面三幅图分别是不加窗，加hanning窗和加高斯窗的时频图。图3 加hanning窗的时频图。图2 加高斯窗的时频图。图1 不加窗的时频图。

原创 S变换matlab实现

【代码】S变换matlab实现。

原创 Matlab一些使用技巧

两个百分号就可以实现代码的分段，不同段之间会以不同的背景色显示，方便调试。

原创 时频分析之S变换

对应信号xt∈L2Rxt∈L2RL2RL^2(R)L2R为能量有限函数空间，xtx(t)xt的S变换的表达式为Sτf∫−∞∞xt∣f∣2πe−t−τ2f22e−j2πftdtSτf∫−∞∞​xt2π​∣f∣​e−2t−τ2f2​e−j2πftdt式中，xtx(t)xt是关于时间的连续函数，τ\tauτ。

原创 高斯函数Gaussian绘制matlab

约翰·卡尔·弗里德里希·高斯，（德语：Johann Carl Friedrich Gauß，英语：Gauss，拉丁语：Carolus Fridericus Gauss）1777年4月30日–1855年2月23日，德国著名数学家、物理学家、天文学家、几何学家、大地测量学家， 毕业于Carolinum学院（现布伦瑞克工业大学），后进入哥廷根大学深造。的时候，高斯函数的积分为1，在这种情况下，它是正态分布随机变量的概率密度函数，期望值。决定了峰值的位置，也就是对称轴的位置，越大的时候，曲线越胖。

原创 matlab的绘图的标题中(title)添加标量以及格式化输出

例如我们要绘制高斯函数，标题中添加均值和方差，那么可以采用sprintf函数进行格式化输出，其中%f表示当前位置为一个浮点数，%.nf表示保留n为小数进行显示。有时候我们需要在matlab绘制的图像的标题中添加一些变量，这样在修改某些参数后，标题会跟着一块儿变。

原创 kittler算法matlab函数-图像二值化分割

其思想：假设灰度图像由目标和背景组成，且目标和背景满足一混合高斯分布，计算目标和背景的均值、方差，根据最小分类误差思想得到的最小误差目标函数，取目标函数最小时的阈值即为最佳阈值。例如，当图像的灰度分布不满足混合高斯分布假设时，或者当目标和背景的灰度差异不明显时，算法的分割效果可能会受到影响。它试图找到一个阈值，使得图像中的目标和背景被正确分类的概率最大化，或者等价地说，使得分类误差最小化。在此基础上，算法的目标是通过计算找到一个最优的阈值，使得图像根据这个阈值被分割成二值图像时，分类误差达到最小。

原创 matlab一下特殊矩阵的操作（对角矩阵、全零行删除）

【代码】matlab一下特殊矩阵的操作（对角矩阵、全零行删除）

原创 LaTeX排版Elsevier公式缩进及参考文献修改等问题

最近打算投Elsevier(爱思唯尔)下的期刊，下载了他们给出的模板，在使用过程中出现了诸多问题，记录如下，以便给遇到同样问题的朋友提供帮助。

原创 visio导出pdf图片变得不清晰的解决办法

从网上查了好久，发现最有效的办法是在visio中选择“文件”->“打印”->“Microsoft Print to PDF”。我一开始选的是Adobe PDF，结果另存为的pdf文件打开之后就是不清晰，后来按照以上办法就清晰多了。

原创 matlab峰值检测

在信号处理中，我们常常需要寻找信号的峰值，对于全局最大值，寻找起来比较容易，而对于局部峰值的寻找就需要考虑更多的条件，例如峰值的绝对幅度以及相邻峰值之间之间的间距；，我如何在信号中找到峰值？如何测量受趋势影响的信号峰值的振幅？信号的峰值似乎以固定间隔出现。例如我们想要寻找如下信号的峰值，设置阈值后，提取到的峰值相对比较准确，但是在第3条竖线和第4条竖线上出现了两个峰值，如下图所示。图中红色圆圈标记的即峰值，可以看到峰值很多，这是因为它把很多局部的小的峰值也检测了出来。各个峰值的绝对值存在一定的差异。

原创 香橙派 AIpro

原创 香橙派 AIpro初体验

原创 Jupyter Lab 软件安装与使用

Jupyter Lab 软件是一个基于web 的交互式开发环境，集成了代码编辑器、终端、文件管理器等功能，使得开发者可以在一个界面中完成各种任务。

原创 香橙派AIpro使用SSH远程登录

香橙派AIpro可以连接HDMI显示器使用，也可以远程登录。这里采用MobaXterm软件远程登录开发板。

原创 香橙派AIpro（OrangePi AIPro）开发板初测评

最近，我拿到手一款Orange Pi AI Pro 开发板，它是香橙派联合华为精心打造的高性能AI 开发板，其搭载了昇腾AI 处理器，可提供8TOPS INT8 的计算能力，内存提供了8GB 和16GB两种版本。可以实现图像、视频等多种数据分析与推理计算，可广泛用于教育、机器人、无人机等场景。

原创 相似系数计算公式

1]李凤娥.自适应抗干扰及盲源分离应用技术研究[D].西安电子科技大学，2021.DOI:10.27389/d.cnki.gxadu.2021.001863.在采用盲分离算法进行信号分离的时候，需要对分离结果进行定量分析和评估，分离效可用相似系数来评估。仅存在幅度上的不同，分离的信号十分纯净；完全不相同，信号没有被分离。的数值越接近1，分离的效果越好。等于1，则表明分离的信号。等于0，则表明分离的信号。

原创 matlab绘图杂谈-stem函数和plot函数

三条曲线应该是用plot函数绘制的，而target哪个绿色的圆圈，我的理解是用stem函数绘制的。它只是1个点，并且没有竖线，stem绘制的默认的是茎叶图，既有线又有标记，像图中这个它就是只有标记。如果自己随便画一个图，那么直接使用语法2，其参数都采用默认参数，如果要自己控制参数，可以采用语法4。值得注意的是，使用plot函数的时候，这些参数的设置方法是和stem函数一样的。它可以分别控制线条、标记的各自的参数，具体可以查阅matlab官方的资料。它设定的线型是点划线，标记是圆圈，颜色是红色。

原创 matlab窗函数-hann窗和hamming窗函数

汉明窗（hamming）和海宁窗（hann）类似，汉明窗的时域波形两端不能到零，而海宁窗时域信号两端是零。简单点说，hann窗比较稳，短期长期发挥稳定，hamming窗比较激进，短期发挥水平高，长期就没劲了。这些窗函数通常在时域上是有限的宽度，并且具有对称性，如矩形窗、汉宁窗、汉明窗和布莱克曼窗等。例如，汉明窗是一种对称窗函数，它可以用来平滑信号，减少频谱能量的泄漏，从而提高傅里叶变换后的频率分辨率和精度。hann窗的matlab函数为hann，使用方法和hamming是相同的，这两种窗有些差异。

原创 matlab抽取与插值

我们假设一个数字信号xnn12Nxnn12...N共有NNN个点，抽取就是每个几个点抽1个点，比如2倍抽取，那么抽取后的信号为yny1x1y2x3y3x5yN2xN−1yny1x1y2x3y3x5...yN/2xN−1，可以看到抽取之后信号的长度变短了。插值和抽取相反，插值之后信号会变得越来越密集，注意绝对时间长度是不变的。例如数字信号xnn12Nxnn12...N。

原创 线性调频信号的解线调(dechirp，去斜)处理matlab仿真

去斜处理的主要原理为用发射副本信号与回波混频得到基带信号，再用傅里叶变换提取基带信号频率，通过时域到频域的转换将时域脉冲信号压缩为频域上的一个峰。匹配滤波和解线调其实都是一个sinc函数乘以一个复指数项的形式，两者的sinc函数一个位于快时间域一个位于距离频率域。去斜处理在一些特殊场合，不仅运算简单，而且可以简化设备，已经广泛应用于SAR和ISAR中作脉冲压缩。