plotdb.m

细胞生长结构可视化工具箱-MATLAB Toolbox1999.zip

% determine decision boundary in actual ip space for 2D Ripley data
	dbp1=zeros(NoPts+1);
	dbp2=zeros(NoPts+1);
	
	maxxdb=1;
	minxdb=-1;
	maxydb=1.2;
	minydb=-0.2;

	xidb=minxdb:(maxxdb-minxdb)/NoPts:maxxdb;
	yidb=maxydb:-(maxydb-minydb)/NoPts:minydb;

	ii=0;
	for i=minxdb:(maxxdb-minxdb)/NoPts:maxxdb
		ii=ii+1;
		jj=NoPts+2;
		for j=maxydb:-(maxydb-minydb)/NoPts:minydb
			jj=jj-1;
			Ip=[i,j];
			for k = 1:size(wvis,1)
				Act(k) = exp(-((norm(Ip-w(k,:),2)^2)./(sigma(k)).^2));
			end
			dbp1(jj,ii) = Act*z1;
			dbp2(jj,ii) = Act*z2;		
		end
	end

	total = dbp1+dbp2;
%	figure(6);
%	imagesc(xidb,yidb,dbp1./total);
%	figure(7);
%	imagesc(xidb,yidb,dbp2./total);
%	figure(5);
	contour(xidb,[minydb:(maxydb-minydb)/NoPts:maxydb],(dbp2./total)-(dbp1./total),[-1000,0,1000],'k');
	title('GCS network in actual 2D input space with decision boundary');