rjtpspline.m

来自「大家可以参考的例子滤波源码」· M 代码 · 共 19 行

function [y] = rjtpSpline(mu,x);
% 
if nargin < 2, error('Not enough input arguments.'); end

[N,d] = size(x);      % N = number of data, d = dimension of x.
y=zeros(N,1);
for j=1:N,
  z=norm(x(j,:)-mu(1,:));             % Euclidean distance.
  y(j,1)=z*log(z+1e-99);        % Gaussian.
end;