|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
对平稳随机信号的自相关序列做FFT变换,求出的应该是功率谱,可是下面的代码运行出来功率谱中怎么会有负数呢?求教各位老师~
clear all;
clc;
fc=4;
D=randint(1,8);
N=length(D);
t=0:0.01:7.99;
d1=2*D-1;
D1=[];
x=ones(1,100);%100个1
fai=zeros(1,N); % 初始化相位常数
fai(1)=0;
for n=1:N,
D1=[D1,d1(n)*x];
end %产生一个+1,-1的矩形随机方波
fai=zeros(1,N);
fai(1)=0; % 将初始相位定为0
%求相位常数fai
for k=2:N*100
if(D1(k)==D1(k-1))
fai(k)=fai(k-1);
else
fai(k)=fai(k-1)+(k-1)*pi;
end
end
for i=1:N*100
if(D1(i)==1)
f(i)=fc+1/4;
else f(i)=fc-1/4;
end
end
m=cos(2*pi*f.*t).*cos(fai)-sin(2*pi*f.*t).*sin(fai); %MSK表达式三角变化得来的式子
rx = xcorr(m,'unbiased');
plot(real(fft(rx)))
|
|