返回   MATLAB中国论坛|MATLAB爱好者之家—不仅仅是MATLAB! > 人工智能|模式识别|机器学习|数字图像处理|SAR图像处理-不仅仅是MATLAB! > 人工神经网络(ANN)


人工神经网络(ANN) A discussion board for Artificial Neural Network(ANN)



回复
 
LinkBack 主题工具 显示模式
旧 2009-04-28, 08:55 PM   #1
初级会员
 
注册日期: 2009-03-13
年龄: 34
帖子: 10
感谢他人: 2
有 1 帖获得 1 感谢
声望力: 11
moneylancy 正向着好的方向发展
呲牙 关于循环谱

关于循环谱
大家好!!

我最近在做调制信号识别,其中需要用到循环谱。

我用了一篇论文里的程序,但是运行却有问题,还请大家指教。

这篇论文的名字是《Detection and identification of cyclostationary signals》。
程序如下:
function [Sx,alphao,fo]=autofam(x,fs,df,dalpha)

if nargin > 4 | nargin < 4
error('Wrong number of arguments.');
end

% Definition of Parameters %
Np = pow2(nextpow2(fs/df));
L = Np/4;
P = pow2(nextpow2(fs/dalpha/L));
N = P*L;

% Input Channelization %
if length(x) < N
x(N) = 0;
elseif length(x) > N
x = x(1:N);
end

NN = (P-1)*L+Np;
xx = x;
xx(NN) = 0;
xx = xx(;
X = zeros(Np,P);

for k=0-1
X(:,k+1) = xx(k*L+1:k*L+Np); end

% Windowing %
a = hamming(Np);
XW = diag(a)*X; %每段加窗
% XW = X; %I think this sentence should be kicked out.

% First FFT %
XF1 = fft(XW);
XF1 = fftshift(XF1);
XF1 = [XF1(:,P/2+1) XF1(:,1/2)];

% Downconversion %
E = zeros(Np,P);

for k = -Np/2:Np/2-1
for m = 0-1
E(k+Np/2+1,m+1) = exp(-i*2*pi*k*m*L/Np);
end
end

XD = XF1.*E;
XD = conj(XD'); %XD转置的复共轭

% Multiplication %
XM = zeros(P,Np^2);

for k = 1:Np
for q = 1:Np
XM(:,(k-1)*Np+q) = (XD(:,k).*conj(XD(:,q)));
end
end

% Second FFT %
XF2 = fft(XM);
XF2 = fftshift(XF2);
XF2 = [XF2(:,Np^2/2+1:Np^2) XF2(:,1:Np^2/2)];
XF2 = XF2(P/4:3*P/4,;
M = abs(XF2); %Absolute value and complex magnitude
alphao = (-1:1/N:1)*fs;
fo = (-.5:1/Np:.5)*fs;
Sx = zeros(Np+1,2*N+1);

for k1 = 1/2+1
for k2 = 1:Np^2
if rem(k2,Np) == 0
q = Np/2-1; %increase 1
else
q = rem(k2,Np)-Np/2-1; %increase 1
end
k = ceil(k2/Np)-Np/2-1; %increase 1
p = k1-P/4-1;
alpha = (k-q)/Np+(p-1)/L/P;
f = (k+q)/2/Np;
if alpha < -1 | alpha > 1
k2 = k2+1;
elseif f < -.5 | f > .5
k2 = k2+1;
elseif rem(k+q,2)==0
kk = 1+Np*(f+.5);
ll = 1+N*(alpha+1);
Sx(kk,ll) = M(k1,k2); %出现问题处,这里计算出的kk有时为非整数。为什么??M到Sx的对应是怎样的??
end
end
end
moneylancy 当前离线   回复时引用此帖
旧 2019-02-21, 03:17 PM   #2
初级会员
 
注册日期: 2019-02-21
帖子: 1
感谢他人: 0
有 0 帖获得 0 感谢
声望力: 0
wzq123 正向着好的方向发展
默认 回复: 关于循环谱

楼主,循环谱研究理论您还有印象吗?想要求教
wzq123 当前离线   回复时引用此帖
回复

书签

主题工具
显示模式

发帖规则
不可以发表新主题
可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码
Trackbacks are 禁用
Pingbacks are 禁用
Refbacks are 启用



所有时间均为北京时间。现在的时间是 07:15 PM


Powered by vBulletin
版权所有 ©2000 - 2019,Jelsoft Enterprises Ltd.

SEO by vBSEO ©2009, Crawlability, Inc.