我有一组.mat数据，采样率为1000000，时间长度为5秒，想要通过连续小波变换获得小波尺度图，x轴是时间，y轴是频率。请给我一个matlab程序，

要在MATLAB中对1MHz采样率、5秒长的时间序列数据进行连续小波变换并生成尺度图，你可以使用`cwt`函数。下面是一个简单的例子：

% 导入数据
data = load('your_data.mat'); % 将'your_data.mat'替换为你的实际文件名

% 确保数据类型是双精度浮点数
data = double(data);

% 设置参数
Fs = 1e6; % 采样率
t = (0:length(data)-1)'/Fs; % 时间数组，从0到5秒
cwt_matrix = cwt(data, 'mexh', t, 'FrequencyRange', [0 Fs/2]); % 使用MEXH小波，t是时间，频率范围设定为整个可用频域

% 计算尺度
scales = linspace(1/Fs, 1/(2*Fs), length(cwt_matrix)); % 小波尺度范围，从最低到最高频率

% 绘制小波尺度图
figure;
imagesc(t, scales, abs(cwt_matrix));
xlabel('Time (seconds)');
ylabel('Scale');
colorbar;

% 如果你想显示频率轴，可以这样：
freq_axis = linspace(0, Fs/2, size(cwt_matrix, 2)) / Fs;
title(['Continuous Wavelet Transform - Frequency vs Time']);
hold on;
plot(freq_axis, ones(size(freq_axis)), 'r--', 'LineWidth', 2); % 假设红色虚线表示频率轴
hold off;

请确保在运行此代码前，已安装了支持小波分析的MATLAB工具箱，如Wavelet Toolbox。如果遇到错误，可能需要检查数据是否导入正确，以及MATLAB是否支持你的操作系统。