MATLAB小波函数小结(小编整理)

时间:2019-05-13 16:39:38下载本文作者:会员上传
简介:写写帮文库小编为你整理了多篇相关的《MATLAB小波函数小结》,但愿对你工作学习有帮助,当然你在写写帮文库还可以找到更多《MATLAB小波函数小结》。

第一篇:MATLAB小波函数小结

通用函数

函数 含义

Allnodes 计算树结点

appcoef

提取一维小波变换低频系数

appcoef2 提取二维小波分解低频系数

bestlevt

计算完整最佳小波包树

besttree

计算最佳(优)树

* biorfilt

双正交样条小波滤波器组

biorwavf

双正交样条小波滤波器

* centfrq

求小波中心频率

cgauwavf

Complex Gaussian小波

cmorwavf

coiflets小波滤波器

cwt

一维连续小波变换

dbaux

Daubechies小波滤波器计算

Dbwavf

Daubechies小波滤波器

dbwavf(W)W='dbN' N=1,2,3,...,50

ddencmp

获取默认值阈值(软或硬)熵标准

depo2ind

将深度-位置结点形式转化成索引结点形式

detcoef

提取一维小波变换高频系数

detcoef2

提取二维小波分解高频系数

disp

显示文本或矩阵

drawtree

画小波包分解树(GUI)

dtree

构造DTREE类

dwt

单尺度一维离散小波变换

dwt2

单尺度二维离散小波变换

dwtmode

离散小波变换拓展模式

* dyaddown

二元取样

* dyadup

二元插值

entrupd

更新小波包的熵值

Fbspwavf

B样条小波

Gauswavf

Gaussian小波

get

获取对象属性值

idwt

单尺度一维离散小波逆变换

idwt2

单尺度二维离散小波逆变换

ind2depo

将索引结点形式转化成深度—位置结点形式

* intwave

积分小波数

isnode

判断结点是否存在istnode

判断结点是否是终结点并返回排列值

iswt

一维逆SWT(Stationary Wavelet Transform)变换

iswt2

二维逆SWT变换

leaves

Determine terminal nodes

mexihat

墨西哥帽小波

Meyer

Meyer小波

meyeraux

Meyer小波辅助函数

morlet

Morlet小波

nodease

计算上溯结点

nodedesc

计算下溯结点(子结点)

nodejoin

重组结点

Nodepar

寻找父结点

nodesplt

分割(分解)结点

noleaves

Determine nonterminal nodes

ntnode

Number of terminal nodes

ntree

Constructor for the class NTREE

* orthfilt

正交小波滤波器组

plot

绘制向量或矩阵的图形

* qmf

镜像二次滤波器

rbiowavf

Reverse biorthogonal spline wavelet filters read

读取二进制数据

Readtree

读取小波包分解树

* scal2frq

Scale to frequency

set

shanwavf

Shannon wavelets

swt

一维SWT(Stationary Wavelet Transform)变换

swt2

二维SWT变换

symaux

Symlet wavelet filter computation.symwavf Symlets小波滤波器

thselect

信号消噪的阈值选择

thodes

References

treedpth

求树的深度

treeord

求树结构的叉数

upcoef

一维小波分解系数的直接重构

upcoef2

二维小波分解系数的直接重构

upwlev

单尺度一维小波分解的重构

upwlev2

单尺度二维小波分解的重构

Wavedec

单尺度一维小波分解

wavedec2

多尺度二维小波分解

wavedemo

小波工具箱函数demo

* wavefun

小波函数和尺度函数

* wavefun2

二维小波函数和尺度函数

wavemenu

小波工具箱函数menu图形界面调用函数

* wavemngr

小波管理函数

Waverec

多尺度一维小波重构

waverec2

多尺度二维小波重构

wbmpen

Penalized threshold for wavelet 1-D or 2-D de-noising wcodemat

对矩阵进行量化编码

wdcbm

Thresholds for wavelet 1-D using Birge-Massart strategy wdcbm2

Thresholds for wavelet 2-D using Birge-Massart strategy Wden

用小波进行一维信号的消噪或压缩

wdencmp

De-noising or compression using wavelets wentropy

计算小波包的熵

wextend

Extend a vector or a matrix

* wfilters

小波滤波器

wkeep

提取向量或矩阵中的一部分

* wmaxlev

计算小波分解的最大尺度

wnoise

产生含噪声的测试函数数据

wnoisest

估计一维小波的系数的标准偏差

wp2wtree

从小波包树中提取小波树

wpcoef

计算小波包系数

Wpcutree

剪切小波包分解树

wpdec

一维小波包的分解

wpdec2

二维小波包的分解

wpdencmp wpfun

wpjoin

wprcoef

wprec

wprec2

wpsplt

wpthcoef

wptree

wpviewcf wrcoef

wrcoef2

Wrev

Write

wtbo

wthcoef

wthcoef2

wthresh

wthrmngr wtreemgr

用小波包进行信号的消噪或压缩

小波包函数

重组小波包

小波包分解系数的重构

一维小波包分解的重构

二维小波包分解的重构

分割(分解)小波包

进行小波包分解系数的阈值处理

显示小波包树结构

Plot the colored wavelet packet coefficients.对一维小波系数进行单支重构

对二维小波系数进行单支重构

向量逆序

向缓冲区内存写进数据

Constructor for the class WTBO

一维信号的小波系数阈值处理

二维信号的小波系数阈值处理

进行软阈值或硬阈值处理

阈值设置管理

管理树结构

第二篇:总结MATLAB中涉及到的小波函数

收集和总结MATLAB中涉及到的小波函数

(1)plot函数:绘制向量或矩阵的图形 用法: plot(Y)plot(X1,Y1,...)plot(X1,Y1,LineSpec,...)plot(...,'PropertyName',PropertyValue,...)plot(axes_handle,...)(2)cwt函数:实现一维连续小波变换的函数。用法:

COEFS=cwt(S, SCALES, 'wname')COEFS=cwt(S, SCALES, 'wname', 'plot')COEFS=cwt(S, SCALES, 'wname', 'PLOTMODE')COEFS=cwt(S, SCALES, 'wname', 'PLOTMODE', XLIM)(3)dwt函数:单尺度一维离散小波变换 用法:

[cA,cD] = dwt(X,'wname')[cA,cD] = dwt(X,'wname','mode',MODE)[cA,cD] = dwt(X,Lo_D,Hi_D)[cA,cD] = dwt(X,Lo_D,Hi_D,'mode',MODE)(4)iswt函数:一维逆SWT(Stationary Wavelet Transform)变换 用法:

X = iswt(SWC,'wname')X = iswt(SWA,SWD,'wname')X = iswt(SWC,Lo_R,Hi_R)X = iswt(SWA,SWD,Lo_R,Hi_R)(5)upcoef函数:一维小波分解系数的直接重构用法:

Y = upcoef(O,X,'wname',N)Y = upcoef(O,X,'wname',N,L)Y = upcoef(O,X,Lo_R,Hi_R,N)Y = upcoef(O,X,Lo_R,Hi_R,N,L)Y = upcoef(O,X,'wname')Y = upcoef(O,X,Lo_R,Hi_R)

(6)upwlev函数:单尺度一维小波分解的重构 用法:

[NC,NL,cA] = upwlev(C,L,'wname')[NC,NL,cA] = upwlev(C,L,Lo_R,Hi_R)

(7)wavedec函数:单尺度一维小波分解

用法:

[C,L] = wavedec(X,N,'wname')[C,L] = wavedec(X,N,Lo_D,Hi_D)

(8)wavefun函数:小波函数和尺度函数

用法:

[PHI,PSI,XVAL] = wavefun('wname',ITER)[PHI1,PSI1,PHI2,PSI2,XVAL] = wavefun('wname',ITER)

[PSI,XVAL] = wavefun('wname',ITER)[...] = wavefun('wname',A,B)

(9)wpdec函数:一维小波包的分解

用法:

T = wpdec(X,N,'wname',E,P)T = wpdec(X,N,'wname')

(10)wpdencmp函数:用小波包进行信号的消噪或压缩

用法:

[XD,TREED,PERF0,PERFL2] =

wpdencmp(X,SORH,N,'wname',CRIT,PAR,KEEPAPP)

[XD,TREED,PERF0,PERFL2] =

wpdencmp(TREE,SORH,CRIT,PAR,KEEPAPP)

(11)wpfun函数:小波包函数 [用法:

WPWS,X] = wpfun('wname',NUM,PREC)

[WPWS,X] = wpfun('wname',NUM)

(12)wpjoin函数:重组小波包 用法:

T = wpjoin(T,N)[T,X] = wpjoin(T,N)T = wpjoin(T)

[T,X] = wpjoin(T)

(13)dbaux函数:Daubechies小波滤波器计算 用法:

W = dbaux(N,SUMW)W = dbaux(N)

(14)meyer函数:Meyer小波 用法:

[PHI,PSI,T] = meyer(LB,UB,N)[PHI,T] = meyer(LB,UB,N,'phi')[PSI,T] = meyer(LB,UB,N,'psi')(15)meyer函数:Meyer小波 用法:

[PHI,PSI,T] = meyer(LB,UB,N)[PHI,T] = meyer(LB,UB,N,'phi')[PSI,T] = meyer(LB,UB,N,'psi')

第三篇:收集和总结MATLAB中涉及到的小波函数

一、收集和总结MATLAB中涉及到的小波函数 1.cwt函数功能:实现一维连续小波变换的函数。11.upwlev函数功能:单尺度一维小波分解的重cwt函数语法格式:

COEFS=cwt(S, SCALES, 'wname')COEFS=cwt(S, SCALES, 'wname', 'plot')COEFS=cwt(S, SCALES, 'wname', 'PLOTMODE')2.dwt函数功能:单尺度一维离散小波变换 函数语法格式:

[cA,cD] = dwt(X,'wname')[cA,cD] = dwt(X,'wname','mode',MODE)[cA,cD] = dwt(X,Lo_D,Hi_D)3.meyer函数功能:Meyer小波 函数语法格式:

[PHI,PSI,T] = meyer(LB,UB,N)[PHI,T] = meyer(LB,UB,N,'phi')[PSI,T] = meyer(LB,UB,N,'psi')4.plot函数功能:绘制向量或矩阵的图形 函数语法格式: plot(Y)plot(X1,Y1,...)plot(X1,Y1,LineSpec,...)5.cgauwavf函数功能:Complex Gaussian小波 函数语法格式:

[PSI,X] = cgauwavf(LB,UB,N,P)6.iswt函数功能:一维逆SWT(Stationary Wavelet Transform)变换 函数语法格式:

X = iswt(SWC,'wname')X = iswt(SWA,SWD,'wname')X = iswt(SWC,Lo_R,Hi_R)7.mexihat函数功能:墨西哥帽小波

函数语法格式:

[PSI,X] = mexihat(LB,UB,N)

8.morlet函数功能:Morlet小波

函数语法格式:

[PSI,X] = morlet(LB,UB,N)

9.symwavf函数功能:Symlets小波滤波器 函数语法格式: F = symwavf(W)10.upcoef函数功能:一维小波分解系数的直接重构

函数语法格式:

Y = upcoef(O,X,'wname',N)Y = upcoef(O,X,'wname',N,L)Y = upcoef(O,X,Lo_R,Hi_R,N)Y = upcoef(O,X,Lo_R,Hi_R,N,L)Y = upcoef(O,X,'wname')Y = upcoef(O,X,Lo_R,Hi_R)

函数语法格式:

[NC,NL,cA] = upwlev(C,L,'wname')[NC,NL,cA] = upwlev(C,L,Lo_R,Hi_R)

12.wavedec函数功能:单尺度一维小波分解

函数语法格式:

[C,L] = wavedec(X,N,'wname')[C,L] = wavedec(X,N,Lo_D,Hi_D)

13. wavefun函数功能:小波函数和尺度函数 函数语法格式:

[PHI,PSI,XVAL] = wavefun('wname',ITER)

14.waverec函数功能:多尺度一维小波重构

函数语法格式:

X = waverec(C,L,'wname')

X = waverec(C,L,Lo_R,Hi_R)

15.wpcoef函数功能:计算小波包系数

函数语法格式: X = wpcoef(T,N)

X = wpcoef(T)

16.wpdec函数功能:一维小波包的分解

函数语法格式:

T = wpdec(X,N,'wname',E,P)T = wpdec(X,N,'wname')

17.wpfun函数功能:小波包函数 [函数语法格式:

WPWS,X] = wpfun('wname',NUM,PREC)

[WPWS,X] = wpfun('wname',NUM)

18.wprcoef函数功能:小波包分解系数的重构

函数语法格式:

X = wprcoef(T,N)

19.wprec函数功能:一维小波包分解的重构

函数语法格式:

X = wprec(T)

20.wrcoef函数功能:对一维小波系数进行单支重构

函数语法格式:

X = wrcoef('type',C,L,'wname',N)X = wrcoef('type',C,L,Lo_R,Hi_R,N)X = wrcoef('type',C,L,'wname')X = wrcoef('type',C,L,Lo_R,Hi_R)

第四篇:MatLab 知识小结

MatLab 知识小结

matlab常用到的永久变量。ans:计算结果的默认变量名。i j:基本虚数单位。

eps:系统的浮点(F10a9Bg个oht): inf: 无限大,例1/0 nan NaN:非数值(N航a nmnb谢)pi:圆周率n(n=3.1415926..)。realmax:系统所能表示的最大数值。realmin: 系统所能表示的最小数值,nargin: 函数的输入参数个数: nargout:函数的输出多数个数

①matlab的所有运算都定义在复数城上。对于方根问题运算只返回处于第一象限的解。

⑦matlab分别用左斜/和右来表示“左除和“右除”运算。对于标量运算而言,这两者的作用没有区别:但对于矩阵运算来说,二者将产生不同的结果。

多项式的表示方法和运算

p(x)=x^3-3x-5 可以表示为p=[1 0 –3 5],求x=5时的值用plotval(p,5)也可以求向量:a=[3 4 5],plotval(p,a)函数roots求多项式的根 roots(p)p=[1 0-3 5];r=roots(p)由根重组多项式poly(根)q=poly(r)

real(q)有时会产生虚根,这时用real抽取实根即可

conv(a,b)函数 多项式乘法(执行两个数组的卷积)a=[1 2 3 4];b=[1 4 9 16];c=conv(a,b)多项式的加减法,低阶的多项式必须用首零填补,使其与高阶多项式有同样的阶次

多项式除法 [q , r]=deconv(c , b)表示b/c q为商多项式,r为余数 多项式的导数 polyder(f)f=[ 2 4 5 6 2 1];s=polyder(f)

多项式的曲线拟合

x=[1 2 3 4 5];

y=[5.6 40 150 250 498.9];

p=polyfit(x,y,n)数据的n次多项式拟合 poly:矩阵的特征多项式、根集对应的多项式

x2=1:0.1:5;n取1时,即为最小二乘法

y2=polyval(p,x2);计算多项式的值(polyvalm计算矩阵多项式)plot(x,y,'*',x2,y2);grid on 最小二乘法 x=[1 2 3 4 5];

y=[5.6 40 150 250 498.9];plot(x,y,’*’),lsline

多项式插值(p158)

YI=interp1(x,y,XI,’method’)一维插值

(XI为插值点的自变量坐标向量,可以为数组或单个数。

method为选择插值算法的方法,包括:

linear(线性插值)cubic(立方插值)spline(三次样条插值)nearst(最近临插值)

例如:人口预测 year=1900:10:1900;

number=[78 91 105 „.每十年的人口数];

x=1900:1:2000;

y=interp1(year,number,x,’spline’);plot(year,numeber,’*’,x,y);grid on

一维博里叶变换插值使用函数interpft实现,计算含有周期函数值的矢量的傅里叶变换

然后使用更多的点进行傅里叶变换的逆变换,函数的使用格式如下:y=interpft(x,n)其中x是含有周期函数值的矢量,并为等距的点,n为返同等间距点的个数。

求解一元函数的最小值

y=fminbnd('humps',0.3,1)humps为一内置函数

求解多元函数的最小值

函数fminserch用于求多元函数的最小值。它可以指定一个开始的矢量,并非指定一个区间。此函数返回一个矢量为此多元函数局部最小函数值对应的自变量

纹理成图功能

由warp函数的纹理成图功能实现平面图像在空间三维曲面上的显示。将文件名为flowers.tif的图像分别投影到圆柱形和球星表面上 i=imread('flowers.tif');[x,y,z]=cylinder;

subplot(1,2,1),warp(x,y,z,i);[x,y,z]=sphere(50);subplot(1,2,2),warp(x,y,z,i);warp(x,y,z,i);

求函数的零点

求函数humps在[1,2]区间上的零点 fzero(‘humps’,[1,2]);

也可以给一个初始值 fzero(‘humps’,0.9);

对于多项式可直接由roots求其根 roots(‘4*x^3+……’);也可以用solve c=sym('c','real');x=sym('x','real');s=solve(x^3-x+c)

函数定积分

q=quadl(‘humps’,0,1)求humps函数在0 1区间上的定积分,也可以用quad语句

二重积分 首先计算内积分,然后借助内积分的中间结果再求出二重积分的值,类似于积分中的分步积分法。Result=dblquad(‘integrnd’,xin,xmax.,ymin,ymax)integrnd为被积函数的名称字符串

符号积分运算int(f)最精确的是符号积分法 计算s=∫12[∫01xydx]dy syms x y 中间为空格,不能为逗号 s=int(int(‘x^y’,’x’,0,1),’y’,1,2)引号可省略 vpa(s)显示s的值 内积分限为函数的二重积分 I=∫14[∫√y2(x2+y2)dx]dy 符号法I=vpa(int(int(‘x^2+y^2’,’x’,sqrt(y),2),’y’,1,4)

微分运算(diff)

微分是描述一个函数在一点处的斜率,是函数的微观性质、因此积分对函数的形状在小范围内的改变不敏感,而微分很敏感。—个函数的小的变化,容易产生相邻点的斜率的大的改变。由干微分这个固有的困难.所以尽可能避免数值微分.特别是对实验获得的数据进行微分。在这种情况,最好用最小二乘曲线拟合这种数据,然后对所得到的多项式进行微分;或用另一种方法对点数据进行三次样条拟合,然后寻找样条微分,但是,有时微分运算是不能避免的,在MATLAB中.用函数diff汁算一个矢量或者矩阵的微分(也可以理解为差分)。a=[1 2 3 3 3 7 8 9];b=diff(a)一次微分 bb=diff(a,2)二次微分 实际上diff(a)=[a(2)-a(1),a(3)-a(2),„„,a(n)-a(n-1)] 对于求矩阵的微分,即为求各列矢量的微分,从矢量的微分值可以判断矢量的单调性、是否等间距以及是否有重复的元素。符号微分运算(diff)syms x t a f =cos(a*x)df =diff(f)由findsym的规则,隐式的指定对x进行微分

dfa=diff(f,'a')指定对变量a进行微分 dfa=diff(f,'a',3)三次微分

diff函数不仅作用在标量上,还可以在矩阵上,运算规则就是按矩阵的元素分别进行微分 syms a x A=[cos(a*x),sin(a*x),-sin(a*x),cos(a*x)];dA=diff(A)微分方程dsolve

在matlab中,符号表达式中包含字母D用来表示微分运算,D2,D3分别对应第二,第三阶导数,D2y表示d2y/dt2 把t缺省了

y=dsolve(‘Dy=f(y)’)单个方程,单个输出

[u,v]=dsolve(‘Du=f(u,v)’,’Dv=g(u,v)’)2个方程,2个输出

s=dsolve(‘Dx=f(x,y,z)’,’Dy=g(x,y,z)’,’Dz=k(x,y,z)’)

s.x s.y s.z 3个方程,架构数组

dsolve('Dx=-a*x')结果:C1*exp(-a*t)没给定初值,所以结果中含参变量 x=dsolve('Dx=-a*x','x(0)=1','s')结果exp(-a*s)给定了初值,独立变量设为s

计算多元函数的梯度

fx=gradient(f)f是一个矢量返回f的一维数值梯度,fx对应于x方向的微分。

[x,y]=meshgrid(-2:.2:2,-2:.2:2);z=x.*exp(-x.^2-y.^2);[px,py]=gradient(z,.2,.2);contour(z),hold on 画等值线 quiver(px,py)

matlab字符串运算 利用sym命令创建表达式

f=sym(‘cos(x)+sin(x)’)或 syms x , f=cos(x)+sin(x)diff(f)求其导数(也

f=diff(‘cos(x)+cos(y)’)

当字符表达式中含有多于一个的变量时,只有—个变量是独立变量。如果不告诉matlab哪一个变量是独立变量,则可以通过findsym命令询问 利用findsym命令查询独立变量 f=sym('sin(a*x)+b')

findsym(f,1)给出独立变量(一个变量,如果为2则给出2个变量)findsym(f)给出所有变量

符号表达式的化简和替换

collect函数 collect(f,v)表示将f表示为关于符号变量v的多项式形式,即关于v合并同类项,v缺省,则用findsym确定的缺省变量 syms x y

f=x^2*y+y*x-x^2-2*x+1 collect(f)得到(-1+y)*x^2+(y-2)*x+1 collect(f,y)

(x+x^2)*y+1-x^2-2*x

expand函数 expand(f)将f展开,写成和的形式 syms x

expand((x-1)^3)得

x^3-3*x^2+3*x-1

horner函数 horner(f)将f写成镶嵌套形式 syms x

horner(x^3-6*x^2)得

(-6+x)*x^2

factor函数 factor(f)将f转换成低阶有理多项式的乘积 syms x

f=x^3-6*x^2+11*x-6

factor(f)得到(x-1)*(x-2)*(x-3)simplify(f)函数 综合化简 simple(f)函数的最简形式 syms x

f=2*sin(x^2)+cos(3*x)

simple(f)如果不想看到中间过程,可z=simple(f)有时使用两次simple命令可以得到最简式

如果想知道哪个简化命令得到最后结果,可以加一个参数how [z,how]=simple(f)

符号表达式的替换 subs(f,new,old)f='a*x^2+b*x+c'

subs(f,'t','x')得到a*(t)^2+b*(t)+c subs是一个符号函数,返回一个符号变量

subexpr函数 有时matlab返回的符号表达式难以理解,用subexpr函数,可以将表达式中重复出现的子式用一个符号表示,从而简化表达形式 c=sym('c','real');x=sym('x','real');s=solve(x^3-x+c)a=subexpr(s)

得到

sigma

=

-108*c+12*(-12+81*c^2)^(1/2)a =

[ 1/6*sigma^(1/3)+2/sigma^(1/3)] [-1/12*sigma^(1/3)-1/sigma^(1/3)+1/2*i*3^(1/2)*(1/6*sigma^(1/3)-2/sigma^(1/3))] [-1/12*sigma^(1/3)-1/sigma^(1/3)-1/2*i*3^(1/2)*(1/6*sigma^(1/3)-2/sigma^(1/3))]

pretty函数有时也能起到同样的作用。Pretty(f)显示函数的习惯书写形式

线性方程组的求解

求解线性方程组,用反斜杠 a=hilb(3)b=[1 2 3]' ab

矩阵的特征值和特征向量

用eig(v,d)函数,[v,d]=eig(A);其中d将返回特征值,v返回相应的特征向量,缺省第二个参数将只返回特征值 syms a b c real A=[a b c;b c a;c a b];[v,d]=eig(A);

为了观察更清楚,使用以前学过的替换函数,这里不用默认的sigma,而改用M,显式的代替繁琐的表达子式 vv=subexpr(v);vs=subs(vv,'m','sigma')运行结果为

vs = [ 1, 1, 1] [-(c+(m)-a)/(c-b),-(c-(m)-a)/(c-b), 1] [-(a-(m)-b)/(c-b),-(a+(m)-b)/(c-b), 1]

再用m替换d中的表达子式 dd=subexpr(d);ds=subs(dd,’m’,’sigma’)运行结果为ds =

[(m), 0, 0] [ 0,-(m), 0] [ 0, 0, c+a+b] note 求特征值也可用以下命令

f=poly(A)poly函数 用来求A的特征多项式

d=solve(f)solve(f)函数用来求多项式的解

svd()函数 求矩阵的奇异值分解,将矩阵分解为两个正交矩阵和对角矩阵的乘积 a=sym(hilb(2))[u,s,v]=svd(a)

代数方程和方程组

代数方程的求解可用solve(f)命令,如果f不含=,matlab将给表达式置零。方程的未知量在默认的情况下由findsym决定或显式指出 syms a b c x

solve(a*x^2+b*x+c)以x为默认变量

solve(a*x^2+b*x+c,a)指定对a为变量

求含有等号的方程的解(一定要加单引号)

f=solve(‘cos(x)=sin(x)’)

x=solve('exp(x)=tan(x)')如果不能求得符号解,就计算可变精度解。求解方程组与单方程类似 解一个三元一次方程

v=solve('a*u^2+v^2','u-v=1','a^2-5*a+6')结果为v =

a: [4x1 sym] u: [4x1 sym] v: [4x1 sym]

一些常用的符号运算 极限运算limit

limit(f)求x到0的极限

limit(f,x,a)或limit(f,a)求x到a的极限

limit(f,a,’left’)limit(f,a,’right’)求x到a的左极限和右极限 limit(f,inf)求x趋于无穷的极限 符号求和symsum(s)

symsum(s)以默认的findsym决定的变量求和

symsum(s,v)以s中指定的变量v求和

symsum(s,a,b)symsum(s,v,a,b)从a到b的有限项求和

syms k n

symsum(k)从0到k求和

symsum(k,0,n-1)从0到n-1求和 symsum(1/k^2,1,inf)无限项求和 泰勒级数taylor(f)

taylor(f)表示求f的5阶talor展开,可以增加参数指定展开的阶数(默认式5),也可以对于多元函数指定展开的变量,还可以指定在哪个点展开 syms x t taylor(exp(-x))

taylor(log(x),6,1)在1点的6阶taylor展开

taylor(x^t,3,t)对t的3阶taylor展开 积分变换

fourier变换和逆变换fourier(f)fourier分析可以将信号转换为不同频率的正弦曲线。可对离散数据进行分析,也可对连续时间系统进行分析,特别在信号和图形处理领域。离散变换(DFT)作用于有限数据的采集,最有效的是快速fourier变换(FFT)F=fourier(f)独立变量x,返回关于参数w的函数

F=fourier(f,v)返回函数F关于符号对象v的函数

F=fourier(f,u,v)对关于u的函数f进行变换,而不是缺省的w,返回函数F是关于v的函数 syms t v w x fourier(1/t)

fourier(exp(-t)*sym('Heaviside(t)'),v)fourier(diff(sym('F(x)')),x,w)Fourier逆变换

f=ifourier(F)缺省独立变量w,返回关于x的函数对w进行积分 f=ifourier(F,v)返回函数f是关于符号对象v的函数,而不是缺省的x f=ifourier(F,u,v)是关于u的函数f进行变换,而不是缺省的x,返回函数f是关于v的函数

Laplace变换和逆变换laplace(f)应用于连续系统(微分方程)中,可以用来求解微分方程的初值问题 laplace(F)缺省独立变量t,缺省返回关于s的函数L

laplace(F,t)返回关于t的函数L,而不是缺省的s

laplace(F,w,z)对函数F的自变量w积分,返回关于z的函数L 逆变换

F=ilaplace(L)缺省独立变量s,返回关于t的函数F F=ilaplace(L,y)返回关于y的函数F,而不是缺省的t F=ilaplace(L,y,x)对函数L的自变量y积分,返回关于x的函数F Z-变换和逆变换ztrans(f)标量符号f的Z-变换

F=ztrans(f)缺省独立变量n,返回关于z的函数

F=ztrans(f,w)返回关于符号变量w的函数F,而不是缺省的z F=ztrans(f,k,w)关于k的符号变量作Z-变换返回关于符号变量w的函数 逆变换iztrans(F)f=iztrans(F)或(F,k)或(F,w,k)

符号绘图函数

符号函数简易绘图函数ezplot(f)f可以包含单个符号变量x的字符串或表达式,默认画图区间(-2pi,2pi),如果f包含x和y,画出的图像是f(x,y)=0的图像,缺省区间是-2pi

syms x t ezplot('t*cos(t)','t*sin(t)',[0,4*pi])绘制符号图像函数fplot(fun,lims,tol,’linespec’,n)其中lims=[xmin,xmax]或[xmin,xmax,ymin,ymax] tol为指定相对误差,默认0.001 ‘linespec’指定绘图的线型 n指定最少以n+1个点绘图

[x,y]=fplot(fun,lims,…)只返回用来绘图的点,并不绘图,可以自己调用plot(x,y)来绘制图形。syms x subplot(2,2,1),fplot('humps',[0,1])f='abs(exp(x*(0:9))*ones(10,1))' subplot(2,2,2),fplot(f,[0,2*pi])subplot(2,2,3),fplot('sin(1./x)',[0.01,0.1],1e-3)matlab绘图 二维图形的绘制

plot 在(x,y)坐标下绘制二维图像 支持多个x-y二元结构

plot3 在(x,y,z)坐标下绘制三维图形 loglog 在(x,y)对数坐标下绘制二维图形

semilogx 在x为对数坐标,y为线性坐标的二维坐标中绘图

semilogy 在x为线性坐标,y为对数坐标的二维坐标中绘图

plotyy 在有两个y轴的坐标下绘图

plot用法

plot(x,y,'--rs','linewidth',2,'markeredgecolor','k',...'markerfacecolor','g','markersize',10)plotyy用法

plotyy(x1,y1,x2,y2)以x1为标准,左轴为y轴绘制y1向量,x2为基准,右轴为y轴,绘制y2向量

plotyy(x1,y1,x2,y2,fun)用字符串fun指

定的绘

(plot ,semilogx,semilogy,loglog,stem)plotyy((x1,y1,x2,y2,fun1,fun2)t=0:pi/20:2*pi;y=exp(sin(t));

plotyy(t,y,t,y,'plot','stem')stem为二维杆图

[ax,h1,h2]=plotyy(„)返回左右两y轴的句柄(分别为ax(1)ax(2),以及在两坐标轴中生成的图形对象的句柄,分别为h1 h2 t=0:900;A=1000;a=0.005;b=0.005;z2=cos(b*t);z1=A*exp(-a*t);

[haxes,hline1,hline2]=plotyy(t,z1,t,z2,'semilogy','plot');axes(haxes(1))

ylabel('semilog plot')对数坐标 axes(haxes(2))ylabel('linear plot')set(hline2,'linestyle','--')其他二维图形绘图指令 bar(x,y)二维条形图 hist(y,n)直方图

histfit(y,n)带拟和线的直方图,n为

直方的个数 stem(x,y)火柴杆图 comet(x,y)彗星状轨迹图 compass(x,y)罗盘图 errorbar(x,y,l,u)误差限图 feather(x,y)羽毛状图

fill(x,y,’r’)二维填充函数 以红色填充

pie(x)饼图

polar(t,r)极坐标图 r为幅值向量,t为角度向量 t=0:0.1:8*pi;r=cos(3*t/2)+1/2;

polar(t,r),xlabel('polar 指令')quiver(x,y)磁力线图 stairs(x,y)阶梯图 loglog(x,y)对数图

semilogx semilogy 半对数图

matlab三维作图 plot3(x,y,z)三维线条图 t=0:pi/50:15*pi;

plot3(sin(t),cos(t),t,'r*')与plot相似 v=axis 返回各个轴的范围

text(0,0,0,'origin')在某个坐标点加入文字

plot3 增加维数可以一次画多个图,使所个二维图形眼一个轴排列

三维网线图的绘制 mesh(x,y,z)网格图

mesh(x,y,z,c)四维作图,(x,y,z)代表空间三维,c代表颜色维 mesh(…,’property name’,property

value,…)设置曲面各属性的值

[x,y,z]=sphere(12);

mesh(x,y,z),hidden off 曲面设置为透明

meshc(x,y,z)画网格图和基本的等值线图

meshz(x,y,z)画包含零平面的网格图 waterfall(x,y,z)与mesh一样,只是在效果上它的网格线只在x轴一个方向出现,呈瀑布状水线

两个变量的标量指令meshgrid(x)或meshgrid(x,y)(p179)

将两个一维向量生成两个二维向量,以便进行z=f(x,y)运算,算出z的所有值,z为x y的标量指令 [X,Y]=meshgrid(x)meshgrid(x,x)的简略式

[X,Y]=meshgrid(x,y)[X,Y,Z]=meshgrid(x,y,z)用于三维图形的绘制

[x,y]=meshgrid([-2:0.1:2]);z=x.*exp(-x.^2-y.^2);plot3(x,y,z)surf(x,y,z,c)着色表面图 surf(x,y,z)隐含着c=z surf(z)隐含着x,y的值为surf指令根据z的尺寸自动生成

surfc 画出具有基本等值线的曲面图 surfl 画出一个具有亮度的曲面图 shading flat 网线图的某整条线段或曲面图的某个贴片都着一种颜色 shading interp 某一线段或贴片上各点的颜色由线或片的顶端颜色经线性插值而得

曲面图不能设成网格图那样透明,但需要时,可以在孔洞处将数据设成nun

等高线的绘制

在二维空间绘制等高线contour contour(x,y,z,n)绘制n条等值线(n可省略)

contour(x,y,z,v)在向量v所指定的高度上绘制等高线(可省)

c=contour(x,y,z)计算等值线的高度值

c=contourc(x,y,z,n)计算n条等高线的x-y坐标数据

c=contourc(x,y,z,v)计算向量v所指定的等高线的x-y坐标数据 clabel(c)给c阵所表示的等高线加注高度标识

clabel(c,v)给向量v所指定的等高线加注高度标识

clabel(c,’manual’)借助鼠标给点中的等高线加注高度标识

三维空间绘制等高线contour3(x,y,z)[x,y,z]=peaks(30);contour3(x,y,z,16,'g')二元函数的伪彩图pcolor(x,y,z)是指令surf的二维等效指令,代表伪彩色,可与contour单色等值线结合画彩色等值线图 [x,y,z]=peaks(30);

pcolor(x,y,z);伪彩色

shading interp 颜色插值,使颜色平均渐变

hold on,contour(x,y,z,20,'k')...画等值线

colorbar('horiz')水平颜色标尺 c=contour(x,y,z,8);clabel(c)标注等高线

矢量场图(速度图)quiver

用于描述函数z=f(x,y)在点(x,y)的梯度大小和方向

[X,Y]=meshgrid(x,y)X,Y为Z阵元素的坐标矩阵

[U,V]=gradient(Z,dx,dy)U,V分别为Z对x对y的导数,dx dy是x y方向上的计算步长

quiver(X,Y,U,V,s,’linespec’,’filled’)U,V为必选项,决定矢量场图中各矢量的大小和方向,s为指定所画箭头的大小,缺省时取1,linespec为字符串,指定合法的线形和彩色,filled用于填充定义的绘图标识符

[x,y]=meshgrid(-2:.2:2,-1:.15:1);z=x.*exp(-y.^2);

[px,py]=gradient(z,.2,.15);contour(x,y,z);

hold on,quiver(x,y,px,py),axis image 多边形的填色fill(x,y,c)

c定义颜色字符串,可以是’r’,’b’等,也可以用RGB三色表示[r,g,b]值为0-1

图形的四维表现

第五篇:matlab图像处理小结

1.function [center, r] = solve_circle(pt1, pt2, pt3)

2.%Effect: solve the circle which across points 'pt1', 'pt2' and 'pt3' 3.%Inputs:

4.%pt1, pt2, pt3: [x, y]

5.%center: the circle center [x0;y0] 6.%r: the radius of the circle 7.%Author: Su dongcai at 2012/1/2 8.A = zeros(2, 2);B = zeros(2, 1);9.[A(1, :), B(1)] = circle2line(pt1, pt2);10.[A(2, :), B(2)] = circle2line(pt2, pt3);11.center = AB;

12.r = norm(pt1'(y2^2 + y2^2)18.%(a-x2)^2 +(b-y2)^2 = r^2 | 19.%Inputs:

20.%pt1, pt2: [x1, y1], [x2, y2] 21.%Outputs:

22.%A: 2[x1-x2, y1-y2]

23.%B:(x1^2 + y1^2)pt2);

26.B = norm(pt1)^2-norm(pt2)^2;

close all;clear;clc;>> i=imread('rice.png');%>> imshow(i);>> background=imopen(i,strel('disk',15));>> i2=imsubtract(i,background);%>> figure,imshow(i2);>> i3=imadjust(i2,stretchlim(i2),[0 1]);%>> figure,imshow(i3);>> level=graythresh(i3);>> bw=im2bw(i3,level);%>> figure,imshow(bw);>> [labeled,numobjects]=bwlabel(bw,4);graindata=regionprops(labeled,'all');

close all;clear;clc;i=imread('rice.png');background=imopen(i,strel('disk',15));i2=imsubtract(i,background);i3=imadjust(i2,stretchlim(i2),[0 1]);level=graythresh(i3);bw=im2bw(i3,level);[labeled,numobjects]=bwlabel(bw,4);data=regionprops(labeled,'all');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %2006.6.2 close all;clear;clc;>> i=imread('r.jpg');%>> figure,imshow(i);>> imgray=rgb2gray(i);>> figure,imshow(imgray)>> background=imopen(imgray,strel('disk',15));>> i2=imsubtract(imgray,background);%>> figure,imshow(i2);>> i3=imadjust(i2,stretchlim(i2),[0 1]);%>> figure,imshow(i3);>> level=graythresh(i3);>> bw=im2bw(i3,level);%>> figure,imshow(bw);>> imnobord=imclearborder(bw,4);%>> figure,imshow(imnobord);>> [labeled,numobjects]=bwlabel(bw,4);>> rgb_label=label2rgb(labeled,@spring,'c','shuffle');>> figure,imshow(rgb_label);>> graindata=regionprops(labeled,'all');hold on;for k=1:numobjects lab=sprintf('%d',k);text(graindata(k).Centroid(1),graindata(k).Centroid(2),lab,'Color','k');end hold off;%剔除碎米粒

>> idxdown=find([graindata.Area]<150);%剔除碎米粒 little=ismember(labeled,idxdown);figure,imshow(little);

[lab_little,num_little]=bwlabel(little,4);rgb_little=label2rgb(lab_little,@spring,'c','shuffle');figure,imshow(rgb_little);

little_data=regionprops(lab_little,'all');hold on;for k=1:num_little lab=sprintf('%d',k);text(little_data(k).Centroid(1),little_data(k).Centroid(2),lab,'Color','k');end hold off;%>> graindata(idxdown,:)=[];%剔除碎米粒 %剔除连接米粒

>> idxup=find([graindata.Area]>250);%剔除连接米粒 big=ismember(labeled,idxup);figure,imshow(big);

[lab_big,num_big]=bwlabel(big,4);rgb_big=label2rgb(lab_big,@spring,'c','shuffle');figure,imshow(rgb_big);

big_data=regionprops(lab_big,'all');hold on;for k=1:num_big lab=sprintf('%d',k);text(big_data(k).Centroid(1),big_data(k).Centroid(2),lab,'Color','k');end hold off;%>> graindata(numup,:)=[];%剔除连接米粒 %获取完整米粒

idxsuit=find([graindata.Area]>=150&[graindata.Area]<=250);suit=ismember(labeled,idxsuit);figure,imshow(suit);%获取完整米粒 [lab_suit,num_suit]=bwlabel(suit,4);suit_data=regionprops(lab_suit,'all');hold on;for k=1:num_suit signature=sprintf('%d',k);text(suit_data(k).Centroid(1),suit_data(k).Centroid(2),signature,'Color','r');end hold off;%获取完整米粒 whos graindata whos little_data whos big_data whos suit_data

>> graindata >> mean([graindata.Area])>> mean([graindata.Eccentricity])>> mean([graindata.MajorAxisLength])>> mean([graindata.MinorAxisLength])>> mean([graindata.EquivDiameter])>> figure,hist([graindata.Area],20);>> figure,hist([graindata.Eccentricity],20);>> figure,hist([graindata.MajorAxisLength],20);>> figure,hist([graindata.MinorAxisLength],20);>> figure,hist([graindata.EquivDiameter],20);

data=[graindata.Area] data=[graindata.Centroid] data=[graindata.BoundingBox] data=[graindata.SubarrayIdx] data=[graindata.MajorAxisLength] data=[graindata.MinorAxisLength] data=[graindata.Eccentricity] data=[graindata.Orientation] data=[graindata.ConvexHull] data=[graindata.ConvexImage] data=[graindata.ConvexArea] data=[graindata.Image] data=[graindata.FilledImage] data=[graindata.FilledArea] data=[graindata.EulerNumber] data=[graindata.Extrema] data=[graindata.EquivDiameter] data=[graindata.Solidity] data=[graindata.Extent] data=[graindata.PixelIdxList] data=[graindata.PixelList]

Area 计算各个连通区域中的象素总数 BoundingBox 包含相应区域的最小矩形 Centroid 给出每个区域的质心

MajorAxisLength 与区域具有相同标准二阶中心矩(又叫标准差)的椭圆的长轴长度 MinorAxisLength 与区域具有相同标准二阶中心矩的椭圆的短轴长度 Eccentricity 与区域具有相同标准二阶中心矩的椭圆的离心率

Orientation 与区域具有相同标准二阶中心矩的椭圆的长轴与x轴的交角 Image 二值图像,与某区域具有相同大小的逻辑矩阵。

FilledImage 与上相同,唯一区别是这是个做了填充的逻辑矩阵!本例中和上面的没有区别,只有 区域有空洞时才有明显差别。

FilledArea 是标量,填充区域图像中的 on 像素个数

ConvexHull 是p行2列的矩阵,包含某区域的最小凸多边形 ConvexImage 二值图像,用来画出上述的区域最小凸多边形 ConvexArea 是标量,填充区域凸多边形图像中的 on 像素个数 EulerNumber 等于图像中目标个数减去这些目标中空洞的个数 Extrema 8行2列矩阵,八方向区域极值点

EquivDiameter 是标量,等价直径:与区域具有相同面积的圆的直径.计算公式为:sqrt(4*Area/pi)

Solidity 是标量,同时在区域和其最小凸多边形中的像素比例。计算公式为: Area/ConvexArea,这也是个仿射特征,实际上反映出区域的固靠性程度。

Extent 是标量,同时在区域和其最小边界矩形中的像素比例。计算公式为:Area除以边界矩 形面积,这也是个仿射特征,实际上反映出区域的扩展范围程度。

PixelIdxList p元向量,存储区域像素的索引下标

PixelList p行ndims(L)列矩阵,存储上述索引对应的像素坐标 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 基于特定原则的区域选择

当你要基于特定准则条件选择某个区域时,将函数 ismember 和 regionprops 联合使用是很有用处的。例如:创建一个只包含面积大于80的二值图像,用以下命令

idx = find([stats.Area] > 80);BW2 = ismember(L,idx);regionprops函数的扩展思路

在regionprops函数的基础上,你可以使用它提供的基本数据来扩展它的功能,比如我就将区域的曲率数据和骨架数据作为它的另外属性值来开发,从而希望它能用来做更细致的特征提取。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %2006.6.2 P221图像粒度测定(雪花)>> i=imread('snowflakes.png');>> figure,imshow(i);>> %(2)>> clahei=adapthisteq(i,'numtiles',[10 10]);>> clahei=imadjust(clahei);>> imshow(clahei);>> gi=imadjust(im2double(i),[],[0 1]);>> figure,imshow(gi),title('adjusted grayscale image');>> %(3)>> se=strel('disk',10);>> topi=imtophat(gi,se);>> figure,imshow(topi),title('top-hat image');>> %(4)>> for counter=0:22 remain=imopen(clahei,strel('disk',counter));intensity_area(counter+1)=sum(remain(:));end >> figure,plot(intensity_area,'m-*'),grid on;>> title('sum of opening(pixels)');>> title('sum of opening values in opened image as a function of radius');>> xlabel('radius of opening(pixels)');>> ylabel('pixel value sum of opened objects(intensity)');>> >> >> >> for counter=0:20 remain=imopen(topi,strel('disk',counter));surfarea(counter+1)=sum(remain(:));end >> figure,plot(surfarea,'m-*'),grid on;>> set(gca,'xtick',[0 2 4 6 8 10 12 14 16 18 20]);>> title('surface area of opened objects as a function of radius');>> xlabel('radius of opening(pixels)');>> ylabel('surface area of opened objects(pixels)');>> %(5)>> intensity_area_prime=diff(intensity_area);>> figure,plot(intensity_area_prime,'m-*'),grid on;>> title('Granulometry(size distrubution)of snowflakes');>> set(gca,'xtick',[0 2 4 6 8 10 12 14 16 18 20 22]);>> xlabel('radius of snowflakes(pixels)');>> ylabel('sum of pixel values in snowflakes as a function of radius');>> derivsurfarea=diff(surfarea);>> figure,plot(derivsurfarea,'m-*'),grid on;>> title('granulometry(size distribution)of stars');>> xlabel('radius of stars(pixels)');>> ylabel('loss of pixels between two successive openings');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %2006.6.2 花椒检测 clc;clear;close all;i=imread('gj.jpg');imshow(i);icanny=edge(i,'canny');imshow(icanny);se90=strel('line',2,90);se0=strel('line',2,0);bwsdil=imdilate(icanny,[se90 se0]);figure,imshow(bwsdil),title('dilated');ifill=imfill(bwsdil,'holes');figure,imshow(ifill);

%bwero=imerode(bwsdil,[se90 se0]);%figure,imshow(bwero);%i2fill=imfill(bwero,'holes');%figure,imshow(bwero);%imshow(i2fill);

%bwnobord=imclearborder(bwsdil,4);%figure,imshow(bwnobord);bwnobord=imclearborder(ifill,4);figure,imshow(bwnobord);se=strel('disk',5);bwc=imclose(bwnobord,se);bwco=imopen(bwnobord,se);figure,imshow(bwc);figure,imshow(bwco);%mask=bwsdil&bwco;%figure,imshow(mask);clc [labeled,numobjects]=bwlabel(bwco);numobjects

jdata=regionprops(labeled,'all');%jdata

jarea=[jdata.Area];mean(jarea)max(jarea)min(jarea)hist(jarea,255)jdata.Eccentricity %std([jdata.Eccentricity])/(Mean([jdata.Eccentricity])jstd=std([jdata.Eccentricity])jmean=Mean([jdata.Eccentricity])jcv=jstd/jmean

>> std([jdata.Area])/ mean([jdata.Area])%面积的变异系数

>> std([jdata.Eccentricity])/ mean([jdata.Eccentricity])%椭圆的变异系数 >> std([jdata.MajorAxisLength])/ mean([jdata.MajorAxisLength])>> std([jdata.MinorAxisLength])/ mean([jdata.MinorAxisLength])>> std([jdata.EquivDiameter])/ mean([jdata.EquivDiameter])

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %2006.06.06 rice.png close all;clear;clc >> i=imread('rice.png');imshow(i);background=imopen(i,strel('disk',15));figure,imshow(background);i2=imsubtract(i,background);figure,imshow(i2);i3=imadjust(i2,stretchlim(i2),[0 1]);figure,imshow(i3);level=graythresh(i3);bw=im2bw(i3,level);figure,imshow(bw);imnobord=imclearborder(bw);[label,numobjects]=bwlabel(imnobord,4);numobjects rgb_label=label2rgb(label,@spring,'c','shuffle');figure,imshow(rgb_label);graindata=regionprops(label,'all');graindata

>> numdown=find([graindata.Area]<150);>> graindata(numdown,:)=[];>> numup=find([graindata.Area]>250);>> graindata(numup,:)=[];>> graindata

>> std([graindata.Area])/ mean([graindata.Area])%面积的变异系数

>> std([graindata.Eccentricity])/ mean([graindata.Eccentricity])%椭圆的变异系数

>> std([graindata.MajorAxisLength])/ mean([graindata.MajorAxisLength])>> std([graindata.MinorAxisLength])/ mean([graindata.MinorAxisLength])>> std([graindata.EquivDiameter])/ mean([graindata.EquivDiameter])%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %2006.06.06 rice的垩白度检测 >> clear;close all;clc;>> rgb=imread('r.jpg');>> close all;>> imshow(rgb);>> i=rgb2gray(rgb);>> j=medfilt2(i,[5 5]);>> figure,imshow(i);>> figure,imshow(j);>> imhist(j,256);>> t=0.3;>> v=imadjust(j,[t 1],[],1);>> imhist(v,256);>> t_c=0.6;>> bw_v=im2bw(v,0.01);>> chalk=imadjust(v,[t_c 1],[],1);>> bw_chalk=im2bw(chalk,0.01);>> figure,imshow(v);>> figure,imshow(bw_v);>> figure,imshow(chalk);>> figure,imshow(bw_chalk);>> degree_chalkness=bwarea(bw_chalk)/bwarea(bw_v)*100 >> bw=im2bw(j,t);>> figure,imshow(bw);>> se=(ones(3,3));>> bw1=imerode(bw,se);%两次腐蚀 >> figure,imshow(bw1);>> bw2=imerode(bw1,se);>> figure,imshow(bw2);

>> [l,num]=bwlabel(bw2);%标记腐蚀后的大米图像 >> t_chalk=100;%设置垩白面积的下限 >> compare=(l)&(chalk>t_chalk);%>> compare=(bw2)&(bw_chalk>t_chalk);>> [r,c]=find(compare);%标记垩白米粒的位置 >> result=bwselect(l,c,r);%显示只含有垩白米粒的图像 >> figure,imshow(result);

>> [l_chalk,num_chalk]=bwlabel(result);%标记垩白米粒图像,便于计数 >> rate_chalky_grains=num_chalk/num*100;>> rate_chalky_grains

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %2006.6.17 bwmorph函数 >> help bwmorph BWMORPH Perform morphological operations on binary image.BW2 = BWMORPH(BW1,OPERATION)applies a specific morphological operation to the binary image BW1.BW2 = BWMORPH(BW1,OPERATION,N)applies the operation N times.N can be Inf, in which case the operation is repeated until the image no longer changes.OPERATION is a string that can have one of these values: 'bothat' Subtract the input image from its closing 'bridge' Bridge previously unconnected pixels 'clean' Remove isolated pixels(1's surrounded by 0's)'close' Perform binary closure(dilation followed by erosion)'diag' Diagonal fill to eliminate 8-connectivity of background 'dilate' Perform dilation using the structuring element ones(3)'erode' Perform erosion using the structuring element ones(3)'fill' Fill isolated interior pixels(0's surrounded by 1's)'hbreak' Remove H-connected pixels 'majority' Set a pixel to 1 if five or more pixels in its 3-by-3 neighborhood are 1's 'open' Perform binary opening(erosion followed by dilation)'remove' Set a pixel to 0 if its 4-connected neighbors are all 1's, thus leaving only boundary pixels 'shrink' With N = Inf, shrink objects to points;shrink objects with holes to connected rings 'skel' With N = Inf, remove pixels on the boundaries of objects without allowing objects to break apart 'spur' Remove end points of lines without removing small objects completely.'thicken' With N = Inf, thicken objects by adding pixels to the exterior of objects without connected previously unconnected objects 'thin' With N = Inf, remove pixels so that an object without holes shrinks to a minimally connected stroke, and an object with holes shrinks to a ring halfway between the hold and outer boundary 'tophat' Subtract the opening from the input image

Class Support-------------The input image BW1 can be numeric or logical.It must be 2-D, real and nonsparse.The output image BW2 is logical.Examples--------BW1 = imread('circles.png');imview(BW1)BW2 = bwmorph(BW1,'remove');BW3 = bwmorph(BW1,'skel',Inf);imview(BW2)imview(BW3)

See also erode, dilate, bweuler, bwperim.Reference page in Help browser doc bwmorph

BW1 = imread('circles.png');figure,imshow(BW1)BW2 = bwmorph(BW1,'erode');figure,imshow(BW2)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %边界提取 b=bwmorph(bw,'remove');b=bwperim(bw,8);%又叫边界象素测定 b=edge(bw,'canny');%又叫边界提取 %去除孤立象素点

nosinglepixel=bwmorph(bw,'clean');%去除小面积物体

nosmall=bwareaopen(bw,CNN);%阈值处理再取反

bw=~im2bw(i,graythresh(i));

%开运算(消除小物体)与闭运算(填充物体内细小空洞)se=strel('disk',6);iopen=imopen(bw,se);iclose=imclose(bw,se);%腐蚀与膨胀联合操作 %(1)创建结构元素 se=strel('rectangle',[40 30]);%(2)使用结构元素腐蚀图像 bw1=imread('circbw.tif');bw2=imerode(bw1,se);imshow(bw2);%(3)逆操作,回复矩形原来大小 bw3=imdilate(bw2,se);figure,imshow(bw3);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %2006.6.18花椒子

%直接对灰度图进行canny运算 >> i=imread('nut.bmp');>> figure,imshow(i);>> ig=rgb2gray(i);>> figure,imshow(ig);%igcanny=edge(ig,'canny');%igcfill=imfill(igcanny,'hole');igcanny_thresh=edge(ig,'canny',(graythresh(ig)*.1));igcfill=imfill(igcanny_thresh,'hole');>> figure,imshow(igcfill);

%先对灰度图滤波,再进行canny运算

>> imed=medfilt2(ig);%中值滤波后对图像边界有一定的损伤!!>> imedcanny=edge(imed,'canny');>> imedfill=imfill(imedcanny,'hole');>> figure,imshow(imedfill);>> nosmall=bwareaopen(imedfill,150);>> figure,imshow(nosmall);

%注意:若对灰度图像先拉氏锐化,在canny提取边界,效果不大好!!%结论:无需拉氏锐化,也不必中值滤波,可直接canny提取边界!!>> ifill=igcfill|imedfill;>> figure,imshow(ifill);>> nosmall=bwareaopen(ifill,150);>> figure,imshow(nosmall);

%当t=0.55时,阈值处理再canny运算的效果 >> imhist(ig);>> t=0.55;>> v=imadjust(ig,[0 t],[],1);>> vcanny=edge(v,'canny');>> vfill=imfill(vcanny,'hole');>> figure,imshow(vfill);>> ifill=igcfill|vfill;>> figure,imshow(ifill);>> nosmall=bwareaopen(ifill,150);>> figure,imshow(nosmall);

%当t=0.6时,阈值处理再canny运算的效果的效果 >> t=0.6;>> v=imadjust(ig,[0 t],[],1);>> vcanny=edge(v,'canny');>> vfill=imfill(vcanny,'hole');>> figure,imshow(vfill);>> ifill=igcfill|vfill;>> figure,imshow(ifill);>> nosmall=bwareaopen(ifill,150);>> figure,imshow(nosmall);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %处理花椒子

>> i=imread('nut.bmp');%figure,imshow(i);ig=rgb2gray(i);figure,imshow(ig);>> imed=medfilt2(ig);imedcanny=edge(imed,'canny');imedfill=imfill(imedcanny,'hole');%figure,imshow(imedfill);nosmall=bwareaopen(imedfill,150);>> figure,imshow(nosmall);>> [labeled,numobjects]=bwlabel(nosmall,4);>> rgb_label=label2rgb(labeled,@spring,'c','shuffle');%>> figure,imshow(rgb_label);>> nutdata=regionprops(labeled,'all');>> min([nutdata.Solidity])

>> rectangle('Position', [253.5000 207.5000 26.0000 28.0000])%画矩形

>> rectangle('Position', [250.5000 50.5000 27.0000 26.0000])>> figure,imshow(nutdata(1).Image)%只显示1号物体的图像

>> figure,imshow(nutdata(1).ConvexImage)%画出1号物体的凸多边形 >> std([nutdata.Eccentricity])/ mean([nutdata.Eccentricity])std([nutdata.Area])/ mean([nutdata.Area])std([nutdata.Solidity])/ mean([nutdata.Solidity])>> std([nutdata.Centroid])/ mean([nutdata.Centroid])std([nutdata.MajorAxisLength])/ mean([nutdata.MajorAxisLength])std([nutdata.MinorAxisLength])/ mean([nutdata.MinorAxisLength])std([nutdata.Orientation])/ mean([nutdata.Orientation])std([nutdata.EquivDiameter])/ mean([nutdata.EquivDiameter])std([nutdata.Extent])/ mean([nutdata.Extent])std([nutdata.Extrema])/ mean([nutdata.Extrema])

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %处理花椒皮 close all;clc;clear;>> i=imread('p.bmp');imshow(i);ig=rgb2gray(i);figure,imshow(ig);imed=medfilt2(ig);imedcanny=edge(imed,'canny');figure,imshow(imedcanny);>> se90=strel('line',2,90);se0=strel('line',2,0);bwsdil=imdilate(imedcanny,[se90 se0]);figure,imshow(bwsdil),title('dilated');ifill=imfill(bwsdil,'holes');figure,imshow(ifill);>> bwero=imerode(ifill,[se90 se0]);>> figure,imshow(bwero);>> nosmall=bwareaopen(bwero,150,4);>> figure,imshow(nosmall);>> nobord=imclearborder(nosmall,4);>> figure,imshow(nobord);>> [labeled,numobjects]=bwlabel(nobord,4);>> numobjects >> pdata=regionprops(labeled,'all');>> max([pdata.Solidity])>> std([pdata.Solidity])/mean([pdata.Solidity])

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %处理混合图像 >> clear;clc;close all;>> i=imread('m.bmp');%>> figure,imshow(i);>> ig=rgb2gray(i);figure,imshow(ig);imed=medfilt2(ig);%>> figure,imshow(imed);imedcanny=edge(imed,'canny');%>> figure,imshow(imedcanny);>> se90=strel('line',2,90);se0=strel('line',2,0);bwsdil=imdilate(imedcanny,[se90 se0]);%figure,imshow(bwsdil),title('dilated');ifill=imfill(bwsdil,'holes');%figure,imshow(ifill);bwero=imerode(ifill,[se90 se0]);%figure,imshow(bwero);>> nosmall=bwareaopen(bwero,150,4);%figure,imshow(nosmall);nobord=imclearborder(nosmall,4);figure,imshow(nobord);>> [labeled,numobjects]=bwlabel(nobord,4);>> numobjects >> rgb_label=label2rgb(labeled,@spring,'c','shuffle');figure,imshow(rgb_label);>> mexdata=regionprops(labeled,'all');hold on;%以下内容画在同一figure中 centr=[mexdata.Centroid];%寻找重心位置 nums=1:numobjects;for k = 1:numobjects soli=mexdata(k).Solidity;soli_string=sprintf('%2.2f',soli);%等价于转字符串 % signal=num2str(nums(k));signal=sprintf('%d',k);%直接使用打印语句打印序号 text(centr(2*k-1),centr(2*k),signal)%按序标记物体

text(centr(2*k-1)-30,centr(2*k)-30,soli_string)%标注每个Solidity值 end

for k=1:numobjects plot(mexdata(k).ConvexHull(:,1),mexdata(k).ConvexHull(:,2),...'b','Linewidth',2)end

%画出1和2号物体的外接矩形

%>> rectangle('position',[9.5000 224.5000 62.0000 63.0000])%>> rectangle('position',[65.5000 141.5000 34.0000 39.0000])%画出每个物体的外接矩形 bb=[mexdata.BoundingBox];for k=1:numobjects rectangle('position',[bb(4*k-3)bb(4*k-2)bb(4*k-1)bb(4*k)])end

%>> figure,imshow(mexdata(1).Image)%只显示1号物体的图像

%>> figure,imshow(mexdata(1).ConvexImage)%画出1号物体的凸多边形 %>> figure,imshow(mexdata(2).Image)%只显示2号物体的图像

%>> figure,imshow(mexdata(2).ConvexImage)%画出2号物体的凸多边形 %画出单个物体的凸多边形的填充图形 for k=1:numobjects figure,imshow(mexdata(k).ConvexImage)end

%只显示Solidity>0.92的物体的图像 >> idx = find([mexdata.Solidity] > 0.92);>> BW2 = ismember(labeled,idx);>> figure,imshow(BW2)

>> mexdata=regionprops(labeled,'all');>> %只显示Solidity<0.92的物体的图像 idx = find([mexdata.Solidity] < 0.92);bw2 = ismember(labeled,idx);figure,imshow(bw2)%mexdata.Solidity;

>> numdown=find([mexdata.Solidity]<0.92);mexdata(numdown,:)=[];>> mexdata

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %2006.6.19 %roipoly函数的用法 I = imread('eight.tif');c = [222 272 300 270 221 194];r = [21 21 75 121 121 75];BW = roipoly(I,c,r);imview(I), imview(BW)

%可以使用下面的方法创建相应的向量: regionprops(L,'Area');allArea = [stats.Area];

%创建一个只包含面积大于80的二值图像 idx = find([stats.Area] > 80);BW2 = ismember(L,idx);

%只显示某个下标所对应的物体图像 bw2=ismember(L,N);figure,imshow(bw2);

%在调用regionprops之前必须将二值图像转变为标注矩阵 L = bwlabel(BW);%或者

L = double(BW);

%将matlab数据写到excel中 a=ones(3);success = xlswrite('c:/matlab/work/myworkbook.xls',a,'A2:C4')%将行矩阵转换为列矩阵 a=[1 2 3 4 5 6];b=transpose(a);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %2006.6.22球形物体的检测和标识(循环检测和标识算法)clc;clear;close all;%Step 1: Read image %Step 2: Threshold the image %Step 3: Remove the noise %Step 4: Find the boundaries %Step 5: Determine which objects are round >> RGB = imread('pillsetc.png');imshow(RGB)>> I = rgb2gray(RGB);threshold = graythresh(I);bw = im2bw(I,threshold);imshow(bw)>> % remove all object containing fewer than 30 pixels bw = bwareaopen(bw,30);>> figure,imshow(bw)>> % fill a gap in the pen's cap se = strel('disk',2);bw = imclose(bw,se);>> figure,imshow(bw)>> % fill any holes, so that regionprops can be used to estimate % the area enclosed by each of the boundaries bw = imfill(bw,'holes');>> figure,imshow(bw)>> [B,L] = bwboundaries(bw,'noholes');>> % Display the label matrix and draw each boundary figure,imshow(label2rgb(L, @jet, [.5.5.5]))>> hold on for k = 1:length(B)boundary = B{k};plot(boundary(:,2), boundary(:,1), 'w', 'LineWidth', 2)end >> stats = regionprops(L,'Area','Centroid');>> stats = regionprops(L,'Area','Centroid');threshold = 0.94;% loop over the boundaries for k = 1:length(B)% obtain(X,Y)boundary coordinates corresponding to label 'k' boundary = B{k};% compute a simple estimate of the object's perimeter delta_sq = diff(boundary).^2;perimeter = sum(sqrt(sum(delta_sq,2)));

% obtain the area calculation corresponding to label 'k' area = stats(k).Area;

% compute the roundness metric metric = 4*pi*area/perimeter^2;

% display the results metric_string = sprintf('%2.2f',metric);% mark objects above the threshold with a black circle if metric > threshold centroid = stats(k).Centroid;plot(centroid(1),centroid(2),'ko');end

text(boundary(1,2)-35,boundary(1,1)+13,metric_string,'Color','y',...'FontSize',14,'FontWeight','bold');end >> title(['Metrics closer to 1 indicate that ',...'the object is approximately round']);

下载MATLAB小波函数小结(小编整理)word格式文档
下载MATLAB小波函数小结(小编整理).doc
将本文档下载到自己电脑,方便修改和收藏,请勿使用迅雷等下载。
点此处下载文档

文档为doc格式


声明:本文内容由互联网用户自发贡献自行上传,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任。如果您发现有涉嫌版权的内容,欢迎发送邮件至:645879355@qq.com 进行举报,并提供相关证据,工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。

相关范文推荐

    MATLAB实训小结

    实训小结 为期一周的MATLAB实训在学习与忙碌中度过了,时间虽短,但我们却真真切切的学到了知识,在现实工作中可以运用的知识。 在第一节课,我们便了解到MATLAB是世界上最流行的、......

    MatLab基础知识小结来源

    MatLab基础知识小结来源: 吴学箫的日志 matlab常用到的永久变量。 ans:计算结果的默认变量名。 i j:基本虚数单位。 eps:系统的浮点(F10a9Bg个oht): inf: 无限大,例1/0 nan NaN:非......

    matlab的nntool工具箱小结

    拟合以及插值还有逼近是数值分析的三大基础工具,通俗意义上它们的区别在于:拟合是已知点列,从整体上靠近它们;插值是已知点列并且完全经过点列;逼近是已知曲线,或者点列,通过逼近使......

    MATLAB实验小结论文 数学建模

    数学建模论文题 目 求π的近似值的数学建模问题 学 院 材料科学与工程 专业班级 学生姓名 成 绩 年 05 月 20 MATLAB 2010 日 摘要 这个学期,我们开了MATLAB的课程,因为是......

    matlab小学期学习感想

    小学期学习感想 大一小学期我们学习了matlab软件,这是一个十分实用和重要的软件。学习MATLAB,感觉这是一个很好的软件,语言简便,实用性强。作为一个做新手,想要学习好这门语言,可......

    小波变换 matlab 总结范文

    小波变换matlab总结 目录 一、预置工具 .................................................................................................................................

    运筹学单纯形法matlab程序(小编整理)

    function [xx,fm]=myprgmh(m,n,A,b,c) B0=A(:,1:m); cb=c(:,1:m); xx=1:n; sgm=c-cb*B0^-1*A; h=-1; sta=ones(m,1); for i=m+1:n if sgm(i)>0 h=1; end end while h>0 [ms......

    matlab作业题

    第一章 MATLAB环境 1、MATLAB通用操作界面窗口包括哪些?命令窗口、历史命令窗口、当前目录窗口、工作空间窗口各有哪些功能? 答:MATLAB通用操作界面窗口包括:命令窗口、历史命令......