askmod.m

无线通信FPGA设计 书的附带代码

clear all
global dt t f df N T %全局变量
close all 
N=2^20; 
dt=0.001; %ms 
df=1/(N*dt); %KHz 
T=N*dt; %截短时间 
Bs=N*df/2; %系统带宽 
t=linspace(-T/2,T/2,N); %时域横坐标 
f=linspace(-Bs,Bs,N)+eps; %频域横坐标 
fm=1; %模拟基带信号频率为1kHz 
fc=6; %正弦载波信号频率为6kHz 
mt=cos(2*pi*fm*t); %模拟基带信号 
m(mt>0)=1;  %发送比特
m(mt<0)=0;
c=cos(2*pi*fc*t); %正弦载波信号 
s=m.*c; %DSB-SC AM已调信号 
S=T2F(s); 
n=awgn(s,20); %AWGN信道噪声
r=s+n; %接受信号 
y=r.*c; 
Lpfil=LPF(fc); %低通滤波器 
Y=T2F(y); 
yrt=real(F2T(Lpfil.*Y)); 
yr=yrt-mean(yrt);
yr(yr>0)=1;  %发送比特
yr(yr<0)=0;
subplot(3,1,1) 
plot(t,m,'LineWidth',1.2)
title('输入信号');
axis([-1,+1,-0.2,1.2*max(m)]) 
xlabel('t (ms)')
ylabel('s(t) (V)') 
subplot(3,1,2) 
plot(t,s,'LineWidth',1.2) 
title('2ASK调制后波形');
axis([-1,+1,1.2*min(c),1.2*max(c)]) 
xlabel('t (ms)')
ylabel('s(t) (V)') 
subplot(3,1,3)
plot(t,yr,'LineWidth',1.2) 
axis([-1,+1,-0.2,1.2*max(yr)])
title(['2ASK解调后波形']) 
xlabel('t (ms)') 
ylabel('s(t) (V)')