clear
clc
theta=5;%1.1,2.1,5
s=1;
sigma=linspace(0,10,100);
randn('state',0)
n=randn(1,1000);%正態(tài)分布的噪音
x=zeros(1,1000);

for jj=1:100
		for ii=1:1000
            x(ii)=s+sigma(jj)*n(ii);
            
            if x(ii)>theta %傳遞函數(shù)
                y(ii)=1;
            else
                y(ii)=0;
            end
		end

S(jj)=sum(y==s)/1000*100%;/100;
end;

plot(sigma,S)
hold on