新建 360压缩 ZIP 文件.zip_zip_共振_随机共振

function x=sr(a,b,h,x1)

x=zeros(1,length(x1));

% for i=1:length(x1)-1

% k1=h*(a*x(i)-b*x(i).^3+x1(i));

% k2=h*(a*(x(i)+k1/2)-b*(x(i)+k1/2).^3+x1(i));

% k3=h*(a*(x(i)+k2/2)-b*(x(i)+k1*(sqrt(2)-1)/2+k2*(2-sqrt(2))/2).^3+x1(i+1));

% k4=h*(a*(x(i)+k3)-b*(x(i)-k2*sqrt(2)/2+k3*(2+sqrt(2))/2).^3+x1(i+1));

% x(i+1)=x(i)+(1/6)*(k1+(2-sqrt(2))*k2+(2+sqrt(2))*k3+k4);

% end

End

clear all

clc

A=0.3; %信号幅度

fs=10; %采样频率

f=0.01; %信号频率

Ts=1/fs;%采样时间

h=1/fs; %时间步长

t=0:Ts:4095*Ts;

%双稳态系统参数

g=SNR(kk);

x1=AWGN(s,g); %噪声

%输入无噪信号傅立叶变换

y=fft(s,4096);

pyy=y.*conj(y)/4096;

ff=fs*(0:2048)/4096;

max=1;

for i=1:2049

if(pyy(i)>pyy(max))

max=i;

%输入信噪比计算

Bwin=0;

for i=1:2049

Bwin=pyy(i)+Bwin;

end

Bwin=(Bwin-pyy(max));

SNRin(kk,1)=10*log(pyy(max)/Bwin)/log(10);

%四阶龙格库塔法对双稳态输出信号求解

x=sr(a,b,h,x1);

%输出信号求傅立叶变换

%输出信噪比计算

Bwout=0;

for i=1:2049

Bwout=py(i)+Bwout;

end

Bwout=(Bwout-py(max));

SNRout(kk,1)=10*log(py(max)/Bwout)/log(10);

Q(kk,1)=SNRout(kk,1)-SNRin(kk,1);

end

clear

clc

theta=5;%1.1,2.1,5

s=1;

sigma=linspace(0,10,100);

randn('state',0)

if x(ii)>theta %传递函数

y(ii)=1;

else

y(ii)=0;

end

end

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

end;

title('加入周期信号较小，不足以翻转')

xlabel('x');

ylabel('U(x)');

v=moviein(20);

subplot(2,1,2)%周期信号较大时

A=3;

uu=yy(i,2)*yy(i,2)/2-b.*yy(i,1).^2/2+c.*yy(i,1).^4/4-A.*yy(i,1).*sin(w*t(i));

uu=-b.*yy(i,1).^2/2+c.*yy(i,1).^4/4-A.*yy(i,1).*sin(w*t(i));

plot(yy(i,1),uu,'.r','markersize',10)

holdoff

v(:,i)=getframe;

%oula.m 子程序:欧拉法解方程

fori=1:(n-1)

x(i+1)=x0+i*h;

y1=y(i)+h*feval(fun,x(i),y(i),r);%调用方程

y2=y(i)+h*feval(fun,x(i+1),y1,r);

y(i+1)=(y1+y2)/2;

end

%f.m 方程文件

%参数 b,c,A,w 为 1,故省略。

%odefilexinhao1.mode 文件

functionydot=xinhao1(t,x,flag,b,z,c,A,w)

ydot=[b*x(1)-c*x(1)^3+A*cos(w*t)+z;...

b*x(2)-3*c*x(1)*x(1)*x(2)-A*w*sin(w*t)];

clear all

clc

s=0.003*cos(2*pi*0.0001*t);

noise=sqrt(2*D)*randn(1,length(t));

x0=s+noise;

figure(1);

plot(t,x0);

f1=abs(fft(x0));

f1=f1(1:round(length(f1)/2));

ff=fs*(0:3999)/8000;

figure(2);

plot(ff,f1);