平行因子分析法

function [Rayleigh]=rzs(X,a,b,c,d,e,f) if nargin==0 disp(' ') disp(' ') disp(' The Primary Rayleigh Zero Setting') disp(' ') disp(' Type <<help rzs>> for more info') disp(' ') disp(' [Rayleigh] = rzs(X,a,b,c,d,e,f);') disp(' Em=[a:b:c]') disp(' Ex=[d:e:f]') disp(' ') disp(' Example: ') disp(' Rayleigh = rzs(X,300,1,700,250,5,600);') disp(' ') disp(' ') disp(' ') elseif nargin<7 error('Input X,Em,Ex must be given'); end if nargin==7 Rayleigh=X; if a>d if (a-d)>10 for i=ceil((a-d-10)/e)+1:ceil((a-d+10)/e); for j=1:fix(((i-(ceil((a-d-10)/e)+1))*e+(ceil((a-d-10)/e)*e-(a-d-10)))/b)+1; Rayleigh(i,j)=0; end end if c>f for i=ceil((a-d+10)/e)+1:fix((f-d)/e)+1 for j=ceil(((i-(ceil((a-d+10)/e)+1))*e+(ceil((a-d-10)/e)*e-(a-d-10)))/b)+1:ceil(((i-(ceil((a-d+10)/e)+1))*e+(ceil((a-d-10)/e)*e-(a-d-10))+21)/b) Rayleigh(i,j)=0; end end elseif c<=f for i=ceil((a-d+10)/e)+1:fix((c-d-10)/e)+1 for j=ceil(((i-(ceil((a-d+10)/e)+1))*e+(ceil((a-d-10)/e)*e-(a-d-10)))/b)+1:ceil(((i-(ceil((a-d+10)/e)+1))*e+(ceil((a-d-10)/e)*e-(a-d-10))+21)/b) Rayleigh(i,j)=0; end end end end if (a-d)<=10 for i=1:ceil((a-d+10)/e); for j=1:fix(((i-1)*e+d+10-a)/b)+1; Rayleigh(i,j)=0; end end if c>f for i=ceil((a-d+10)/e)+1:fix((f-d)/e)+1 for j=ceil(((i-(ceil((a-d+10)/e)+1))*e+d+10-a)/b)+1:ceil(((i-(ceil((a-d+10)/e)+1))*e+d+10-a+21)/b) Rayleigh(i,j)=0; end end elseif c<=f for i=ceil((a-d+10)/e)+1:fix((c-d-10)/e)+1 for j=ceil(((i-(ceil((a-d+10)/e)+1))*e+d+10-a)/b)+1:ceil(((i-(ceil((a-d+10)/e)+1))*e+d+10-a+21)/b) Rayleigh(i,j)=0; end end end end end if a<d if c>f for i=1:fix((f-d)/e)+1 for j=ceil(((i-1)*e+d-a-10)/b)+1:ceil(((i-1)*e+d-a-10+21)/b) Rayleigh(i,j)=0; end end elseif c<=f for i=1:fix((c-d-10)/e)+1 for j=ceil(((i-1)*e+d-a-10)/b)+1:ceil(((i-1)*e+d-a-10+21)/b) Rayleigh(i,j)=0; end end end end if a==d for i=1:fix(10/e)+1 for j=1:fix(((i-1)*e+10)/b)+1 Rayleigh(i,j)=0; end end if c>f for i=ceil(10/e)+1:fix((f-d)/e)+1 for j=ceil(((i-(ceil(10/e)+1))*e)/b)+1:ceil(((i-(ceil(10/e)+1))*e+21)/b) Rayleigh(i,j)=0; end end elseif c<=f for i=ceil(10/e)+1:fix((c-d-10)/e)+1 for j=ceil(((i-(ceil(10/e)+1))*e+(ceil(10/e)*e-10))/b)+1:ceil(((i-(ceil(10/e)+1))*e+(ceil(10/e)*e-10+21))/b) Rayleigh(i,j)=0; end end end end end end