冬瓜班小朋友的博客

原创 LaTex常用技巧

罗马字母加粗\usepackage{bm}\bm{\theta}或\boldsymbol{\theta}强制换页\clearpage一行放置两幅图（实例为2*3）\begin{figure}[!ht]\centering\subfigure[The 1st iteration.]{\includegraphics[width=5.5cm]{u1.eps}}\quad...

原创 2019年新学期随笔

2019年发生了两件大事1.我结婚了，因为爱情，恋爱207天我们火速领证了2.我和我的先生在北欧的一个国家里开始了留学之旅关于爱情，和所有俗套的爱情一样，我很幸运找到了一个可以包容我所有缺点的爱人，而且在我看来他既优秀又高颜值，特别有我喜欢的大双眼皮高鼻梁，可以为我的后代提供优质的基因。关于留学，我是通过国家留学基金委（CSC)的审核，拿到了两年的联合培养的博士的资助，工资大概是国内的4-...

转载 雷达通信术语中英文对照

A/BADC Analog to Digital Converter模数转换器，将模拟(连续)信号转换成数字(数字)信号的电子设备。AESA Active Electronically Scanned Array有源电子扫描阵列，每个辐射阵元都配装有一个发射/接收组件，每一个组件都能发射和接收电磁波。AEWR Airborne Early Warning Radar机载预警雷达，预...

原创 Matlab常用绘图命令

选择图形窗口命令1.打开不同的图形窗口命令 figurefigure(1);figure(2);figure(n);用来打开不同的图形窗口，以便绘制不同的图形2.图形窗口拆分命令 subplotsubplot(m,n,p):分割图形显示窗口，m表示上下分割个数，n表示左右分割个数，p表示子图编号坐标轴相关命令1.axis([xmin xmax ymin ymax])表示x轴最小值...

原创 中文核心的套路模板

论文大纲一、摘要本文就……问题探讨了……的解决方案。（……问题略写，……解决方案用30个字左右写完）。不同于之前在解决该问题的（A.B.C……）方法，在解决该问题时遇到的……困难，本文提出的方法通过……的步骤（算法里的步骤都要简洁明确的写出来），即使在遇到……（之前说过的困难），也会取得……的效果，达到……的指标/目的（重点突出）。最后套话：仿真实验……，验证了……。二、引言（第一自然段）...

原创 压缩感知里的稀疏为什么是非线性模型

Sparsity is a highly nonlinear model，since the choice of which dictionary elements are used can change from signal to signal[1].[1]Devore R A. Nonlinear approximation[J]. Acta Numerica, 1998, 7(7):51...

原创 CS采样与传统采样方式的对比

从信号类型上看传统采样：长度不限，连续时间信号CS：有限长度信号从采样方式上看传统采样：直接在时间点上采样CS：在经过函数变换后的（内积形式上）采样，是一种线性测量采样从信号重构方式上看传统采样：插值，是一种线性的处理CS: 非线性方法待补充...

原创 宽带信号的DOA估计学习笔记（四）：空间分辨率

时域/频域傅里叶限制指的是傅里叶变换的积分区间不是[−∞,+∞][-\infty,+\infty][−∞,+∞]，而是[−t,+t]/[−ω,+ω][-t,+t]/[-\omega,+\omega][−t,+t]/[−ω,+ω]空域傅里叶限阵列的物理孔径限，线阵就是阵列长度，面阵就是阵列（长*宽中的长，圆阵中的直径）。提高空域处理精度的办法：增加阵列的物理孔径（等效于减小波束），也就是增加...

原创 一篇闲扯——500毫升的水为什么老是想找1000毫升的容器装自己

这就是一篇闲扯，因为我没有其他地方发泄这个情绪我想写的就是最近频繁有人让我给他/她介绍对象的事我现在最害怕听到别人让我给他们介绍对象了。先说男生，男生比较简单。从始至终要求就很简单——长的漂亮。长的漂不漂亮不是我说的算的，他们觉得漂亮就行了（至少要跟我差不多水平的，但是还不能像我这么高以及这么高学历的）。我听过最多的话就是（我怎么也得找个你这样的），这句话听起来甚是奇怪，我这样的怎么了，我虽...

原创 宽带信号的DOA估计学习笔记（三):宽带信号阵列接收模型

在宽带阵列信号处理中，阵列接收到的信号包络不再恒定，相位延迟和时延不再是简单的线性关系。宽带时域模型MMM个阵元组成的阵列，接收空间中PPP个宽带信号，入射角度分别为θ1,θ2,⋯ ,θP\theta_1,\theta_2,\cdots,\theta_Pθ1​,θ2​,⋯,θP​，则第kkk个阵元接收到的信号可表示为xk(t)=∑i=1Psi[t+τk(θi)]+n...

原创 宽带信号的DOA估计学习笔记（二): 信源数估计

一.窄带方法（高白噪声）信息论方法只能对独立信号源的总数作出估计统一表达式J(k)=L(k)+P(k)J(k)=L(k)+P(k)J(k)=L(k)+P(k)，其中L(k)L(k)L(k)是对数似然函数，P(kP(kP(k)是罚函数nnn为待估计的信号源数，LLL为采样数，Λ(n)\Lambda(n)Λ(n)为似然函数AIC准则（Akaike Information Theore...

原创 宽带信号的DOA估计学习笔记（一）

1.宽带信号的定义：https://blog.csdn.net/weixin_38452468/article/details/888198782.宽带信号与窄带信号在公式表达上的区别：窄带信号的导向矢量：角度宽带信号的导向矢量：角度+频率分布加粗样式→ 窄带信号的DOA估计方法无法完全应用在宽带信号处理上。3.宽带信号的DOA估计方法的研究现状参考文献：【1】 G. Su, ...

转载 窄带信号的三种定义

原创 我的第一篇SCI——辛酸过程的回忆录和给各位的一点SCI写作小预告

写在最前面：这篇论文被接收的时候，已经大修了无数次（自己修改）拒稿三次（两次杂志社认为和他们的杂志不是很相符，一次是审稿人觉得我没有数学证明（在下真是无能），写的也不高大上），去年12月初完稿到现在已经过去了快一年，写一个回忆录和一个SCI论文写作流程的总结。现在看来，我这篇文章很垃圾，真的很垃圾，无论是方法还是语言，我甚至一度认为，如果有一天我成为学术界的牛人，那么这篇论文就会成为我的屈辱史...

原创 稀疏表示和DOA估计的联系与区别(the link and the gap)

DOA估计跨过了漫长的二阶统计量算法（MUSIC,ESPRIT）来到了稀疏表示的新时代，现在总结下二者的联系与区别。1.稀疏表示用的是很多个向量组成的矩阵当作字典（我就把他们理解成矩阵论里的基，然而基是必须独立的，字典里的每一列可以是相关的）来线性的表示测量值，比如我现在有10个原子组成的字典[a1,a2,⋯ ,a10][\mathbf{a}_1,\mathbf{a}...

原创 2018新学期随笔

2018.9.3，博士第二年的开学第一天。 其实我已经在实验室呆了三周了，但是生活需要一点仪式感所以写篇博客记录下。1992年的我已经27岁了，未婚未育，独生女，长的还行，学习也就那样，老师抬举硕博连读走到了今天。回头想当年收到研究生录取结果的激动心里暗暗发誓一定要好好努力的想法现在被减了大半，很多时候都只想偷懒。虽然已经对DOA有个一个相对广泛的了解但是很多细节都不是很清楚，或者说...

原创 解惑——为什么导向矢量和噪声子空间正交

相信所有研究DOA方向的同学一定看过一个理论，就是子空间算法中导向矢量张成的子空间与噪声子空间正交。由此，若干个子空间类算法也因此得到了各种求得DOA估计值的方法，那么你是不是跟我一样，迷惑为什么导向矢量和噪声子空间就正交了呢？为了更好地帮助理解，我这篇博客也必须使用一系列的公式推导，首先简单说一下各个符号代表的含义。设有PPP个阵元（天线），MMM个信号，NNN快拍，信号来波方向为{θ1...

原创 矩阵求逆引理（matrix inversion lemma)

论文：Li J, Stoica P. An adaptive filtering approach to spectral estimation and SAR imaging[J]. IEEE Transactions on Signal Processing, 2002, 44(6):1469-1484.问题描述： 已知 求Q^(ω)−1Q^(ω)−1\hat{Q}(\omega)^{...

原创 Capon filter的推导

最近在看的论文总是绕不过这个概念 但是公式（2）到底是怎样来的呢？1.不得不提的波束形成 这里引用一篇硕士论文。 [1]喻敏. 声矢量传感器的Capon方位估计[D].哈尔滨工程大学,2004. 那么看到这里，相信公式2-18前都是相当容易理解的，那么公式2-18是如何推导到2-19这一步的呢？可以看到现在的问题是解决 minω ωHR^ωs.t...

原创 0-1规划问题

0代表放弃，1代表选取，可用于决策问题%2018.5.23%intlinprog函数vsbintprog%求解0，1规划问题%max z=x1+2x2+2x3-6x4-4x5%s.t 3x1+2x2-x3+x4+2x5<=5% 2x1+4x2-2x3-x4-2x5<=5% xi=0或1（i=1,2,...,5)f=-[1,2,2,-6,-4]';a=...

原创 linprog函数

%2018.5.23%linprog函数使用练习%用函数linprog求线性规划%使用要点：%1.目标函数f要写成列矩阵的形式%2.约束要写成<=的形式%3.函数是为求最小值设置的%[x,fv,exitflag,output,lambda]=linprog(f,a,b,aeq,beq,lb,ub)%x表示最优解，以列矩阵表示%fval为最优值%exitflg=1有最优...

原创 CVX error:Conversion to double from cvx is not possible解决方案

昨天在做MATLAB仿真的时候遇到了这个问题，真的是解决了好几个小时。背景：用CVX工具箱做蒙特卡洛实验（就是多次重复实验）报错：如题百思不得其解。解决过程1.CVX USER'S GUIDE;他说要把变量声明一下，写成expression variables的形式，于是我就把提示报错行的那块在cvx_begin后声明了一下，但是依然报错；2.继续搜百度，在cv

原创 阵元校正

参考文献只贴幅度误差校正，相位误差校正类似。%2017.22.12%code author：WYN%function:ULA gain error calibrating%ref:%paper author:A. Paulraj and T. Kailath,%title:Direction of arrival estimation by eigenstruct

原创 unique函数去除重复行不排序

unique函数在使用时会自动升序排序，但是在实际应用中，我们并不希望改变原有矩阵的顺序，只是希望去除重复行。上段代码可以实现去除重复行而不对原有矩阵进行排序。

原创 Khatri-Rao积（KR积）

输入：A,B两个矩阵，要求两矩阵列数一致输出：维度为（A行数*B行数）*列数的矩阵。计算过程：A,B矩阵相同列做kron积运算，并逐列摆放组成结果矩阵。举例：A=[1 2;3 4],B=[5 6;7 8]其中AB均为2*2的矩阵，得到（2*2）*2=4*2的矩阵。第一列：用A的第一列和B的第一列做kron积运算[1*5=51*7=73*5=153*7=21]

原创 L1-SRACV

%2017.7.11%L1-SRACVclcclear allclose allC = 3e8;Fc = 100e6; %信号频率为100MHzlamda = C/Fc; %波长L = 200;

原创 施密特正交化

[row,col] = size(X);if row < col fprintf('行数小于列数，不能规范正交化！\n'); return;endY = zeros(row,col);Y(:,1) = X(:,1);for i = 1:col Y(:,i) = X(:,i); for j = 1:i - 1 Y(:,i) = Y(:,i)

原创 Joint SVD

%2017.6.24%JSVD%Joint SVD of 2 cross_correlation matrices to achieve automatic paring in%2-D angle estimation problemsclcclose all;clear allMx=10;

原创 CCM-ESPRIT，L型阵列二维角度配对算法

论文：贴代码%2017.6.21%L型阵列配对%Pair-matching method for estimating 2d angle of arrival with a cross%correlation matrixclcclose all;clear allMx=10;

原创 PM算法，不需要奇异值分解，快速获得噪声子空间

%2017.6.22%PM算法，不需要svd，快速获得噪声子空间。clcclose all;clear allM=10; %x轴阵元数d_lamda=1/2; %阵元间距

原创 root-MUSIC算法

%2017.6.15%ROOT_MUSIC ALOGRITHM%DOA ESTIMATION BY ROOT_MUSICclc;clear all;close all;K=2; %信源数M=8;

原创 L1_SVD方法的稀疏矩阵的DOA估计

写点题外话，最近，一直在研究Malioutouv的这是他毕业论文摘出来的一篇论文，他的毕业论文下载地址：D. M. Malioutov. (2003, Jul.) A Sparse Signal Reconstruction Perspective for Source Localization With Sensor Arrays. M.S. thesis, Mass. Inst. Techno...

原创 LP范数球

一维情况%2016.6.5%author:Lola%description:p范数clcclose allclear allx1= -2:0.01:2;x=abs(x1);p=[0.5 1 2];y1=x.^(p(1));y2=x.^(p(2));y3=x.^(p(3));plot(x1,y1,'c:','linewidth',2);hold onplot(x1

原创 等距网格细化的单快拍稀疏矩阵的DOA估计

方法：这里采用的是等距划分。想了很久也试了很多，最后是二分法给我的启发二分法的方法参照：function [i,x,fx]=erfenfa_new_use1(f,bot,top,err) %输入：f为要求解的方程函数表达式，bot为求解区间的下界，top为求解区间的上界，err为所要求的误差范围 %输出：i为二分法求解的次数，x为最终的根，fx为方程的根x对应的函数值（精度

原创 利用CVX工具箱实现多快拍采样的稀疏矩阵DOA估计（多快拍视为独立的若干个单快拍，结果求均值）

参考文献：同上篇方法：论文中给出的仿真结果：实验条件：入射角度：60°，70°，108°。贴代码：%date:2017.6.2%author:Lola%处理多点采样，将多点采样中的每一点看成是独立的单点采样，从而进行稀疏矩阵的DOA估计clc clear allclose allM = 8;

原创 利用CVX工具箱实现单快拍的稀疏矩阵DOA估计

最近在学习压缩感知方面的知识，正好自己的研究方向是DOA估计，看到相关方法就拿来试一试，现在把过程记录下来吧，作为一个小菜鸡，成功的那一秒钟真想给自己鼓鼓掌。参考文献：D. Malioutov, M. Cetin and A. S. Willsky, "A sparse signal reconstruction perspective for source localization with

原创 MATLAB logspace( )

一直以来都是看别人写的博客，觉得能把自己所学所知以文字的形式表达清楚，是一件并不容易的事情，作为一个初入博士生涯的小菜鸡，也要开始自己的博客之路。今天在学习CVX工具箱手册的时候，看到了一个全新的logspace函数，现将此函数的说明及使用记述如下。首先在MATLA重点内容B中 help logspace； 得到 翻译过来就是“对数间隔矢量” 它有两种使用形式： 1.logspace(