mbops.m

来自「人体脑信号波形检测与分析」· M 代码 · 共 30 行

nfft=9192;
    noverlap=512;
    dflag='mean';
    window=hamming(1024);
    [pxxb,fb]=psd(chb,nfft,sf,window,noverlap,dflag);
    [pxxo,fo]=psd(cho,nfft,sf,window,noverlap,dflag);
    loo=length(pxxo);
    fb=fb(1:loo);
    fo=fo(1:loo);
    pxb=pxxb(1:loo);
    pxo=pxxo(1:loo);
    
    set(gcf,'Name','Power Spectral');
    pxb=pxb/max(pxb);    
    pxo=pxo/max(pxo);
%    pxb=RDIVIDE(pxb,max(pxb));
 %   pxo=RDIVIDE(pxo,max(pxo));
    fb=fb(1:loo*0.2);%????????
    pxb=pxxb(1:loo*0.2);%????????
    fo=fo(1:loo*0.2);%????????    
    pxo=pxxo(1:loo*0.2);%????????
    plot(fb,pxb,'r');
    hold on;
    plot(fo,pxo,'b');
    hold off;
    title('Power Spectral');
    ylabel('Power Spectral');
    xlabel('frequency (Hz)');
    grid on;
    set(gca,'NextPlot','replacechildren');