ASK调制与解调代码.docVIP

ASK调制与解调代码 2ASK代码 close all; clear all; clc; t0=0.049999; % 0.004999 50个码元 0.049999 500个码元 0.099999 1000个码元 % 码元宽度 0.0001s， 每个码元100个采样点 ts=0.000001; fc=50000; %snr=0:1:20; snr=6; %信噪比，以dB来表示 fs=1/ts; df=0.3; %频率分辨率 t=[0:ts:t0]; %计算仿真的数字序列长度 n=size(t); N=n(2)/100; %载波信号生成，频率为fc c=cos(2*pi*fc*t); %本地载波信号 c_local=cos(2*pi*fc.*t); %生成消息信号（二进制数字序列） for i=1:N m2(i)=fix(2*rand(1)); end %将二进制数字序列转换为时域信号 for k=1:N for i=(100*k-99):(100*k) m2_time(i)=m2(k); end end %清零处理 for i=(100*(N/2)+1):(100*N) m2_time(i)=0; end %{ %绘制二进制数字基带序列（不归零码型） plot(t,m2_time) axis([0 0.001 -0.5 1.5]) grid on; title(二进制数字基带序列) %} ask2=m2_time.*c; % figure(1) subplot(3,1,1) plot(t,m2_time) axis([0 0.001 -0.5 3.5]) title(调制信号); grid on; subplot(3,1,2) plot(t,c) axis([0 0.001 -3.5 3.5]) title(载波信号); grid on; subplot(3,1,3) plot(t,ask2) axis([0 0.001 -3.5 3.5]) title(ASK已调信号); grid on; %} ask2_noise=awgn(ask2,snr,measured); %{ for j=1:length(snr) ask2_noise=awgn(ask2,snr(j),measured); %} % %绘制ASK信号波形（受噪声影响） figure(2) %subplot(4,1,4) plot(t,ask2_noise) axis([0 0.001 -2 2]) grid on; title(ask信号波形（加噪）) %} %本地载波与接收已调信号相乘，实现相干解调 ask2_dem1=ask2_noise.*c_local; %将相乘结果转换为频域 [ASK2_DEM1,ask2_dem1,df1]=fftseq(ask2_dem1,ts,df); ASK2_DEM1=ASK2_DEM1/fs; f=[0:df1:df1*(length(ask2_dem1)-1)]-fs/2; % %将已调信号转换为频域 [ASK2_NOISE,ask2_noise,df1]=fftseq(ask2_noise,ts,df); ASK2_NOISE=ASK2_NOISE/fs; %} %低通滤波器设计 f_cutoff=50000; %截止频率 n_cutoff=floor(f_cutoff/df1); H=zeros(size(f)); H(1:n_cutoff)=2*ones(1,n_cutoff); H(length(f)-n_cutoff+1:length(f))=2*ones(1,n_cutoff); %接收信号通过低通滤波器 ASK2_DEM2=H.*ASK2_DEM1; ask2_dem2=real(ifft(ASK2_DEM2))*fs; %{ %绘制抽样前信号的频谱 subplot(3,1,1) plot(t,m2_time(1:length(t))) axis([0 0.001 -1 2]) title(调制信号); grid on; subplot(3,1,2) plot(t,ask2_dem2(1:length(t))) axis([0 0.001 -1 2]) title(抽样前信号波形); grid on; %} figure(3) %subplot(3,1,3) plot(f,abs(fftshift(ASK2_NOISE))) axis([-150000 150000 0 0.001]) xlabel(频率); title(已调信号频谱); grid on; %} %抽样判决 k=1; for i=50:1