clc; clear; %设置八个幅值 A1=3;A2=8;A3=5;A4=2;A5=6;A6=7;A7=4;A8=1; %设置八个频率分量 f1=1;f2=2;f3=3;f4=4;f5=5;f6=6;f7=7;f8=8; %定义采样区间 x=0:0.001:8; %构建八个信号函数 y1=A1*sin(2*pi*f1*x).*(0<x&x<1); y2=A2*sin(2*pi*f2*x).*(1<x&x<3); y3=A3*sin(2*pi*f3*x).*(2<x&x<3); y4=A4*sin(2*pi*f4*x).*(2<x&x<4); y5=A5*sin(2*pi*f5*x).*(4<x&x<5); y6=A6*sin(2*pi*f6*x).*(5<x&x<6); y7=A7*sin(2*pi*f7*x).*(6<x&x<7); y8=A8*sin(2*pi*f8*x).*(0<x&x<8); %八个信号叠加求和构建非平稳信号函数 y=y1+y2+y3+y4+y5+y6+y7+y8;继续编写matlab代码，选用一个合适的小波画出y的时频能量图

可以使用小波变换来绘制非平稳信号的时频能量图。下面是使用MATLAB编写的代码：

% 导入小波变换工具箱
import wavelet.*

% 小波变换参数设置
wname = 'db4'; % 选用 Daubechies 4 小波
level = 5;     % 小波分解的层数

% 进行小波分解
[c, l] = wavedec(y, level, wname);

% 计算小波系数的能量
energy = zeros(level+1, length(y));
for i = 1:level+1
    energy(i, :) = wrcoef('a', c, l, wname, i-1);
end

% 绘制时频能量图
t = 0:0.001:8;
figure;
imagesc(t, 1:level+1, energy);
colormap(jet);
xlabel('时间（秒）');
ylabel('尺度');
title('非平稳信号的时频能量图');

这段代码首先导入小波变换工具箱，然后选用 Daubechies 4 小波对信号进行五层小波分解，然后计算小波系数的能量，并绘制时频能量图。运行这段代码，就可以得到非平稳信号的时频能量图。