相同的噪声信号，Matlab与Python

Question

我试图用Matlab中的代码从python中获取类似的信号。

MATLAB中的代码：

Fs = 1e3;               % sampling frequency (in Hz)
L = 1e5;                % signal length (number of samples)
f0 = 0.1*Fs;            % cycle frequency (in Hz)   

x = rand(L, 1);
a = [1 -2*cos(2*pi*.2)*.9 .9^2];
x = filter(1,a,x);
x = x.*(1 + sin(2*pi*(0:L-1)'*f0/Fs));
x = x +std(x)*rand(L,1);
histogram(x,100);

Python中的代码：

Fs=10**3
L=10**5
f0=0.1*Fs
x=np.random.normal(0,1,L)
a=[1,-2*np.cos(2*np.pi*.2)*.9,.9**2]
x=sps.lfilter([1],a,x)
Random_modulated_signal=x*(1+np.sin(2*np.pi*np.arange(0,L)*f0/Fs))
Rms_Whitenoise=x+np.std(x)*np.random.normal(0,1,L)
plt.hist(Rms_Whitenoise,bins=100,edgecolor='k')
plt.show()

如果我绘制两个信号的直方图，它们不会共享相同的行为。添加std(x)*rnd(L,1)后，信号会有所不同。

奇怪的是，我使用了与前面所示相同的代码，但是我使用rand('twister', 2)和np.random.seed(2)代替使用正态分布，在两个程序中生成相同的值并且它有效。< / p>

请有人澄清我做错了什么？

两个信号的直方图

Answer 1

Matlab＆＃39; s rand是一个统一的随机分布，而不是正态分布。 matlab中numpy.random.normal的等效值为randn。 matlab rand的等价物是numpy.random.random。