kthspec.m

完整的幅相调制信号可k阶非线性功率谱分析

function [P,f]=kthspec(x,k,fs)
%周期图法计算信号k阶功率谱
%参数说明
%    P     输出功率谱幅值，行向量
%    f     输出频率，行向量
%    x     输入调制信号
%    k     非线性变换阶数
%    fs    信号采样率

%函数体
x=interp(x,k);  %k倍插值，调整采样率
x=x.^k;  %k次方
window=hanning(length(x));
nfft=4096;
[P,f]=periodogram(x,window,nfft,k*fs);
P=P';%变换为行向量
f=f';