lc21.asv

在matlab下输入guide命令后打开main_main.gig 实现卷积编码

%本程序中时间单位是微秒
%频率单位为MHz
%码元速率单位是Mb/s

global dt df t f N
close all

N=2^13;   %采样点数
L=64;    %每码元的采样点数
M=N/L     %码元数
Rb=2;
Ts=0.5		%码元宽度是0.5us
dt=Ts/L;        
df=1/(N*dt)     %MHz
T=N*dt          %截短时间
Bs=N*df/2       %系统带宽
n0=0.001;

t=[-T/2+dt/2:dt:T/2]; %时域横坐标
f=[-Bs+df/2:df:Bs];   %频域横坐标

figure(1)
set(1,'Position',[10,50,400,400])
%设定窗口位置及大小 
figure(2)
set(2,'Position',[410,50,400,400])
%设定窗口位置及大小 
EP=zeros(size(f))+eps;
EPR=zeros(size(f))+eps;
for ii=1:30               
    a=round(rand(1,M));
    s=zeros(1,N);  
    for ii=1:0.5*L,
        s(ii+[0:M-1]*L)=a; ,
    end%产生单极性归零码
    
    sita=n0*Bs; %信道中噪声功率
    n_ch=sqrt(sita)*randn(size(t)); %信道噪声
    sr=s+n_ch;       %接收信号
    
    S=t2f(s);    %S是s的傅氏变换
    P=S.*conj(S)/T;    %功率谱
    SR=t2f(sr);    %接收信号傅氏变换
    PR=SR.*conj(SR)/T;    %功率谱
    EP=(EP*(ii-1)+P+eps)/ii;
    EPR=(EPR*(ii-1)+PR+eps)/ii;
end


figure(1)
subplot(2,1,1)
plot(t,s,'b')
grid
axis([-3,+3,-1.5,1.5'])
xlabel('t (ms)')
ylabel('s(t) (V)')

figure(1)
subplot(2,1,2)
aa=30+10*log10(EP+eps);    %加eps以避免除以零
plot(f,aa,'b');
grid
axis([-10,+10,-50,50])
xlabel('f  (MHz)')
ylabel('Ps(f)   (dBm/MHz)')


figure(2)
subplot(2,1,1)
plot(t,sr,'b')
grid
axis([-3,+3,-1.5,1.5'])
xlabel('t (ms)')
ylabel('s(t) (V)')

figure(2)
subplot(2,1,2)
aa=30+10*log10(EPR+eps);    %加eps以避免除以零
plot(f,aa,'b');
grid
axis([-10,+10,-50,50])
xlabel('f  (MHz)')
ylabel('Ps(f)   (dBm/MHz)')