xiaobo

clear all I=imread( lena.bmp ) figure imshow(I) I2=imrotate(I,-4, bilinear ) 逆时针旋转4度 figure imshow(I2) I3=fliplr(I) 垂直镜像 figure imshow(I3) I4=imresize(I,0.5, bilinear ) 缩小为原图的1/2 figure imshow(I4) A=double(I) 转换为double类型 计算7个不变矩 [nc,nr]=size(A) [x,y]=meshgrid(1:nr,1:nc) 得到网格 x=x(:) y=y(:) A=A(:) m.m00=sum(A) if m.m00==0 m.m00=eps end m.m10=sum(x.*A) m.m01=sum(y.*A) 计算均值 xmean=m.m10/m.m00 ymean=m.m01/m.m00 计算中心矩(li9_23.rar cm.cm00=m.m00 cm.cm02=(sum((y-ymean).^2.*A))/(m.m00^2) cm.cm03=(sum((y-ymean).^3.*A))/(m.m00^2.5) cm.cm11=(sum((x-xmean).*(y-ymean).*A))/(m.m00^2) cm.cm12=(sum((x-xmean).*(y-ymean).^2.*A))/(m.m00^2.5) cm.cm20=(sum((x-xmean).^2.*A))/(m.m00^2) cm.cm21=(sum((x-xmean).^2.*(y-ymean).*A))/(m.m00^2.5) cm.cm30=(sum((x-xmean).^3.*A))/(m.m00^2.5) im(1)=cm.cm20+cm.cm02 im(2)=(cm.cm20-cm.cm02)^2+4*cm.cm11^2 im(3)=(cm.cm30-3*cm.cm12)^2+(3*cm.cm21-cm.cm03)^2 im(4)=(cm.cm30+cm.cm12)^2+(cm.cm21+cm.cm03)^2 im(5)=(cm.cm30-3*cm.cm12)*(cm.cm30+cm.cm12)... *((cm.cm30+cm.cm12)^2-3*(cm.cm21+cm.cm03)^2)... +(3*cm.cm21-cm.cm03)*(cm.cm21+cm.cm03)... +(3*(cm.cm30+cm.cm12)^2-(cm.cm21+cm.cm03)^2) im(6)=(cm.cm20-cm.cm02)*((cm.cm30+cm.cm12)^2-(cm.cm21+cm.cm03)^2)... +4*cm.cm11*(cm.cm30+cm.cm12)*(cm.cm21+cm.cm03) im(7)=(3*cm.cm21-cm.cm03)*(cm.cm30+cm.cm12)... *((cm.cm30+cm.cm12)^2-3*(cm.cm21+cm.cm03)^2)... +(3*cm.cm12-cm.cm30)*(cm.cm21+cm.cm03)... *()

使用差分法解一类Bessel方程在不同边界条件下的离散解(the use of differential method for a type Bessel equation in different boundary conditions of the Discrete Solutions)

Generalized Nyquist diagram for stability analysis of linear time invariant systems (Nyquist on arbitrary contour)

许多程序还提供了图形用户界面.可直接用于各种智能算法设计. (many procedures also provide a graphical user interface. Can be directly applied to various smart algorithm design.)

粗糙集理论的属性约简方法及应用研究，很有用的哈！！！！(Rough set attribute reduction method and application, very useful, ha! ! ! !)

Example of key performance indicators applied on a real problem

This code is only intended to explain the basics behind the bandelet transform, and to be able to play with the major concepts and properties of these orthogonal basis. It does not contain any coding scheme, nor a fast implementation(This code is only intended to explain the basics behind the bandelet transform, and to be able to play with the major concepts and properties of these orthogonal basis. It does not contain any coding scheme, nor a fast implementation)

合成地震记录是用声波测井或垂直地震剖面资料经过人工合成转换成的地震记录（地震道） 合成记录的制作是一个简化的一维正演的过程，合成记录F(t)是地震子波S(t)与反射系数R(t)褶积的结果。(Synthetic seismogram is sonic or vertical seismic profile data after conversion into synthetic seismograms (seismic trace) the production of synthetic seismogram is a simplified one-dimensional forward modeling process, synthetic seismogram F (t) is the seismic wavelet S (t) and the reflection coefficient R (t) the convolution results.)

设计采用模糊控制的直接转矩控制，和传统PI调节器相比较(Design of fuzzy control of direct torque control, and compared to the traditional PI regulator)

线性调频脉冲雷达中脉冲压缩与相参积累程序，经过验证(Linear FM pulse radar pulse compression and coherent integration program.)