给我一个设置40dB信噪比的程序

以下是一个设置40dB信噪比的MATLAB程序示例：

% 生成信号
fs = 1000; % 采样率
t = 0:1/fs:1; % 时间序列
f = 100; % 信号频率
x = sin(2*pi*f*t); % 信号

% 生成噪声
noise_power = 10^(-40/10); % 噪声功率
noise = randn(size(x)) * sqrt(noise_power); % 均值为0，方差为noise_power的高斯白噪声

% 加噪声
snr = 40; % 信噪比
y = x + noise * sqrt((x'*x)/(noise'*noise)*10^(-snr/10)); % 加噪声

% 绘制信号和噪声
figure;
plot(t, x, 'b', t, noise, 'r');
xlabel('时间（s）');
ylabel('幅度');
legend('信号', '噪声');

% 绘制加噪声后的信号
figure;
plot(t, x, 'b', t, y, 'r');
xlabel('时间（s）');
ylabel('幅度');
legend('信号', '加噪声后的信号');

该程序首先生成一个频率为100Hz的正弦信号，然后生成一个方差为10^(-40/10)的高斯白噪声，并将其加到信号上，从而得到信噪比为40dB的加噪声信号。程序最后绘制了信号、噪声以及加噪声后的信号的波形图。