不借助 MATLAB 函数库，自行实现连续小波变换（CWT）及注意要点

以下是一个不使用 MATLAB 内置 `cwt` 函数，自行实现连续小波变换（CWT）的代码： ```matlab % 生成信号 t = 0:0.01:1; % 时间向量，从0到1，采样间隔为0.01 x = sin(2*pi*10*t) + sin(2*pi*20*t); % 生成包含两个正弦波的信号 % 定义尺度范围 scales = 1:100; % 定义小波函数（这里以 Morlet 小波为例） function psi = morlet_wavelet(t) psi = pi^(-0.25) * exp(1i * 2 * pi * t) .* exp(-t.^2 / 2); end % 初始化小波系数矩阵 coefs = zeros(length(scales), length(x)); % 进行连续小波变换 for i = 1:length(scales) scale = scales(i); % 生成不同尺度下的小波 tau = (-5*scale:0.01:5*scale); % 时间平移范围 psi = morlet_wavelet(tau / scale) / sqrt(scale); % 进行卷积操作 for j = 1:length(x) % 对信号进行平移 shifted_x = [zeros(1, length(tau)), x, zeros(1, length(tau))]; shifted_x = shifted_x((length(tau)+1 - j):(length(tau)+length(x) - j)); % 计算小波系数 coefs(i, j) = sum(shifted_x .* conj(psi)) * 0.01; end end % 绘制原始信号 figure; subplot(2,1,1); plot(t, x); xlabel('Time (s)'); ylabel('Amplitude'); title('Original Signal'); % 绘制小波变换结果 subplot(2,1,2); imagesc(t, scales, abs(coefs)); set(gca, 'YDir', 'normal'); % 反转y轴方向 xlabel('Time (s)'); ylabel('Scale'); title('Continuous Wavelet Transform'); colorbar; % 显示颜色条 ``` 代码注释解释： 1. 首先生成一个包含 10Hz 和 20Hz 正弦波的信号。 2. 定义了 `morlet_wavelet` 函数，用于生成 Morlet 小波。 3. 初始化一个矩阵 `coefs` 来存储小波系数。 4. 通过两层循环实现连续小波变换：外层循环遍历不同的尺度，内层循环遍历信号的每个时间点。 5. 在每个尺度下，生成相应的小波，并对信号进行平移，然后计算小波系数。 6. 最后将原始信号和小波变换结果绘制在同一个图中。 我在以往的开发过程中，碰到过以下几点需要格外小心处理的地方： 1. 小波函数的定义要准确，不同的小波函数有不同的表达式，需要根据具体需求进行选择和实现。 2. 尺度和时间平移范围的选择会影响小波变换的精度和计算量，需要合理设置。 3. 在进行卷积操作时，要注意信号的平移和边界处理，避免出现索引越界等错误。 ######[AI写代码神器 | 757点数解答 | 2025-10-16 10:56:43]