
[cory,lag]=xcorr(f,'unbiased');
plot(lag/fs,cory); %自相关函数(无偏差的),其中,f为原函数,cory为要求的自相关函数,lag为自相关函数的长度。fs为地函数的取样频率。
fc=fft(cory);
cm=abs(fc);
fl=(0:length(fc)-1)'fs/length(fc);
plot(fl,cm);%自相关函数的傅里叶变换:即功率谱密度。其中。cm是cory付里叶变换后的幅值。fl为fc的长度。
bnet2 = mk_bnet(dag, node_sizes);
seed = 0;
rand('state', seed);
%随机初始化网络参数
bnet2CPD{C} = tabular_CPD(bnet2, C);
bnet2CPD{R} = tabular_CPD(bnet2, R);
bnet2CPD{S} = tabular_CPD(bnet2, S);
bnet2CPD{W} = tabular_CPD(bnet2, W);
bnet3 = learn_params(bnet2, samples);
CPT3 = cell(1,N);
for i=1:N
s=struct(bnet3CPD{i}); % violate object privacy
欢迎分享,转载请注明来源:内存溢出
微信扫一扫
支付宝扫一扫
评论列表(0条)