matlab自带函数mafdr,当ttest数较多时,可直接用[FDR, Q]=mafdr(P):但是Storey procedure在p值少于1000个时会崩溃,此时应改用BH FDR方法:mafdr(P,'BHFDR', true).该方法对于少量ttest更稳健,但是更保守. 参考链接: https://stat.ethz.ch/pipermail/bioconductor/2014-January/056992.html http://home.52brain.com/forum.php

MATLAB中“fitgmdist”的用法及其GMM聚类算法 作者:凯鲁嘎吉 - 博客园http://www.cnblogs.com/kailugaji/ 高斯混合模型的基本原理:聚类——GMM,MATLAB官方文档中有关于fitgmdist的介绍:fitgmdist.我之前写过有关GMM聚类的算法:GMM算法的matlab程序.这篇文章主要应用MATLAB自带的函数来进行聚类. 1. fitgmdist函数介绍 fitgmdist的使用形式:gmm = fitgmdist(X,k,Name,V

MATLAB中冒号的用法 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ matlab中冒号代表步长,用实例来说明 >> A=[1 2 3 4;5 6 7 8;9 10 11 12;13 14 15 16;17 18 19 20] A = 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 >> B=A(1:2:5,2:3) B = 2 3 10 11 18 19 >> C

MATLAB中冒号的用法解析 1.: 表示所有的意思. (1)如:a(1,:) 表示a的第1行,示例: 结果: 同样的如果a(2,:)表示a的第2行 (2)反过来,a(:,2) 表示a的第3列,示例: 结果: (3)如果不写数字即a( : ),则表示全部元素,并且返回的结果a变为列向量.示例: 结果: 2.: 表示到,表示一个范围. (1)如a(1:2,3) 表示第1行到第2行的第3列元素.示例: 结果: (2)a(1:2,:) 表示第1行到第2行的所有元素.示例: 结果: (3)同样有a(:,

https://blog.csdn.net/wangyang20170901/article/details/78745587 MATLAB中mean的用法 转载仙女阳 最后发布于2017-12-07 20:39:27 阅读数 23473  收藏 展开 1.函数功能     求数组的平均数或者均值2.使用方法    ① M = mean(A)     返回沿数组中不同维的元素的平均值.     如果A是一个向量,mean(A)返回A中元素的平均值.     如果A是一个矩阵,mean(A)将其中

对于一个变量a,matlab中定义需要这样 syms a: 定义之后就可以写由变量组成的式子,比如 c=(1+a^-1)^5; 可以用expand(c) 就能把c展开成多项式,每一项的系数就可以看到. freqz的用法 MATLAB提供了专门用于求离散系统频响特性的函数freqz(),调用freqz()的格式有以下两种: l        [H,w]=freqz(B,A,N) B和A分别为离散系统的系统函数分子.分母多项式的系数向量,N为正整数,返回量H则包含了离散系统频响 在 0——pi范围内

struct在matlab中是用来建立结构体数组的.通常有两种用法: s = struct('field1',{},'field2',{},...)  这是建立一个空的结构体,field1,field2是域名称 s = struct('field1',values1,'field2',values2,...)  这是在建立结构体的时候就赋值,values1,values2就是针对field1,field2所赋的值. 例如:s = struct('type',{'big','little'},'c

@是Matlab中的句柄函数的标志符,即间接的函数调用方法. 1 句柄函数 主要有两种语法: handle = @functionname handle = @(arglist)anonymous_function handle = @functionname:返回一个特别的Matlab函数句柄,它提供了一种间接访问函数的方式,也被成为函数的函数(function functions),是一种标准的Matlab数据类型.在C/C++中,有个类似的用法称为引用(使用标识符&),引用只是它绑定的对象

Matlab中max函数在矩阵中求函数大小的实例如下:(1)C = max(A)返回一个数组各不同维中的最大元素.如果A是一个向量,max(A)返回A中的最大元素.如果A是一个矩阵,max(A)将A的每一列作为一个向量,返回一个行向量,向量的第i个元素是矩阵A的第i列上的最大值. 如果A是多维数组,max(A) treats the values along the first non-singleton dimension as vectors, returning the maximum v

1.在matlab中 ,infro.mat中存有很多子矩阵(比如:mean_FA.mat, mean_e1.mat和 mean_e2.mat),调出某一个矩阵时,命令行为:load([path,'\',infro.mat],'mean_FA'); 2.load(path) 导出的形式是结构体:而importdata(path)导出的形式是矩阵: 3.选择性读入temp: temp = dir([path,'\','ROICorrelation_FisherZ_*.mat']):

Fspecial函数用于创建预定义的滤波算子,其语法格式为:h = fspecial(type)h = fspecial(type,parameters,sigma) 参数type制定算子类型,parameters指定相应的参数,具体格式为: type='average',为均值滤波,参数为n,代表模版尺寸,用向量表示,默认值为[3,3]. type= 'gaussian',为高斯低通滤波器,参数有两个,n表示模版尺寸,默认值为[3,3], sigma表示滤波器的标准差,单位为像素,默认值为 0

cell元包是matlab中提供的一种数据类型,功能强大. 关于cell的创建: 1.跟一般创建举证一样,直接使用C = {A B D E}这种形式,不过这里把"[]"改成了}"{}" 2.使用cell创建cell,C= {C1 C2},这里C1和C2都是cell类型,直接把C1 C2作为一个cell得到的C3是一个1x2的cell数组:使用C3 = [C1 C2],是使用C1和C2的内容来创建cell数组 3.使用cell函数,A=cell(2,2),创建一个空的

来源:https://blog.csdn.net/qq_15971883/article/details/78334380 转载自:http://blog.csdn.net/u013066730/article/details/56554332 以及http://blog.csdn.net/u010740725/article/details/51557202 功能:对任意类型数组或多维图像进行滤波.用法:B = imfilter(A,H) B = imfilter(A,H,option1,op

一. 读写图像文件 1. imread imread函数用于读入各种图像文件,如:a=imread('e:/w01.tif') 注:计算机E盘上要有w01相应的.tif文件. 2. imwrite imwrite函数用于写入图像文件,如:imwrite(a,'e:/w02.tif',’tif’) 3. imfinfo imfinfo函数用于读取图像文件的有关信息,如:imfinfo('e:/w01.tif') 二. 图像的显示 1. image image函数是MATLAB提供的最原始的图像显示

一.理论准备 matlab函数randn:产生均值为0,方差 σ^2 = 1,标准差σ = 1的正态分布的随机数或矩阵的函数. 用法:Y = randn(n),返回一个n*n的随机项的矩阵.如果n不是个数量,将返回错误信息.       Y = randn(m,n) 或 Y = randn([m n]),返回一个m*n的随机项矩阵. Y = randn(m,n,p,...) 或 Y = randn([m n p...]),产生随机数组(感觉就是三维数组,请看如下例子). 1: >> rand(

一.调用方法X=FFT(x):X=FFT(x,N):x=IFFT(X);x=IFFT(X,N) 用MATLAB进行谱分析时注意: (1)函数FFT返回值的数据结构具有对称性. 例:N=8;n=0:N-1;xn=[4 3 2 6 7 8 9 0];Xk=fft(xn) →Xk = 39.0000           -10.7782 + 6.2929i        0 - 5.0000i   4.7782 - 7.7071i   5.0000             4.7782 + 7.707

一. 功能: 寻找非零元素的索引和值 二.相关函数语法: ind = find(X) ind = find(X, k) ind = find(X, k, 'first') ind = find(X, k, 'last') [row,col] = find(X, ...) [row,col,v] = find(X, ...) 三.说明: index= find(X) 找出矩阵X中的所有非零元素,并将这些元素的线性索引值(linear indices:按列)返回到向量index中. 如果X是一个行向

来源:https://ww2.mathworks.cn/help/matlab/ref/axis.html?searchHighlight=axis&s_tid=doc_srchtitle axis 设置坐标轴范围和纵横比 全页折叠 语法 axis(limits) axis style axis mode axis ydirection axis visibility lim = axis [m,v,d] = axis('state') ___ = axis(ax,___)   说明 示例 ax

参考:https://ww2.mathworks.cn/help/matlab/ref/uigetfile.html?searchHighlight=uigetfile&s_tid=doc_srchtitle uigetfile 打开文件选择对话框 全页折叠 语法 file = uigetfile [file,path] = uigetfile [file,path,indx] = uigetfile ___ = uigetfile(filter) ___ = uigetfile(filter,

写成subplot(m,n,p)或者subplot(mnp). subplot是将多个图画到一个平面上的工具.其中,m表示是图排成m行,n表示图排成n列,也就是整个figure中有n个图是排成一行的,一共m行,如果第一个数字是2就是表示2行图.p是指你现在要把曲线画到figure中哪个图上,最后一个如果是1表示是从左到右第一个位置.下面是两个例子,可加深理解, >> t=0:0.001:1;>> y1=sin(10*t);>> y2=sin(15*t);>>

C = max(A) 返回一个数组各不同维中的最大元素.如果A是一个向量,max(A)返回A中的最大元素.如果A是一个矩阵,max(A)将A的每一列作为一个向量,返回一行向量包含了每一列的最大元素. C = max(A,B) 返回一个和A和B同大小的数组,其中的元素是从A或B中取出的最大元素. 当Amxn是一个矩阵的时候,有以下几种情况:A = fix (rand (5,3)*50)%生成5X3随机矩阵B = A;%矩阵备份一次C = max(max(A))%矩阵最值D = max(A,[],1