sym4小波滤波

clear all;
%设置采样点数为1000
N=1000;
t=1:N;
pi=3.14;
%前250个点为幅值220，中间250点为幅值110，后500个点为幅值220
%构成电压暂降模型
sig（1:250）=220*sin（100*pi*（1:250））;
sig（251:500）=110*sin（100*pi*（251:500））;
sig（501:1000）=220*sin（100*pi*（501:1000））;
%生成随机白噪声
y=20*randn（1N）;
%给电压暂降模型添加噪声，模拟实际电网含噪声模型
x=sig+y;
figure（1）;
%绘制电压暂降信号
subplot（211）;plot（sig）;title（‘原始信号‘）;
xlabel（‘采样点‘）;
ylabel（‘幅值/V‘）;
%绘制含噪声电压暂降信号
subplot（212）;
plot（x）;
title（‘含噪声信号‘）;
xlabel（‘采样点‘）;
ylabel（‘幅值/V‘）;
%用sym4小波对含噪声电压暂降信号进行3层分解并提取系数
%aj为低