Question

matlab人眼虹膜识别程序调试

Y =imread('test1.jpg');
Y1=rgb2gray(Y);
I= im2bw(Y1);
[sy,sx]=size(I);
[y,x]=find(I);
totalpix = length(x);
HM_tmp = zeros(sy*sx,1);
b=1:sy;
a=zeros(sy,totalpix);
if nargin==1
R_min=1;
R_max=max(max(x),max(y));
else
R_min=Rp(1);
R_max=Rp(2);
end
y=repmat(y',[sy,1]);
x=repmat(x',[sy,1]);
HPN=0;
for R=R_min:R_max
R2=R^2;
b1=repmat(b',[1,totalopix]);
b2=b1;
a1=(round(x - sqrt(R2-(y-b1).^2)));
a2=(round(x + sqrt(R2-(y-b2).^2)));
b1=b1(imag(a1)==0 & a1>0 & a1<sx);
a1=a1(imag(a1)==0 & a1>0 & a1<sx);
b2=b2(imag(a2)==0 & a2>0 & a2<sx);
a2=a2(imag(a2)==0 & a2>0 & a2<sx);
ind1=sub2ind([sy,sx],b1,a1);
ind2=sub2ind([sy,sx],b2,a2);
ind=[ind1;ind2];
val=ones(length(ind),1);
data=accumarray(ind,val);
HM_tmp(1:length(data))=data;
HM2_tmp=reshape(HM_tmp,[sy,sx]);
imshow(HM2,[]);
maxval=max(max(HM2_tmp));
if maxval>HPN
HPN=maxval;
HM=HM2_tmp;
Rc=R;
end
end
[B,A]=find(HM==HPN);
C=[mean(A),mean(B),Rc];

if nargin==1这一行报错，错误使用 nargin
您只能从 MATLAB 函数中调用 nargin/nargout。
出错 Untitled (line 10)
if nargin ==1

Answer 1

pdf函数 function a=pdf_8psk(x,y,SNR) a=1/16/pi*(exp(-((x-sqrt(SNR))).^2/2)+exp(-(x+sqrt(SNR)).^2/2)+exp(-(y-sqrt(SNR)).^2/2)+exp(-(y+sqrt(SNR)).^2/2)+exp(-(x-sqrt(SNR/2)).^2/2-(y-sqrt(SNR/2)).^2/2)... +exp(-(x+sqrt(SNR/2)).^2/2-(y-sqrt(SNR/2)).^2/2)+exp(-(x-sqrt(SNR/2)).^2/2-(y+sqrt(SNR/2)).^2/2)+exp(-(x+sqrt(SNR/2)).^2/2-(y+sqrt(SNR/2)).^2/2));