a=load( q.txt )

width=input( 请输入窗宽 )

step=input( 请输入步长 )

m=length(a)

for i=1:floor((m-width)/step+1)

s(i,1:width)=a(1,1+(i-1)*step:width+(i-1)*step)

end

b=max(max(s))

M=s./b

Rs=input( 请输入本体集半径 )

nd=input( 请输入检测器数模 )

k=1

while k<=nd

d{k}=rand(1,width)

for i=1:floor((m-width)/step+1)

if (norm(d{k}-s(i,:)))^2<Rs

break

end

end

k=k+1

end

for j=1:nd

for i=1:floor((m-width)/step+1)

R(i)=(norm(d{j}-s(i,:)))^2-a = load (' q.txt' ) width = input (' Please enter the window width' ) step = input (' Enter step size' ) m = length (a) for i = 1: floor ((m-width)/step+1) s (i, 1: width) = a (1,1+ (i-1)* step: width+ (i-1)* step)　end b = max (max (s)) M = s./b Rs = input (' Please enter the body set radius' ) nd = input (' Please enter the number of detector modules' ) k = 1 while k < = nd d {k} = rand (1, width) for i = 1: floor ((m-width)/step+1) if (norm (d {k}-s (i ,:)))^ 2 <Rs

　　　　　　break

　　　　end

　　end

　　k=k+1

end

for j=1:nd

　　 for i=1:floor((m-width)/step+1)

　　R(i)=(norm(d{j}-s(i,:)))^2