ask.m

对ASK调制解调进行了仿真,包括程序仿真和模型仿真

n=1:8192;
m=1:128; x(n)=randint(1,8192,2); x=[x(n)]';
y(n)=zeros(1,8192); z(m)=zeros(1,128);
for n=1:8192 
    for m=1:128
        if n==64*m-63
            z(m)=x(n);
            if m==ceil(n/64)
                y([(64*m-63):(64*m)]')=z(m);
            end
        end
    end
end
n=1:8192; rm2=y(n);
x2=rm2;
n=[1:(2^13)];
x1=cos(n.*1e9*2*pi/4e9);
x=x1.* x2;
b=blackman(2^13);
x=b'.* x;
x3=[ones(1,64) zeros(1,8128)];
y1=x(1:(2^13)); y4=x1.*x3;
Y1=fft(y1,(2^13));magY1=abs(Y1(1:1:(2^12)+1))/(200);
Y4=fft(y4,(2^13));magY4=abs(Y4(1:1:(2^12)+1))/(37);
k1=0:(2^12);w1=(2*pi/(2^13))*k1;
u=(2*w1/pi)*1e9 ;
figure(1)
subplot(2,1,1);
plot(u,magY1,'b',u,magY4,'r');grid
title('ASK'); axis([4e8,1.6e9,0,1.1])
x2=b'.* x2;
y2=x2(1:(2^13));
Y2=fft(y2,(2^13)); magY2=abs(Y2(1:1:(2^12)+1))/(200)+eps;
k1=0:(2^12); w1=(2*pi/(2^13))*k1;
u=(2*w1/pi)*1e9;
Y3=fft(x3,(2^13)); magY3=abs(Y3(1:1:(2^12)+1))/(35)+eps;
subplot(2,1,2); semilogy(u,magY2,'b',u,magY3,'r'); grid
title('ASKr-modulation'); axis([0,1.2e9,3e-2,3])
figure(2)
subplot(2,1,1); plot(n,x2); title('ASKr');
axis([300,800,-0.02,0.05]); grid
subplot(2,1,2); plot(n,x);
axis([300,800,-0.1,0.1]); grid