通信原理上机指导_实现一个基本的AM 或DSB 调制解调系统.docxVIP

实验内容及步骤3.1 实验内容3.1.1 基本内容：实现一个基本的AM 或DSB 调制解调系统3.1.2 拓展内容：实现一个SSB 调制解调系统3.2 实验步骤Step1. 打开matlab新建一个M 文件Step2. 编写幅度调制代码（1）产生一个信息信号（2）产生一个载波信号（3）调制Step3. 添加噪声Step4. 带通滤波Step5. 相干解调（1）经过乘法器（2）低通滤波主要函数介绍：? 添加噪声：y = awgn(x,SNR) 在信号x中加入高斯白噪声。信噪比SNR以dB为单位。x的强度假定为0dBW。如果x是复数，就加入复噪声。 y = awgn(x,SNR,SIGPOWER) 如果SIGPOWER是数值，则其代表以dBW为单位的信号强度；如果SIGPOWER为measured，则函数将在加入噪声之前测定信号强度。 y = awgn(…,POWERTYPE) 指定SNR和SIGPOWER的单位。POWERTYPE可以是dB或linear。如果POWERTYPE是dB，那么SNR以dB为单位，而SIGPOWER以dBW为单位。如果POWERTYPE是linear，那么SNR作为比值来度量，而SIGPOWER以瓦特为单位。? 滤波器设计：在MATLAB下设计IIR滤波器可使用Butterworth函数设计出巴特沃斯滤波器，使用Cheby1函数设计出契比雪夫I型滤波器，使用Cheby2设计出契比雪夫II型滤波器，使用ellipord函数设计出椭圆滤波器。下面主要介绍第一个方法。主要函数用到的函数有：[B , A] = butter( n, Wn,‘ftype’)其中n代表滤波器阶数，Wn代表滤波器的截止频率，这两个参数可使用buttord函数来确定。buttord函数可在给定滤波器性能的情况下，求出巴特沃斯滤波器的最小阶数n，同时给出对应的截止频率Wn。buttord函数的用法为：[n ,Wn] = buttord( Wp, Ws, Rp, Rs)其中Wp和Ws分别是通带和阻带的拐角频率（截止频率），取值范围为0至1之间（用fs/2归一化）。如设置40-60Hz的带通滤波器就设置下面代码中的a为[40 60], b只要大于这个范围如[30 75]都可以，与Rp, Rs (Rp和Rs分别是通带和阻带区的波纹系数)一起调整，在设计可实现的前提下，Rp的取值应越小越好，Rs的取值则越大越好。当其值为1时代表采样频率的一半。不同类型（高通、低通、带通和带阻）滤波器对应的Wp和Ws值遵循以下规则：? 高通滤波器：Wp和Ws为一元矢量且WpWs；? 低通滤波器：Wp和Ws为一元矢量且WpWs；? 带通滤波器：Wp和Ws为二元矢量且WpWs，如Wp=[0.2,0.7],Ws=[0.1,0.8];? 带阻滤波器：Wp和Ws为二元矢量且WpWs，如Wp=[0.1,0.8],Ws=[0.2,0.7]。下面给出设计一个Butterworth型IIR滤波器的主要程序：Wp=a/(fs/2);若为低通滤波器a,b分别为一个固定的值，且ab；若为高通滤波器a,b分别为一个固定的值，且ab；若为带通滤波器a,b分别为一个区间如[a1 a2]，[b1 b2];且a1b1,a2b2。Ws=b/(fs/2;Rp=c; %% c为通带最大衰减分贝Rs=d; %% d为阻带最小衰减分贝[N,Wn]= Buttord(Wp,Ws,Rp,Rs) ; %%%求滤波器的阶数和截止频率[B,A]=Butter(N,Wn,stype); %%stype为滤波器类型可选：bandpass、low、high等q=filtfilt(B,A,p); 将信号p通过所设计的滤波器。(3) 频谱绘制第一步：Y = fft(X); 将信号X进行傅里叶变换，点数为信号长度。第二步：f=(0:40000)*fs/40001-fs/2;第三步：plot(f,fftshift(abs(Y))); fftshift是将FFT的直流分量移到频谱中心。4. 实验报告内容及要求(1) 选择AM或DSB中的任何一个写出：实验原理框图，实验结果及分析，程序源代码;(2) 要求画出信号频率为10Hz，载波频率为50Hz,采样率为1000Hz,信噪比为5时，对应框图每一点的波形及频谱图.下面是解答：程序源代码%f=10;t=-0.2:0.001:0.2;m=cos(2*pi*f*t);figure;plot(t,m);fs=1000;M=fft(m);f2=(0:400)*fs/401-fs/2;figure;plot(f2,fftshift(abs(M)));%F=50;A=1;t=-0.2:0.001:0.2 s=A.*cos(2*pi*F*t);figure;plot(t,s);S=fft(s);f3