apaam/mlpack
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

polished niche/fastica

Browse Source
This commit is contained in:
tekhnofiend
2008-01-17 18:04:45 +00:00
parent c18ef480c6
commit 44d80efb58
2 changed files with 12 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files
fastlib/u/niche/functional/prelim_funcica.m
+8 -2
View File
@@ -30,15 +30,21 @@ x = x - repmat(mean(x')', 1, N);
% generate b-spline basis curves
t = linspace(0,1,3000);
t = linspace(0,1,1000);
load s1s2_10;
s = [s1(t); s2(t)]';
z = normrnd(zeros(length(t), N), 1);
z = normrnd(zeros(length(t), N), .2);
data = s * x;
data = data + z;
indices = 1:10:1000;
t = t(indices);
data = data(indices,:);
%data = noisy_data;
fastlib/u/niche/functional/test_smooth_funcica.m
+4 -4
View File
@@ -7,13 +7,13 @@
N = size(data, 2);
p = 30; % hardcoded for now
mybasis = create_bspline_basis([0 1], p, 4);
mybasis = create_bspline_basis([0 .9910], 100, 4);
basis_curves = eval_basis(t, mybasis);
basis_inner_products = full(eval_penalty(mybasis, int2Lfd(0)));
myfd_data = data2fd(data, t, mybasis);
cut_fraction = .5;
cut_fraction = .01;
cut = round(cut_fraction * N);
data_coef = getcoef(myfd_data);
@@ -26,8 +26,8 @@ num_tests = 1;
%data = data(indices,:);
%lambda_set = 0;
lambda_set = [0 1e-6 1e-5 1e-4 1e-3 5e-3 1e-2];
lambda_set = 0;
%lambda_set = [0 1e-6 1e-5 1e-4 1e-3 5e-3 1e-2];
myfdPar_set = cell(1,length(lambda_set));
for lambda_i = 1:length(lambda_set)