You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
38 lines
848 B
38 lines
848 B
function smooth_compare_DH()
|
|
load DH_afterSmooth25.mat DH
|
|
DH_after = DH;
|
|
|
|
load DH_beforeSmooth25.mat DH
|
|
DH_before= DH;
|
|
|
|
nx=60;ny=1;nz=20;
|
|
i=1;j=1; epss=0.014;
|
|
% i=1;j=2; epss=0.20;
|
|
% epss=0.5;
|
|
Dij_after = zeros(nx,ny,nz);
|
|
Dij_before= zeros(nx,ny,nz);
|
|
for ix=1:nx
|
|
for iy=1:ny
|
|
for iz=1:nz
|
|
Dij_after(ix,iy,iz) =DH_after{ix,iy,iz}(i,j);
|
|
Dij_before(ix,iy,iz)=DH_before{ix,iy,iz}(i,j);
|
|
% Dij_after(ix,iy,iz) =trace(DH_after{ix,iy,iz});
|
|
% Dij_before(ix,iy,iz)=trace(DH_before{ix,iy,iz});
|
|
end
|
|
end
|
|
end
|
|
|
|
figure()
|
|
imagesc(squeeze(Dij_before(:,1,:)).')
|
|
axis equal; axis tight; axis off;colorbar;
|
|
set(gca,'YDir','normal')
|
|
|
|
|
|
figure()
|
|
error=abs(Dij_after-Dij_before)./Dij_before;
|
|
error=error.*(abs(Dij_before)>epss);
|
|
imagesc(squeeze(error(:,1,:)).')
|
|
axis equal; axis tight; axis off;colorbar;
|
|
set(gca,'YDir','normal')
|
|
mean(error,'all')
|
|
end
|