基于matlab的am调制系统仿真

基于matlab的am调制系统仿真内容摘要：

（ 2）将所需类型数字滤波器的边界频率转换成相应的模拟滤波器的边界频率，转换公式为Ω =2/T tan( ) (3)将相应类型的模拟滤波器技术指标转换成模拟低通滤波器技术指标。 （ 4）设计模拟低通滤波器。 （ 5）通过频率变换将模拟低通转换成相应类型的过渡模拟滤波器。 （ 6）采用双线性变换法 将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。 我们知道，脉冲响应不变法的主要缺点是会产生频谱混叠现象，使数字滤波器的频响偏离模拟滤波器的频响特性。 为了克服之一缺点，可以采用双线性变换法。 下面我们介绍用窗函数法设计 FIR 滤波器的步骤。 如下： （ 1）根据对阻带衰减及过渡带的指标要求，选择串窗数类型（矩形窗、三角窗、汉宁窗、哈明窗、凯塞窗等），并估计窗口长度 N。 先按照阻带衰减选择窗函数类型。 原则是在保证阻带衰减满足要求的情况下，尽量选择主瓣的窗函数。 （ 2）构造希望逼近的频率响应函数。 （ 3）计算 h(n).。 （ 4）加窗得到设计结果。 接下来，我们根据语音信号的特点给出有关滤波器的技术指标： 低通滤波器的性能指标： 通带边界频率 fp=300Hz，阻带截止频率 fc=320Hz，阻带最小衰减 As=100db , 通带最大衰减 Ap=1dB 在 Matlab中，可以利用函数 fir1设计 FIR滤波器，利用函数 butter,cheby1和 ellip 设计 IIR 滤波器，利用 Matlab 中的函数 freqz 画出各步步器的频率响应。 hn=fir1(M， wc， window)，可以指定窗函数向量 window。 如果缺省 window参 数，则 fir1默认为哈明窗。 其中可选的窗函数有 Rectangular Barlrtt Hamming Hann Blackman 窗，其相应的都有实现函数。 MATLAB 信号处理工具箱函数 buttp buttor butter 是巴特沃斯滤波器设计函数，其有 5 种调用格式，本课程设计中用到的是基于 matlab的 AM— DSB 调制系统仿真 8 [N,wc]=butter(N,wc,Rp,As,’ s’ ),该格式用于计算巴特沃斯模拟滤波器的阶数 N 和 3dB 截止频率 wc。 MATLAB 信号处理工具箱函数 cheblap,cheblord和 cheeby1 是切比雪夫 I型滤波器设计函数。 我们用到的是 cheeby1 函数，其调用格式如下： [B,A]=cheby1(N,Rp,wpo,’ ftypr’ ) [B,A]=cheby1(N,Rp,wpo,’ ftypr’ ,’ s’ ) 函数 butter,cheby1 和 ellip 设计 IIR 滤波器时都是默认的双线性变换法，所以在设计滤波器时只需要代入相应的实现函数即可。 下面我们将给出 FIR 数字滤波器的主要程序。 %=========================FIR低通滤波器 ======================= Ft=2020。 %采样频率 fpts=[100 120]。 %通带边界频率 fp=100Hz，阻带截止频率 fs=120Hz mag=[1 0]。 dev=[ ]。 %通带波动 1%，阻带波动 5% [n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev,Ft)。 %kaiserord 估计采用凯塞窗设计的 FIR滤 波器的参数 b21=fir1(n21,wn21,Kaiser(n21+1,beta))。 %由 fir1 设计滤波器 [h,w]=freqz(b21,1)。 %得到频率响应 plot(w/pi,abs(h))。 grid on title(39。 FIR 低通滤波器 39。 )。 基于 matlab的 AM— DSB 调制系统仿真 9 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 100 . 20 . 40 . 60 . 811 . 21 . 4F I R 低通滤波器 如图 FIR 低通滤波器 三 AM 解调 %=========================AM信号解调 ======================= t=1::1。 A0=10。 %载波信号振幅 A1=5。 %调制信号振幅 A2=3。 %已调信号振幅 f=3000。 %载波信号频率 w0=2*f*pi。 m=。 %调制度 k=。 %DSB 前面的系数 mes=A1*cos(*w0*t)。 %调制信号 Uam=A2*(1+m*mes).*cos((w0).*t)。 %AM 已调 信号 Dam=Uam.*cos(w0*t)。 %对 AM 调制信号进行解调 subplot(2,1,1)。 plot(t,Dam)。 grid on。 title(39。 滤波前 AM 解调信号波形 39。 )。 subplot(2,1,2)。 Y5=fft(Dam)。 % 对 AM 解调信号进行傅里叶变换 plot(abs(Y5)),grid。 基于 matlab的 AM— DSB 调制系统仿真 10 title(39。 滤波前 AM 解调信号频谱 39。 )。 axis([187960,188040,0,202000])。 1 0 . 8 0 . 6 0 . 4 0 . 2 0 0 . 2 0 . 4 0 . 6 0 . 8 10246滤波前 AM 解调信号波形1 . 8 7 9 6 1 . 8 7 9 7 1 . 8 7 9 8 1 . 8 7 9 9 1 . 8 8 1 . 8 8 0 1 1 . 8 8 0 2 1 . 8 8 0 3 1 . 8 8 0 4x 1 0500 . 511 . 52x 1 05滤波前 AM 解调信号频谱 如图 AM 解调信号 %=======================AM解调信号 FIR滤波 ================= t=1::1。 A0=10。 %载波信号振幅 A1=5。 %调制信号振幅 A2=3。 %已调信号振幅 ************************************* f=6000。 %载波信号频率 w0=f*pi。 m=。 %调制度 ************************************** Uc=A0.*cos(w0*t)。