利用matlab对am，dsb，ssb，ask，fsk，bpsk信号进行正交调制解调仿真

clear echo on t0=2; ts=0.001; fc=50;%载波频率 fs=1/ts;%采样频率 df=0.3; t=[-t0/2:ts:t0/2]; m=cos(2*pi*2*t); c=cos(2*pi*fc*t); b=sin(2*pi*fc*t); Il=m;Ql=imag(hilbert(m));%下边带信号的I路和Q路 p=Il.*c+Ql.*b; % 下边带信号的表达式 u=awgn(p,30); %下边带信号加入噪声 Iu=m;Qu=-imag(hilbert(m));%上边带信号的I路和Q路 q=Iu.*c+Qu.*b; % 上边带信号的表达式 v=awgn(q,30); %上边带信号加入噪声 [M,m,df1]=fftseq(m,ts,df); M=M/fs; [U,u,df1]=fftseq(u,ts,df); U=U/fs; [V,v,df1]=fftseq(v,ts,df); V=V/fs; f=[0:df1:df1*(length(m)-1)]-fs/2; figure(1) subplot(2,2,1) plot(t,m(1:length(t))); xlabel('时间');title('未调信号') subplot(2,2,2) plot(t,c(1:length(t)));axis([-0.1,0.1,-1.5,1.5]) xlabel('时间');title('载波') subplot(2,2,3) plot(t,u(1:length(t))); axis([-0.2,0.2,-1,1.2]);xlabel('时间'); title('上边带已调信号') subplot(2,2,4) plot(t,v(1:length(t))); axis([-0.2,0.2,-1,1.2]);xlabel('时间'); title('下边带已调信号'); figure(2) subplot(2,1,1); plot(f,abs(fftshift(M))) xlabel('频率');title('未调信号的频谱') subplot(2,1,1) plot(f,abs(fftshift(U))) title('上边带已调信号的频谱');xlabel('频率'); subplot(2,1,2); plot(f,abs(fftshift(V))) title('下边带已调信号的频谱');xlabel('频率'); %上边带信号的解调 figure(3) [ar,br]=butter(3,[48 52]/500); r=filter(ar,br,v); r=r(1:length(c)); xui=r.*cos(2*pi*fc*t+pi);%上边带信号I路信号 [aui,bui]=butter(3,2/500);%经过低通滤波器 yui=filter(aui,bui,xui); xuq=-r.*sin(2*pi*fc*t);%上边带信号Q路信号 [auq,buq]=butter(3,2/500);%经过低通滤波器 yuq=filter(auq,buq,xuq); %根据算法I路输出即为调制信号 subplot(2,1,1) plot(t,m(1:length(t))); title('基带信号') subplot(212) plot(t,yui);title('解调信号') %下边带信号的解调 figure(4) [au,bu]=butter(3,[48 52]/500); s=filter(ar,br,u); s=s(1:length(c)); xli=s.*cos(2*pi*fc*t+pi);%下边带信号I路信号 [ali,bli]=butter(3,2/500);%经过低通滤波器 yli=filter(ali,bli,xli); xlq=-s.*sin(2*pi*fc*t);%下边带信号Q路信号 [alq,blq]=butter(3,2/500);%经过低通滤波器 ylq=filter(alq,blq,xlq); %根据算法I路输出即为调制信号 subplot(2,1,1) plot(t,m(1:length(t))); title('基带信号') subplot(212) plot(t,yli);title('解调信号')