基于小波变换的电能质量扰动时间定位 构造理想的暂降信号波形函数(观测时间长度取0.6Su(t) ui(t) + uz(t) + uz(t)其中: 0<t≤0.2u(t)= Asin(w(t)) uz(t) = A,sin(wz(t)) uz(t) = Azsin(w;(t)) 0.2 <t≤0.30.3<t≤0.5 令A=A2= 8500VA =5200V.步长取0.001，仿真出理想的暂降波形，并对其进行小波分解，根据信号奇异性表现，计算出暂降持续时间。请给出Matlab程序。

好的，下面是基于小波变换的电能质量扰动时间定位的Matlab程序，包括构造理想的暂降信号波形函数、进行小波分解和计算暂降持续时间三个部分：

% 构造理想的暂降信号波形函数
t = 0:0.001:0.5; % 时间步长为0.001，总时间为0.5
A = 8500; % A的值为8500VA
A2 = 5200; % A2的值为5200V
w = 2*pi*50; % 频率参数w
wz = 2*pi*200; % 频率参数wz
ui = A*sin(w*t); % 计算ui(t)
uz1 = A*sin(wz*t(1:round(length(t)*0.2))); % 计算uz(t)在0<t≤0.2的时间段内的值
uz2 = A*sin(wz*t(round(length(t)*0.2)+1:round(length(t)*0.3))); % 计算uz(t)在0.2<t≤0.3的时间段内的值
uz3 = A2*sin(wz*t(round(length(t)*0.3)+1:round(length(t)*0.5))); % 计算uz(t)在0.3<t≤0.5的时间段内的值
uz = [uz1, uz2, uz3]; % 拼接uz(t)
s = ui + uz; % 拼接ui(t)和uz(t)得到s(t)

% 进行小波分解
level = 5; % 小波分解的层数为5层
wname = 'db4'; % 选择使用db4小波
[c, l] = wavedec(s, level, wname); % 进行小波分解
A5 = appcoef(c, l, wname); % 提取第5层低频系数
D5 = detcoef(c, l, 5); % 提取第5层高频系数
D4 = detcoef(c, l, 4); % 提取第4层高频系数
D3 = detcoef(c, l, 3); % 提取第3层高频系数
D2 = detcoef(c, l, 2); % 提取第2层高频系数
D1 = detcoef(c, l, 1); % 提取第1层高频系数

% 计算暂降持续时间
energy_ratio = cumsum(D5.^2)./sum(D5.^2); % 计算能量比
duration_index = find(energy_ratio >= 0.95, 1); % 寻找能量比达到95%的位置
duration = (duration_index-1)*0.001; % 计算暂降持续时间（单位为秒）

% 绘图
figure;
subplot(2, 1, 1);
plot(t, s);
title('理想的暂降波形');
xlabel('时间（秒）');
ylabel('电压（伏特）');
subplot(2, 1, 2);
plot(0:0.001:duration, D5(1:duration_index));
title('第5层高频系数随时间的变化');
xlabel('时间（秒）');
ylabel('第5层高频系数');

以上程序可以绘制出理想的暂降波形和第5层高频系数随时间的变化曲线，并且计算出暂降持续时间。需要注意的是，在计算暂降持续时间时，我们选择能量比达到95%的位置作为暂降持续时间的结束位置。