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

在Matlab中，你可以使用`cwt`函数来执行Morlet小波变换，并创建小波尺度图。首先，你需要安装Wavelet Toolbox，如果你还没有安装，可以考虑安装。以下是基本的步骤：

% 假设你的.mat文件名为data.mat并且已经读取进来了，变量名可能是'your_data'
load('data.mat', 'your_data');

% 设置小波参数
fs = 1000000; % 采样率
time_length = 5; % 时间长度 (秒)
scale_range = logspace(-3, 3, 100); % 小波尺度范围

% 使用Morlet小波
[wavelet, scales, freqs] = cwt(your_data, morl, scales);

% 频谱计算和可视化
figure;
surf(scales, freqs, abs(wavelet).^2);
title('Morlet小波变换 - 尺度图');
xlabel('Scale');
ylabel('Frequency (Hz)');
zlabel('Amplitude Squared');

% 创建尺度图
scale_profile = squeeze(mean(abs(wavelet), 2)); % 平均幅度
figure;
plot(scales, scale_profile);
title('小波尺度图');
xlabel('Scale');
ylabel('Average Amplitude');

在这个代码中，`morl`通常代表 Morlet 小波函数，需要你自己导入或者生成。`mean(abs(wavelet), 2)` 是为了得到每个尺度的平均幅度。