バンドパスフィルターで雷鳴を合成する方法

Question

MATLABで雷鳴を合成する作業を進めています。 私はピンクノイズを生成していたし、今私は30Hzの上のバンドパスフィルタを追加しようとしているが、私はエラーを取得しています：

The frequency specifications Fpass, Fstop must have increasing values. 

ここに私のコードです：たぶん

dur=[0:1/Fs:time-1/Fs]; 
sz=length(dur); 
g=randn(1, sz); % 
Y=fft(g); 

NumUniquePts = sz/2 + 1; 
n = 1:NumUniquePts; 
n = sqrt(n); 
Y(1:NumUniquePts) = Y(1:NumUniquePts)./n; 

Y(NumUniquePts+1:sz) = real(Y(sz/2:-1:2)) -1i*imag(Y(sz/2:-1:2)); 

y = ifft(Y); 

y = real(y(1, 1:sz)); 

y = y - mean(y); 
yrms = sqrt(mean(y.^2)); 
y = y/yrms; 

d = fdesign.lowpass('Fp,Fst,Ap,Ast',30,10,0.5,50,48e4); 
Hd1 = design(d,'equiripple'); 
Hd2 = design(d,'butter'); 

out = filter(Hd1,y); 

そこに他の方法それを合成する？

Answer 1

I am getting an error. "The frequency specifications Fpass, Fstop must have increasing values."

このエラーはFstopがFpassより大きくなければならないことを意味しています。あなたはFpass=30Hzが必要だと述べました。したがって、Fstop > 30が必要ですが、10Hzの値を指定しました。このエラーを修正するには、Fstopの値を増やす必要があります（たとえば、最初は100Hz、好きなところまで調整してください）。