clear;
clc;
t=1:2000;
s1=zeros(1,2000);
s2=exp(-t./30);
s2=s2./max(abs(s2));
s3=(exp(-t./30)-exp(-3*t./30));
s3=s3./max(abs(s3));
s4=exp(-t./30).*sin(2*pi*0.02*t);
s4=s4./max(abs(s4));
s5=((exp(-t./30)-exp(-3*t./30)).*sin(2*pi*0.02*(t)));
s5=s5./max(abs(s5));
s0=[s1,s2,s3,s4,s5];%原信号
wgns=awgn(s0,0,'measured');%加入高斯白噪
%mixs=wgns;
x=linspace(0,0.0002,1e4);%正弦干扰
noise_sin_f1=100e4;%设置正选噪声频率
noise_sin_f2=300e4;
noise_sin_f3=500e4;
noise_sin=sin(2*pi*noise_sin_f1*x)+sin(2*pi*noise_sin_f2*x)+sin(2*pi*noise_sin_f3*x);
noise_sin=noise_sin./max(abs(noise_sin));
% mixs=s0+0.2*noise_sin;
mixs=wgns+0.2*noise_sin;
Lbp=10;
g=1/30*ones(1,Lbp);%结构元素:扁平,元素宽度Lbp
% Lcos=13;
% y=0:Lcos-1;
% g=1/30*sin(pi/(Lcos-1)*y);%结构元素:正弦,元素宽度Lcos
%
% Lcir=35;
% y=0:Lcir-1;
% g=sqrt((Lcir-1)^2/4-(y-(Lcir-1)/2).^2);
% g=1/30*g./max(g);%结构元素:椭圆
r=ochunhe(mixs,g);
subplot(2,1,1)
plot(mixs)
xlabel('a)叠加两种干扰的局放信号','Fontsize',14)
subplot(2,1,2)
plot(r)
xlabel('b)交替混合滤波器处理结果','Fontsize',14)
snr=snr(s0,r)
emsr=sqrt(1e-4*sum((s0-r).^2))
cor=sum(s0.*r)/sqrt(sum(s0.^2)*sum(r.^2))
评论1