LucyLiHHU的博客

原创 回归相关系数r

残差平方和y的离差平方和x的离差平方和回归相关系数r

原创 matlab之xlsread读取数据均为0

跑数据的时候，xlsread读取数据均为0，以为是xlsread语法错误原来是数据格式错误。当excel 显示的时候左上角有实心小三角就是代表文本格式解决办法：将“文本格式存储的数字”转换为“数字”...

原创 matlab之xlswrite

原创 Matlab曲线拟合之polyfit及交互式曲线拟合工具

1、polyfit     p=ployfit（x，y，n）；     x——自变量；     y——因变量；     n——多项式阶数；例如：     对于ln（1+x）在【0，1】内的采样数据进行多项式拟合；>>x=0:0.1:1.0;>>y=log(1+x);>>p=ployfit(x,y,3);运行结果为：p=...

原创 求多项式的值之polyval 和polyvalm

1、pv=polyval（p,a);     pv=polyval（p,A);     p——多项式表示；     a——标量输入；     A——矩阵输入；当A为M×N矩阵时，函数返回值pv也是M×N矩阵，且pv（i，j）=polyval（p，A（i，j））。2、pv=polyvalm（p，A）；      A——只能是N阶方阵。 ...

转载 excel之mod、函数

mod函数是一个求余函数，其格式为： mod(nExp1,nExp2)，即是两个数值表达式作除法运算后的余数。特别注意：在EXCEL中，MOD函数是用于返回两数相除的余数，返回结果的符号与除数（divisor）的符号相同。语法：MOD(number,divisor)参数：Number 为被除数。Divisor 为除数。如果 divisor 为零，函数 MOD 返回值 为#DIV/0...

原创 等值线绘制（等高线）

绘制等高线1、Spatial Analyst Tools——Surface——Contour2、Input surface选择要生成等高线的栅格数据集3、Contour interval设置等高距4、Base contour指定等高线基准高程5、Z factor设定高程变换系数6、Output features输入结果文件存储与文件名称...

转载 ArcGIS影像多种裁剪方法

在互联网上下载的遥感影像都进行过分幅处理，下载下来的影像多是规则的四方形，而在进行遥感影像研究时，多是针对特定区域来进行，比如研究北京市的遥感影像，不在北京市范围内的影像对于研究者就没有利用意义，如果不进行裁剪处理，而是对整幅影像进行分析处理，就会增加很多不必要的工作，既浪费时间又浪费资源，所以根据研究区的范围裁剪出研究区的影像显得十分重要。本文的目的在于介绍arcgis中的多样裁剪功能，以北

转载 ArcGIS栅格转矢量

1. 栅格数据向矢量数据的转换　　栅格向矢量转换处理的目的，是为了将栅格数据分析的结果，通过矢量绘图装置输出，或者为了数据压缩的需要，将大量的面状栅格数据转换为由少量数据表示的多边形边界，但是主要目的是为了能将自动扫描仪获取的栅格数据加入矢量形式的数据库。由栅格数据可以转换为 3 种不同的矢量数据，分为点状、线状和面状的矢量数据。下面以栅格数据转换为面状矢量数据为例进行说明，其他两种转换操作大

原创 ArcGIS 图像合并至新栅格图层（Mosaic To New Raster和Mosaic）

1、介绍 Mosaic To New Raster1）Input Raster。输入图层，没什么好说的，把要合并的图层加进来就好了；2）Output Location。这个是结果文件保存的文件夹，注意，是文件夹，我第一次的时候选错了，一直提示我出错，差点气死我；3）Raster dataset name with extension。这里可以写输出文件的文件名了，带后缀名的。只支持ERDAS IM...

转载 深入浅出空间参考——对ArcGIS空间参考的理解-- 转载

摘要：空间参考描述了一个地物在地球上的真实位置。为了正确的对位置进行描述，需要引入一个可供测量和计算的框架，使得大地测量的结果能够在这个框架上进行描述。而地球是一个不规则形状的椭球体，那么使用什么样的方法来模拟地球的形状，又该如何将球面上的坐标投影在平面的地图上？这就需要先了解大地水准面、参考椭球体、基准面的概念，和它们之间的关系。另外，本文还对我国常用的北京54和西安80两种坐标系统进行了详细的...

转载 ArcMap开始编辑遇到一个或多个有错误或警告的图层，并且无法继续

开始编辑后可能出现的错误：如果 ArcMap 在所选数据上启动编辑会话时遇到问题，将弹出一个对话框以提供附加信息。您可能会收到错误、警告或信息消息。出现错误 时用户不可以启动任何编辑会话。只有解决了问题，才能编辑数据。如果没有可编辑的数据源、许可授权方面出现问题（例如，在使用 ArcGIS for Desktop Basic 时尝试编辑某些类型的地理数据库要素）或某些图层存在于底图图层内，就会收到...

转载 ArcGIS栅格数据SetNull IsNull Con函数

当两个栅格进行叠加，有时会有一部分没有数据，即用identify点击该区域，Value为NoData，而不是像其他非空区域一样有值。如果不进行处理的话，nodata数据跟其他图层做运算是会变为nodata。所以nodata区域要赋予0值，因为nodata任何数=nodata，因此当需要对栅格空值数据赋值时，要采用条件查询函数将NoData的地方进行赋值。以下以arcgis10.0为例：（1）方法是...

原创 Matlab的xlsread数据读取问题（line247）

问题原因： 由于你的系统中装了一个福昕阅读器（或其他pdf阅读器），该程序占用了excel的com端口，所以当Matlab再去调用这个端口时就会出现异常。解决方法：打开任意一个Excel文件，点击文件——选项——加载项——点击下方的管理旁边的下拉菜单，选择COM加载项——点击转到——把福昕阅读器的前面的勾去掉，然后确定。你重新回到MATLAB读取Excel文件试试，应该就可以了。如果还不行，那我就

转载 matlab函数汇总

函数表Ⅰ-1  概率密度函数函数名 对应分布的概率密度函数betapdf 贝塔分布的概率密度函数binopdf 二项分布的概率密度函数chi2pdf 卡方分布的概率密度函数exppdf 指数分布的概率密度函数fpdf f分布的概率密度函数gampdf 伽玛分布的概率密度函数geopdf 几何分布的概率密度函数hygepdf 超几何分布的概率密度函数nor

转载 tic/toc/cputime测试时间的区别

cputime测试代码运行时间可能不及tic/toc准确是众所周知的事情。本文并非旧话重提，而是期望起到抛砖引玉的效果，从而找到cputime与tic/toc内在的区别。望不吝赐教！用tic/toc测试如下代码ticpause(10);t=toc复制代码会发现matlab处于busy状态并持续10秒，然后输出:t =   10.0007

转载 Matlab统计工具箱中各类概率分布函数使用方法介绍

1.      Matlab Statistics Toolbox 中概率分布函数介绍Matlab Statistics Toolbox 提供了对几乎所有概率分布的支持，可以方便产生服从各类分布的随机数及其PDF/CDF函数。其中，既提供了针对所有分布的通用分布函数，其参数是指定的特定分布；也提供了针对特定分布的专用分布函数，其函数名直接针对某种特定分布，故该函数的无需再提供分布特性参数。

原创 matlab画图函数汇总（三）

8、fplot在指定的限度内绘制函数的图像函数fplot在指定的限度内绘制函数的图像。该函数的形式必须是y=f（x），其中x是一个向量，它的范围定义绘制的i俺都，y是一个与x相同维数的向量且包含函数在点x处的值。如果函数对于给定的x返回多于一个值，则y是一个矩阵，矩阵的列包含f（x）的每一个分量。【语法】fplot（‘function’，limits）在由limits定义的两个限度

原创 matlab画图函数汇总（二）

4、ezplot容易使用的函数绘制器【语法】ezplot（f）绘制表达式f=f（x）在默认区间-2π对于隐式含义的函数f=f（x，y），绘制在默认区间-2πezplot（f，[min,max]）在区间minezplot（f,[xmin,xmax,ymin,ymax]）绘制区间xminezplot（x,y）绘制在默认区间0ezplot（x,,[tmin,tmax]）在

原创 matlab画图函数汇总（一）

1、ezcontour容易使用的等值线绘制器，会自动添加标题和轴的坐标【语法】ezcontour（f）绘制f（x,y）的等值线，这里f是一个代表两个变量，如x和y的数学函数的字符串。ezcontour（f，domain）在指定的区间上绘制f的曲线，如-4ezcontour（…，n）使用n×n的网格绘制f在默认定义域区间上的图像，如果不指定，n的默认值是60.2、ezmesh

转载 matlab sort函数

结构(struct)数组要在MALTAB中实现比较复杂的编程，就不能不用struct类型。而且在MATLAB中实现struct比C中更为方便。4. 3.1 结构数组的创建MATLAB提供了两种定义结构的方式：直接应用和使用struct函数。1. 使用直接引用方式定义结构与建立数值型数组一样，建立新struct对象不需要事先申明，可以直接引用，而且可以动态扩充。比如建立一

转载 matlab 结构数组

结构(struct)数组要在MALTAB中实现比较复杂的编程，就不能不用struct类型。而且在MATLAB中实现struct比C中更为方便。4. 3.1 结构数组的创建MATLAB提供了两种定义结构的方式：直接应用和使用struct函数。1. 使用直接引用方式定义结构与建立数值型数组一样，建立新struct对象不需要事先申明，可以直接引用，而且可以动态扩充。比如建立一

转载 matlab unifrnd函数

unifrnd可以创建随机的连续均匀分布的数组。1.R = unifrnd(A,B)returns an array of random numbers chosen from the continuous uniform distribution on the interval from A to B. The size of R is the common size of A and

转载 matlab repmat函数

B = repmat(A,m,n)B = repmat(A,[m n])B = repmat(A,[m n p...])这是一个处理大矩阵且内容有重复时使用，其功能是以A的内容堆叠在（MxN）的矩阵B中，B矩阵的大小由MxN及A矩阵的内容决定，如果A是一个3x4x5的矩阵，有B = repmat(A,2,3)则最后的矩阵是6x12x5例如：>>B=repmat( [1 2;3

转载 matlab 匿名函数

匿名函数（anonymous function）是MATLAB 7.0版提出的一种全新的函数描述形式，和内联函数类似，可以让用户编写简单的函数而不需要创建M文件，因此，匿名函数具有inline函数的所有优点，并且效率比inline函数高，同时还具有一些独有的特点。       事实上匿名函数很像C语言中的宏，这个从匿名函数的构造就能够看得出来，但是它比宏的功能更加强大，就是可以支持变量替换，这

转载 matlab中关于函数文件和脚本文件使用实例

matlab中关于函数文件和脚本文件使用实例如下：一般来说分两种情况：如果首行(不包括注释)不是以function开头的，这种叫做脚本或者命令。如果首行(不包括注释)以function开头，这种叫做函数。对于脚本，可以直接复制到命令行运行，也可以保存成m文件执行。而函数必须保存成m文件，然后在脚本里或者其他函数里调用。所以m文件也就分成脚本文件和函数文件。区别就看开头是否为

转载 sort

sort函数功能把数组元素按升序或降序排列 使用方法B = sort(A)把数组A沿不同维度上的元素按照升序排列。 如果A是向量，sort(A) 对A中元素按照升序排列。如果A是矩阵，sort(A) 对A按每一列元素按照升序排列。如果A是多维数组，sort(A) 对沿着第一非单元集的元素像向量一样处理并返回一个升序排列的向量数组。如果A是字符串元胞数

转载 差分进化算法

差分进化算法在 1997 年日本召开的第一届国际进化优化计算竞赛(ICEO)]表现突出，已成为进化算法(EA)的一个重要分支，很多学者开始研究 DE 算法，并取得了大量成果．2006年 CEC 国际会议将其作为专题讨论，由此可见 DE 算法已成为学者的研究热点，具有很大的发展空间.每个个体的优劣程度根据已定义好的适应度函数来评价，这与被解决的问题有关．基本的差分进化算法实现过程如下：

原创 几种求函数最值的算法

1、遗传算法2、粒子群算法3、模拟退火4、蚁群算法

转载 matlab中text 函数在显示字符串时的使用方法

功能 在当前轴中创建text对象。函数text是创建text图形句柄的低级函数。可用该函数在图形中指定的位置上显示字符串。用法 text(x,y,'string')在图形中指定的位置(x,y)上显示字符串stringtext(x,y,z,'string') 在三维图形空间中的指定位置(x,y,z)上显示字符串stringtext(x,y,z,’string’.'PropertyName

转载 matlab的findobj函数

句柄图形是对底层图形函数集合的总称，它实际上进行生成图形的工作。这些函数一般隐藏于 M 文件内部， 但是它们非常地重要， 因为程序员可以利用它对图象或图片的外观进行控制 MATLAB图形系统是建立图形对象的等级系统之上，每一个图形对象都有一个独立的名字，这个名字叫做句柄。每一个图形对象都有它的属性，我们可以通过修改它的属性来修改物体的行为。例如，一条曲线是图形对象的一种。曲线对

转载 matlab的line函数

line([起点横坐标，终点横坐标],[起点纵坐标，终点纵坐标])，例line([1,2],[3,4])将画出（1,3）到(2,4)的一条直线，而不是(1,2)到(3,4)。当line(a,b)中，a b是相同大小的矩阵时，将会在对应的每一列做一条直线。当line(a,b,c)时，相应地会在三维图中画一条线。a b c 均为2XN矩阵。

转载 matlab的size函数

size（）：获取矩阵的行数和列数（1）s=size(A),         当只有一个输出参数时，返回一个行向量，该行向量的第一个元素时矩阵的行数，第二个元素是矩阵的列数。（2）[r,c]=size(A),         当有两个输出参数时，size函数将矩阵的行数返回到第一个输出变量r，将矩阵的列数返回到第二个输出变量c。（3）size(A,n)如果在size函数的输入参

转载 NORMINV

NORMINV是返回指定平均值和标准偏差的正态累积分布函数的反函数计算公式为NORMINV(probability,mean,standard_dev)，其中Probability：正态分布的概率值；0Mean： 分布的算术平均值Standard_dev ：分布的标准偏差；Standard_dev> 0；如果 mean = 0 且 standard_dev = 1，函数 NOR

转载 matlab中subplot函数的功能

subplot 功能分割figure，创建子坐标系语法h = subplot(m,n,p) or subplot(mnp)       subplot(m,n,p,'replace')       subplot(m,n,P)       subplot(h)       subplot('Position',[left bottom width height])

转载 matlab坐标轴函数 semilogx loglog

主要是学习semilogx函数，其中常用的是semilogy函数，即后标为x的是在x轴取对数，为y的是y轴坐标取对数。loglog是x y轴都取对数。例子，clc;clear;close all;x = 0:.1:10;y = 2*x+3;subplot(211);plot(x,y);grid onsubplot(212);semilogy(x,y);grid

转载 Matlab num2cell函数的用法

Matlab num2cell函数的用法num2cell 指令可将一数值阵列转换成异质阵列, 其使用语法为: C = num2cell(A, dim) 其中 dim 是代表 「被切割」的维度，若不指定，则将每个元素视为异质阵列里面的一个 1×1 矩阵例一：将数值阵列转换成异质阵列(cont)>> A = [1 2 3;4 5 6]    % 建立一个数值阵列 AA =