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為低