ndy.m

认知无线电系统中多载波信号带外干扰的抑制

global dt df t f N%全局变量
close all
k=input('取样样点=2^k,k>=[14]')
if k==[],k=14;end
g1=[1,0,0,0,0,1,1;0,1,0,0,1,0,1;0,0,1,0,1,1,1;0,0,0,1,1,1,0]%线组分组码生成矩阵
g2
N=2^k;				%采样点数
L=64;					%每码元的采样点数
M=N/L;				%码元数
Rb=2;					%码速率是2Mb/s
Ts=1/Rb;				%码元间隔
dt=Ts/L;				%时域采样间隔
df=1/(N*dt);		%频域采样间隔
T=N*dt;				%截短时间
Bs=N*df/2;			%系统带宽
t=[-T/2+dt/2:dt:T/2];%时域横坐标
f=[-Bs+df/2:df:Bs];		%频域横坐标
EP=zeros(size(f));
for jj=1:100
a=round(rand(1,M));	%码元矢量
s=zeros(1,N);
s1=zeros(1,N/4*7);
for qq=
for ii=1:L
   s(ii+[0:M-1]*L)=a;
end						%产生NRZ码
S=t2f(s);
P=S.*conj(S)/T;		%功率谱
EP=(EP*(ii-1)+P)/ii; %累计平均
aa=30+10*log10(EP+eps);	%加eps以避免除以零
end

figure(1)
set (1,'position',[10,50,300,200]);%设定窗口位置及大小
figure(2)
set(2,'position',[350,50,300,200]);%设定窗口位置及大小
figure (1)
plot(f,aa,'g');
grid
axis([-10,+10,-50,+50])
xlabel('f(MHZ)')
ylabel('Ps(f)(MHZ)')
figure(2)
plot(t,s,'r')
grid
axis([-10,+10,-1,+2])
xlabel('t(us)')
ylabel('s(t)(V)')