2023年信号与通信仿真实验报告.pdf

《2023年信号与通信仿真实验报告.pdf》由会员分享，可在线阅读，更多相关《2023年信号与通信仿真实验报告.pdf（27页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、现代通信原理与技术信号与通信仿真设计实习报告姓名：班级:学号:一 实 验 目的在本实验中使用的软件工具是M A T L A B。设计本实验的目的是希望在以下几方面有所收获：1 .会M A T L A B软件的最基本运用。M A T L A B是一种很实用的数学软件,它易学易用。M A T L A B对于许多的通信仿真类问题来说是比较合适的。2 .了解计算机仿真的基本原理及方法，知道如何通过仿真的方法去研究通信问题。3 .加深对信号与系统和通信原理及其相关课程内容的理解。二实验特点与硬件实验相比,软件实验具如下一些特点：1 .软件实验具有广泛的实用性和极好的灵活性。在硬件实验中改变系统参数也许意

2、味着要重做硬件,而在软件实验中这只是该一两个数据，或者只是在屏幕上按几下鼠标。2.软件实验更有助于我们较为全面地研究通信系统。有许多问题，通过硬件实验来研究也许非常困难，但在软件实验中却易于解决。3 .硬件实验的精确度取决于元器件及工艺水平，软件实现的精确度取决于CPU的运算速度或者说是程序的运算量。4 .软件实验开发周期短，成本低。三 上机实验规定1 .掌握m a t l a b的基本操作及了解基本的仿真方法，分析运营范例程序。2 .按以下规定编制仿真程序并调试运营(1)基本信号的仿真(2 )模拟调制与解调的仿真(3)数字基带传输码型的仿真(4)数字调制与解调的仿真(5)脉冲编码调制仿真四

3、实 验 内 容1、基本信号的仿真(1)产生并绘出以下信号：a 单位脉冲序列b 单位阶跃序列c 正弦信号及其频谱d 周期锯齿波sawto o th()e 周期方波s q u a re()f 实指数序列y(n)=2ng s i n2 n fit*c o s 2 n f2t f l=50Hz f2=2023H z(2)产生一条一 2 到 2 之间的Sa(200 t)曲线。(3)产生下面信号，并绘出频谱 t 0 tt0/4s(t)=-t+tO/4 tO/4t 3 t0/4 假设 t 0=0.5 st-t 0 3t0/4t 0 x(i )=1；o,20.1e l s e x (i)=0 ；0e n d单

4、位阶跃序列e n d10.9p l o t (i,x);s t e m(x);o s0.7t i t l e(，单位阶跃序列，);n c 6a x i s (-1 0 3 0 0 1 );050 40.3C正弦信号及其频谱 0 20.1%正弦信号及其频谱,s i n e.m o-ii-5 0 5 10 15 20 25 30n=L 1:0 .0 1:1 ；d周期锯齿波s a w t o o t h()%周期锯齿波,j i z h i b o.mn=-3.1 4*4：0.0 1 ：3.1 4 *4 ;x=s a w t o o t h(n)；P l o t(n,x)f l=5 0 ;f 2=2

5、0 2 3 ;y=s i n(2 *p i *c o s(2*pi*f2*n;p l o t (n,y)t i t l e (y =s i n 2 n f 1 n c o s 2 兀 f 2 n (f 1=5 0 H zf 2 =2 0 2 3 H z)t i t l e（周期锯齿波）心得体会:e周期方波s q u a r e。%周期方波，f a n g b o.my=s q u a r e (n );p l o t(n,y)t i t 1 e(周期方波)a x i s(0 1 0 -2 2 )f实指数序歹U y (n )=2 n通过实验熟悉了 m a t l a b的一些基本指令的使用，知道

6、了一些产生基本信号的编程方法,为后续实验打下了基础。%实指数序列 y (n)=2.n,s h i z hi s h u.mend曲线。%产 生 一 条-2到2之 间 的Sa(200t)曲end线.s in c.mY,y,d f l=f f t_seq(y,t s,d f)；t=-2:0.0 1:2;f=O:df l:df 1*(1 eng t h(y)-1)-F s/2;x=s i n c(2 00*t/pi);%频 率 矢 量s u bplot(2,l,1 );p 1 ot(t,y(1:1 e nBs=10;%高斯噪声限带gth(t)no=0.01;%单边带功率谱密度title(St 信号,

7、)；tO=1;%信号连续时间s ubplot(2,1,2);p 1 o t(f,abs(f f t sts=O.0 0 1 1;%时间采样间隔hift(Y)；Fc=250;%载波中心频率t i t i e St信号频谱);Fs=1/ts；%采样频率)；)；%调制信号%载波信号%已调信号-500-400-300-20D-100 0 100 200 300 400 500a xis(0 I-1 1 );t itle(调制信号);%DSB 调制,DSB.ms u bplot(3,2,3)；p 1 o t(t,c(l:length(t)%调 制 信 号 时 域 表 达 式 u(t)=axis(l 0

8、0.1-1 1);tit1e(载波信号)；m(t)c(t)=A c*m(t)c o s (2 fct)subpl o t(3,2,5);p 1 o t(t,u(1 ：length(t);%函数M,m,d f=fft_ s eq(m,tsax i s(0 1 -1 1)；t i t l e 已调信号);df)求取频谱函数sub p lot(3,2,2);p 1 o t(f,abs(f f ts h if t(M)；%高 斯 噪 声 可 表 达 成 x=sqr t(B stitle。调制信号频谱,);*no)*r andn(l,M)su b p 1 ot(3,2,4);p 1 ot(f,abs(f

9、 ft s hif t(C);title d 载波信号频谱）；title，滤波器输入信号频谱）;s u bp 1 ot(3,2,6);plot(f,a b s(ff t s hisubplot(2,2,3);p lot(t,dem(1 :1 e n gth(t);ft(U)；t i t l e-滤波器输出信号;t i t l e d 己调信号频谱）;sub p lot(2,2,4);pl o t(f,a bs(f f t sh i f t(D E M);%相干解调，加入高斯白噪声noisenoi s e=s q r t(Bs*n o)*ran d n(1,leng t h(u);Y,y,df

10、1=fft_ s e q(y,ts,d f);Y=Y/Fs;y=(u+n o i s e).*c;%低通滤波f_ c u t of f=15 0;%滤波器截止频率n _c u toff=floor(150/d f l);f=O:dfl:dfl*(1 e n gth(m)-1)-Fs/2；H=zer o s(size(f);H(l:n_ c utoff)=2*one s(l,n _ cut o ff);H(le n gt h(f)-n_ c u tof f+1:length(f)=2*one s(l,n_cut o f f)；D E M=H.*Y;%滤 波器输出频谱d em=real(ifft(

11、D E M)*Fs；fi g u r e(2);s ubplot(2,2,l);plot(t,y(1 ：length(t);t i t le(,滤波器输入信号)；subpl o t(2,2,2);plot(f,a b s(fftsht i t l e。滤波器输出信号频谱）;调制信号载波信号频谱0.500.5-100-2000200就被器除入信号频谱i ft(Y);s u bp 1 ot(3,2,3);p lot(t,c(1:1 ength(t)%AM调制信号的时域表达式为:(t);u(t)=AC1+a mn(t)c os(2 n f c t)axis(0 0.1-1 1)title(,载波信号

12、)；%a 是调制指数,mn(t)是通过归一化s u b p lo t (3,2,5);p lot(t,u(l:1 ength(t);解决的消息信号，函数M,m,d f a x i s(0 1 -2 2)ti t 1 e。已调信号,);=f f t_seq(m,ts,d f)求取频谱函数,高斯噪声可表达成B s=1 0;%高斯噪声限带no=0.01;%单边带功率谱密度a=0.8；%调制指数t0=1;%信号连续时间ts=0.0 0 1;%时间采样间隔Fc=2 5 0；%载波中心频率F s=l/ts;%采样频率df=0.3;%频率分辨率t=0:ts:t0;%时 间矢量subplot(3,2 ,2);

13、pl o t(f,ab s(f ftshift(M);axis(-10 1001)t i 1调制信号频谱s ubpl o t(3,2,4)；p 1 ot(f,abs(ff t sh i f t(C);axis(-300 300 0 200)title。载波信号频谱；su b p 1 o t(3,2,6);p 1 o t(f,a bs(fftshift(U)；ax i s(-3 0 0 30 0 0 I 1);titl e(，己调信号频谱);%相干解调，加入高斯白噪声noisen o i se=s qr t(Bs*no)*ra n d n(1 ,1 ength(u);m=s in(2*p i.*

14、t)；%调制信号y=(u+n oi s e).*c;c=cos(2*pi*F c.*t)；%载波信Y,y,df 1=fft_s e q(y,ts,d f)；Y=Y/Fs;%低 通滤波f_cu t off=1 5 0;n_ c u t off=flo or(l 5 0/dfl);H=zeros(s i z e(f);H(1：n _ cu t o ff)=2*o n e s(1,n _cuto f f)；H(length(f)-n_ c u tof f+1 :1 e n g th(0)=2*0 n e s(1,n_ c u tof f)；DEM=H.*Y；%滤波器输出频谱d e m=rea 1 (

15、ifft(DE M)*Fs;%滤波器的输出f ig u re(2);sub p lot(2,2,l);pl o t(t,y(klength(t);a x is(0 1 -4 4);title。滤波器输入信号)；s ubp 1 o t(2,2,2);p l o t (f,ab s(ffts hift(Y);ax i s(0 10 0 0.5);title(z 滤波器输入信号频谱)；sub p lot(2,2,3 );p lot(t,d e m(1 :1 e ng t h(t);a x i s(0 1 -1 4);t itle C滤波器输出信号)；subplo t(2,2,4);p I ot(f,

16、abs(fftsaxis(-1 0 10 0 1 )；t i t le(，滤波器输出信号频谱)；载波信号 载波信号频谱心得体会：通 过AM与DSB实验,我对模拟调制解调的方法及原理有了清楚直观的结识，知道了各类调制解调方法的优缺陷和互相之间的区别。3、数字基带信号的码型的仿真(1)单极性非归零码function y=djxbg 0(x)%本函数实现将输入的一段二进制代码编为相应的单极性不归零码输出，编码规则:零电平与正电%平分别相应着二进制代码0 和 1%输入X为二进制码，输出y 为编%好的码，给出计算每一个码元的%点数，由于我们只有用离散的点%来得出连续的函数表达。在命令hift(DEM);

17、%窗口键入如下指令即会出现图形%t=l 0 0 1 1 0 0 0 0 1 0 axis(0,i,-0.1,1.1)；1；%采用了 title命令来实现标记出各码元相应的二元信息%djxb g0(t);title6 00 1 1 0 000 1 O r );ts=0.0 0 1;%时间采样间隔gri d=l/ts;t=O:ts:le n g th(x);fo r i=l：leng t h(x),%计算码元的值if(x(i)=l),%假如信息为1f or j=1：gr i d,%该码元相应的点值取1y(i-l)*gr i d+j)=1；endelsefor j=l:gri d,%反之，信息为0,

18、码元相应点取0y(i-1)*g r id+j)=0；e n de ndendy=y,x(i);x 1 abe IC单极性不归零码）滤波器输入信号檄器输入信号0 0.5 1%为了画图注意序列加上最后一位su b plot(2,l,1);plot(t,y);微器输出信号 蹦器输出信号频谱10.5010011000 0码 不 的 值10 2 4 6 8 10 12单极性不归零码 e n di)=1),%假如信息为1J =1：g r i d/2,%该码元相应的点值取y(i-1)*g r i d+j)=l;(2)双极性归零码，占空比5 0%f u n c t i o n y=s j x g O(x)%本

19、函数实现将输入的一段二进制代码编为相应的双极性归零码输出%编码规则:负电平与正电平分别相应着二进制代码0和1,每个码元内的脉冲都回到零电平%输 入x为二进制码，输 出y为编好的码%给出计算每一个码元的点数，由于我们只有用离散的点来得出连续的函数表达。%在命令窗口键入如下指令即会出现图形%t =L 1 oo 1 1 ooooio u；%s j x gO(t);g r i d=3O O；t=0：1/g r i d:1 e n g t h(x);f o r i=l:1 e n g t h (x),%计算f o r j=g r i d/2:g r i d,%占空比 5 0%,归 0y (i-l)*g

20、r i d+j)=O；ende I s ef o r j=1:g r i d/2,%反之，信 息 为0,码元相应点取-1y (i -l)*g r i d+j)=-1;endf o r j=g r i d/2:g r i d,%占空比 5 0%,归 0y (i -l)*g r i d+j)=O;ende n dendy=y,x ;%为了画图，注意将序列加上最后一位M=m a x (y );m=m i n(y);s u b p 1 o t(2,l,1)；p l o t (t,y);x l a b e l (双极性归零码)a x i s(0,i,m-0.l,M+0.1)%采用了 t i t l e命

21、令来实现标记出各码元相应的二元信息t i t 1 e (1 0 0 1 1 0 0 00 1 0 1 );1 0 0 1 1 0 0 0碣元曲值0 1(3)双相码f unction y,t=shuan g xiang(x)%本函数实现将输入的一段二进制代码编为相应的双向%码 输出，编码规则：“0”码 用“01”表达，“1”码 用“1 0”%表达，输入x 为二进制码，输出y 为编好的码，给出计算，每一个码元的点数，由于我们只有用离散的点来得出连续的函数表达。%在命令窗口键入如下指令即会出现图形%t=ll 00 I 1 0 0 00 I 0 1;%shu a n g xi a n g(t);t s

22、=0.0 0 1;%时间采样间隔g rid=l/ts；t=O:t s：length(x);e nde ls e%假如信息为0,该码元相应的点值取01for j=1:grid/2,y(i-1)*gr i d+j)=0;e n df o r j=grid/2:gr i d,y(i-l)*grid+j)=l;ende n den dy=y,x;%为了画图，注意将序列加上最后一位subp 1 o t(2,1,1 );pl o t(t,y);a x is(0,i,-0.1,1.1)；%采用了 t i t l e 命令来实现标记出各码元相应的二元信for i=l：1 eng t h(x),%计算息X la

23、bel(双向码);%am i(t)；title(4 0 0 1 1 0 0 0 0 1 0r);娴(4)AMI 码f unc t io n y=ami(x)%本函数实现将输入的一段二进制代码编为相应的AMI码输出%编码规则：二 进 制“1 ”交替变换为“+1”，“/，“0”保持不变%输 入 x 为二进制码，输 出 y 为编好的码%给出计算每一个码元的点数，由于我们只有用离散的点来得出连续的函数表达%在命令窗口键入如下指令即会出现图形%t=1 0 0 1 1 0 0 0 0 1 0gri d=30 0;s=0;%s为符号标志t=0:1 /gri d：length(x);for i=1 ：leng

24、th(x),%计算码元的值if(x=1),%假如信息为1i f(s=0),f o r j=1:grid,%s=0 时，该码元相应的点值取1y(i-l)*grid+j)=1；s=l;endelsef o r j=1:gr i d,%s=l 时,该码元相应的点值取-1y(i-l)*gri d+j)=-l;s=0;ende nde 1 sefor j=l:g rid,%反之,信息为0,码元相应点取0y(i-l)*gr i d+j)=O；en den d1;endy=y,x；%为了画图，注意将序列Y2(i),y2(i)=Hh(f(i),0.5,1);%加上最后一位a=0.5s u bplo t(2,1

25、,1);Y 3(i),y3(i)=H h(f(i),P lo t(t,y)；1,1);%a=1axis(0 i 1.1 1.11)En dxlabel(AM I 码);for i=(N+1)：(2*N 1)e n d%采用了 t i t l e 命令来实现标记出f(i)=(i-N)/100；各码元相应的二元信息Y1(i),yl(i)=Hh(f(i),0,1);%a=0ti t 1 e(z 1 0 0 1 1 0 0 0 0 1 0IY2(i),y2(i)=Hh(f(i),0.5,1)；%a=0.51)；Y3(i),y3(i)=Hh(f(i),1,1)；%a=l(3)a=0,0.5,1 的升余弦

26、滚降系统频谱及其各自相应的时域波形%g u njiang.m%取 a=0,0.5,1的升余弦滚降系统频谱及各自相应的时域波形%函数Y,yJ=H(f,a,IT)为升余弦滚降系统频谱及各自相应的时域表达式N=400;for i=l:Nf(i)=(i-N)/100；%采样间隔Y1(i),yl(i)=Hh(f(i),0,1);%a=0figure(l);p l o t (f,Yl);hold onplot(f,Y2/r);hold onplo t(f,Y3；k);le g e n d(a=0r a=0.5 Va=1 1)t i t l 贝升余弦滚降频谱，);ax i s(2 2 0 2);figu r

27、 e(2);p lot(f,yl);ho 1 d o nA M贿H=(T s/2)*(1 +sin(piabs(f)*2*pi*Ts)/(2*a);升余弦滚降时域波形p 1 o t(f,y2/r);hold on%成 间隔为N-1个 0 的序列P 1 o t(f,y3,k)；Ts=l；%连续时间title C升余弦滚降时域波形);N_sample=1 7;legend(z a=0；a=0.5；a=r)e ye_n u m=7;a xis(-4 4 1 1);alph a=1;N_data=l 0 0 0;%Hh.mdt=Ts/N_sample；%函数Y,y=H(U 1 T)为升余弦滚降t=-3

28、*Ts:dt：3*T s;%时间矢量%产生双极性数字信号d=si g n(ra ndn(L N_d a ta)；d d=si g expan d(d,N _ s a mpie);%系统输入信号%基带系统冲击响应(升余弦)h t=sin c(t/Ts).*(cos(alpha*pi*t/Ts)./(1 -4*al p haA2*t/2/TsA2)；st=con v(d d,ht);%系统输出信号tt=-3*T s：dt:(N_dat a+3)*N_ samp 1 e*dtdt;%时间矢量 ttfigu re(l)subplot(41 1 )；plo t(d d);axis(0 10 0-2 21

29、);t i tie(双极性数字信号)；sub p lot(412)p 1 o t(t,ht)；x 1 a bel(*t/T s )；tit1 e(，系统函数波形(a=l);su b plot(4 13)plot(t t,st);ax i s(0 2 0-1.21.2);xl a be 1 (tt/Ts1);t itle(,基带信s ub p lot(414);xl a b el(ttt7T s);ti t 1 e(基带信号眼图 )；%画眼图ss=z e ros(1，eye_num*N_ s ample);t 11=0:dt:ey e _num*N_ s a m p le*d t dt；for

30、k=3：50s s=st(k*N_samp 1 e+l:(k+eye_ num)*N_ samp ie);d raw n ow;plot(t tt,ss);hold on;endf u nctio n o u t =sig e x pand(d,M)%将输入的序列扩展成间隔为N 1个0的序列N=1 e n gth(d);ou t=zeros(M,N);o ut(1,：)=d;out=res h a pe(o u t,M*N);心得体会通过数字基带信号的码型的仿真实验，我熟知了各种码型的编程规则及各自的特点，并且对实现无码间串扰传输与通信系统的关系有了清楚的结识。4、数字调制与解调的仿真(1)2

31、FSK调制和解调%2FSK 调制，FSKnoise.m号 );df=0.3;%频率分辨率ts=O.0005；%时间采样间隔t i tle(l 0 0 1 1 0 1 00 1 0 1);Fs=1/ts;%采样频率x=1 0 0 1 1 0 1 0 0 1 0U；n=1 en g t h(x);sL t=djx b g0(x);%基带信号为单极性不归零码s2=1 si;%单 极性不归零码的反码cl=sin(2*pi*2 0.*t)；%载 波信号c l,f=20hzsubplot。,2,2);p lo t(t,s 2);axi s(0,n,-0.1y label。基带信号反码s2);ti t le

32、fl 0 0 1 1 0 1 0 0 1 0 I1)；s u b plot(3,2,3);pl o t(t,c 1);a xis(0,1,-1 eC载波信号cl);sub p lot(3,2,4);pl o t(t,c2)；a x is(0,l,-1,it l e d 载波信号c 2、);s ubplot(3,2,5)；p l o t (t,e i 1);ax is(0,n,-l,l)；tit 1 e(，已调信号 eil);subplot(3,2,6);p lo t(t,ei 2);c2=sin(2*p i*1 20.*t);%载波信号c2,f=120hzeil=cl.*s 1；%已调信号e

33、i 1e i 2=c2.*s2；%已 调 信 号 e i 2e=ei 1+ei2;%2FSK 信号noise=0.1 *rand n(l,le n g th(t)；%加入噪声e=e+noise;f igure(1);subp 1 ot(3,2,l);双极性数字信号1 A _A _A _i i V0 1 0 20 30 40 50 60 70 80 90 1 0 0系统函嬲形（利-2-10123t/T s基制资基带信布图。k x x 3 X X XX X2厂 一 一 一 厂 一 厂 一 厂 一 厂 一*01234567tttffsplot(t,s 1);a xis(0,n,-0.1,1.1)；y

34、labe 1（基带信号s）；f i gur e(2)；11001011010L 1wU pIlla 0.5iW o01001101001015 1 0载 姊 号cl0 0.5-1-10subplol(32t i tsubpl(1)；N o t(t,e(1:1 e ngth(t);2F5I【言号e)；22);)o t(f,a b s(fft s h i f t (E)axis(f-2 0 0 Q)O 0 1 0 0 0 0);鼬 信 号c 2频 谱E);p 1 o t(t,eol(1:leng t h(t);）；己调信号的0 5 1 0,信 号eo l(f =2 0 h z),);51112,4

35、)；plot(f,abs(ffts h i ft(E 1);1,二o1o2sN出胆修构o：(%相 干 解 调bl=f i rl(l 0 1,10/8 0 02 0/800);b2=f i rl(101,9 0/8 00110/800);e o l=filter(b 1 ,l,e);eo 2=fi 1 ter(b2,1,e);%通过带通滤波器后的信号E,e,df 1 =fft_ seq(e,t s,d f);E 1,eol,d f 1 =fft_ s eq(eol,ts,df)；E2,e o 2,d f 1 J=fft _ s eq(eo2,ts,df);subplot(3,2,6);pl o

36、t(f,ab s(f ftshift(E2);axis(-2 00 20 0 0 1 00 0 0);t itle C带通滤波器输出信号频谱E2);swl=e ol.*e ol；sw2=e o 2.*e o 2；bn=fir 1 (101,12/8 0010/8 00J);stl=filter(b n,1,sw 1);s t 2=filter(bn,l,sw 2);fig u re(3)s u bp 1 ot(3 1 1);plo t(t,st 1 (1:1 e ng t h(t);ylabeK7 幅度上x 1 a bel(r)tit le(通过低通滤波器后的波形str)s ubplot(3

37、12);plo t(t,st 2(klength(t);f=0:d fl:df 1*(length(e)-l)-F s/2;title（7通过低通滤波器后的波形st2）;ylabelf幅度);xla b e 1 (t )2 F S K信 号e10X 05X01 0 1 50-2 0 0带 通 滤 波 器 输 出 信 号eo2(H 2O hz)带 是1 0 0 0 0 L0-1IIIL IJ0 55 0 0 0o 1-2 0 02 F S K信 号e2F S唯号频耻带 通 海 被 器 输 出 信 号e o 1(f=2 0 h z)20魏 解 牌 出 信 号 解E 11 0 0 0 05 0 0

38、01 0 0 2 0 000 5 1 0 1 5-2 0 0 -1 0 0带 通 源 液 晶 输 出 信 号eo2 20hz)带通凝牌出信号解E 2%判决f o r i=l:l e ngt h(t)i f(s t2(i)=stl(i)s t(i)=O;else st(i)=l;endendsubp 1 o t(3 1 3);plo t(t,s t(1:1 en g t h(t);titled通过抽样判决器后的波形st);ylabe 1 (幅度);xlab e K7 f)(1)2P S K 调制和解调%2PSK 调制，PS K n o i se.m ylabe 1 (基带信号 s)t s=0.0

39、 01；s ubp 1 ot(3,2,2);p l o t (t,s2)；a x i s(0 1,1.1)；F s=1/t s;出le(l 0 0 1 1 0 1 0 0 1 0 D ;df=0.1;x=1 0 0 1 1 0 1 0 0 1 01；i=1 ength(x);si,t=d j xbgO(x);s 2=1 -s 1;cl=s i n(2*pi*l 0,*t);%载波信号c l,初相位位0c2=sin(2*pi*1 0 p i)；%载波信号c 2,初相位位180el=cl.*sl;%已调信号 e 1e2=c 2.*s2;%已 调信号经过低通毓器削脂的破物1e2e i=e l+e 2

40、；%2PSK 信号noise=0.1 *rand n(l,len g th(t);%噪声ei=e i+n o is e;figure(1);subpl o t(3,2,l)：plot(t,sl)；axis(0,i,0.1,1.11);title(10 0 1 1 0 1 00 1 0 1 9 ;p 1 ot(t,e2);a x is(0,i,-l,1);经过低遨融器州佰的破形stlloS0246 3国 油 1 2loS%频率矢量ylabel(基带信号反码s 2);subplot(3,2,3);plot(t,c 1);axis(0,1,1,1);t i t Ie。载波信号c 1(初相位位0),)

41、；subplot(3,2,4)；p lot(t,c 2);axis(0,1,-1,1)；t i tl e(7载 波 信 号 c2(初相位位1 8 0)；s u bplot(3,2,5);p lot(t,e l)；axi s(0,i,1,1)；tille(已调信号e)；subpl o t(3,2,6);H=zer o s(s i ze(f)；H(1 ：n_cutoff)=on e s(l,n_ c utolf);H(length(f)-n_ c u toff+1:length(f)=ones(1,n _ cut o ff)；D EM=H.*E;%滤波器输出频谱dem=r eal(i f ft(D

42、 EM);%滤波器的输出f i gure(2)；su b p 1 o t(3,2,1);pl o t(t,eo(1：1 eng t h(t);axis(0,i,-2,2);tit 1 e(滤波器输入信号eo1);s u bplo t(3,2,2);plot(f,abs(f ftsh i ft(E);dem(n)=1 ;els edem(n)=0；ende ndsubpl o t(3,2,5);p lot(t,dem(l:leng t h(t);a x is(0,i,-0.1,1.11);t itle(,恢复信号 d e m)；D EM=r eal(i f f t(dem);subp 1 ot(

43、3,2,6);p 1 o t(f,a bs(fft s h ift(D E M);ax i s(l0,2,-0.1,0.5)；title,恢复信号频谱DEM);f i g ur e(3);subp 1 ot(2,1,1);p 1 ot(t,e i(1：len g t h(t)；axi s(0,5,-1,1);title-2PSK 信号);sub p 1 ot(2,1 ,2);p l o t (f,ab s(fftsh i ft(E I);心得体会通 过 2FSK与 2PSK实验，我对数字调制解调的实现方法及原理有了清楚直观的结识，知道了各类调制解调方法的优缺陷和互相之间的区别。axis(L-1 0,1 0 0.1,20 2 3);t i tl e(2 P SK 信号频谱，)；总结通过信号与通信仿真实验，我对通信技术中的一些知识和原理有了更深的结识，对通信系统的总体把握更加全面和准确，为此后的进一步学习扎实了基础。在此感谢狄卫国老师的悉心指导。