基于MATLAB实现的rfid（2ASK信号的调制与解调）可以直接出仿真结果+使用说明文档.rar

%2ASK信号的调制与解调源代码 close all clear all clc bit=[1 1 0 1 0 1 0 0];%g为比特流 N=100; %码元速率 fc=1000;%载波频率 swr=20;%信噪比 M=8;%码元个数 bitrate=100; fs=N*bitrate; t=linspace(0,1/bitrate,N); ts=1/fs; %载波信号 c=cos(2*pi*fc*t); figure(1) plot(c,'r'); title('Carrier signal');%载波信号 y1=fftshift(abs(fft(c))); figure(2) plot(y1,'r'); title('Carrier spectrum');%载波信号频谱 axis([-100,200,0,100]); %数字基带信号产生 bits=[]; sig=[]; s=[]; for i=1:length(bit)%%length(m)表示信号m的长度 if bit(i)==0 bits=zeros(1,N); else bits=ones(1,N); end sig=[sig,bits]; %数字基带信号 s=[s,cos(2*pi*fc*t).*bits];%已调信号 end figure(3) plot(sig,'r'); title('Digital baseband signal');%数字基带信号波形 axis([0,900,-0.5,2]); y2=fftshift(abs(fft(sig))); figure(4) plot(y2,'r'); title('Digital baseband spectrum');%数字基带信号频谱 axis([0,900,0,450]); %基带信号功率谱 nfft=2048; cxn=xcorr(sig,'unbiased'); %计算序列的自相关函数 CXk=fft(cxn,nfft); Pxx=abs(CXk); index=0:round(nfft/2-1); k=index*fs/nfft; plot_Pxx=10*log10(Pxx(index+1)); figure(5) plot(k,plot_Pxx); title('Baseband signal power spectrum');%数字基带信号功率谱 %已调信号的波形以及频谱 figure(6) plot(s,'k'); title('Modulated signal');%已调信号波形 y3=fftshift(abs(fft(s))); figure(7) plot(y3,'k'); title('Modulated signal spectrum');%已调信号频谱 axis([-100,1000,0,250]) %加噪声 figure(8) y4= awgn(s,swr); plot(y4,'r'); title('Modulated signal with noise');%加噪声后的已调信号 %求加噪声后的功率谱密度 nfft=2048; cxn=xcorr(y4,'unbiased'); %计算序列的自相关函数 CXk=fft(cxn,nfft); Pxx=abs(CXk); index=0:round(nfft/2-1); k=index*fs/nfft; plot_Pxxn=10*log10(Pxx(index+1)); figure(9) plot(k,plot_Pxxn); title('Baseband signal power spectrum with noise');%加噪声后数字基带信号功率谱 %带通滤波 a=[700,1300];c=0.1;d=35; Wp=a/(fs/2); %%若为低通，a 为一个频率值，若为带通a 为一个频率范围[a1 a2] Rp=c; %% c 为通带最大衰减分贝 Rs=d; %% d 为阻带最小衰减分贝 [b,a]=ellip(4,Rp,Rs,Wp); s1=filter(b,a,s); figure(10) plot(s1,'r'); title('Digital baseband signal after bandpass filtering');%带通滤波后的数字基带信号 %带通滤波后信号功率谱密度 nfft=2048; cxn=xcorr(s1,'unbiased'); %计算序列的自相关函数 CXk=fft(cxn,nfft); Pxx=abs(CXk); index=0:round(nfft/2-1); k=index*fs/nfft; plot_Pxx=10*log10(Pxx(index+1)); figure(11) plot(k,plot_Pxx); title('Power spectrum after bandpass filter');%带通滤波后信号功率谱 %抽样判决（解调过程） sigout=[]; for i=1:M if s1(N*(i-1)+N/2)>0.5 bits=ones(1,N); else bits=zeros(1,N); end sigout=[sigout,bits]; end figure(12) plot(sigout,'b'); title('Demodulated signal after sampling and decision');%抽样判决后的解调信号波形 axis([0,900,-0.5,1.5]); %解调信号功率谱密度 nfft=2048; cxn=xcorr(sigout,'unbiased'); %计算序列的自相关函数 CXk=fft(cxn,nfft); Pxx=abs(CXk); index=0:round(nfft/2-1); k=index*fs/nfft; plot_Pxx=10*log10(Pxx(index+1)); figure(13) plot(k,plot_Pxx); title('Power spectrum of demodulated signal after filtering');%滤波后的解调信号功率谱 %最后将原基带信号和解调后的基带信号统一绘图作比较 figure(14) subplot(2,1,1); plot(sig,'k'); title('The source digital baseband signal');%原数字基带信号 axis([0,1000,-0.5,1.5]); subplot(2,1,2); plot(sigout,'k'); title('Demodulated signal');%解调信号 axis([0,1000,-0.5,1.5]);