m序列编码成的BPSK信号仿真

MTALB对BPSK信号的仿真程序！编码是m序列来的，频谱仿真！

clc;
clear all;
close all;

connections1=[1 0 1 0];
registers1=[1 1 1 1];
sequence1=ssmlsrs1_Fab（connections1registers1）; %（1/0）

Tb1=0.3e-3; %编码脉冲宽度，0.3ms

Tf=6e-3;     %发射脉冲宽度，6ms
Tz=0.02e-3;  %子码宽度，0.02ms
Fz=1/Tz; 
f1=300000;   %载频300KHz

fs=12e6;     %采样率12M%%%%%%%%%%%%%%07
timestep=1/fs; %采样间隔
t=0:timestep:Tb1-timestep;

fd=0;
fc=f1+fd; 
y=sin（2*pi*fc*t）; %基本码型 sin函数作为载波

N2=fs/fc;       %fd=0时一个sin周期上的采样点数
N1=N2*6;     %一个码元填充6个周期的sin，故采样点数为6倍

for ii=1:length（sequence1） 
    for jj=1:N1
        seq1（（ii-1）*N1+jj）=sequence1（ii）;      %seq1:采样后的子码（M序列）
    end
end

N3=N1*15;  %0.3ms采样出3600个点，40*6*15

for ii=1:N3
     if seq1（ii）==1
             y1（ii）=y（ii）;
     else  y1（ii）=-y（ii）;    %y1:基本码型经m序列调制之后所需的宽带信号！
     end
end

figure（1）;
subplot（311）;
plot（seq1）;            %M序列
subplot（312）;
plot（y）;                   %y为sin数据
subplot（313）;
plot（y1）;                %y1为0.3ms数据

y11=fft（y1）;
% y11=fftshift（y11）;
figure（2）;
plot（abs（y11））;       %宽带信号的频谱图

% [YI]=max（abs（y11））;
% ff=fs*I/N3               %载波频率估计


Z1=0;
Y=y1‘;
PP=1024;%FFT点数
TT=20;%谱的平均次数
a3=-1024000:4000:1024000;%循环频率
ff=-（1024000-2000）:2000:1024000-2000;%频率
Z=zeros（length（ff）length（a3））;%幅度
% DD=10*log10（var（Y））;%信号功率
% noise=wgn（100000015）;%产生噪声
% noise=randn（64001）;%产生噪声
for T1=1:TT
    u=fft（Y（（351+（T1-1）*350）:（PP+T1*350）））;%信号频谱
%     u=zeros（PP1）;
%     v=fft（noise（（151+（T1-1）*150）:（PP+T1*150）））;%噪声频谱
    v=zeros（PP1）;
    UV=u+v;uv=zeros（PP1）;
    uv（1:（PP/2）1）=UV（（PP/2+1）:PP）;uv（（PP/2+1）:PP1）=UV（1:（PP/2））;%平移到正负频率点上
    uv=repmat（uv21）;
    z=zeros（length（ff）length（a3））;
    for k=1:length（a3）
        for i=1:length（ff）
             if （（（0.0005*（ff（i）-a3（k）/2）+513）<=0）|（（0.0005*（ff（i）+a3（k）/2）+513）<=0））%考虑超过范围时，利用FFT频谱的周期性平移
                if （（0.0005*（ff（i）-a3（k）/2）+513）<=0）
                    z（ik）=uv（0.0005*（ff（i）-a3（k）/2）+513+1024）*conj（uv（0.0005*（ff（i）+a3（k）/2）+513））;
                else 
                    z（ik）=uv（0.0005*（ff（i）-a3（k）/2）+513）*conj（uv（0.0005*（ff（i）+a3（k）/2）+513+1024））;
                end
                Z（ik）=Z（ik）+（1/TT）*z（ik）;
            else
                z（ik）=uv（0.0005*（ff（i）-a3（k）/2）+513）*conj（uv（0.0005*（ff（i）+a3（k）/2）+513））;
                Z（ik）=Z（ik）+（1/TT）*z（ik）;
            end
        end  
    end
 Z=abs（Z）;
Z1=Z1+Z/TT;
end

[xy]=meshgrid（a3ff）;

mesh（xyabs（Z1））;  %画三维图
xlabel（‘循环频率（Hz）‘）;
ylabel（‘频率（Hz）‘）;
zlabel（‘幅度‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        469  2012-08-16 10:12  ssmlsrs1_Fab.m

     文件       2745  2013-03-08 20:54  kuandaixinhao_sin_xunhuanpuxiangguan.m

----------- ---------  ---------- -----  ----

                 3214                    2