关于伪随机信号功率谱的问题
在沈善德书《电力系统辨识》中p39有伪随机序列的功率密度谱曲线，我是竭尽我所能了，就是得不到他的那种形状啊！在此请教此方面的高手指点。伪随机序列我是由程序genepncode生成的，程序如下：
function =genpncode(sm,a,ts,num)
% 产生长度等于 pow2(sm)-1 m序列，序列的幅值为a，Y为输出的M序列，m为M序列的周期等于 pow2(sm)-1
l=pow2(sm)-1;
num=l*num;
n=;
if sm<3
errordlg('错误:sm不应小于3');
return;
end
if sm>10
errordlg('错误:sm不应大于10');
return;
end
for i=1:num
temp=0;
% 选择多项式
if sm==3
temp=n(1)+n(3);
end
if sm==4
temp=n(1)+n(4);
end
if sm==5
temp=n(2)+n(5);
end
if sm==6
temp=n(1)+n(6);
end
if sm==7
temp=n(3)+n(7);
end
if sm==8
temp=n(2)+n(3)+n(4)+n(8);
end
if sm==9
temp=n(9)+n(4);
end
if sm==10
temp=n(10)+n(3);
end
% mod(2)
temp=mod(temp,2);
for j=sm:-1:2
n(j)=n(j-1);
end
n(1)=temp;
y(i)=n(sm);
if y(i)==0
y(i)=-a;
else
y(i)=a;
end
end
Y=y;
m=l;
请高人对我的分析谱的程序指点下：
clc;
clear all;
sm=7; %寄存器个数
a=1; %幅值
ts=1/12.7; %间隔时间
%n=2; %产生n个周波的伪随机信号
=genpncode(sm,a,ts,1); %得到M序列为u，序列的周期为m，M序列的幅值为a，寄存器的个数为sm
%l=length(u)
%自己加的几个语句，作用是将每个采样点重复K次再频谱分析
x=u;
w1=u;
k=10;
n=length(u);
for i=1:n
for j=0:(k-1)
x(1+k*(i-1)+j)=w1(i);
end
end
%间接法，现求取自相关函数，再对自相关函数fft，求取功率谱===========================
fs=127;
N=length(u);
Nfft=256;
maxlags=100;
=xcorr(u,maxlags,'unbiased'); %计算序列的自相关函数
figure(1)
plot(lags/fs,cxn)
CXk=fft(cxn,Nfft);
pxx4=abs(CXk);
f=(0:(Nfft/2-1))*fs/Nfft;
plot_Pxx=10*log10(pxx4(1:(Nfft/2)));
figure(2)
plot(f,plot_Pxx);
谢谢了先！！！
页:
[1]