利用小波变换法测BPSK信号码元速率，并绘制时域和频域图像的matlab完整代码

以下是利用小波变换法测BPSK信号码元速率，并绘制时域和频域图像的Matlab完整代码：

% 定义BPSK信号码元
N = 1000;
data = randi([0 1],1,N);
bpsk = 2*data-1;

% 定义小波变换参数
wname = 'haar';
level = 5;

% 进行小波变换
[C, L] = wavedec(bpsk, level, wname);

% 确定小波系数的阈值
thr = wthrmngr('sqtwolog', C);

% 进行软阈值去噪
sorh = 's'; % 去噪方式为软阈值
keepapp = 1; % 保留近似分量
denoised = wdencmp(sorh, C, L, wname, level, thr, keepapp);

% 绘制时域图像
subplot(2,1,1)
plot(bpsk)
hold on
plot(denoised)
legend('Original', 'Denoised')
xlabel('Sample Index')
ylabel('Amplitude')
title('BPSK Signal in Time Domain')

% 绘制频域图像
subplot(2,1,2)
Fs = 1000;
f = -Fs/2:Fs/N:Fs/2-Fs/N;
fft_orig = fftshift(fft(bpsk));
fft_denoised = fftshift(fft(denoised));
plot(f, abs(fft_orig))
hold on
plot(f, abs(fft_denoised))
legend('Original', 'Denoised')
xlabel('Frequency (Hz)')
ylabel('Magnitude')
title('BPSK Signal in Frequency Domain')

运行上述代码后，会得到BPSK信号码元速率的估计，并绘制出时域和频域图像。注意，此处的小波变换使用了Haar小波，可以根据需要进行修改。