局放源模拟_40Khz_窄带信号_局放_局部放电仿真_四种衰减函数

%%%%%%%%使用db5小波进行5层的小波变换%%%%%%%%%%%%%%%%%%%% %分解得到近似和细节系数使用db5小波进行5层的小波变换%%%%%% %%%%%%%%%%%%%%%%%%分解得到近似和细节系数%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%未进行消噪处理%%%%%%%%%%%%%%%%%%% clc; clear all; close all; %%%%%%%%%%%%%%%%%%%%根据放电模型生成数据%%%%%%%%%%%%%% %生成目标信号 t=0:0.0000001:0.001; data_t=zeros(10001,1); %目标信号 data_n=zeros(10001,1); %加噪信号 data1=zeros(10001,1);%单指数衰减信号 T1=20 * 10^(-6); data2=zeros(10001,1);%单指数衰减振荡信号 T2=32 * 10^(-6); data3=zeros(10001,1);%双指数衰减信号 T3=T1; data4=zeros(10001,1);%单指数衰减信号 T4=64 * 10^(-6); data5=zeros(10001,1);%单指数衰减振荡信号 T5=T2; data6=zeros(10001,1);%双指数衰减信号 T6=T1; data7=zeros(10001,1);%白噪声 f=10000 * 10^3; %采样频率 T=1/f; fc=40 * 10^3; %信号频率 for i=201:1:10001 %单指数衰减信号 data1(i,1)=exp(-1/T1*(i-201)*T); end for i=2001:1:10001 %单指数衰减振荡信号 data2(i,1)=exp(-1/T2*(i-2001)*T)*sin(2*pi*fc*(i-2001)*T+pi/2); end for i=4001:1:10001 %双指数衰减信号 data3(i,1)=exp(-1.3/T3*(i-4001)*T)-exp(-2.2/T3*(i-4001)*T); end for i=6001:1:10001 %双指数衰减振荡信号 data4(i,1)=(exp(-1.3/T4*(i-6001)*T)-exp(-2.2/T4*(i-6001)*T))* sin(2 * pi * fc * (i-6001) * T +pi/2); end % for i=1003:1:2001 %单指数衰减信号 % data4(i,1)=exp(-1/T4*(i-1003)*T); % end % % for i=1207:1:2001 %单指数衰减振荡信号 % data5(i,1)=exp(-1/T5*(i-1207)*T)*sin(2*pi*fc*(i-1207)*T+pi/2); % end % % for i=1607:1:2001 %双指数衰减信号 % data6(i,1)=exp(-1.3/T6*(i-1607)*T)-exp(-2.2/T6*(i-1607)*T); % end data7=0.05*randn(10001,1); figure(1); subplot(6,2,1); plot(t,data1) title('单指数衰减信号'); subplot(6,2,2); plot(data1) subplot(6,2,3); plot(t,data2) title('双指数衰减信号'); subplot(6,2,4); plot(data2) subplot(6,2,5); plot(t,data3) title('单指数衰减振荡信号'); subplot(6,2,6); plot(data3) subplot(6,2,7); plot(t,data4) title('双指数衰减振荡信号'); subplot(6,2,8); plot(data4) subplot(6,2,9); plot(data7) title('噪声信号'); figure(2); data_T=data1+data2+data3+data4+data5+data6; plot(t,data_T); title('目标信号'); xlabel('时间t/s'); S1_fft_abs = abs(fft(data4 - mean(data4))); figure(4); plot(S1_fft_abs) figure(5); ST_fft_abs = abs(fft(data_T - mean(data_T))); plot(ST_fft_abs) figure(3); data_N=data1+data2+data3+data4+data5+data6+data7; plot(t,data_N); title('加噪信号'); xlabel('时间t/s'); % subplot(3,3,5); % plot(t,data5) % subplot(3,3,6); % plot(t,data6) % subplot(3,3,7); % plot(t,data7)