123456789101112131415 |
- function [sample_pnt,sampleDen] = sphDensity(sphpoint,sample_pntNum,kernel_sig)
- if ~exist('kernel_sig','var'), kernel_sig = [];end
- if isempty(kernel_sig),kernel_sig = .01;end
- gaussFunc = @(x,sig) 1./(sig*sqrt(2*pi)).*exp(-x.^2./(2*sig^2));
- if ~isa(sample_pntNum,'floQ')
- sample_pnt = CPS_sampling('pnt_num',sample_pntNum);
- sample_pnt = floQ(sample_pnt);
- else
- sample_pnt = sample_pntNum;
- end
- pnt_dist = real(asin(norm(sample_pnt+conj(sphpoint))./2)).*2;
- weightMat = gaussFunc(pnt_dist,kernel_sig);
- sampleDen = mean(weightMat,2);
- end
|