Homec4science
Diffusion Slepians (master) Browse

Slepi/Code/Script_Figure1.m
3fc67bd6510bmaster

Last Change
Raw File

Script_Figure1.m
View Options

% Construction of the noise and adding it
BB_noise1 = zeros(n_ROI,n_TP);
BB_noise2 = zeros(n_ROI,n_TP);
for t = 1:n_TP
rvec = randperm(n_ROI);
for i = 1:50
BB_noise1(rvec(i),t) = 1;
BB_noise2(rvec(i),t) = 1;
end
for i = 51:150
BB_noise2(rvec(i),t) = 1;
end
end
X_noisy_11 = X_GT + 0.5 * BB_noise1 .* randn(n_ROI,n_TP);
X_noisy_12 = X_GT + 3.5 * BB_noise1 .* randn(n_ROI,n_TP);
X_noisy_21 = X_GT + 0.5 * BB_noise2 .* randn(n_ROI,n_TP);
X_noisy_22 = X_GT + 3.5 * BB_noise2 .* randn(n_ROI,n_TP);
figure; gsp_plot_signal(G,X_noisy_11(:,100)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_11(:,300)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_11(:,500)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_11(:,700)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_11(:,900)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_12(:,100)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_12(:,300)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_12(:,500)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_12(:,700)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_12(:,900)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_21(:,100)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_21(:,300)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_21(:,500)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_21(:,700)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_21(:,900)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_22(:,100)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_22(:,300)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_22(:,500)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_22(:,700)); colormap(jet); caxis([-0.25 1.25]);
figure; gsp_plot_signal(G,X_noisy_22(:,900)); colormap(jet); caxis([-0.25 1.25]);
CM = brewermap(100,'RdBu');
figure; gsp_plot_signal(G,X_noisy_11(:,100)); colormap(flipud(CM)); caxis([-4 4]);
figure; gsp_plot_signal(G,X_noisy_11(:,300)); colormap(flipud(CM)); caxis([-4 4]);
figure; gsp_plot_signal(G,X_noisy_11(:,500)); colormap(flipud(CM)); caxis([-4 4]);
figure; gsp_plot_signal(G,X_noisy_11(:,700)); colormap(flipud(CM)); caxis([-4 4]);
figure; gsp_plot_signal(G,X_noisy_11(:,900)); colormap(flipud(CM)); caxis([-4 4]);
figure; gsp_plot_signal(G,X_noisy_22(:,100)); colormap(flipud(CM)); caxis([-4 4]);
figure; gsp_plot_signal(G,X_noisy_22(:,300)); colormap(flipud(CM)); caxis([-4 4]);
figure; gsp_plot_signal(G,X_noisy_22(:,500)); colormap(flipud(CM)); caxis([-4 4]);
figure; gsp_plot_signal(G,X_noisy_22(:,700)); colormap(flipud(CM)); caxis([-4 4]);
figure; gsp_plot_signal(G,X_noisy_22(:,900)); colormap(flipud(CM)); caxis([-4 4]);
CM = brewermap(100,'RdYlBu');
figure;
set(gcf,'color','w');
surf(0:20:200,0:0.5:4,median(Error_GFT,3));
colormap(flipud(CM));
caxis([0 0.3]);
figure;
set(gcf,'color','w');
surf(0:20:200,0:0.5:4,median(Error_LGFT,3));
colormap(flipud(CM));
caxis([0 0.3]);
figure;
set(gcf,'color','w');
surf(0:20:200,0:0.5:4,median(Error_SLEP,3));
colormap(flipud(CM));
caxis([0 0.3]);
CM = brewermap(100,'RdBu');
figure;
set(gcf,'color','w');
surf(0:20:200,0:0.5:4,median(Error_GFT,3)-median(Error_LGFT,3));
colormap(flipud(CM));
caxis([-0.1 0.1]);
figure;
set(gcf,'color','w');
surf(0:20:200,0:0.5:4,median(Error_GFT,3)-median(Error_SLEP,3));
colormap(flipud(CM));
caxis([-0.1 0.1]);
c4science · Help