基于matlab的阀值系统随机共振现象的模拟仿真,含仿真操作录像

clc; clear; close all; warning off; addpath(genpath(pwd)); rng('default') n=10000; theta=1.1;%阀值 s=randn(1,n); for i=1:n if(s(i)>=0) s(i)=1; else s(i)=0; end end for kk=1:101 %产生高斯噪声 p=(kk-1)/10; P(kk)=p; y=randn(1,n); y=y/std(y); y=y-mean(y); b=sqrt(p); y=b*y; sadd=s+y; for i=1:n if(sadd(i)>=theta) Y(i)=1; else Y(i)=0; end end S=0; %相似度 for i=1:n if(Y(i)==s(i)) S=S+1; end end Se(kk)=S/n; end theta=2.1;%阀值 for kk=1:101 %产生高斯噪声 p=(kk-1)/10; P(kk)=p; y=randn(1,n); y=y/std(y); y=y-mean(y); b=sqrt(p); y=b*y; sadd=s+y; for i=1:n if(sadd(i)>=theta) Y(i)=1; else Y(i)=0; end end S=0; %相似度 for i=1:n if(Y(i)==s(i)) S=S+1; end end Se2(kk)=S/n; end theta=5;%阀值 for kk=1:101 %产生高斯噪声 p=(kk-1)/10; P(kk)=p; y=randn(1,n); y=y/std(y); y=y-mean(y); b=sqrt(p); y=b*y; sadd=s+y; for i=1:n if(sadd(i)>=theta) Y(i)=1; else Y(i)=0; end end S=0; %相似度 for i=1:n if(Y(i)==s(i)) S=S+1; end end Se3(kk)=S/n; end figure(1); plot(P,Se);grid on;xlabel('噪声强度D');ylabel('输入－输出信号相似度S');title('阀值系统中的随机共振现象'); hold on; plot(P,Se2,'red'); plot(P,Se3,'green');