How to insert noise in a signal using matlab?

12 次查看（过去 30 天）

Nur Fauzira Saidin 2015-11-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/253187-how-to-insert-noise-in-a-signal-using-matlab

评论： Star Strider 2015-11-5

Hi. I want to create a signal with noise. I found this coding but I don't know how it works. Can somebody explain this to me?

x= cos(2*pi*12*[0:0.001:1.23]);
x(end) = [];
[b a] = butter(2,[0.6 0.7],'bandpass');
filtered_noise = filter(b,a,randn(1, length(x)*2));
x = (x + 0.5*filtered_noise(500:500+length(x)-1))/length(x)*2;

This is the result when I plot the x.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Star Strider 2015-11-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/253187-how-to-insert-noise-in-a-signal-using-matlab#answer_198620

It creates a signal in ‘x’, then creates band-limited noise by bandpass-filtering a Gaussian random vector (created by the randn call), the same length as ‘x’, and stores it in the ‘filtered_noise’ variable. It then adds that to the original ‘x’ and scales it to create the final value of ‘x’. (There is a bit more involved in the code than I described, but that is the essence of it.)

6 个评论
显示 4更早的评论隐藏 4更早的评论

Nur Fauzira Saidin 2015-11-5

在 MATLAB Online 中打开

The length of 'x' is 1230, so the final amplitude is 1/(1230/2) = 1.63x10-3. Okay I understand it now. Thanks! Now if you don't mind, can you help me with the next step? I want to plot the frequency spectrum of signal 'x'.

%plot magnitude spectrum of the signal
X_mags = abs(fft(x));
figure,plot(X_mags);
xlabel('DFT Bins');
ylabel('Magnitude');

Then, this code is to plot the first half of DFT. I don't quite understand this part. Can you explain this to me?

%plot first half of DFT (normalised frequency)
num_bins = length(X_mags);
figure,plot([0:1/(num_bins/2 -1):1], X_mags(1:num_bins/2));
xlabel('Normalised frequency (\pi rads/sample)');
ylabel('Magnitude');