失效有关曲率处理matlab源码程序 - matlab算法设计

下载频道> 资源分类> matlab源码> 算法设计> 失效有关曲率处理matlab源码程序

标题：失效有关曲率处理matlab源码程序
分享到：

所属分类: 算法设计资源类型：程序源码文件大小: 4.34 MB 上传时间: 2019-06-17 22:49:19 下载次数: 7 资源积分：1分提供者: zhangsan456 code

内容：
失效有关曲率处理matlab源码程序

%This program calculate the wave guide by a photonic crystal.

%For two dimension case.

%Only for TM case, in this version.

%Trial Version

clear;

%tic

%Initial parameters and other things.

W=0.36; %Normalized frequency

%The following parameters are control parameters.

WaveGuide=1; %If this program is for Wave Guide? If so, please specify 1.

IsMovie=1; %If you want to play movie, please use 1.

IsFigure=1; %If it will plot the figures? If so, please specify 1.

WantToSeeEp=1; %Do you want to see the distrubution of Ep? If so, please specify 1.

%End of defining control parameters

MLatx=11; %How many Lattice cell in x direction.

MLaty=11; %How many Lattice cell in y direction.

NMlat=21; %The gird number in each Lattice Cell.

%SHOULD BE ODD INTEGER!!!!

if mod(NMlat,2)==0

NMlat=NMlat+1;

end %Force it to be a odd integer!

NTx=MLatx*NMlat+1; %It is the number of the Grid along x axis.

NTy=MLaty*NMlat+1; %It is the number of the Grid along y axis.

if WaveGuide==1

Nrow=4; %The row number of coloumns between the PML boundary and the waveguide.

end

NPML=12; %How many PML s will be used in our computation.

NTimeSteps=2000; %Total number of Time Steps

Meach=20; %Define the interval for plot figures if IsFigure==1.

%This also works for saving intervals.

R=0.2; %The radius of dielectric columns

ea=11.4; %The dielectric constant of these columns.

Zmax=0.6; %The maximum value for z axis when plotting figures.

Colormax=0.6; %The maximum value for colormap when plotting figures.

%Some constants

mu0=4*pi*1.0e-7; %Epsilon Zero, if using Gauss Unit, it equals to 1.

e0=8.85*1e-12; %Mu Zero, if using Gauss Unit, it equals to 1.

c=1/sqrt(mu0*e0); %The light speed.

factor=mu0/e0; %The factor between conductivity and permeability.

%Permeability=Conductivity*factor, in PML.

a=1;%e-6; %The lattice constant.

W=W*(2*pi*c/a); %frequency

Dx=a/NMlat; %Delta x.

Dy=Dx; %Delta y.

Dt=1/sqrt(1/(Dx*Dx)+1/(Dy*Dy))/c; %Time interval

%tic

%In the following partm we define the dielectric constants:

Ep=ones(NTx-1,NTy-1)*e0;

Ep_cell=ones(NMlat,NMlat)*e0;

x=-(NMlat-1)/2*Dx:Dx:(NMlat-1)/2*Dx;

[X,Y]=meshgrid(x);

X=X';

Y=Y';

flag=find(sqrt(X.^2+Y.^2)<R);

Ep_cell(flag)=e0*ea;

Ep=repmat(Ep_cell,MLatx,MLaty);

if WaveGuide==1

eb=1; %The dielectric constant in the waveguide

Ep(1:6*NMlat,5*NMlat+1:(5+1)*NMlat)=e0*eb;

%Ep(1:(MLatx)*NMlat,4*NMlat+1:(4+1)*NMlat)=e0*eb;

for ii=6:10

Ep(ii*NMlat+1:(ii+1)*NMlat,ii*NMlat+1:(ii+1)*NMlat)=e0*eb;

Ep((ii)*NMlat+1:(ii+1)*NMlat,(10-ii)*NMlat+1:(10-ii+1)*NMlat)=e0*eb;

end

%Ep((MLatx-Nrow-1)*NMlat+1:(MLatx-Nrow)*NMlat,Nrow*NMlat+1:(Nrow+1)*NMlat)=Ep_cell;

%Ep((MLatx-Nrow-2)*NMlat+1:(MLatx-Nrow-1)*NMlat,(Nrow+1)*NMlat+1:(Nrow+2)*NMlat)=e0*eb;

end

%toc

if WantToSeeEp==1

x=0:Dx:(NMlat*MLatx-1)*Dx;

x=x-(NMlat*MLatx-1)*Dx/2;

y=0:Dy:(NMlat*MLaty-1)*Dy;

y=y-(NMlat*MLaty-1)*Dy/2;

[X,Y]=meshgrid(x,y);

X=X';

Y=Y';

surf(X,Y,Ep/e0);

shading interp;

view(0,90);

axis([min(x), max(x),min(y), max(y)])

axis off;

disp('Press any key to continue...');

pause

end

%End of defining the Ep.

if IsFigure==1

%Define the X Y coordinate for figures

x=0:Dx:(NMlat*MLatx-1+2*NPML)*Dx;

x=x-(NMlat*MLatx-1+2*NPML)*Dx/2;

y=0:Dy:(NMlat*MLaty-1+2*NPML)*Dy;

y=y-(NMlat*MLaty-1+2*NPML)*Dy/2;

[X,Y]=meshgrid(x,y);

X=X';

Y=Y';

end

%Define the Ez, Hx, Hy, which are in the inside region.

Ez=zeros(NTx-1,NTy-1);

Hx=zeros(NTx-1,NTy);

Hy=zeros(NTx,NTy-1);

%Parameters about PML:

n=4; %The order of the polynomial that decribes the conductivity profile.

R=1e-10;

Delta=NPML*Dx;

SigmaMax=-(n+1)*e0*c*log(R)/(Delta*2);

NUM=NPML*2:-1:1;

Sigmax=SigmaMax*((NUM*Dx/2+Dx/2).^(n+1)-(NUM*Dx/2-Dx/2).^(n+1))/(Delta^n*Dx*(n+1));

Sigmay=Sigmax;

SigmaBound=SigmaMax*(Dx/2).^(n+1)/(Delta^n*Dx*(n+1));

%Sigmax=SigmaMax*(NUM/6).^(n+1); %Another way, the definition.

EzxPML1=zeros(NPML,NPML);

EzyPML1=zeros(NPML,NPML);

HxPML1=zeros(NPML,NPML);

HyPML1=zeros(NPML,NPML); %Zone 1

Sigmax_z1=repmat(Sigmax(2:2:NPML*2)',1,NPML);

Sigmax_x1=repmat(Sigmax(2:2:NPML*2)',1,NPML);

Sigmax_y1=repmat(Sigmax(1:2:NPML*2-1)',1,NPML);

Sigmay_z1=fliplr(repmat(Sigmax(2:2:NPML*2),NPML,1));

Sigmay_x1=fliplr(repmat(Sigmax(1:2:NPML*2-1),NPML,1));

Sigmay_y1=fliplr(repmat(Sigmax(2:2:NPML*2),NPML,1)); %Zone 1

EzxPML2=zeros(NPML,NPML);

EzyPML2=zeros(NPML,NPML);

HxPML2=zeros(NPML,NPML);

HyPML2=zeros(NPML,NPML); %Zone 2

Sigmax_z2=flipud(Sigmax_z1);

Sigmax_x2=flipud(Sigmax_x1);

Sigmax_y2=flipud(Sigmax_y1);

Sigmay_z2=Sigmay_z1;

Sigmay_x2=Sigmay_x1;

Sigmay_y2=Sigmay_y1; %Zone 2

EzxPML3=zeros(NPML,NPML);

EzyPML3=zeros(NPML,NPML);

HxPML3=zeros(NPML,NPML);

HyPML3=zeros(NPML,NPML); %Zone 3

Sigmax_z3=Sigmax_z1;

Sigmax_x3=Sigmax_x1;

Sigmax_y3=Sigmax_y1;

Sigmay_z3=fliplr(Sigmay_z1);

Sigmay_x3=fliplr(Sigmay_x1);

Sigmay_y3=fliplr(Sigmay_y1); %Zone 3

EzxPML4=zeros(NPML,NPML);

EzyPML4=zeros(NPML,NPML);

HxPML4=zeros(NPML,NPML);

HyPML4=zeros(NPML,NPML); %Zone 4

Sigmax_z4=flipud(Sigmax_z1);

Sigmax_x4=flipud(Sigmax_x1);

Sigmax_y4=flipud(Sigmax_y1);

Sigmay_z4=fliplr(Sigmay_z1);

Sigmay_x4=fliplr(Sigmay_x1);

Sigmay_y4=fliplr(Sigmay_y1); %Zone 4

EzxPMLA=zeros(NTx-1,NPML);

EzyPMLA=zeros(NTx-1,NPML);

HxPMLA=zeros(NTx-1,NPML);

HyPMLA=zeros(NTx,NPML); %Zone A

Sigmay_zA=repmat(Sigmay_z1(1,:),NTx-1,1);

Sigmay_xA=repmat(Sigmay_x1(1,:),NTx-1,1);

Sigmay_yA=repmat(Sigmay_y1(1,:),NTx,1); %Zone A

EzxPMLB=zeros(NTx-1,NPML);

EzyPMLB=zeros(NTx-1,NPML);

HxPMLB=zeros(NTx-1,NPML);

HyPMLB=zeros(NTx,NPML); %Zone B

Sigmay_zB=fliplr(Sigmay_zA);

Sigmay_xB=fliplr(Sigmay_xA);

Sigmay_yB=fliplr(Sigmay_yA); %Zone B

EzxPMLC=zeros(NPML,NTy-1);

EzyPMLC=zeros(NPML,NTy-1);

HxPMLC=zeros(NPML,NTy);

HyPMLC=zeros(NPML,NTy-1);%Zone C

Sigmax_zC=repmat(flipud(Sigmay_z1(1,:)'),1,NTy-1);

Sigmax_xC=repmat(flipud(Sigmay_x1(1,:)'),1,NTy);

Sigmax_yC=repmat(flipud(Sigmay_y1(1,:)'),1,NTy-1); %Zone C

EzxPMLD=zeros(NPML,NTy-1);

EzyPMLD=zeros(NPML,NTy-1);

HxPMLD=zeros(NPML,NTy);

HyPMLD=zeros(NPML,NTy-1);%Zone D

Sigmax_zD=flipud(Sigmax_zC);

Sigmax_xD=flipud(Sigmax_xC);

Sigmax_yD=flipud(Sigmax_yC); %Zone D

if IsMovie==1

Movie=moviein(NTimeSteps/Meach+1);

Mnum=1;

EzAll=[EzxPML3+EzyPML3,EzxPMLC+EzyPMLC,EzxPML1+EzyPML1;

EzxPMLB+EzyPMLB,Ez,EzxPMLA+EzyPMLA;

EzxPML4+EzyPML4,EzxPMLD+EzyPMLD,EzxPML2+EzyPML2];

HxAll=[HxPML3,HxPMLC,HxPML1;HxPMLB,Hx,HxPMLA;HxPML4,HxPMLD,HxPML2];

HyAll=[HyPML3,HyPMLC,HyPML1;HyPMLB,Hy,HyPMLA;HyPML4,HyPMLD,HyPML2];

figure(2)

surf(X,Y,real(EzAll))

shading interp;

axis([-0.5*a*MLatx-Dx*NPML 0.5*a*MLatx+Dx*NPML -0.5*a*MLaty-Dx*NPML 0.5*a*MLaty+Dx*NPML -Zmax Zmax])

zlabel('Ez');

Movie(:,1)=get ;

end

if WaveGuide==1

Npy=NMlat*Nrow+round((NMlat+1)/2);

Npx=NMlat*Nrow+round((NMlat+1)/2);

else

Npx=round(NTx/2);

Npy=round(NTy/2);

end

expboundary=exp(-SigmaBound*Dt/e0);

tic

for m=1:NTimeSteps

%tic

%The following part is for the inside region.

%H components.

Hx(:,2:NTy-1)=Hx(:,2:NTy-1)-Dt*(Ez(:,2:NTy-1)-Ez(:,1:NTy-2))/(Dy*mu0);

Hy(2:NTx-1,:)=Hy(2:NTx-1,:)+Dt*(Ez(2:NTx-1,:)-Ez(1:NTx-2,:))/(Dx*mu0);

Hx(:,NTy)=expboundary*Hx(:,NTy)-(1-expboundary)*...

(EzxPMLA(:,1)+EzyPMLA(:,1)-Ez(:,NTy-1))/(SigmaBound*factor*Dy); %Boundary A

Hx(:,1)=expboundary*Hx(:,1)-(1-expboundary)*...

(Ez(:,1)-EzxPMLB(:,NPML)-EzyPMLB(:,NPML))/(SigmaBound*factor*Dy); %Boundary B

Hy(1,:)=expboundary*Hy(1,:)+(1-expboundary)*...

(Ez(1,:)-EzxPMLC(NPML,:)-EzyPMLC(NPML,:))/(SigmaBound*factor*Dx); %Boundary C

Hy(NTx,:)=expboundary*Hy(NTx,:)+(1-expboundary)*...

(EzxPMLD(1,:)+EzyPMLD(1,:)-Ez(NTx-1,:))/(SigmaBound*factor*Dx); %Boundary D

%The following part is for ZONE A. H components.

HxPMLA(:,1:NPML-1)=exp(-Sigmay_xA(:,1:NPML-1)*Dt/e0).*HxPMLA(:,1:NPML-1)-...

(1-exp(-Sigmay_xA(:,1:NPML-1)*Dt/e0))./(Sigmay_xA(:,1:NPML-1)*factor*Dy).*...

(EzxPMLA(:,2:NPML)+EzyPMLA(:,2:NPML)-EzxPMLA(:,1:NPML-1)-EzyPMLA(:,1:NPML-1));

HyPMLA(2:NTx-1,:)=HyPMLA(2:NTx-1,:)+...

Dt*(EzxPMLA(2:NTx-1,:)+EzyPMLA(2:NTx-1,:)-EzxPMLA(1:NTx-2,:)-EzyPMLA(1:NTx-2,:))/(mu0*Dx);

HyPMLA(1,:)=expboundary*HyPMLA(1,:)+(1-expboundary)*...

(EzxPMLA(1,:)+EzyPMLA(1,:)-EzxPML1(NPML,:)-EzyPML1(NPML,:))/(SigmaBound*factor*Dx); %Boundary Left

HyPMLA(NTx,:)=expboundary*HyPMLA(NTx,:)+(1-expboundary)*...

(EzxPML2(1,:)+EzyPML2(1,:)-EzxPMLA(NTx-1,:)-EzyPMLA(NTx-1,:))/(SigmaBound*factor*Dx); %Boundary Right

HxPMLA(:,NPML)=HxPMLA(:,NPML-1); %Boundary Upper;

%The following part is for ZONE B.

%H components.

HxPMLB(:,2:NPML)=exp(-Sigmay_xB(:,2:NPML)*Dt/e0).*HxPMLB(:,2:NPML)-...

(1-exp(-Sigmay_xB(:,2:NPML)*Dt/e0))./(Sigmay_xB(:,2:NPML)*factor*Dy).*...

(EzxPMLB(:,2:NPML)+EzyPMLB(:,2:NPML)-EzxPMLB(:,1:NPML-1)-EzyPMLB(:,1:NPML-1));

HyPMLB(2:NTx-1,:)=HyPMLB(2:NTx-1,:)+...

Dt*(EzxPMLB(2:NTx-1,:)+EzyPMLB(2:NTx-1,:)-EzxPMLB(1:NTx-2,:)-EzyPMLB(1:NTx-2,:))/(mu0*Dx);

HyPMLB(1,:)=expboundary*HyPMLB(1,:)+(1-expboundary)*...

(EzxPMLB(1,:)+EzyPMLB(1,:)-EzxPML3(NPML,:)-EzyPML3(NPML,:))/(SigmaBound*factor*Dx); %Boundary Left

HyPMLB(NTx,:)=expboundary*HyPMLB(NTx,:)+(1-expboundary)*...

(EzxPML4(1,:)+EzyPML4(1,:)-EzxPMLB(NTx-1,:)-EzyPMLB(NTx-1,:))/(SigmaBound*factor*Dx); %Boundary Right

HxPMLB(:,1)=HxPMLB(:,2); %Boundary Bottom;

%The following part is for ZONE C.

%H components.

HxPMLC(:,2:NTy-1)=HxPMLC(:,2:NTy-1)-...

Dt*(EzxPMLC(:,2:NTy-1)+EzyPMLC(:,2:NTy-1)-EzxPMLC(:,1:NTy-2)-EzyPMLC(:,1:NTy-2))/(mu0*Dy);

HyPMLC(2:NPML,:)=exp(-Sigmax_yC(2:NPML,:)*Dt/e0).*HyPMLC(2:NPML,:)+...

(1-exp(-Sigmax_yC(2:NPML,:)*Dt/e0))./(Sigmax_yC(2:NPML,:)*factor*Dx).*...

(EzxPMLC(2:NPML,:)+EzyPMLC(2:NPML,:)-EzxPMLC(1:NPML-1,:)-EzyPMLC(1:NPML-1,:));

HxPMLC(:,1)=expboundary*HxPMLC(:,1)-(1-expboundary)*...

(EzxPMLC(:,1)+EzyPMLC(:,1)-EzxPML3(:,NPML)-EzyPML3(:,NPML))/(SigmaBound*factor*Dy); %Boundary Down

HxPMLC(:,NTy)=expboundary*HxPMLC(:,NTy)-(1-expboundary)*...

(EzxPML1(:,1)+EzyPML1(:,1)-EzxPMLC(:,NTy-1)-EzyPMLC(:,NTy-1))/(SigmaBound*factor*Dy); %Boundary Upper

HyPMLC(1,:)=HyPMLC(2,:); %Boundary Left;

%The following part is for ZONE D.

%H components.

HxPMLD(:,2:NTy-1)=HxPMLD(:,2:NTy-1)-...

Dt*(EzxPMLD(:,2:NTy-1)+EzyPMLD(:,2:NTy-1)-EzxPMLD(:,1:NTy-2)-EzyPMLD(:,1:NTy-2))/(mu0*Dy);

HyPMLD(1:NPML-1,:)=exp(-Sigmax_yD(1:NPML-1,:)*Dt/e0).*HyPMLD(1:NPML-1,:)+...

(1-exp(-Sigmax_yD(1:NPML-1,:)*Dt/e0))./(Sigmax_yD(1:NPML-1,:)*factor*Dx).*...

(EzxPMLD(2:NPML,:)+EzyPMLD(2:NPML,:)-EzxPMLD(1:NPML-1,:)-EzyPMLD(1:NPML-1,:));

HxPMLD(:,1)=expboundary*HxPMLD(:,1)-(1-expboundary)*...

(EzxPMLD(:,1)+EzyPMLD(:,1)-EzxPML4(:,NPML)-EzyPML4(:,NPML))/(SigmaBound*factor*Dy); %Boundary Down

HxPMLD(:,NTy)=expboundary*HxPMLD(:,NTy)-(1-expboundary)*...

(EzxPML2(:,1)+EzyPML2(:,1)-EzxPMLD(:,NTy-1)-EzyPMLD(:,NTy-1))/(SigmaBound*factor*Dy); %Boundary Upper

HyPMLD(NPML,:)=HyPMLD(NPML-1,:); %Boundary Right;

%The following part is for ZONE 1.

%H components.

HxPML1(:,1:NPML-1)=exp(-Sigmay_x1(:,1:NPML-1)*Dt/e0).*HxPML1(:,1:NPML-1)-...

(1-exp(-Sigmay_x1(:,1:NPML-1)*Dt/e0))./(Sigmay_x1(:,1:NPML-1)*factor*Dy).*...

(EzxPML1(:,2:NPML)+EzyPML1(:,2:NPML)-EzxPML1(:,1:NPML-1)-EzyPML1(:,1:NPML-1));

HxPML1(:,NPML)=HxPML1(:,NPML-1);

HyPML1(2:NPML,:)=exp(-Sigmax_y1(2:NPML,:)*Dt/e0).*HyPML1(2:NPML,:)+...

(1-exp(-Sigmax_y1(2:NPML,:)*Dt/e0))./(Sigmax_y1(2:NPML,:)*factor*Dx).*...

(EzxPML1(2:NPML,:)+EzyPML1(2:NPML,:)-EzxPML1(1:NPML-1,:)-EzyPML1(1:NPML-1,:));

HyPML1(1,:)=HyPML1(2,:);

%The following part is for ZONE 2.

%H components.

HxPML2(:,1:NPML-1)=exp(-Sigmay_x2(:,1:NPML-1)*Dt/e0).*HxPML2(:,1:NPML-1)-...

(1-exp(-Sigmay_x2(:,1:NPML-1)*Dt/e0))./(Sigmay_x2(:,1:NPML-1)*factor*Dy).*...

(EzxPML2(:,2:NPML)+EzyPML2(:,2:NPML)-EzxPML2(:,1:NPML-1)-EzyPML2(:,1:NPML-1));

HxPML2(:,NPML)=HxPML2(:,NPML-1);

HyPML2(1:NPML-1,:)=exp(-Sigmax_y2(1:NPML-1,:)*Dt/e0).*HyPML2(1:NPML-1,:)+...

(1-exp(-Sigmax_y2(1:NPML-1,:)*Dt/e0))./(Sigmax_y2(1:NPML-1,:)*factor*Dx).*...

(EzxPML2(2:NPML,:)+EzyPML2(2:NPML,:)-EzxPML2(1:NPML-1,:)-EzyPML2(1:NPML-1,:));

HyPML2(NPML,:)=HyPML2(NPML-1,:); %Boundary Right;

%The following part is for ZONE 3.

%H components.

HyPML3(2:NPML,:)=exp(-Sigmax_y3(2:NPML,:)*Dt/e0).*HyPML3(2:NPML,:)+...

(1-exp(-Sigmax_y3(2:NPML,:)*Dt/e0))./(Sigmax_y3(2:NPML,:)*factor*Dx).*...

(EzxPML3(2:NPML,:)+EzyPML3(2:NPML,:)-EzxPML3(1:NPML-1,:)-EzyPML3(1:NPML-1,:));

HyPML3(1,:)=HyPML3(2,:); %Boundary Left;

HxPML3(:,2:NPML)=exp(-Sigmay_x3(:,2:NPML)*Dt/e0).*HxPML3(:,2:NPML)-...

(1-exp(-Sigmay_x3(:,2:NPML)*Dt/e0))./(Sigmay_x3(:,2:NPML)*factor*Dy).*...

(EzxPML3(:,2:NPML)+EzyPML3(:,2:NPML)-EzxPML3(:,1:NPML-1)-EzyPML3(:,1:NPML-1));

HxPML3(:,1)=HxPML3(:,2); %Boundary Bottom;

%The following part is for ZONE 4.

%H components.

HxPML4(:,2:NPML)=exp(-Sigmay_x4(:,2:NPML)*Dt/e0).*HxPML4(:,2:NPML)-...

(1-exp(-Sigmay_x4(:,2:NPML)*Dt/e0))./(Sigmay_x4(:,2:NPML)*factor*Dy).*...

(EzxPML4(:,2:NPML)+EzyPML4(:,2:NPML)-EzxPML4(:,1:NPML-1)-EzyPML4(:,1:NPML-1));

HxPML4(:,1)=HxPML4(:,2); %Boundary Bottom;

HyPML4(1:NPML-1,:)=exp(-Sigmax_y4(1:NPML-1,:)*Dt/e0).*HyPML4(1:NPML-1,:)+...

(1-exp(-Sigmax_y4(1:NPML-1,:)*Dt/e0))./(Sigmax_y4(1:NPML-1,:)*factor*Dx).*...

(EzxPML4(2:NPML,:)+EzyPML4(2:NPML,:)-EzxPML4(1:NPML-1,:)-EzyPML4(1:NPML-1,:));

HyPML4(NPML,:)=HyPML4(NPML-1,:); %Boundary Right;

%The following part is for inside region.

%Ez component.

Ez=Ez+Dt*((Hy(2:NTx,1:NTy-1)-Hy(1:NTx-1,1:NTy-1))/Dx-...

(Hx(1:NTx-1,2:NTy)-Hx(1:NTx-1,1:NTy-1))/Dy)./Ep; %Inside region.

%The following part is for source

Ez(1,5*NMlat+1:5*NMlat+21)=Ez(1,5*NMlat+1:5*NMlat+21)+sin(W*m*Dt);%*exp(-(m*W*Dt-3)^2);

%End of Source

%The following part is for ZONE A.

%Ez component.

EzxPMLA=EzxPMLA+Dt*(HyPMLA(2:NTx,:)-HyPMLA(1:NTx-1,:))/(e0*Dx);

EzyPMLA(:,2:NPML)=exp(-Sigmay_zA(:,2:NPML)*Dt/e0).*EzyPMLA(:,2:NPML)-...

(1-exp(-Sigmay_zA(:,2:NPML)*Dt/e0))./(Sigmay_zA(:,2:NPML)*Dy).*...

(HxPMLA(:,2:NPML)-HxPMLA(:,1:NPML-1));

EzyPMLA(:,1)=exp(-Sigmay_zA(:,1)*Dt/e0).*EzyPMLA(:,1)-...

(1-exp(-Sigmay_zA(:,1)*Dt/e0))./(Sigmay_zA(:,1)*Dy).*...

(HxPMLA(:,1)-Hx(:,NTy));

%The following part is for ZONE B.

%Ez component.

EzxPMLB=EzxPMLB+Dt*(HyPMLB(2:NTx,:)-HyPMLB(1:NTx-1,:))/(e0*Dx);

EzyPMLB(:,1:NPML-1)=exp(-Sigmay_zB(:,1:NPML-1)*Dt/e0).*EzyPMLB(:,1:NPML-1)-...

(1-exp(-Sigmay_zB(:,1:NPML-1)*Dt/e0))./(Sigmay_zB(:,1:NPML-1)*Dy).*...

(HxPMLB(:,2:NPML)-HxPMLB(:,1:NPML-1));

EzyPMLB(:,NPML)=exp(-Sigmay_zB(:,NPML)*Dt/e0).*EzyPMLB(:,NPML)-...

(1-exp(-Sigmay_zB(:,NPML)*Dt/e0))./(Sigmay_zB(:,NPML)*Dy).*...

(Hx(:,1)-HxPMLB(:,NPML));

%The following part is for ZONE C.

%Ez component.

EzxPMLC(1:NPML-1,:)=exp(-Sigmax_zC(1:NPML-1,:)*Dt/e0).*EzxPMLC(1:NPML-1,:)+...

(1-exp(-Sigmax_zC(1:NPML-1,:)*Dt/e0))./(Sigmax_zC(1:NPML-1,:)*Dx).*...

(HyPMLC(2:NPML,:)-HyPMLC(1:NPML-1,:));

EzxPMLC(NPML,:)=exp(-Sigmax_zC(NPML,:)*Dt/e0).*EzxPMLC(NPML,:)+...

(1-exp(-Sigmax_zC(NPML,:)*Dt/e0))./(Sigmax_zC(NPML,:)*Dx).*...

(Hy(1,:)-HyPMLC(NPML,:));

EzyPMLC=EzyPMLC-Dt*(HxPMLC(:,2:NTy)-HxPMLC(:,1:NTy-1))/(e0*Dy);

%The following part is for ZONE D.

%Ez component.

EzxPMLD(2:NPML,:)=exp(-Sigmax_zD(2:NPML,:)*Dt/e0).*EzxPMLD(2:NPML,:)+...

(1-exp(-Sigmax_zD(2:NPML,:)*Dt/e0))./(Sigmax_zD(2:NPML,:)*Dx).*...

(HyPMLD(2:NPML,:)-HyPMLD(1:NPML-1,:));

EzxPMLD(1,:)=exp(-Sigmax_zD(1,:)*Dt/e0).*EzxPMLD(1,:)+...

(1-exp(-Sigmax_zD(1,:)*Dt/e0))./(Sigmax_zD(1,:)*Dx).*...

(HyPMLD(1,:)-Hy(NTx,:));

EzyPMLD=EzyPMLD-Dt*(HxPMLD(:,2:NTy)-HxPMLD(:,1:NTy-1))/(e0*Dy);

%The following part is for ZONE 1.

%Ez component.

EzxPML1(1:NPML-1,:)=exp(-Sigmax_z1(1:NPML-1,:)*Dt/e0).*EzxPML1(1:NPML-1,:)+...

(1-exp(-Sigmax_z1(1:NPML-1,:)*Dt/e0))./(Sigmax_z1(1:NPML-1,:)*Dx).*...

(HyPML1(2:NPML,:)-HyPML1(1:NPML-1,:));

EzxPML1(NPML,:)=exp(-Sigmax_z1(NPML,:)*Dt/e0).*EzxPML1(NPML,:)+...

(1-exp(-Sigmax_z1(NPML,:)*Dt/e0))./(Sigmax_z1(NPML,:)*Dx).*...

(HyPMLA(1,:)-HyPML1(NPML,:));

EzyPML1(:,2:NPML)=exp(-Sigmay_z1(:,2:NPML)*Dt/e0).*EzyPML1(:,2:NPML)-...

(1-exp(-Sigmay_z1(:,2:NPML)*Dt/e0))./(Sigmay_z1(:,2:NPML)*Dy).*...

(HxPML1(:,2:NPML)-HxPML1(:,1:NPML-1));

EzyPML1(:,1)=exp(-Sigmay_z1(:,1)*Dt/e0).*EzyPML1(:,1)-...

(1-exp(-Sigmay_z1(:,1)*Dt/e0))./(Sigmay_z1(:,1)*Dy).*...

(HxPML1(:,1)-HxPMLC(:,NTy));

%The following part is for ZONE 2.

%Ez component.

EzxPML2(2:NPML,:)=exp(-Sigmax_z2(2:NPML,:)*Dt/e0).*EzxPML2(2:NPML,:)+...

(1-exp(-Sigmax_z2(2:NPML,:)*Dt/e0))./(Sigmax_z2(2:NPML,:)*Dx).*...

(HyPML2(2:NPML,:)-HyPML2(1:NPML-1,:));

EzxPML2(1,:)=exp(-Sigmax_z2(1,:)*Dt/e0).*EzxPML2(1,:)+...

(1-exp(-Sigmax_z2(1,:)*Dt/e0))./(Sigmax_z2(1,:)*Dx).*...

(HyPML2(1,:)-HyPMLA(NTx,:));

EzyPML2(:,2:NPML)=exp(-Sigmay_z2(:,2:NPML)*Dt/e0).*EzyPML2(:,2:NPML)-...

(1-exp(-Sigmay_z2(:,2:NPML)*Dt/e0))./(Sigmay_z2(:,2:NPML)*Dy).*...

(HxPML2(:,2:NPML)-HxPML2(:,1:NPML-1));

EzyPML2(:,1)=exp(-Sigmay_z2(:,1)*Dt/e0).*EzyPML2(:,1)-...

(1-exp(-Sigmay_z2(:,1)*Dt/e0))./(Sigmay_z2(:,1)*Dy).*...

(HxPML2(:,1)-HxPMLD(:,NTy));

%The following part is for ZONE 3.

%Ez component.

EzxPML3(1:NPML-1,:)=exp(-Sigmax_z3(1:NPML-1,:)*Dt/e0).*EzxPML3(1:NPML-1,:)+...

(1-exp(-Sigmax_z3(1:NPML-1,:)*Dt/e0))./(Sigmax_z3(1:NPML-1,:)*Dx).*...

(HyPML3(2:NPML,:)-HyPML3(1:NPML-1,:));

EzxPML3(NPML,:)=exp(-Sigmax_z3(NPML,:)*Dt/e0).*EzxPML3(NPML,:)+...

(1-exp(-Sigmax_z3(NPML,:)*Dt/e0))./(Sigmax_z3(NPML,:)*Dx).*...

(HyPMLB(1,:)-HyPML3(NPML,:));

EzyPML3(:,1:NPML-1)=exp(-Sigmay_z3(:,1:NPML-1)*Dt/e0).*EzyPML3(:,1:NPML-1)-...

(1-exp(-Sigmay_z3(:,1:NPML-1)*Dt/e0))./(Sigmay_z3(:,1:NPML-1)*Dy).*...

(HxPML3(:,2:NPML)-HxPML3(:,1:NPML-1));

EzyPML3(:,NPML)=exp(-Sigmay_z3(:,NPML)*Dt/e0).*EzyPML3(:,NPML)-...

(1-exp(-Sigmay_z3(:,NPML)*Dt/e0))./(Sigmay_z3(:,NPML)*Dy).*...

(HxPMLC(:,1)-HxPML3(:,NPML));

%The following part is for ZONE 4.

%Ez component.

EzxPML4(2:NPML,:)=exp(-Sigmax_z4(2:NPML,:)*Dt/e0).*EzxPML4(2:NPML,:)+...

(1-exp(-Sigmax_z4(2:NPML,:)*Dt/e0))./(Sigmax_z4(2:NPML,:)*Dx).*...

(HyPML4(2:NPML,:)-HyPML4(1:NPML-1,:));

EzxPML4(1,:)=exp(-Sigmax_z4(1,:)*Dt/e0).*EzxPML4(1,:)+...

(1-exp(-Sigmax_z4(1,:)*Dt/e0))./(Sigmax_z4(1,:)*Dx).*...

(HyPML4(1,:)-HyPMLB(NTx,:));

EzyPML4(:,1:NPML-1)=exp(-Sigmay_z4(:,1:NPML-1)*Dt/e0).*EzyPML4(:,1:NPML-1)-...

(1-exp(-Sigmay_z4(:,1:NPML-1)*Dt/e0))./(Sigmay_z4(:,1:NPML-1)*Dy).*...

(HxPML4(:,2:NPML)-HxPML4(:,1:NPML-1));

EzyPML4(:,NPML)=exp(-Sigmay_z4(:,NPML)*Dt/e0).*EzyPML4(:,NPML)-...

(1-exp(-Sigmay_z4(:,NPML)*Dt/e0))./(Sigmay_z4(:,NPML)*Dy).*...

(HxPMLD(:,1)-HxPML4(:,NPML));

if mod(m,Meach)==0

if IsFigure==1

%The following part is for visual figures.

EzAll=[EzxPML3+EzyPML3,EzxPMLC+EzyPMLC,EzxPML1+EzyPML1;

EzxPMLB+EzyPMLB,Ez,EzxPMLA+EzyPMLA;

EzxPML4+EzyPML4,EzxPMLD+EzyPMLD,EzxPML2+EzyPML2];

HxAll=[HxPML3,HxPMLC,HxPML1;HxPMLB,Hx,HxPMLA;HxPML4,HxPMLD,HxPML2];

HyAll=[HyPML3,HyPMLC,HyPML1;HyPMLB,Hy,HyPMLA;HyPML4,HyPMLD,HyPML2];

figure(1);

clf

surf(X,Y,real(EzAll));

caxis([-Colormax Colormax]);

shading interp;

axis([-0.5*a*MLatx-Dx*NPML 0.5*a*MLatx+Dx*NPML -0.5*a*MLaty-Dx*NPML 0.5*a*MLaty+Dx*NPML])

zlabel('Ez');

if IsMovie==1

Mnum=Mnum+1;

Movie(:,Mnum)=get ;

end

view(0,90);

axis off;

pause(0.2)

end

%toc

end

toc

文件列表（点击上边下载按钮，如果是垃圾文件请在下面评价差评或者投诉）:

失效有关曲率处理
失效有关曲率处理\Face Segmentation
失效有关曲率处理\Face Segmentation\cropped images
失效有关曲率处理\Face Segmentation\cropped images\forehead
失效有关曲率处理\Face Segmentation\cropped images\left cheek
失效有关曲率处理\Face Segmentation\cropped images\left eye
失效有关曲率处理\Face Segmentation\cropped images\lips
失效有关曲率处理\Face Segmentation\cropped images\nose
失效有关曲率处理\Face Segmentation\cropped images\right cheek
失效有关曲率处理\Face Segmentation\cropped images\right eye
失效有关曲率处理\Face Segmentation\cropped_images
失效有关曲率处理\Face Segmentation\fuzzy
失效有关曲率处理\Face Segmentation\hk
失效有关曲率处理\Face Segmentation\hk\h
失效有关曲率处理\Face Segmentation\hk\k
失效有关曲率处理\Face Segmentation\hk\k0
失效有关曲率处理\Face Segmentation\hk\k135
失效有关曲率处理\Face Segmentation\hk\k45
失效有关曲率处理\Face Segmentation\hk\k90
失效有关曲率处理\Face Segmentation\photos
失效有关曲率处理\Face Segmentation\results
失效有关曲率处理\Face Segmentation.zip
失效有关曲率处理\Face Segmentation\code.m
失效有关曲率处理\Face Segmentation\cropped images\forehead\forehead1.jpeg
失效有关曲率处理\Face Segmentation\cropped images\forehead\forehead2.jpeg
失效有关曲率处理\Face Segmentation\cropped images\forehead\forehead3.jpeg
失效有关曲率处理\Face Segmentation\cropped images\forehead\forehead4.jpeg
失效有关曲率处理\Face Segmentation\cropped images\forehead\forehead5.jpeg
失效有关曲率处理\Face Segmentation\cropped images\left cheek\left cheek1.jpeg
失效有关曲率处理\Face Segmentation\cropped images\left cheek\left cheek2.jpeg
失效有关曲率处理\Face Segmentation\cropped images\left cheek\left cheek3.jpeg
失效有关曲率处理\Face Segmentation\cropped images\left cheek\left cheek4.jpeg
失效有关曲率处理\Face Segmentation\cropped images\left cheek\left cheek5.jpeg
失效有关曲率处理\Face Segmentation\cropped images\left eye\left eye1.jpeg
失效有关曲率处理\Face Segmentation\cropped images\left eye\left eye2.jpeg
失效有关曲率处理\Face Segmentation\cropped images\left eye\left eye3.jpeg
失效有关曲率处理\Face Segmentation\cropped images\left eye\left eye4.jpeg
失效有关曲率处理\Face Segmentation\cropped images\left eye\left eye5.jpeg
失效有关曲率处理\Face Segmentation\cropped images\lips\lips1.jpeg
失效有关曲率处理\Face Segmentation\cropped images\lips\lips2.jpeg
失效有关曲率处理\Face Segmentation\cropped images\lips\lips3.jpeg
失效有关曲率处理\Face Segmentation\cropped images\lips\lips4.jpeg
失效有关曲率处理\Face Segmentation\cropped images\lips\lips5.jpeg
失效有关曲率处理\Face Segmentation\cropped images\nose\nose1.jpeg
失效有关曲率处理\Face Segmentation\cropped images\nose\nose2.jpeg
失效有关曲率处理\Face Segmentation\cropped images\nose\nose3.jpeg
失效有关曲率处理\Face Segmentation\cropped images\nose\nose4.jpeg
失效有关曲率处理\Face Segmentation\cropped images\nose\nose5.jpeg
失效有关曲率处理\Face Segmentation\cropped images\right cheek\right cheek1.jpeg
失效有关曲率处理\Face Segmentation\cropped images\right cheek\right cheek2.jpeg
失效有关曲率处理\Face Segmentation\cropped images\right cheek\right cheek3.jpeg
失效有关曲率处理\Face Segmentation\cropped images\right cheek\right cheek4.jpeg
失效有关曲率处理\Face Segmentation\cropped images\right cheek\right cheek5.jpeg
失效有关曲率处理\Face Segmentation\cropped images\right eye\right eye1.jpeg
失效有关曲率处理\Face Segmentation\cropped images\right eye\right eye2.jpeg
失效有关曲率处理\Face Segmentation\cropped images\right eye\right eye3.jpeg
失效有关曲率处理\Face Segmentation\cropped images\right eye\right eye4.jpeg
失效有关曲率处理\Face Segmentation\cropped images\right eye\right eye5.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\forehead1.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\forehead2.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\forehead3.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\forehead4.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\forehead5.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\left cheek1.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\left cheek2.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\left cheek3.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\left cheek4.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\left cheek5.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\left eye1.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\left eye2.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\left eye3.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\left eye4.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\left eye5.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\lips1.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\lips2.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\lips3.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\lips4.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\lips5.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\nose1.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\nose2.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\nose3.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\nose4.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\nose5.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\right cheek1.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\right cheek2.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\right cheek3.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\right cheek4.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\right cheek5.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\right eye1.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\right eye2.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\right eye3.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\right eye4.jpeg
失效有关曲率处理\Face Segmentation\cropped_images\right eye5.jpeg
失效有关曲率处理\Face Segmentation\fuzzy\FUZZY.fis
失效有关曲率处理\Face Segmentation\hk\h\h1.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h10.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h11.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h12.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h13.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h14.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h15.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h16.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h17.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h18.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h19.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h2.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h20.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h21.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h22.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h23.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h24.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h25.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h26.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h27.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h28.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h29.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h3.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h30.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h31.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h32.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h33.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h34.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h35.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h4.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h5.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h6.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h7.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h8.jpeg
失效有关曲率处理\Face Segmentation\hk\h\h9.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k01.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k010.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k011.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k012.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k013.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k014.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k015.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k016.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k017.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k018.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k019.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k02.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k020.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k021.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k022.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k023.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k024.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k025.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k026.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k027.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k028.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k029.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k03.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k030.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k031.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k032.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k033.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k034.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k035.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k04.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k05.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k06.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k07.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k08.jpeg
失效有关曲率处理\Face Segmentation\hk\k0\k09.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k1351.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13510.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13511.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13512.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13513.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13514.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13515.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13516.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13517.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13518.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13519.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k1352.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13520.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13521.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13522.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13523.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13524.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13525.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13526.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13527.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13528.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13529.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k1353.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13530.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13531.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13532.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13533.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13534.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k13535.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k1354.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k1355.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k1356.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k1357.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k1358.jpeg
失效有关曲率处理\Face Segmentation\hk\k135\k1359.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k451.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4510.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4511.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4512.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4513.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4514.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4515.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4516.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4517.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4518.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4519.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k452.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4520.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4521.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4522.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4523.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4524.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4525.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4526.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4527.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4528.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4529.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k453.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4530.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4531.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4532.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4533.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4534.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k4535.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k454.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k455.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k456.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k457.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k458.jpeg
失效有关曲率处理\Face Segmentation\hk\k45\k459.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k901.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9010.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9011.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9012.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9013.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9014.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9015.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9016.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9017.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9018.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9019.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k902.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9020.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9021.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9022.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9023.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9024.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9025.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9026.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9027.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9028.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9029.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k903.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9030.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9031.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9032.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9033.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9034.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k9035.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k904.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k905.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k906.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k907.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k908.jpeg
失效有关曲率处理\Face Segmentation\hk\k90\k909.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k1.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k10.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k11.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k12.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k13.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k14.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k15.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k16.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k17.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k18.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k19.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k2.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k20.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k21.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k22.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k23.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k24.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k25.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k26.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k27.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k28.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k29.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k3.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k30.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k31.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k32.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k33.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k34.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k35.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k4.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k5.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k6.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k7.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k8.jpeg
失效有关曲率处理\Face Segmentation\hk\k\k9.jpeg
失效有关曲率处理\Face Segmentation\loop.asv
失效有关曲率处理\Face Segmentation\loop.m
失效有关曲率处理\Face Segmentation\photos\c1.jpeg
失效有关曲率处理\Face Segmentation\photos\c2.jpeg
失效有关曲率处理\Face Segmentation\photos\c3.jpeg
失效有关曲率处理\Face Segmentation\photos\c4.jpeg
失效有关曲率处理\Face Segmentation\photos\c5.jpeg
失效有关曲率处理\Face Segmentation\project presentation.pptx
失效有关曲率处理\Face Segmentation\read me.txt
失效有关曲率处理\Face Segmentation\results\forehead1.jpg
失效有关曲率处理\Face Segmentation\results\forehead2.jpg
失效有关曲率处理\Face Segmentation\results\forehead3.jpg
失效有关曲率处理\Face Segmentation\results\forehead4.jpg
失效有关曲率处理\Face Segmentation\results\forehead5.jpg
失效有关曲率处理\Face Segmentation\results\left cheek1.jpg
失效有关曲率处理\Face Segmentation\results\left cheek2.jpg
失效有关曲率处理\Face Segmentation\results\left cheek3.jpg
失效有关曲率处理\Face Segmentation\results\left cheek4.jpg
失效有关曲率处理\Face Segmentation\results\left cheek5.jpg
失效有关曲率处理\Face Segmentation\results\left eye1.jpg
失效有关曲率处理\Face Segmentation\results\left eye2.jpg
失效有关曲率处理\Face Segmentation\results\left eye3.jpg
失效有关曲率处理\Face Segmentation\results\left eye4.jpg
失效有关曲率处理\Face Segmentation\results\left eye5.jpg
失效有关曲率处理\Face Segmentation\results\lips1.jpg
失效有关曲率处理\Face Segmentation\results\lips2.jpg
失效有关曲率处理\Face Segmentation\results\lips3.jpg
失效有关曲率处理\Face Segmentation\results\lips4.jpg
失效有关曲率处理\Face Segmentation\results\lips5.jpg
失效有关曲率处理\Face Segmentation\results\nose1.jpg
失效有关曲率处理\Face Segmentation\results\nose2.jpg
失效有关曲率处理\Face Segmentation\results\nose3.jpg
失效有关曲率处理\Face Segmentation\results\nose4.jpg
失效有关曲率处理\Face Segmentation\results\nose5.jpg
失效有关曲率处理\Face Segmentation\results\right cheek1.jpg
失效有关曲率处理\Face Segmentation\results\right cheek2.jpg
失效有关曲率处理\Face Segmentation\results\right cheek3.jpg
失效有关曲率处理\Face Segmentation\results\right cheek4.jpg
失效有关曲率处理\Face Segmentation\results\right cheek5.jpg
失效有关曲率处理\Face Segmentation\results\right eye1.jpg
失效有关曲率处理\Face Segmentation\results\right eye2.jpg
失效有关曲率处理\Face Segmentation\results\right eye3.jpg
失效有关曲率处理\Face Segmentation\results\right eye4.jpg
失效有关曲率处理\Face Segmentation\results\right eye5.jpg
失效有关曲率处理\IEEEFuzzyJournal_11July.pdf