wulimmya的博客

前段时间磕盐接触到了希尔伯特频谱，它是一种信号分解方法，1998年提出来的，主旨是把复杂信号分解为简单信号的加权和，就像傅里叶变换小波变换一样，但是他和傅里叶变换等方法的区别是他是纯粹时间域的分解，但是每个子信号却可以表示不同的频率成分，于是可以得到像小波变换那样的时频平面，但是这个方法明显比小波分解冷门的多，而且在我的实验结果里确实远远远远弱于小波分解，不过也算是自己辛苦几天看论文和写代码的成果...

gpuDeviceans = CUDADevice (具有属性): Name: 'GeForce RTX 2060 SUPER' Index: 1 ComputeCapability: '7.5' SupportsDouble: 1 ...

clf; % 清除当前图形窗口所有可见的图形对象，不关闭图形窗口！！！close%关闭当前窗口，等效于close(gcf)

poolobj = gcp('nocreate');% If no pool, create new one.if isempty(poolobj) poolsize = 0; CoreNum=6; %设定机器CPU核心数量 % 在本地虚拟出6台可以运行matlab的工作站，这样用分布式计算工具箱可以进行并行计算（matlabpool这个命令是并行计算工具箱里的）。观察...

CoreNum=6; %设定机器CPU核心数量parpool(CoreNum);>> parallelemotionStarting parallel pool (parpool) using the 'local' profile ... connected to 6 workers.测试：竟然快了接近20倍！！！ticparfor i=1:3, c(:,i) =...

自己总结的，不一定准确或者符号所有期刊的要求，我用来统一规范自己的图，仅供参考不加title，图注代替title的职能横纵坐标标签等文本 20号字体，保证看起来清晰可以axis square就square单位用圆括号括起来，用简写，比如秒就是s，不写seconds网格线用虚线，图像线条用实线线宽尽量最小1磅，不然看不见...

很多论文中的图像都是用的虚线网格哦，图像用实线，比较好看分明set(gca,'gridlinestyle',':');

shading interp; % 对曲面或图形对象的颜色着色进行色彩的插值处理，使色彩平滑过渡三维时频谱图（STFT）， 不用超级丑用了顺眼的多二维时频谱图（STFT），用了才不好看...

左边是imagesc个人觉得pcolor（伪彩色图）的分辨率更高更好看一些具体原因暂时没时间追究，看到网上imagesc处理图片倒是特别好看的

求一维数组的最大值及其索引[value,ind]=max(a)

好看[X,Y] = meshgrid(-2:.2:2, -2:.2:2); Z = X .* exp(-X.^2 - Y.^2); [c h] = contourf(Z);

直接全选可以选到所有的边框，可以同时修改更好看的颜色，我比较喜欢浅灰色的边框比黑色好看选中一条线后按住shift可以同时多选很多其他线条或者边框，同时修改颜色，不需要一个一个的来区域图线形图阶梯图线宽任何在程序中不好设置的比如颜色（代码可控颜色较少），图的类型（区域图）··都可以编辑坐标轴是否加框等matlab画出的图形的所有的东西都是可以直接编辑的...

把三维矩阵的所有页竖着拼接起来，不能用reshapea = rand(2,3,4);b = [];for i=1:4b = cat(1,b,a(:,:,i));end% b = reshape(a,8,3); % 效果不对a=val(:,:,1) = 0.8147 0.1270 0.6324 0.9058 0.9134 0.0975...

矩阵就是一个m行n列的数表说明标量（不管实数还是复数），在matlab中都是以矩阵形式存储的。判断两个数是否相等空矩阵不占存储空间For nonempty arrays, cat(2, A, B) is the same as [A, B], and cat(1, A, B) is the same as [A; B].cat(3, A, B)得到两页，cat(4, A, B)...

markersize,用圆圈画每个点，离散的点之间不连线，不用圈圈的话plot默认把每个离散点用直线连起来linewidth，线宽MATLAB是不可以画出连续函数的图像的，只能用一群离散点（划分更细）去在视觉上像连续的。x1=(0:12)*pi/6;y1=cos(3*x1);x2=(0:360)*pi/180;y2=cos(3*x2);figure(1)subplot(221);pl...

t=linspace(0,10,100);s=sin(2*pi/5*t);noise=0.2*rand(size(t));%xnoise = x + 2.5*gallery('normaldata',size(t),4);%标准正态分布的随机数据矩阵xs=s+noise;y=zeros(size(t));A=[1 -0.9];B=[.05 .06];y=filter(B,A,x...

隐枚举法就是在 枚举过程中引入过滤条件。0-1整数规划是整数规划的特殊情形，部分或全部变量只能取0或者1, 这个约束表示为0≤xj≤1,xj∈N,j=1,2,⋯n0\leq x_j\leq1, x_j \in N, j=1,2,\cdots n0≤xj​≤1,xj​∈N,j=1,2,⋯n或者表示为xj(1−xj)=0,j=1,2,⋯nx_j(1-x_j)=0, j=1,2,\cdots ...

本文参考了维基和minitab网站，只是汇总梳理和强调了相关知识点。Spearman’s rank correlation coefficient，是衡量两个变量的依赖性的 非参数 指标（衡量分级定序变量之间的相关程度），利用单调方程评价两个统计变量的相关性。 如果数据中没有重复值， 并且当两个变量完全单调相关时，斯皮尔曼相关系数则为+1或−1。用途是描述数据间的相关性。可以看到斯皮尔...

matlab画图subplot(211)表示2行1列第一个图，结果今天程序需要画400多个图，到第210个图画完就出错了，错误使用 subplot (line 171)需要使用有效的轴句柄作为输入。出错 gps_interp (line 37) subplot(211)原来是因为211！！！！这时候就必须写成subplot(2，1，1)拉！！哈哈哈有意思...

文章目录先说变换是什么再说傅里叶变换有什么局限简述FTFT因无法处理时变信号而裹足不前STFT为了克服FT的缺点（不能处理时变信号）诞生了小波变换为了克服STFT的缺点（单一分辨率）诞生了仔细说说小波变换的真正含义小波是什么（主要讲母小波）母小波伸缩平移得到子小波连续小波变换先说变换是什么变换就是定义域和值域都是函数集而非数值的函数。 它仍然是一种映射，只是从函数集映射到函数集。比如傅里叶变...

clearclc[x,y]=meshgrid(1:10);% 测量数据h=[0 .02 -.12 0 -2.09 0 -.58 -.08 0 0; .02 0 0 -2.38 0 -4.96 0 0 0 -.1; 0 .1 1 0 -3.04 0 -.53 0 .1 0 ; 0 0 0 3.52 0 0 0 0 0 0; -.43 -1.98 0 0 ...

a是数组逻辑输出0表示假，输出1表示真any（x）返回0，如果x是零向量或者零矩阵一些用于测试的函数

字符串类型可以和double类型运算whos查看工作区所有变量，包括详细信息who只给所有变量，不给信息eps是一个数和它后面的紧接着的那个浮点数的距离，总是正的，也用于测量浮点数的精度。比如双精度的精度就在2−522^{-52}2−52左右。angle函数求得的是弧度制,单位是radNaN not a number非数值量...

matlab是交互式环境用于算法开发，数据可视化MathWorks公式成立于1984!命令clc是清除命令行的内容clear清除工作区的变量命令cd,dir查看matlab安装版本的license，但我是破解的，所以这个license根本没用，是假的。path命令查看matlab的所有搜索路径（一般只有文件区和toolbox）添加搜索路径lookfor命令...

input()函数：clear allclcprompt='Do u want more? Y/N [Y]: ';reply=input(prompt,'s');if isempty(reply) reply='Y';end prompt：提示backslash:’'反斜线clear allclca=magic(5);expression=input('...

0、运算优先级一、算术运算符二、关系运算符format compact可以使得输出比较紧凑，想恢复稀松的输出形式用format looseeps的强大作用：避免发生除零清除工作区某个变量清除两个变量，空格隔开三、逻辑运算符xor(x,y)，x,y均是零或者均是非零则输出逻辑0；有一个为零一个非零则输出逻辑1...

一、单元（cell）数组类型单元数组是一种无所不包的广义矩阵。单元数组的每一个元素叫一个单元。建立空的元胞数组，作用是预先分配到连续的内存空间，以节省内存和提高效率用class()查看元素的数据类型元胞数组的合并删除某个单元reshape函数的使用二、map容器类型容器：container一个map容器是map类的一个对象，里面有多个键值对，keys/value...

一、结构体类型matlab中的结构体和C语言的结构体相似，一个结构体通过字段存储多个不同类型的数据。因此结构体就像是一个数据容器，把多个相关联的不同类型的数据封装到一个结构体对象中。1.直接对结构体的四个字段赋值，创建了一个结构体对象Student。Student是个1*1的结构体数组可以用struct函数创建结构体对象，或者任意尺寸的结构体数组（数组中所有结构体对象的字段相同）。...

matlab中，一个m文件中只能有一个主函数，即与m文件同名的函数，可以有其他多个函数，但都不是主函数。函数的调用分为直接调用和间接调用两种，直接调用只允许同一个m文件中的函相互调用；但函数句柄（间接调用）允许跨m文件的调用，也可以调用库函数等。...

我可能真的太闲了吧···在网上看到这图觉得挺好看：t=0:pi/360:2*pi;x=sin(t);y=cos(t);z=2*x.^2+y.^2;figureplot3(x,y,z,'Color','r','LineWidth',2);xlabel('x');ylabel('y');zlabel('z');title('z=(cos t)^2+2*(sin t)^2');ax...

数组对于MATLAB是非常重要的数据结构，数组运算是MATLAB完成计算的基础因为在MATLAB中，所有数据！！都要以数组的形式被存储！！数组的元素可以是任何数据类型，比如数值，字符串矩阵是特殊的数组，一定要明白矩阵在MATLAB中也是以数组形式存储的由于MATLAB是面向对象的，所以数值数组是最重要的内建数据类型创建数组：clear allformat compact % 紧凑显...

clear all clcformat compacta=hadamard(2)b=rosserc=hilb(4)d=invhilb(4)e=c*d % cd可逆f=toeplitz(2:5,2:1:6) % 输入列和第一个元素必须和输入行的第一个元素相同，否则会报错g=toeplitz(2:5) % 对称a = 1 1 1 -1b ...

文章目录矩阵合并删除行、列：（赋予空矩阵/数组）共轭转置矩阵的加减乘除幂运算矩阵乘法除法矩阵元素的查找，排序，求和，求积，差分矩阵元素的查找矩阵元素排序：向量元素排序矩阵元素求和矩阵元素求积矩阵元素差分矩阵合并a=zeros(2,3)b=ones(2,4)c=7:9d=[a b] % 水平合并，行数相同e=[a;c] % 垂直合并，列数相同a = 0 0 ...

文章目录一、基础概念插值是什么拟合是什么插值和拟合的相同点插值和拟合的不同点二、常用的基本插值方法高次多项式插值拉格朗日多项式插值牛顿插值差商矩阵低次多项式插值（不易震荡）分段线性插值Hermite插值三次样条插值一、基础概念插值是什么求出过已知个有限数据点的近似函数。拟合是什么求出一个不要求过已知数据点的近似函数，不要求过那些数据点，只要求在这些点上的总偏差最小。插值和拟合的相同点...

多项式拟合：x=1:10;y=randn(size(x));p=polyfit(x,y,7); % 7阶多项式，返回拟合多项式的降幂系数x1=0:.1:10;y_hat=polyval(p,x1);plot(x1,y_hat,x,y,'r*')或者直接使用matlab的拟合工具箱写好x, y值x=1:10;y=randn(size(x));打开工具箱或者在命令行输入...

文章目录顺序结构if-else-end分支结构switch-casetry-catch循环结构for循环while循环顺序结构if-else-end分支结构switch-casetry-catch循环结构for循环while循环

clear all clcfid=fopen('test.txt','wt'); % 打开一个文件，返回一个指标fidfprintf(fid,'This is the database of class 1. \n'); % 直接写入一句话% 每赋值一次，写入一次name='Alice';type=1;x=3.4;y=32;ans='yes';fprintf(fid,'%s Type ...

文章目录变量本质上，变量代表一段可操作的内存，是内存的符号化表示局部变量全局变量永久变量特殊变量关键字变量程序中定义变量用于保存和处理数据。我们需要使用内存时，就定义某种数据类型的变量。编译器会根据数据类型分配一定的内存空间，程序通过变量名访问内存。本质上，变量代表一段可操作的内存，是内存的符号化表示注意matlab的变量只能以字母开头，不可以下划线开头，这是和其他语言的区别。&gt...

文章目录文件操作指令文件流控制指令matlab的工具库里有大量预定义的m文件，如addpath.m, newplot.m等，他们在安装matlab的时候就直接被存放在安装目录里。文件操作指令对于m文件，有大量的操作和控制指令：what 列出安装目录下所有的文件>> what当前文件夹 C:\Program Files\MATLAB\R2016b\bin 中的 MATL...

定义： 含有大量0元素的矩阵不是稀疏矩阵，则称满矩阵。full matrixmatlab对稀疏矩阵有专门的的运算命令，为了节省存储空间和运算时间。怎么节省存储空间？一个m×nm\times nm×n的满矩阵，元素的数据类型是double，则总共需要m×n×8m\times n\times 8m×n×8个Bytes。如果它满足稀疏特性，即大量元素为0，则可以存储为稀疏矩阵，只需要存储所...