通信工程专业 2ASK 2PSK 课程设计报告(1)

1、通信工程专业 2ask 2psk 课程设计报告(1) 2ask仿真程序 clc; clear all; close all; %信源 a=randint(1,15,2); t=0:0.001:0.999; m=a(ceil(15*t+0.01); subplot(511) plot(t,m); axis(0 1.2 -0.2 1.2); title(信源); %载波 f=150; carry=cos(2*pi*f*t); *sk调制 st=m.*carry; subplot(512); plot(t,st) axis(0 1.2 -1.2 1.2) title(2ask信号) %加高斯噪声 n

2、st=awgn(st,70); %解调部分 nst=nst.*carry; subplot(513) plot(t,nst) axis(0 1.2 -0.2 1.2); title(乘以相干载波后的信号) %低通滤波器设计 wp=2*pi*2*f*0.5; ws=2*pi*2*f*0.9; rp=2; as=45; n,wc=buttord(wp,ws,rp,as,s); b,a=butter(n,wc,s); %低通滤波 h=tf(b,a); %转换为传输函数 dst=lsim(h,nst,t); subplot(514) plot(t,dst) axis(0 1.2 -0.2 1.2);

3、title(经过低通滤波器后的信号); %判决器 k=0.25; pdst=1*(dst0.25); subplot(515) plot(t,pdst) axis(0 1.2 -0.2 1.2); title(经过抽样判决后的信号) %频谱观察 %调制信号频谱 t=t(end); df=1/t; 成绩： 批阅老师 ： n=length(st); f=(-n/2:n/2-1)*df; sf=fftshift(abs(fft(st); figure(2) subplot(411) plot(f,sf) title(调制信号频谱) %信源频谱 mf=fftshift(abs(fft(m); subp

4、lot(412) plot(f,mf) title(信源频谱) % 乘以相干载波后的频谱 mmf=fftshift(abs(fft(nst); subplot(413) plot(f,mmf) title(乘以相干载波后的频谱) %经过低通滤波后的频谱 dmf=fftshift(abs(fft(pdst); subplot(414) plot(f,dmf) title(经过低通滤波后的频谱); 2psk仿真程序 clear all; close all; clc; max=10 g=zeros(1,max); g=randint(1,max);%长度为max的随机二进制序列 cp=;mod1=

5、;f=2*2*pi;t=0:2*pi/199:2*pi; for n=1:length(g); if g(n)=0; a=zeros(1,200);%每个值200个点 else g(n)=1; a=ones(1,200); end cp=cp a; %s(t),码元宽度200 c=cos(f*t);%载波信号 mod1=mod1 c;%与s(t)等长的载波信号,变为矩阵形式 end figure(1);subplot(4,2,1);plot(cp);grid on; axis(0 200*length(g) -2 2);title(二进制信号序列); cm=;mod=; for n=1:len

6、gth(g); if g(n)=0; b=ones(1,200);%每个值200个点 c=cos(f*t); %载波信号 else g(n)=1; b=ones(1,200); c=cos(f*t+pi); %载波信号 end cm=cm b; %s(t),码元宽度200 成绩： 批阅老师 ： mod=mod c; %与s(t)等长的载波信号 end tiaoz=cm.*mod;%e(t)调制 figure(1);subplot(4,2,2);plot(tiaoz);grid on; axis(0 200*length(g) -2 2);title(2psk调制信号); figure(2);s

7、ubplot(4,2,1);plot(abs(fft(cp); axis(0 200*length(g) 0 400);title(原始信号频谱); figure(2);subplot(4,2,2);plot(abs(fft(tiaoz); axis(0 200*length(g) 0 400);title(2psk信号频谱); %带有高斯白噪声的信道 tz=awgn(tiaoz,10);%信号tiaoz中加入白噪声，信噪比为10 figure(1);subplot(4,2,3);plot(tz);grid on axis(0 200*length(g) -2 2);title(通过高斯白噪声

8、信道后的信号); figure(2);subplot(4,2,3);plot(abs(fft(tz); axis(0 200*length(g) 0 400);title(加入白噪声的2psk信号频谱); jiet=2*mod1.*tz;%同步解调 figure(1);subplot(4,2,4);plot(jiet);grid on axis(0 200*length(g) -2 2);title(相乘后信号波形) figure(2);subplot(4,2,4);plot(abs(fft(jiet); axis(0 200*length(g) 0 400);title(相乘后信号频谱);

9、%低通滤波器 fp=500;fs=700;rp=3;rs=20;fn=11025; ws=fs/(fn/2); wp=fp/(fn/2);%计算归一化角频率 n,wn=buttord(wp,ws,rp,rs);%计算阶数和截止频率 b,a=butter(n,wn);%计算h(z） jt=filter(b,a,jiet); figure(1);subplot(4,2,5);plot(jt);grid on axis(0 200*length(g) -2 2);title(经低通滤波器后信号波形) figure(2);subplot(4,2,5);plot(abs(fft(jt); axis(0

10、200*length(g) 0 400);title(经低通滤波器后信号频谱); %抽样判决 for m=1:200*length(g); if jt(m) jt(m)=1; else jt(m)=0; jt(m)=0; end end figure(1);subplot(4,2,6);plot(jt);grid on axis(0 200*length 信源10.50010-1010.50010.50010.5000.20.40.60.811.20.20.40.60.8经过抽样判决后的信号11.20.20.40.60.8经过低通滤波器后的信号11.20.20.40.60.8乘以相干载波后的信

11、号11.20.20.40.62ask信号0.811.2调制信号频谱4002000-600500-400-2000信源频谱2004006000-6004002000-6004002000-600-400-2000200乘以相干载波后的频谱400600-400-2000200经过低通滤波后的频谱400600-400-2000200400600 2psk仿真结果 成绩： 批阅老师 ： 二进制信号序列20-20500100015002000通过高斯白噪声信道后的信号20-20500100015002000经低通滤波器后信号波形20-20500100015002000原始信号频谱40020000500100015002000加入白噪声的2psk信号频谱40020000500100015002000经低通滤波器后信号频谱40020000500100015002000 2psk调制信号20-205001