标准实验报告(四)

1、实 验 报 告（二）学生姓名：王超楠学号： 2013020904011指导教师：廖红舒/ 张花国实验地点：科研二号楼b453 实验时间：周二晚一、实验室名称：信息对抗系统专业实验室二、实验项目名称： 典型数字通信信号调制识别三、 实验学时： 4 学时实验原理：bpsk、 qpsk、 fsk 、 ask 数字通信信号具有不同的频谱特征，如 bpsk 平方和四次方后的傅立叶变换出现单根离散谱线，qpsk 四次方后才有单根离散谱线等。 因此针对这些信号特征， 可通过设置特征门限区分不同通信信号， 达到 信号调制识别目的。因此可让学生通过实际上机matlab 编程实验，对上述通信信号的特征进行仿真验证

2、，加深理解不同通信信号的调制识别方法。五、实验目的：利用 matlab 软件编程提取数字通信信号的频谱、二次方谱、四次方谱等特征， 同时使用决策树的分类方法进行调制信号的识别。 让学生通过实际上机实验，加深理解不同数字通信信号的特点。六、实验内容：（ 1） 把上次实验产生的 bpsk、 qpsk、 2ask 、 2fsk 信号分别画出它们的频谱、二次方谱和四次方谱，描述不同信号谱特征之间的差异。（ 2） 分别用介绍的两种方法编程提取这四种信号的频谱、二次方谱、四次方谱的谱峰个数（根据每种信号的特点分别提取频谱或二次方谱或四次方谱的谱峰个数）。信噪比从0db变化到10db （即0: 2: 10）

3、,分 别画出bpsk、qpsk、2ask、2fsk用两种方法获得谱峰个数特征随 每个信噪比的变化曲线，然后设定区分门限。比较这两种方法的优缺点。（ 3） 选择其中一种方法提取谱峰特征， 设计识别决策树， 并编写程序实现。画出bpsk、qpsk、2ask、2fsk信号的正确识另率随信噪比 (0: 2:10)的变化曲线。说明在多高信噪比下能保证每种信号的正确识别率 均在90%以上。七、实验器材(设备、元器件)：计算机、matlab计算机仿真软件八、实验步骤：1、在编写的信号源基础上，根据实验内容提取信号特征并进行调制识别。(a) bpsk、qpsk、2ask、2fsk频谱特征提取程序(b) bps

4、k、qpsk、2ask、2fsk数字通信信号调制识另程序九、实验数据及结果分析(1) ask信号的频谱、二次方谱、四次方谱bpsk信号的频谱、二次方谱、四次方谱豆 j % e iqpsk信号的频谱、二次方谱、四次方谱fsk信号的频谱、二次方谱、四次方谱程序如下：1 .askfc=300;%载波频率fs=2000;痛样频率fb=200;%tb比特周期符号速度fd=1/2*fbfd=200;t=2;a=randi(0,1,1,t*fs);%产生二进制码元%根升余弦成型滤波器nt=60;r=0.3;rate=fs/fb;%ii采卞率=fs/ 码速率delay=nt/(2*rate);t=0:1/fs

5、:t-1/fs;x=exp(j*2*pi*fc*t).*a;ask=rcosflt(x,fd,fs,sqrt ,r,delay);%m 形成型figure(1)subplot(411)plot(real(ask);axis(0 2000 -1 1);grid on;title( ask 时域波形)；subplot(412)plot(fftshift(abs(fft(ask);grid on;title( ask的频域波形)；subplot(413)plot(fftshift(abs(fft(ask.a2);grid on ;title( ask 的二次方谱 );subplot(414)plot

6、(fftshift(abs(fft(ask.a4);grid on ;title( ask 的四次方谱 );2.bpskfc=200;%载波频率fs=2000;%采样频率fb=200; %tb 比特周期fd=200;%符号速度t=2; %时间长2sm=2; %二元%产生产生 0 ， 1 序列或符号序列m=randi(0,1,1,fs);%相位映射为复基带信号a=pskmod(m,m);%基带消息t=0:length(a)-1;x=cos(2*pi*fc*t+angle(a);%生成bpsk%根升余弦成型滤波器nt=60;r=0.3;rate=fs/fb;%过采样率=fs/ 码速率delay=n

7、t/(2*rate);xbpsk=rcosflt(x,fd,fs,sqrt ,r,delay);%波形成型%画图figure(1)subplot(411)plot(xbpsk);axis(0 2000 -1 1);grid on ;title( bpsk 时域波形图 );subplot(412)plot(fftshift(abs(fft(xbpsk);grid on ;title( bpsk 频域波形图 );subplot(413)plot(fftshift(abs(fft(xbpsk.a2);title( bpsk 信号二次方谱 );subplot(414)plot(fftshift(abs

8、(fft(xbpsk.a4);title( bpsk 信号四次方谱 );3.qpskclear all ;close all ;clc;fc=200;%载波频率fs=2000;%采样频率fb=100; %tb 比特周期fd=200;%符号速度t=2; %时间长2sm=4; %四元n=t*fs;r=fs/fd;filtorder = 60;%滤波器阶数60a=0.3;%滚降因子0.3delay = filtorder/(r*2);h=rcosfir(a, delay, r,1,sqrt );h=sqrt(r)*h/norm(h);s=randsrc(n,1,0:m-1);p=pskmod(s,m

9、);x_base=zeros(r,n);x_base(1,:)=p.;x_base=x_base(:);ss=conv(h.,x_base);x_qpsk=x_base.*exp(j*(2*pi*fc/fs*0:length(x_base)-1.+rand*2*pi);% figure(1)subplot(411) plot(real(x_qpsk);axis(0 500 -1 1);title( qpsk 时域波形图 ); grid on ;subplot(412)plot(fftshift(abs(fft(x_qpsk);title( qpsk 频域波形图 ); grid on ;subp

10、lot(413)plot(fftshift(abs(fft(x_qpsk.a2);title( qpsk 二次方谱 );grid on ;subplot(414)plot(fftshift(abs(fft(x_qpsk.a4);title( qpsk 四次方谱 );grid on ;4.fsk% 因为调频指数为0.6400hz ，可以计算出 f1 f2; fc=(f1+f2)/2h=(f2-f1)/fb%载波频率%采样频率%tb 比特周期fc=400;fs=1200;fb=200;fd=200;f2=460;f1=340;t=2;% 生成 an bn 随机序列%产生二进制码元a=randi(0

11、,1,1,t*fs);n=length(a);b=zeros(1,n);for i=1:nb(i)=1-a(i);end%根升余弦成型滤波器nt=60;r=0.3;rate=fs/fb;%过采样率=fs/ 码速率delay=nt/(2*rate);t=0:1/fs:t-1/fs;c=a.*exp(j*2*pi*f1*t)+b.*exp(j*2*pi*f2*t);fskk=rcosflt(c,fd,fs,sqrt ,r,delay);%波形成型 figure(1)subplot(411)plot(real(fskk);axis(0 2000 -1 1);grid on ;title( fsk 时

12、域波形 );subplot(412)plot(fftshift(abs(fft(fskk);grid on ;title( fsk 的频域波形 );subplot(413)plot(fftshift(abs(fft(fskk.a2);grid on ;title( fsk 二次方谱)；subplot(414)plot(fftshift(abs(fft(fskk.a4);grid on;title( fsk 四次方谱);(2)谱峰个数的提取方法(1)以频谱或二次方谱或四次方谱中的最大值为参考值，设置一个门限，搜索大于该门限以上的谱峰个数方法(2)以频谱或二次方谱或四次方谱 p中的最大值为参考值，

13、记录最大 值的位置为pos,选取最大值左右两侧第 n个点的谱的幅度值，计算比值p p(pos)p(pos- n) p(pos - n)设置门限，搜索 r大于门限以上的谱峰个数；如果有多个谱峰同时存在，则把幅度最高的谱峰去掉后再依次搜索幅度次高的谱峰。a)2ask信号频谱的谱峰个数提取方法一：ask的频谱图q figure 2,q i回131文件的编辑查看(v)括入(1 tmcd 桌囿营口 w超助(h)不口国 e07000i60005000(400030002000woo0 cask的频域波形* * i i l臼l- -, -* *匕= w. j -)0.511.622.533.544x 105

14、分析可得，该频谱图只有一根离散谱线，在方法一中，可以设置门限r=5000提取程序如下：%提取谱峰%方法一a=zeros(1,length(abs(fft(ask);a=abs(fft(ask);num=0;for i=1:length(abs(fft(ask)if a(i)5000num=num+1;fprintf(方法一：2ask信号频谱的谱峰数目为num=%dn ,num);endend ask方法一：2ask信号频谱的谱博数目为5=1al加噪后识别的谱峰数随信噪比的变化:igure 1方法一：信噪比为-1 odb时2ask信号频谱的诺雌数目为nu际1方法一：信噪比为-9db时2ask信号

15、频谱的语蟾数目为皿:1 方法一：信噪比为-3db时2a5k信号频谱的道幅数目为mm= 1 方法一二信噪比为一：db时2ask信号撅谐的谐嶂数目为皿m 1 方法一二信噪比为-6此时2ask信号频谱的诸博数目为euh二i 方法一二信噪比为-5db时2ask信号频谱的谱峰数目为nun=l 方法一：信嗓比为t讪时2ask信号频谱的谱峰数目为讥5= 1 方法一：信噪比为-3db时2ask信号频谱的造睚数目为皿=1 方法一：信噪比为-edb时2ask信言频谱的造蛭数目为皿1 方法一二信噪比为-1 db时2ask信号缴谱的借嵯数目为皿目1 方法一二信噪比为0北时窈3k信号频谙的谱喧数目为片5=1 方法一二信

16、噪比为1超时2a5k信号频谱的谱幄数目为num=1 方法一：信噪比为2db时2ask信号频谱的诸龌勤目为mm=1 方法一：信噪比为3db时2a5k信覆频谱的谱醺数目为huh 1 方法一：信噪比为4db时2a1信号频谱的谱醺翱目为近皿二1 方法一二信噪比为5db时方sk信号频谱的谱嵯题目为mm二i 方法一：信噪比为6db时裔$长信号频谱的谱峨数目为皿 方法一:信噪比为；df时2a5k信号频谱的谱摩数目为num= 1 方法一：信噪比为8db时2a5k信号频谱的谱嘘却目为mm=l 方法一：信噪比为9db时插5k信号粉谱的谱雌题目为inm=1 fx i程序：% 加噪处理后 识别谱峰数随信噪比的变化曲线

17、 num=zeros(1,20);a=zeros(1,length(abs(fft(ask) j=1;for snr=-10:9askk=awgn(ask,snr, measured );a=abs(fft(askk);n=0;for i=1:length(abs(fft(askk)if a(i)5000n=n+1;end endfprintf( 方法一：信噪比为 ddb时2ask信号频谱的谱峰数目为num=%dn ,snr,n);num(j)=n;j=j+1;snr=snr+1;endsnr=-10:9;plot(snr,num);axis(-10 9 0 10);xlabel( 信噪比 /

18、db);ylabel（谱峰数目/个）；title（谱峰数随信噪比变化曲线 ）；grid on;方法二：老师，第二种算法学生实在没有看懂，自认为那个算法有很大歧义, 因此只做了一种算法提取，忘老师见谅。b）2fsk信号频谱的谱峰个数提取fsk加噪频谱图：q figure 1i 0 .回 汉 ta谱峰提取程序：根据频谱图可以设置门限 r=2000;%谱峰提取程序a=zeros(1,length(abs(fft(fsk);a=abs(fft(fsk);num=0;for i=1:length(abs(fft(fsk)if a(i)2000num=num+1;endfprintf(end方法一：2fs

19、k信号频谱的谱峰数目为num=%dn ,num）;1命令行窗口方法一二2fek信号频谱的谱幢数目为皿2谱峰数随信噪比变化曲线：文件时 (ej葺看(v)插入虚面窗口 (w)帮助(h)国。结w 0婆/图i 口图i 国0987654321 ojh02468信噪比ae方法一二信噪比为-10db时2fsk信号频谱的谱龌数目为115=2 方法一二信噪比为-9db时2fsk信号频谱的谱雌数目为n32 方法一二信噪比为-8de时2f5k信号频谱的谱嘘数目为mm=2 方法一：信嗪比为-;db时2f5k信号频谱的谱崂数目为nuiz 方法一：信噪比为-6db时2fsk信号频谱的谱蟾数目为mm=2 方法一:信i曼比为

20、-5db时2fsk信号频谱的谱嵯数目为nu&2 方法一二信噪比为7db时2fsk信号频谱的谱解数目为 方法一二信噪比为-3db时2fskft号撅谓的谱蜂数目为n皿二2 方法一:信嗪比为-2db时2fsk信3频谱的谱嘘数目为3=2 方法一：信i桑比为-idb时2fsk信号频谱的谱蟾数目为nuw2 方法一：信噪比为ode时2fsk信号频谱的谱陲数目为理皿二2 方法一二信噪比为idb时2fsk信言频谱的谓嵯数目为mm二2 方法一:信噪比为2db时2fsk信号频谱的谓帷数目为nun=2 方法一：信噪比为3db时2fsk信号频谱的谱峥数目为mm=2 方法一：信噪比为4db时芬sk信号频谱的谱瞋数目为nu

21、=2 方法一：信噪比为5de时2fek信号版谱的诸嵯数目为皿=2 方法一二信噪比为6de时2fek信号r频谱的谱嵯数目为mm二2 方法一：信噪比为7db时2fek信号频谱的谱帷数目为nun=2 方法一：信噪比为时2f3k信号频谱的谱嘘数目为mm=2 方法一：信嗓比为9de时2f5k信号频谱的诺峨数目为mi游2程序：clear all close all ; clc;fl f2; fc=(f1+f2)/2%因为调频指数为 0.6 中心频率为400hz ,可以计算出h=(f2-f1)/fb fd=200 f2=460fc=400;%载波频率fs=1200;fb=200;%采样频率%tb比特周期f1

22、=340t=2;%生成an bn随机序列a=randi(0,1,1,t*fs);%产生二进制码元n=length(a);b=zeros(1,n);for i=1:nb(i)=1-a(i);end%根升余弦成型滤波器nt=60;r=0.3;rate=fs/fb;%ii采卞率=fs/ 码速率 delay=nt/(2*rate);t=0:1/fs:t-1/fs;c=a.*exp(j*2*pi*f1*t)+b.*exp(j*2*pi*f2*t);fsk=rcosflt(c,fd,fs,sqrt ,r,delay);%波形成型%谱峰提取程序num=zeros(1,20);a=zeros(1,length

23、(abs(fft(fsk)j=1;for snr=-10:9fskk=awgn(fsk,snr,measured );a=abs(fft(fskk);n=0;for i=1:length(abs(fft(fskk) if a(i)2000n=n+1;endendfprintf( 方法一：信噪比为 ddb时2fsk信号频谱的谱峰数目为num=%dn ,snr,n);num(j)=n;j=j+1;snr=snr+1;endsnr=-10:9;plot(snr,num);axis(-10 9 0 10);xlabel( 信噪比/db );ylabel( 谱峰数目 / 个 );title( 谱峰数随信

24、噪比变化曲线 );grid on ;c)bpsk 信号频谱的谱峰个数提取bpsk 加噪频谱图(信噪比为-5db) ：bpsk二次方谱谱峰提取程序：根据二次方频谱图可以设置门限r=5000帝管仃匐口方法一：efsk信号二次方谓的谱瞋数目为a ”%二次方谱谱峰数目提取a=zeros(1,length(abs(fft(bpsk.a2);a=abs(fft(bpsk.a2);num=0;for i=1:length(abs(fft(bpsk.a2) if a(i)5000 num=num+1;endfprintf(end方法一：bpsk信号二次方谱的谱峰数目为num=%dn ,num);谱峰数随信噪比

25、变化曲线:文件附编款e）|日口金| 4isl3-30 5工具（d莫面（口）窗口 （wd要/一国| 目| 帮勖(h)o0987654 3 21 045000n=n+1;endendfprintf( 方法一：信噪比为 ddb时bpsk信号频谱的谱峰数目为num=%dn ,snr,n);num(j)=n;j=j+1;snr=snr+1;endsnr=-10:9;plot(snr,num);axis(-10 9 0 10);xlabel( 信噪比 /db);ylabel(谱峰数目/个,);title(谱峰数随信噪比变化曲线);grid on;d)qpsk信号频谱的谱峰个数提取qpsk加噪频谱图:谱峰数

26、随信噪比变化曲线:qpsk四次方谱谱峰提取程序：根据频谱图设置门限r=2000“菱口方法一:qpsk信号频谱的谱峥数目为n皿=1i%四次方谱的提取a=zeros(1,length(abs(fft(qpsk);a=abs(fft(qpsk.a4);num=0;for i=1:length(abs(fft(qpsk.a4)if a(i)2000 num=num+1;endendfprintf( 方法一：qpsk信号频谱的谱峰数目为num=%dn ,num);plot(fftshift(abs(fft(qpsk.a4);title( qpsk 四次方谱)；4-20246310(8噪比*b2000150010000-108h figure 1文件的编瑁(e)查看(v)插入(i) zmcd 桌面(5 窗口 (w)帮助(hjn昌u jtq式耳/包卜口百彳回谱峰数随信噪比变化曲线 4000 3600 3000 2500方法一：信噪比为todb时qp5k信号频谱的谱瞳数目为皿3942 方法一：信噪比为-9dfi