想请问一下这个代码哪里有问题，求助各位大佬

[Echoheaven]

clc
clear

x = xlsread('D:\掘进数据\研究环原始掘进参数.xlsx', 1, 'B2:B501');
signal = x';

noise = 0.5 * randn(size(signal));
noisy_signal = signal + noise;


level = 4; 
wname = 'db8'; 
[c, l] = wpdec(noisy_signal, level, wname); 

threshold = thselect(noisy_signal, 'heursure'); 

denoised_coeffs = wthresh(c, 's', threshold); 
denoised_signal = wprec(denoised_coeffs, l, wname); 

subplot(3,1,1);
plot(1:length(signal), signal);
title('原始信号');
subplot(3,1,2);
plot(1:length(noisy_signal), noisy_signal);
title('含噪信号');
subplot(3,1,3);
plot(1:length(denoised