Nobuhito Mori
Navigation :
開始行:
*フーリエスペクトルの計算 [#bc0c7995]
関数fftでは,データ長が値にかかってくる.
f: データ
fの振幅とスペクトルの振幅を合わせたいときは
P = 2*abs(fft(f))/length(f)
Powersectrum
P = c.*conj(c)/length(c);
とする.
*** Sample of cpsd [#eb73fa00]
Fs = 1/dt;
filt = hanning(nfft/16);
n_overlap = 15;
[P_cpsd, f] = cpsd(x,x,filt, n_overlap, nfft, Fs);
*** Sample of csd [#eb73fa00]
Fs = 1/dt;
filt = hanning(nfft/16);
n_overlap = 15;
[S,freq] = csd( eta2, eta2, nfft, Fs, filt, n_overlap, '...
* パワースペクトルの計算方法 [#l26be0c9]
パワースペクトルは,複素フーリエ係数cにより定義するため,...
c*conj(c)で計算できる.
しかし,パワースペクトルは確率密度関数であり,適切なフィ...
MATLABでは,フィルター付パワースペクトル推定法としてcsd, ...
*** データの用意 [#w563a5ef]
dt = 0.001;
t = 0:dt:0.6;
tmax = max(t);
df = 1/tmax;
Fs = 1/dt;
x = 5*sin(2*pi*50*t)+sin(2*pi*120*t);
nfft = length(x);
plot(t,x)
disp('mean energy of x')
var(x)
*** fftで直接計算 [#naf72ac4]
c = fft(x,nfft);
S_fft = c.*conj(c)/length(x);
P_fft = S_fft;
disp('mean energy of x by fft')
sum(P_fft)/length(x)
*** cpsdで計算 [#k9fe556e]
[P_cpsd, f] = cpsd(x,x,[],[],nfft,Fs);
disp('mean energy of x by cpsd')
sum(P_cpsd)/length(x)*Fs
*** csdで計算 [#p6586c74]
[P_csd, f] = csd(x,x,nfft,Fs);
disp('mean energy of x by csd')
sum(P_csd)/length(f)
* ETC [#nf9204d4]
** cpsdとcsd [#p6a18c8e]
csdはそのうち無くなり,cpsdが今後サポートされる予定.
しかし,cpsdはcpdとちょっと仕様が異なるので移行に注意が必...
cpsdはパワースペクトルが周波数で割られており,単位周波数...
csdはそうでなく,fftベースで求められたパワー量をベースと...
信号のパワーの平均を得るためにはcsdの方は,fftベースと同...
終了行:
*フーリエスペクトルの計算 [#bc0c7995]
関数fftでは,データ長が値にかかってくる.
f: データ
fの振幅とスペクトルの振幅を合わせたいときは
P = 2*abs(fft(f))/length(f)
Powersectrum
P = c.*conj(c)/length(c);
とする.
*** Sample of cpsd [#eb73fa00]
Fs = 1/dt;
filt = hanning(nfft/16);
n_overlap = 15;
[P_cpsd, f] = cpsd(x,x,filt, n_overlap, nfft, Fs);
*** Sample of csd [#eb73fa00]
Fs = 1/dt;
filt = hanning(nfft/16);
n_overlap = 15;
[S,freq] = csd( eta2, eta2, nfft, Fs, filt, n_overlap, '...
* パワースペクトルの計算方法 [#l26be0c9]
パワースペクトルは,複素フーリエ係数cにより定義するため,...
c*conj(c)で計算できる.
しかし,パワースペクトルは確率密度関数であり,適切なフィ...
MATLABでは,フィルター付パワースペクトル推定法としてcsd, ...
*** データの用意 [#w563a5ef]
dt = 0.001;
t = 0:dt:0.6;
tmax = max(t);
df = 1/tmax;
Fs = 1/dt;
x = 5*sin(2*pi*50*t)+sin(2*pi*120*t);
nfft = length(x);
plot(t,x)
disp('mean energy of x')
var(x)
*** fftで直接計算 [#naf72ac4]
c = fft(x,nfft);
S_fft = c.*conj(c)/length(x);
P_fft = S_fft;
disp('mean energy of x by fft')
sum(P_fft)/length(x)
*** cpsdで計算 [#k9fe556e]
[P_cpsd, f] = cpsd(x,x,[],[],nfft,Fs);
disp('mean energy of x by cpsd')
sum(P_cpsd)/length(x)*Fs
*** csdで計算 [#p6586c74]
[P_csd, f] = csd(x,x,nfft,Fs);
disp('mean energy of x by csd')
sum(P_csd)/length(f)
* ETC [#nf9204d4]
** cpsdとcsd [#p6a18c8e]
csdはそのうち無くなり,cpsdが今後サポートされる予定.
しかし,cpsdはcpdとちょっと仕様が異なるので移行に注意が必...
cpsdはパワースペクトルが周波数で割られており,単位周波数...
csdはそうでなく,fftベースで求められたパワー量をベースと...
信号のパワーの平均を得るためにはcsdの方は,fftベースと同...
ページ名:
MENU
INTRODUCTION
RESEARCH
EDUCATION
MEMO
COMPUTER
RESEARCH
ETC
LINK
最新の5件
2024-03-23
Publications/Journal Papers
2024-03-17
Introduction/CV in English
Introduction/CV in Japanese
2024-02-01
Publications/学会発表講演集・アブストラクト
2024-01-26
Publications/メディア等
カスタム検索
Counter
Total : 0
Today : 0
Yesterday : 0
Edit Menu
[
新規
] [
ヘルプ
]
Site Search
AND
OR