信号检测基于matlab双稳随机共振微弱信号检测含Matlab源码 1701期

Posted 2022-02-17 紫极神光

一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【信号检测】基于matlab双稳随机共振微弱信号检测【含Matlab源码 1701期】

获取代码方式2：
通过订阅紫极神光博客付费专栏，凭支付凭证，私信博主，可获得此代码。

备注：
订阅紫极神光博客付费专栏，可免费获得1份代码（有效期为订阅日起，三天内有效）；

二、部分源代码

 function [xu]=xulietu(alpha,beta,deta,mu)
clc;
% clear  all;
% close all;
fs=100;%采样频率
Ts=1/fs;%采样时间
% h=1/fs;%时间步长
% t=0:Ts:4095*Ts;
t=0:Ts:409500*Ts;
% t=0:Ts:511*Ts;
n=length(t);

% alpha=1.2;
% beta=0;
% deta=10;
%  mu=0;

for i=1:n
  v=unifrnd(-pi/2,pi/2);
  w=exprnd(1);
%   if 0<alpha<1
    if (0<alpha)&&(alpha<1)
     deta0=deta*(1+beta.^2*tan(pi*alpha/2).^2).^(1/(2*alpha));
     beta1=2*atan(beta*tan(pi*alpha/2))/(pi*alpha);
     v0=-atan(beta*tan(pi*alpha/2))/alpha;
     x1=sin(alpha*(v-v0))/(cos(v)).^(1/alpha)*(cos(v-alpha*(v-v0))/w).^((1-alpha)/alpha);
     xu(i)=deta0*x1+mu;
   elseif alpha==1
       beta1=beta;
       deta0=2/pi*deta;
       x1=(pi/2+beta1*v)*tan(v)-beta1*log10(w*cos(v)/(pi/2+beta*v));
       xu(i)=deta0*x1+mu+2/pi*deta0*beta*log(deta0);
    
%    elseif 1<alpha<2
   elseif (1<alpha)&&(alpha<2)
          beta1=2*atan(beta*tan(pi*alpha/2))/(pi*(alpha-2));
           deta0=deta*(1+beta.^2*tan(pi*alpha/2).^2).^(1/(2*alpha));
          v0=-atan(beta*tan(pi*alpha/2))/alpha;
          x1=sin(alpha*(v-v0))/(cos(v)).^(1/alpha)*(cos(v-alpha*(v-v0))/w).^((1-alpha)/alpha);
          xu(i)=deta0*x1+mu;
  
    elseif alpha==2
        v0=-atan(beta*tan(pi*alpha/2))/alpha;
        x1=sin(alpha*(v-v0))/(cos(v)).^(1/alpha)*(cos(v-alpha*(v-v0))/w).^((1-alpha)/alpha);
        xu(i)=deta0*x1+mu;
     else
       xu(i)=randn;
    end
end



% plot(t,xu);xlabel('采样点');
% % ylim([-100,100]);
% ylabel('幅值');
% xlim([0,2000]);
% ylim([-4000,4000]);
%title('a稳定噪声序列');

三、运行结果

四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社，2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社，2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社，2018.