?%%初始化程序
clearclc
t1=clock;
?%%?載入7.6日00:15到7.21日24:00的數據
load（‘SJ76721.mat‘）;%matrix
YSJ=?SJ76721;
disp（‘------7.6日00:15到7.21日24:00的降噪處理結果-------‘）;?
?%?載入7.6日00:15到7.24日24:00的數據
%?load（‘SJ76724.mat‘）;
%?YSJ=?SJ76724;
%?disp（‘------7.6日00:15到7.24日24:00的降噪處理結果-------‘）;
?%%?數據預處理
[cl]=size（YSJ）;
Y=[];
for?i=1:c
????Y=[YYSJ（i:）];
end
[c1l1]=size（Y）;
X=[1:l1];
?%%?繪制原始信號圖像
figure（1）;
plot（XY）;
xlabel（‘負荷個數/（個）‘）;
ylabel（‘kW·h‘）;
title（‘原始信號‘）;
?%%?硬閾值處理
lev=3;
xd=wden（Y‘heursure‘‘h‘‘one‘lev‘db4‘）;%硬閾值去噪處理后的信號序列
figure（2）
plot（Xxd）
xlabel（‘負荷個數/（個）‘）;
ylabel（‘kW·h‘）;
title（‘硬閾值去噪處理‘）
set（gcf‘Color‘[1?1?1]）
?%%?軟閾值處理
lev=3;
xs=wden（Y‘heursure‘‘s‘‘one‘lev‘db4‘）;%軟閾值去噪處理后的信號序列
figure（3）
plot（Xxs）
xlabel（‘負荷個數/（個）‘）;
ylabel（‘kW·h‘）;
title（‘軟閾值去噪處理‘）
set（gcf‘Color‘[1?1?1]）
%%?改進閾值后的去噪處理
lev=3;
xz=wden（Y‘sqtwolog‘‘s‘‘sln‘lev‘db4‘）;%改進閾值去噪處理后的信號序列
figure（4）
plot（Xxz）;
xlabel（‘負荷個數/（個）‘）;
ylabel（‘kW·h‘）;
title（‘改進閾值后的去噪處理‘）
set（gcf‘Color‘[1?1?1]）
%%?計算信噪比SNR
Psig=sum（Y*Y‘）/l1;
Pnoi1=sum（（Y-xd）*（Y-xd）‘）/l1;
Pnoi2=sum（（Y-xs）*（Y-xs）‘）/l1;
Pnoi3=sum（（Y-xz）*（Y-xz）‘）/l1;
SNR1=10*log10（Psig/Pnoi1）;
SNR2=10*log10（Psig/Pnoi2）;
SNR3=10*log10（Psig/Pnoi3）;
%%?計算均方根誤差RMSE
RMSE1=sqrt（Pnoi1）;
RMSE2=sqrt（Pnoi2）;
RMSE3=sqrt（Pnoi3）;
%%?計算均方根誤差RMSE
disp（[‘硬閾值去噪處理的SNR=‘num2str（SNR1）‘，RMSE=‘num2str（RMSE1）]）;
disp（[‘軟閾值去噪處理的SNR=‘num2str（SNR2）‘，RMSE=‘num2str（RMSE2）]）;
disp（[‘改進閾值后的去噪處理SNR=‘num2str（SNR3）‘，RMSE=‘num2str（RMSE3）]）;
t2=clock;
tim=etime（t2t1）;
disp（[‘------------------運行耗時‘num2str（tim）‘秒-------------------‘]）

?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1855??2020-03-13?15:26??小波降噪matlab代碼\main.m
?????文件????????4733??2020-03-12?18:09??小波降噪matlab代碼\SJ76721.mat
?????文件????????5590??2020-03-12?18:10??小波降噪matlab代碼\SJ76724.mat