我有一个只有一列的.mat数据，采样率为1000000，共计5000000个点，想要通过连续小波变换获得小波尺度图，x轴是时间，y轴是频率。请给我一个matlab程序，

您可以使用MATLAB的内置函数`cwt`来进行连续小波变换（CWT）。假设您的数据存储在一个名为`data`的一维数组中，以下是创建尺度图的简单程序：

% 加载数据
load('your_datafile.mat') % 替换为您的数据文件名
data = data_column; % 假设数据已加载到data_column，现在直接使用data

% 确定采样率
fs = 1e6; % 采样率是1000000 Hz

% 创建时间向量
t = (0:size(data, 1))/fs; % 时间从0开始，单位是秒

% 选择合适的Ricker小波（常用的选择之一）
wavelet = ricker(64, 1/fs); % 64点的小波，中心频率对应于采样率的1/64

% 使用cwt函数进行连续小波变换
[cwt_matrix, frequencies] = cwt(data, wavelet, logspace(-3, 2)); % 调整尺度范围

% 对连续小波矩阵取绝对值并绘制图像
figure;
surf(t, frequencies, abs(cwt_matrix));
xlabel('Time (s)');
ylabel('Frequency (Normalized to Fs/2)');
zlabel('Magnitude');
title('Continuous Wavelet Transform Scaleogram');

% 保存图像
saveas(gcf, 'scaleogram.jpg', 'jpg'); % 以jpg格式保存结果