彬彬有礼的专栏

也就是说，矩阵Z中的第(r,s)个元素等于X的第r行之和乘以Y的第s行之和。因此，可以先分别将X和Y所有列对应相加，得到长为d的列向量Xsum和Ysum，则Z等于Xsum乘以Ysum的转置。注意到第一个括号是将X所有列对应相加，第二个括号是将Y所有列对应相加，因此Z也就是前面提到的Xsum乘以Ysum的转置。其中X(:,i), Y(:,j)分别表示矩阵X, Y的第i列和第j列，易知Z为d×h矩阵。注意到第一个括号是X的第r行之和，第二个括号是Y的第s行之和。这两天一直在琢磨如下矩阵计算问题。

题目：Friedman test以及后续检验Nemenyi test和Bonferroni-Dunn test在做算法对比时，往往需要对实验结果进行统计检验。Friedman test是一种常用的检验，用来比较k个算法在N个数据集上的整体表现性能。但Friedman test只能给出k个算法的性能之间是否存在差异，如果存在差异，还需要进行“后续检验”(post-hoc test)，常用的后续检验方法包括Nemenyi test和Bonferroni-Dunn test。Nemenyi test适用于

题目：一句话记住Bregman距离Bregman距离，又称Bregman散度(Bregman Divergence)，两个点与之间的Bregman距离定义如下： 其中为某凸函数，表示求导，表示求内积；可以看出Bregman不满足对称性。回忆一下Taylor公式，可以看出实际上是函数的一阶Taylor近似，因此用一句话概括就是：【与之间的Bregman距离实际上可以理解为函数与其一阶Taylor近似之差】。...

线性时不变(LTI)系统能实现频谱搬迁么？作者：jbb0523(彬彬有礼)前言  曾经，老师问过我一个问题：线性时不变(Linear Time Invariant, LTI)系统能实现频谱搬迁么？当时我真不知如何回答，因为实在不知道从何处下手去分析这个问题。今天也算对此问题进行一个尝试性的回答，因为我始终感觉有很多细节还没明白，LTI这个概念在我脑子里变得扑朔迷离……已知某LTI系统，

雷达的波段    雷达波段(radar frequency band)是指雷达发射电波的频率范围。其度量单位是赫兹(Hz)或周/秒(C／S)。大多数雷达工作在超短波及微波波段，其频率范围在30－300000兆赫，相应波长为10米至1毫米，包括甚高频(VHF)、特高频(UHF)、超高频(SHF)、极高频(EHF)4个波段。第二次世界大战期间，为了保密，用大写英文字母表示雷达波段。将230—1

===========================第一篇===========================旧的无线电波段划分中L、S、C、X、Ku、Ka、W波段频率分为分别是多少？这种划分方式是雷达业内的通俗叫法，没有一个严格、统一的标准。通常的划分是：L波段 1~2GHz；S波段 2~4GHz；C波段 4~8GHz；X波段 8~12GHz；Ku波段 12~18GHz；K波

冲激响应不变法或双线性变换法中的参数T为什么是一个无关紧要的参数？在数字信号处理课程中，IIR滤波器设计方法一般只讲间接法，即首先根据实际要求设计一个模拟滤波器，然后再将这个模拟滤波器转换成数字滤波器。从模拟滤波器转换成数字滤波器的方法常见的有三种：冲激响应不变法、阶跃响应不变法、双线性变换法，其中以冲激响应不变法和双线性变换法更为常见；从模拟滤波器转换成数字滤波器的本质是从拉普拉斯(Laplace)域(或称为S域)映射到Z域。详细设计方法可以参见《姚天任，江太辉. 数字信号处理（第二版）[M]. 华中科技

离散正(余)弦信号的时域与FFT变换后所得频域之间的关系作者：jbb0523(彬彬有礼)        正弦信号在信号处理中是很常见的，比如通信领域的载波。由于正弦与余弦只是相差π/2的初相，因此这里统称正弦信号。给出连续正弦信号的表达式：式中，A为振幅，Ω为模拟角频率(rad/s)，φ为初相，f为模拟频率(Hz)，Ω=2πf 。        在满足奈奎斯特采样定理条

题目：对称矩阵、Hermite矩阵、正交矩阵、酉矩阵、奇异矩阵、正规矩阵、幂等矩阵        看文献的时候，经常见到各种各样矩阵，本篇总结了常见的对称矩阵、Hermite矩阵、正交矩阵、酉矩阵、奇异矩阵、正规矩阵、幂等矩阵七种矩阵的定义，作为概念备忘录吧，忘了可以随时查一下。1、对称矩阵（文献【1】第40页）其中上标T表示求矩阵的转置（文献【1】第38-39页）

题目：二次规划基础：二次型、正定矩阵、海塞矩阵        学习二次规划时，需要知道一些基础知识方能看懂，二次型、正定矩阵和半正定矩阵是绕不开的概念，海塞矩阵也是经常要用的名词，所以本篇首先学习一下这三个概念。1、二次型（文献【1】第127-128页）主对角线上的元素分别为x2,y2,z2的系数（分别为1，0，3），第1行第2列和第2行第1列两

前些时间需要衡量多个分类数据之间两两相关程度，想找出最相关的一对分类数据；于是想到了曾经看过的克莱姆相关系数，但在网上搜了好久之后，即没发现Matlab现成的built-in函数，也没找到别人分享的Matlab代码，于是决定自己动手写一个~ 克莱姆V（Cramer’s V），又称为克莱姆相关系数、克莱姆关联系数、独立系数等，是双变量相关分析的一种方法，专门用于衡量...

前面发过三篇有关dds/nco方面的博文：1）有关DDS的一些初步理解：相位噪声和无杂散动态范围（http://blog.csdn.net/jbb0523/article/details/7879999）2）DDS的对称性质——参考时钟是100M，则产生的90MHz与10MHz的正弦信号频率相同，相位相反（http://blog.csdn.net/jbb0523/article/d

DDS，Direct Digital Frequency Synthesis，直接数字频率合成1）原理两个重要的组成：1.1）相位累加器，例如范围从0~2^N-11.2）只读存储器，例如存储着正弦波的波形，此处假设里面正常存储，即存储0~2pi相位，而并没有采用其它的一些技术，如只存0~pi/2范围，这里有两个重要的概念，1.2.1）共存储了0~2pi范围内N1个点，此处假如存

FPGA信号截位策略研究作者：jbb0523（彬彬有礼）        在FPGA中，随着信号处理的层次加深，对信号进行乘、累加、滤波等运算后，可能输入时仅为8位位宽的信号会扩展成几十位位宽，位宽越宽，占用的硬件资源就越多，但位宽超过一定范围后，位宽的增宽并不会对处理精度带来显著的改善，这时就需要对信号进行截位。写过FPGA HDL代码的童鞋都应该知道，截位是最为经常的一种操作。

先在matlab中运行以一命令，观察一下结果： >> b=[1 2 3 4 5 6 7 8]b =     1     2     3     4     5     6     7     8>> fft(b)ans =  Columns 1 throu

我们经常利用FFT来计算相关运算，这样速度会快很多……例如两个长为N的序列A和B，可以按如下方法计算其相关值： ifft(fft(A).*conj(fft(B)))，conj为取共轭。利用这种方法得到的相关结果是实现的循环相关，并且是取共轭的序列相对于未取共轭的序列向右移

前些天写了如下日志：如何计算一个有符号数的补码表示？http://blog.csdn.net/jbb0523/article/details/6668237 现在补充一点：按照【谭浩强.C程序设计(第三版).北京：清华大学出版社，2005】的第40页至第4

特性阻抗      根据传输线理论和信号的传输理论，信号不仅仅是时间变量的函数，同时还是距离变量的函数，所以信号在连线上的每一点都有可能变化。因此定义连线的交流阻抗，即变化的电压和变化的电流之比为传输线的特性阻 抗。Z(w):理想传输线的特性阻抗，单位Ω；L: 理想传输线的电感，H/mm;C:理想传输线的电容，F/mm。      传输线的特性阻抗只

使用matlab进行频谱分析时若干问题解释作者：jbb0523（彬彬有礼）本文共说明了以下问题：一、在matlab中如何表示频率为f1，以采样率f抽样后所得到的数字信号？如此表示的依据是什么？二、使用matlab画出的频谱(一般是幅度谱或称振幅谱)的横坐标轴的意义是什么？如何根据横坐标轴的值得到其所对应的实际频率？三、实数序列的频谱除第零个点和第N/2个(当N为偶数时)点外（从

一、“取反加一”的新理解在前面的两篇博文中分别讨论的“取反加一”的两个作用。（取反加一应该是按位取反加一，下同）第一个应用：在博文《如何计算一个有符号数的补码表示？》（链接http://blog.csdn.net/jbb0523/article/details/6668237）提到了：若想求一个负数的补码，则对其绝对值进行“取反加一”即可；第二个应用：在博文《知道某负数补码后如

题目：DDS的对称性质——参考时钟是100M，则产生的90MHz与10MHz的正弦信号频率相同，相位相反 DDS的原理大家都知道，在FPGA中的可以调用IP核生成，ISE中DDS的IP核界面如下图所示：现在的说明的问题是：若产考时钟为fs，那么由DDS生成的fo与fs-fo的两个正弦信号频率相同，均为fo和fs-fo中小于fs/2的一个，而相位相反，若产生余弦信号，则完全相同。

题目： 数字信号处理中：采样率、采样最小记录时间的一点直观认识 首先上一幅截图，来自姚天任的《数字信号处理（第二版）》，华科出版社，第92页，写的很不错的一本书。1、文中提到了频率分频率等于最小记录长度的倒数，这个可以这样直观的来认识：假设你采样的数据时长为tp，那么按照上述理论，频率分辨率是F=1/tp，这是因为你采样的数据只有tp长，对于频率小于F的信号，它的周期将大于tp，也

求助：关于信噪比的确切含义——Eb/No和SNR的区别在《现代通信系统——使用matlab》一书中，一些调制方式的Monte Carlo仿真性能曲线中，横轴是用Eb/No来表示的，而matlab代码是用的却是SNR，很多人都说这两个信噪比的含义是不一样的，我也曾见过其他的系统性能曲线，仿真条件相同，而信噪比分别为Eb/No和SNR时，曲线是不同的。谁能解释一下这两个信噪比的区别？他们之

在某些时候，由于某些数据符号位太多，截位时我们需要截掉高几位的符号位，但截高位时就需要注意一个问题：硬限幅！硬限幅的概念可以直观的用如图看出：如图所示，下面的正弦波就是上面正弦波的一个硬限幅。这样做有如下好处：虽然你认为你截掉的都是符号位，但难免还是会有一些数较大，比如上图中正弦的被截掉的那些值，硬限幅可以使这些过大的数不改变自己的符号位，只是在幅度上有一些失真。例如：对于

按照【谭浩强.C程序设计(第三版).北京：清华大学出版社，2005】的第40页至第41页所述，求一个负数的补码的方法如下：【例】求-10的补码的方法如下：1）取-10的绝对值10；2）10的绝对值的二进制形式为1010；3）对1010取反得1111 1111 11