DUAN的专栏

假设圆半径为r，则[x, y] = meshgrid( -r : r );circle_mat = x.^2 + y.^2;figure, imshow( circle_mat );title( Circle Mask );   

UNIX/LINUX环境中由于dir函数无法忽略大小写，为此可做一简单改造。function file_list = dir_unix( path, ext_type ) %------------------------------------------------------ % This function lists files with specified extension % ign

Armijo算法实现： function mk = armijo( fun, xk, rho, sigma, gk ) assert( rho > 0 && rho < 1 ); assert( sigma > 0 && sigma < 0.5 ); mk = 0; max_mk = 100; while mk <= max_mk x = xk - rho^mk * gk;

function [label, model, llh] = emgm(X, init) % EM algorithm for Gaussian mixture model %% initialization fprintf('EM for Gaussian mixture: running ... '); R = initialization(X,init); tol = 1e-6; maxi