第1页共53页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第1页共53页通信仿真实验学号：10071213姓名：胡帆目录1实验二滤波器的设计和仿真实现··········22实验四随机过程的计算机仿真··········123实验七数字基带调制··········164信号的功率谱与带限信道的滤波器等效··········205实验：基带传输最佳接收机仿真··········2664psk传输的误码率仿真··········3474PAM系统的MonteCarlo仿真··········378信号的功率谱与带限信道的滤波器等效·········41实验二滤波器的设计和仿真实现第2页共53页第1页共53页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第2页共53页————实验日期：2013.10.27实验目的：各种滤波器的设计与仿真实现。实验内容：1.题目一：试设计一个模拟低通滤波器，fp=3500Hz，fs=4500Hz，αp=3dB，αs=25dB。分别用巴特沃斯和椭圆滤波器原型，求出其3dB截止频率和滤波器阶数，传递函数，并作出幅频、相频特性曲线。实验程序1巴特沃斯型：clear;f_p=3500;f_s=4500;R_p=3;R_s=25;[n,fn]=buttord(f_p,f_s,R_p,R_s,'s');Wn=2*pi*fn;[b,a]=butter(n,Wn,'s');f=0:100:10000;s=j*2*pi*f;H_s=polyval(b,s)./polyval(a,s);figure(1);subplot(2,1,1);plot(f,20*log10(abs(H_s)));axis([010000-401]);xlabel('频率Hz');ylabel('幅度dB');subplot(2,1,2);plot(f,angle(H_s));xlabel('频率Hz');ylabel('相角rad');figure(2);freqs(b,a);实验结果1：第3页共53页第2页共53页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第3页共53页实验程序2椭圆滤波器型：clear;f_p=3500;f_s=4500;R_p=3;R_s=25;[n,fn]=ellipord(f_p,f_s,R_p,R_s,'s');Wn=2*pi*fn;[b,a]=ellip(n,R_p,R_s,Wn,'s');f=0:100:10000;s=j*2*pi*f;H_s=polyval(b,s)./polyval(a,s);figure(1);subplot(2,1,1);plot(f,20*log10(abs(H_s)));axis([010000-401]);xlabel('频率Hz');ylabel('幅度dB');subplot(2,1,2);plot(f,angle(H_s));xlabel('频率Hz');ylabel('相角rad');figure(2);freqs(b,a);实验结果2：第4页共53页第3页共53页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第4页共53页2.题目二：试设计一个巴特沃斯型数字低通滤波器，设采样率为8000Hz，fp=2100Hz，fs=2500Hz，αp=3dB，αs=25dB。并作出幅频、相频特性曲线。实验程序：f_N=8000;f_p=2100;f_s=2500;R_p=3;R_s=25;Ws=f_s/(f_N/2);Wp=f_p/(f_N/2);第5页共53页第4页共53页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第5页共53页[n,Wn]=buttord(Wp,Ws,R_p,R_s);[b,a]=butter(n,Wn);figure(1);freqz(b,a,1000,8000)subplot(2,1,1);axis([04000-303])figure(2);f=0:40:4000;z=exp(j*2*pi*f./(f_N));H_z=polyval(b,z)./polyval(a,z);subplot(2,1,1);plot(f,20*log10(abs(H_z)));axis([04000-401]);xlabel('频率Hz');ylabel('幅度dB');subplot(2,1,2);plot(f,angle(H_z));xlabel('频率Hz');ylabel('相角rad');实验结果:第6页共53页第5页共53页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第6页共53页3.题目三：试设计一个切比雪夫1型高通数字滤波器，采样率为8000Hz，fp=1000Hz，fs=700Hz，αp=3dB，αs=20dB。并作出幅频、相频特性曲线。实验程序：f_N=8000;f_p=1000;f_s=700;R_p=3;R_s=20;Ws=f_s/(f_N/2);Wp=f_p/(f_N/2);[n,Wn]=cheb1ord(Wp,Ws,R_p,R_s);[b,a]=cheby1(n,R_p,Wn,'high');freqz(b,a,1000,8000)subplot(2,1,1);axis([04000-303])实验结果：4.题目四：试设计一个椭圆型带通数字滤波器。设采样率为10000Hz，fp=[1000,1500]Hz，fs=[600,1900]Hz，αp=3dB，αs=20dB。并作出幅频、相频特性曲线。实验程序：f_N=10000;第7页共53页第6页共53页编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第7页共53页f_p=[1000,1500];f_s=[600,1900];R_p=3;R_s=20;Ws=f_s/(f_N/2);Wp=f_p/(f_N/2);[n,Wn]=ellipord(Wp,Ws,R_p,R_s);[b,a]=ellip(n,R_p,R_...
