有哪些歌可以听出耳机音质好坏？ 第1页

       赫里斯托日沃米罗夫(教授)于1987年出生于保加利亚瓦尔纳。 他分别于2010年和2012年在瓦尔纳技术大学获得通信设备和技术的学士和硕士学位。 2016年，他获得了通信理论的博士学位。 日沃米罗夫先生目前是瓦尔纳技术大学电气工程理论和测量系的副教授，IEEE、保加利亚科学-技术联盟联合会 和保加利亚科学家联盟的成员。 他的研究兴趣包括信号处理、电气和电子测量领域以及Matlab编程（信号处理、数据采集、数据可视化等）。     

       %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Pink Noise Generation with MATLAB Implementation % % % % Author: M.Sc. Eng. Hristo Zhivomirov 07/30/13 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  function y = pinknoise(N)  % function: y = pinknoise(N) % N - number of samples to be returned in row vector % y - row vector of pink (flicker) noise samples  % The function generates a sequence of pink (flicker) noise samples. % Pink noise has equal energy in all octaves (or similar log bundles) of frequency. % In terms of power at a constant bandwidth, pink noise falls off at 3 dB per octave.  % define the length of the vector % ensure that the M is even if rem(N,2) M = N+1; else M = N; end  % generate white noise x = randn(1, M);  % FFT X = fft(x);  % prepare a vector for 1/f multiplication NumUniquePts = M/2 + 1; n = 1:NumUniquePts; n = sqrt(n);  % multiply the left half of the spectrum so the power spectral density % is proportional to the frequency by factor 1/f, i.e. the % amplitudes are proportional to 1/sqrt(f) X(1:NumUniquePts) = X(1:NumUniquePts)./n;  % prepare a right half of the spectrum - a copy of the left one, % except the DC component and Nyquist frequency - they are unique X(NumUniquePts+1:M) = real(X(M/2:-1:2)) -1i*imag(X(M/2:-1:2));  % IFFT y = ifft(X);  % prepare output vector y y = real(y(1, 1:N));  % ensure unity standard deviation and zero mean value y = y - mean(y); yrms = sqrt(mean(y.^2)); y = y/yrms;  end