eps.rar_EPS_EPS滤波_edge preserving_边缘平滑_边缘平滑 matlab

clear clc for i=1:60 if i<=30 A(i)=0; else A(i)=1; end end % plot(A(1,1:60)); % axis([0 60 -0.5 1.5]) jj=awgn(A,10); tianjia1=ones(1,10)*jj(1); tianjia2=ones(1,10)*jj(60); data=[tianjia1,jj,tianjia2]; %---------------------------------------eps---------------------------- for i=11:70 k=1; for j=0:10 c(k,:)=data(i+j-10:i+j); k=k+1; end for z=1:10 ss(z)=std(c(z,:)); end mins=min(ss); d=c(1,:); for q=1:10 if mins==std(c(q,:)); d=c(q,:); end ave=mean(d); E(i-10)=ave; end plot(E); axis([1 60 -0.5 1.5]); end %---------------------------------------2DEPS---------------------------- clear clc % load l2_sc_sc_pc.sgy; % sig=l2_sc_sc_pc.sgy; A=ones(5,5); a1=10*A; b1=8*A; C1=[a1,b1];C2=[b1,a1]; c=[C1;C2]; c1=[c,c,c,c,c,c]; c2=[c1;c1;c1;c1;c1;c1]; c3=[c2,c2(:,1:5)]; c4=[c3;c3(1:5,:)]; %%%求取信号% obj=awgn(c4,0.5); [m,n]=size(obj); E=[]; obj1=obj; for i=3:m-2 for j=3:n-2 d=obj(i-2:i+2,j-2:j+2); z=1; for x=2:4 for y=2:4 d1=d(x-1:x+1,y-1:y+1); std(z)=std2(d1); z=z+1; end end mins=min(std); Z=d(1:3,1:3); for i1=2:4 for j1=2:4 if mins==std2(d(i1-1:i1+1,j1-1:j1+1)); Z=d(i1-1:i1+1,j1-1:j1+1); end end end ave=mean2(Z); E(i-2,j-2)=ave; end end figure(1); imagesc(obj); axis([1 60 1 60]); figure(2); imagesc(E) axis([1 60 1 60]); % z=0; % for k=1:10 % z=z+1; % [fx,fz]=gradient(obj); % [fxx,fxz]=gradient(fx); % [fzx,fzz]=gradient(fz); % [m,n]=size(obj); % dt=0.2; % miu=0.01; % H=fspecial('gaussian',3,1); % for i=1:m % for j=1:n % d=[fx(i,j);fz(i,j)]; % J=d*d'; % Nj=filter2(H,J); % [v,d]=eig(Nj); % alfa=sqrt((Nj(1,1)-Nj(2,2))^2+4*Nj(1,2)^2); % c1=miu; % if d(1,1)==d(2,2); % c2=miu; % else c2=miu+(1-miu)*exp(-1/alfa^2); % end % D=v*[c1,0;c2,0]*v'; % a=D(1,1); % c=D(2,2); % b=D(1,2); % sig1(i,j)=obj(i,j)+dt*(a*fxx(i,j)+c*fzz(i,j)+2*b*fxz(i,j)); % end % end % obj=sig1; % end %