stochastic_resonance.m

来自「随机共振相似度的Matlab计算程序」· M 代码 · 共 26 行

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

for jj=1:100
		for ii=1:1000
            x(ii)=s+sigma(jj)*n(ii);
            
            if x(ii)>theta %传递函数
                y(ii)=1;
            else
                y(ii)=0;
            end
		end

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

plot(sigma,S)
hold on