ssb调制与解调

《ssb调制与解调》由会员分享，可在线阅读，更多相关《ssb调制与解调（23页珍藏版）》请在装配图网上搜索。

1、中用M文件实现SSB解调一、课程设计目的本次课程设计是对通信原理课 程理论教学和实验教学的综合和总结。通过这次课程设计，使同学认识和理解通信系统， 掌握信号是怎样经过发端处理、被送入信道、然后在接收端还原。要求学生掌握通信原理 的基本知识，运用所学的通信仿真的方法实现某种传输系统。能够根据设计任务的具体要 求，掌握软件设计、调试的具体方法、步骤和技巧。对一个实际课题的软件设计有基本了 解，拓展知识面，激发在此领域中继续学习和研究的兴趣，为学习后续课程做准备。二、 课程设计内容（1）熟悉中M文件的使用方法，掌握SSB信号的解调原理， 以此为基础用M文件编程实现SSB信号的解调。（2）绘制出SSB

2、信号解调前后在时 域和频域中的波形，观察两者在解调前后的变化，通过对分析结果来加强对SSB信号解 调原理的理解。（3）对信号分别叠加大小不同的噪声后再进行解调，绘制出解调前后信 号的时域和频域波形，比较未叠加噪声时和分别叠加大小噪声时解调信号的波形有何区别， 借由所得结果来分析噪声对信号解调造成的影响。（4）在老师的指导下，独立完成课程 设计的全部内容，并按要求编写课程设计论文,文中能正确阐述和分析设计和实验结果。三、 设计原理1、SSB解调原理在单边带信号的解调中，只需要对上、下边带的其中一个边 带信号进行解调，就能够恢复原始信号。这是因为双边带调制中上、下两个边带是完全对 称的，它们所携带

3、的信息相同，完全可以用一个边带来传输全部消息。单边带解调通常采 用相干解调的方式，它使用一个同步解调器，即由相乘器和低通滤波器组成。在解调过程 中，输入信号和噪声可以分别单独解调。相干解调的原理框图如图a所示：S SSB (t ) S1 (t )低通滤波器S 2 (t ) c(t)图a相干解调原理框图 此图表 示单边带信号首先乘以一个同频同相的载波，再经过低通滤波器即可还原信号。单边带信 号的时域表达式为S SSB (t ) = 1 1 ? m(t ) cos w c t m m(t ) sin w c t 2 2 m(t )表示基带信 号 其中取“-”时为上边带，取“+时为下边带。乘上同频同

4、相载波后得S1 (t ) = S SSB (t ) cos w c t = 1 1 1 ? m(t ) + m(t ) cos 2w c t m m(t ) sin 2w c t 4 4 4 ? m(t ) 表示 m(t )的希 尔伯特变换经低通滤波器可滤除2 wc的分量，所得解调输出为S 2 (t ) = 1 m(t ) 4由此 便可得到无失真的调制信号。2、SSB解调的实现、SSB信号的产生在本次课程设 计中，我选用频率为50Hz,初相位为0的余弦信号m作为原始基带信号， 在中表示为m=cos(2*pi*50.*t);” ；设置载波信号c为频率500Hz的余弦信号 程序中表示为fc=500

5、; c=cos(2*pi*fc.*t);”。在调制过程的实现中，通过 语句 b=sin(2*pi*fc.*t); lssb=m.*c+imag(hilbert(m).*b;产生下边带信号。接下来的解调过程主 要就针对该下边带信号进行。SSB解调实现根据相干解调原理，在解调的实现中，我先用下边带调制信号乘以本地载波， 再通过自己设计的低通滤波器滤除高频成分，如此便得到了解调后的信号。相关的程序语 句 为：y=lssb.*c; Y,y,dfl=fft_seq(y,ts,df); Y=Y/fs; n_cutoff=floor(fl/dfl); H=zeros(size(f); H(1:n_cutof

6、f)=2*o nes(1,n_cutoff);H(le ngth(f)-n_cutoff+1:le ngth(f)=2*o nes(1,n_cutoff);DEM=H.*Y; dem=real(ifft(DEM)*fs;在此段程序中，y是下边带信号与本地载波相乘的结果, 用数学表达式可以表示为 1 1 1 ? y = lssb(t ) cos w c t = m(t )+ m(t ) cos 2w c t m m(t ) sin 2w c t 4 4 4根据此式可知，此低通滤波器的作用在于去除y中包含的具有频率为2wc的 分量，滤波后就得到了基带信号。Y是对y进行傅立叶变换求得的频谱函数。通过

7、“Y=Y/fs” 对结果进行缩放，用于满足显示要求和计算需要。低通滤波器的设计是由 “n_cutoff=floor(fl/dfl);”到“DEM=H*Y之间的语句实现的，其中“fl”表示低通滤波器的截止频 率，混频信号y中高于此频率的分量将被滤除，得到解调信号；而“DEM=H.*Y;”表示Y与 截止系数H相乘后，其高频部分被滤去，所得结果DEM即为过滤后的频谱函数，就是 解调后信号的频谱。最后再由“dem=real(ifft(DEM)*fs;”对DEM进行傅立叶逆变换，取其 实部，乘以采样频率后可得出滤波后的时域信号函数dem，它就是解调后的信号。另外 要注意对采样频率fs的选取。由于仿真模拟

8、信号是通过抽样来确定此信号的， 所以能否通过抽样完全确定信号直接关系到仿真的成败。我所设的载波频率为fc=500Hz， 根具抽样定理，要使得载波信号c被所得到的抽样值完全确定，必须以T 2fc。故fs最小值应为2fc=1OOOHz,为了确保结果 合理，我取fs=2500Hz。其他一些参数设置为：信号保持时间t0=0.5s，滤波器截止频率fl=700Hz,频率分辨率df=0.1Hz。在求频谱函数时，我调用了一个自己设定的函数 “fft_seq”，它是对快速傅立叶变换函数fft的改进。它通过提高频率分辨率，能更好地区分 两个谱峰。若设信号m的长度为T秒，则其频率分辨率为?f =1/T(Hz)。因为

9、输入计 算机中的连续信号被转换成了离散的序列，而设其序列共有N个点，则两根谱线间的距 离?f =fs/N，它也是频率分辨率的一种度量。如果该式的?f不够小，我们可以增加N来 减小?f o而一种有效的方法是在离散序列后面补零，这不会破坏原来的信号。函数“fft_seq” 就是基于以上思想编成。我在中实现SSB解调的仿真就是基于以上的程序和 参数设定，运行调用了 plot函数的M文件，可以分别得出解调前后的频域与时域对比图， 为了便于参照，将LSSB信号的时域与频域图也一同进行比较。三者时域图比较结果如 图b所示：图b调制信号、LSSB信号、解调信号的时域对比图由图b可以看出，解调信号能够基 本还

10、原原始基带信号，只有少许失真，可见解调过程较为成功。对于失真接问题接下来将 会进行具体分析。为了更好分析解调的结果，还必须将时域图与频域图进行比较，通过观 察频域图能够更为直观的了解信号的频率分布，借此来讨论滤波器的作用特性。三者的频 域图比较如图c所示：图c调制信号、LSSB信号、解调信号的频域对比图由频域图 可以明显地看出，下边带信号将基带信号的频谱向两侧进行了搬移，然后去除了其中上边 带的那一半。而解调后的频谱与基带信号在频率小于700Hz的范围内几乎没有区别，用 肉眼很难观察出不同之处；解调信号在频率大于700Hz范围外的频谱已经被低通滤波器 完全滤除，故没有波形，这就是解调信号与基带

11、信号频谱差异最明显之处。由于在理论上， 解调过程可以完全恢复基带信号，但实际中仍然会产生少许失真，故失真的原因可以从参数的设置上入手找寻。在解调过程中，相乘器是严格按照数学关系 运算的，不会产生误差，失真的产生因该在低通滤波器环节。为了探究低通滤波器的设置 对解调的影响，我将滤波器的截止频率fl分别改为1000Hz和50Hz进行比较，根据结 果来判断原因。首先我选fl=1000Hz，它等于2fc,它应该无法滤除混频信号中的2 ec分 量。实际得出的时域和频域图结果如图d和图e：图d fl=1000Hz时的解调信号时域对 比图图e fl=1000Hz时的解调信号频域对比图 从时域比较图中可以明显

12、看出解调信号严重失 真，还保留着载波信号的频率变化；而频谱图中，解调信号频谱的两侧将近1000Hz处有 两个高峰，可见解调信号频谱与基带信号频谱存在巨大差异。原因很简单，低通滤波器的 截止频率高于2fc，故无法滤除相乘器输出的混频信号y的高频2 ec分量，完全没有起 到低通滤波的作用，所以解调信号保留有高频分量，无法正确还原调制信号。接下来将截 止频率fl设为50Hz，观察在fl较低的情况下的解调输出结果。此时，fl远远小于2fc, 不存在无法滤除高频分量的问题，不会产生上一种情况的错误。通过分析这种情况下的结 果可以得出低截止频率对解调的影响。在fl=50Hz时得出的时域和频域解调结果对比如

13、图 f和图g所示：图f fl=50Hz时解调信号时域对比图图g fl=50Hz时解调信号频域对比图由时域比较图中可以清晰地看出解调信号能够保持 平滑的余弦波形，但幅度小了将近一倍，且两端的幅度明显小于中间。由此可以估计出一 种可能性：基带信号的一部分丢失了，解调信号只是原基带信号的一部分。因此才会造成 幅度下降的结果。频谱对比图很好地验证了这种可能性。观察发现解调信号的频谱被滤去 频率大于50Hz的部分之后，虽然保持了原基带信号的主体部分，但两边的分量已经全部 滤除。这 就是问题所在。两边的分量也占据了一定的比重，滤除的越多，原基带信号损失 的部分也越多，就如时域图中所示。幅度减小的原因就在于

14、此。通过在fl=50Hz时对解 调信号频域和时域的观察，可以得出结论：fl较低时，在滤去高频2 3C分量的同时，也 会连带去除基带信号中的一部分。因此为了保全基带信号的完整性，fl不宜过低。综上所 述，截止频率fl过高时会导致高频2 wc分量无法滤除；而fl较低，又会使得基带信号 有部分被滤去。经过调整分析，我认为fl应该在能够滤去2 wc分量的情况下，尽可能高一点。这样就能保证基带信号被滤去的部分减到最少。所以开始所设的 fl=700Hz是一个较为合适的值，解调结果也很接近基带信号。实现解调还有一种方法，就 是调用专用于SSB信号解调的函数ssbdemod(u,fc,fs),此函数内部集成了

15、解调过程，不 需要自行再设置相乘器与滤波器。为了与前一种方法对比，现简要用此函数来实现SSB的 解调。为了增强对比性，除了 fl无需再设置外，其它的参数设置均与上一种方法相同：采 样频率fs=2500Hz，信号保持时间t0=0.5s，频率分辨率df=0.1Hz,且调制信号不变。在 本程序中，实现调制过 程使用语句u=ssbmod(m,fc,fs);” ，解调单边带信号则用 y=ssbdemod(u,fc,fs);”来实现。运行本程序，所得基带信号、LSSB信号、解调信号的时 域图对比情况如图h所示：图h方法2中基带、LSSB、解调信号的时域对比图解调信号通过中函数ssbdemod()实现了解调

16、并通过函数本身内设的低通滤波器，解调出了基带信号m(t)o解调后的信号基本与基带信号相同，只是边缘处有少许失真。再观察三者的频域比较如图i所示：图i方法2中基带、LSSB、解调信号的频域对比图由 频域图中可以看出，解调信号的频谱与基带信号很接近，不包含相乘器输出的混频信号中 的高频2 wc分量，也没采用第一种方法中设定截止频率，将高于一定频率的部分完全滤 除的做法。但仔细观察可以发现，解调信号中分布在频率较高处的波形比相同频率处的基 带信号波形起伏更大一些，这应该是本解调方式的缺陷，会造成一定误差。在时域图中的 解调信号边缘失真的情况或许来源于此。这种方法不用自己设置滤波器，程序简单，在实 现

17、解调时较为简便。前一种方法可以根据自己的需要调整滤波器的截止频率，从而得到更精确的解调信号。这两种解调方式 可谓各有千秋。3、叠加噪声时的SSB解调在实际信号传输过程中，通信系统不可避免 的会遇到噪声，例如自然界中的各种电磁波噪声和设备本身产生的热噪声、散粒噪声等， 它们很难被预测。而且大部分噪声为随机的高斯白噪声，所以在设计时引入噪声，才能够 真正模拟实际中信号传输所遇到的问题，进而思考怎样才能在接受端更好地恢复基带信号。 信道加性噪声主要取决于起伏噪声，而起伏噪声又可视为高斯白噪声，因此我在此环节将 对单边带信号添加高斯白噪声来观察噪声对解调的影响情况。为了具体而全面地了解噪声 的影响问题

18、，我将分别引入大噪声(信噪比为3dB)与 小噪声(信噪比为30dB)作用于 单边带信号，再分别对它们进行解调，观察解调后的信号受到了怎样的影响。在此过程中， 我用函数“awgn(x,snr,measured来添加噪声，此函数功能为向信号X中添加信噪比为snr 的高斯白噪声。其中的一个参数measured的功能是在叠加噪声前先测量信号x的功率， 再根据所测结果来产生所设信噪比的噪声。具体实现语句为：snr1=3; snr2=30; si =awgn(lssb,sn门,measured); s2 = awgn(lssb,snr2,measured);其中 sn门 表示大信噪比， snr2表示小信噪

19、比，s1表示叠加大噪声后的信号，s2表示 叠加小噪声后的信号。再通过 以下傅立叶变换程序来求出它们的频谱图：S1,s1,dfl=fft_seq(s1,ts,df); S1=S1/fs; S2,s2,dfl=fft_seq(s2,ts,df); S2=S2/fs;为了便于比较，我先显示了下边带信号加入两种噪 声后的时频波形图，再将基带信号、叠加大噪声后的解调信号以及叠加小噪声后的解调信 号放在一起，分别比较三者的时域与频域图。运行程序，求得下边带信号加入两种噪声后的时频波形图如图j所示：图j下边带信号加 入两种噪声后的时频波形比较图基带信号、叠加大噪声后的解调信号以及叠加小噪声后的 解调信号三者

20、之间的时域对比图如图k所示：图k基带信号、叠加大小噪声后解调信号的时域对比图基带信号、叠加大噪声后的解调信 号以及叠加小噪声后的解调信号三者之间的频域对比图如图l所示：图I基带信号、叠加 大小噪声后解调信号的频域对比图由图i可知，首先下边带信号因为叠加噪声而发生了变 化。再进一步观察图1中的时域变化：当叠加大噪声时，解调波形出现明显的紊乱，产生 了很多不规则的毛刺杂波，与基带信号相去甚远；而叠加小噪声时，波形图只是局部稍有 不平滑现象，总体还是与未加噪声时的解调信号相同，出入不大。再比较图k中各频谱之 间的关系，可以清晰地看出，加大噪声后，解调信号的频谱波形两边杂乱无章，起伏远大 于未加噪声时

21、的波形；而加小噪声后，解调信号频谱基本与未加噪声时一致，看不出明显 变化。造成此现象的原因是当信噪比较小时，噪声的功率在解调信号中所占比重较大，所以 会造成杂波较多的情况；而信噪比很大时，噪声的功率在解调信号中所占比重就很小了， 噪声部分造成的杂乱波形相对就不是很明显，甚至可以忽略。综上所述，叠加噪声会造成解调信号的失真，信噪比越小，失真程度越大。所以当信噪比 低于一定大小时，会给解调信号带来严重的失真，导致接收端无法正确地接收有用信号。 所以在解调的实际应用中，应该尽量减少噪声的产生。4、遇到的问题及解决办法 我在 该课程设计中曾遇到过一些问题，具体如下：其一是频谱显示问题。起初我得出的频谱

22、显 示波形并非是频率为0点居中，而是0点在最左端。这样不利于我观察频谱特性，为了能 让频谱图居中显示，我对程序设置进行分析，发现关键在于频率矢量的设置上。起初所设 置的值f = (0:length(U)-1)*fs/length(U),不能将图形显示在0Hz两边，经过改进为f = (0:length(U)-1)*fs/length(U)-fs/2,使得频谱居中显示成功。其中-fs/2就是将图形向右搬移 半个长度，这样0点就能居中了。其二是滤波器截止频率fl的设置问题。开始我将截止频 率设置为fl=50Hz,但仿真结果显示解调信号的幅度明显小于基带信号。这说明基带信号 没有完全被还原，有部分损失

23、。经分析我认为是滤波器的截止频率太低，将基带信号的成 分也滤去了一些。然而滤波器的截止频率也不宜过大，这样会使它无法滤除混频信号中的 2 3C高频分量。对于其合适的参数我进行了多次测试，力求在可以滤去2 3C高频分量的 情况下使解调信号损失最小。最终我发现将fl设为700Hz较为合理，解决了此问题。四、 老师的问题1、SSB信号是如何得到？有两种方法得到SSB信号，分别是相移法和滤波 法，本次设计用的是相移法。2、SSB信号是如何解调的？是用相干解调实现的，在解调 的实现中，我先用下边带调制信号乘以本地载波，再通过低通滤波器滤除高频成分，如此 便得到了解调后的信号。111?用数学表达式可以表示

24、为y = lssb(t ) cos w c t = m(t ) + m(t ) cos 2w c t m m(t ) sin 2w c t 4 4 4 3 、哪些语句是生成噪声的？我用函数 awgn(x,snr,measured)”来添加噪声，此函数功能为向信号x中添加信噪比为snr的高斯 白噪声。其中的一个参数measured的功能是在叠加噪声前先测量信号x的功率，再根据 所测结果来产生所设信噪比的噪声。具体实现语句为：snr1=3; snr2=30; si = awgn(lssb,sn门,measured); s2 = awgn(lssb,snr2,measured);五、 总结在这次课程

25、设计中，我通过多方面地搜集资料，对SSB的解调原理进行深入理解分析, 成功地用7.1编写出M程序完成了对SSB解调过程的仿真实现。在设计过程中，我面对着一系列的难点：将与通信原理结合起来运用，根据书中介绍的解调原理自主设计所要运用的器件并调整其参数，对解调信号结合理论进行分析并总结，掌握产 生具有一定信噪比的高斯白噪声的方法，等等。这些具体的问题在书上只有少量的理论介 绍，要在实际运用中结合这些理论是一件不太容易的事。接到课程设计任务书后，我积极 投入，去图书馆借阅相关资料，上网浏览相关信息，获得了大量的信息，为完成本次课程 设计奠定了基础。对于我开始还有许多未解之处，为了了解各个函数的用法，

26、我多次查阅了其中的help文件，大体了解了所调用函数的使用方法。然后通过不断的实验 与探究总结，终于逐个克服了这些难点，完成了任务。通过这次的实践，我明白了要将理 论与实际相结合的道理，尽管这个过程会有一些辛苦，但通过努力实现后，就能大大深化 我对知识的理解程度，增长实践经验。这表现在我对SSB的解调原理、的功能特性都有了进一步的认识。总之，我在本次课程设计中辛苦不少，获益良多。参考文献1东方渔民无线电广播与收音机发展的历史回眸一一纪念费森登发明无线电广 播100 周 年.广 播 之 家 网 站，. 2樊昌信，张甫翊，徐炳祥 吴成柯国防工业出版社，2007. . . 3孙学军，王秉均.通信原理

27、电子工业出版社，2005.4孙祥，徐流美，吴清.7.0基础教程北京：清华大学出版社2006. . . . 5张辉，曹丽娜现代通信原理与技术西安：电子科技大学出版社，2007. . . . 6唐向宏，岳恒立, 邓雪峰及在电子信息类课程中的应用电子工业出版社， 2006.附录1： SSB解调程序方法1清单程序名称：ssb.m %程序功能：自行设计解调器实 现 SSB 解调 %= = t0=0.05; ts=1/2500; fc=500; fs=1/ts; fl=700; df=0.1; snr1=3; snr2=30; t=O:ts:tO; m=cos(2*pi*50.*t); c=cos(2*p

28、i*fc.*t);b=si n( 2*pi*fc.*t);lssb=m.*c+ imag(hilbert(m).*b;s1 =awgn(lssb,sn门,measured); s2 = awgn(lssb,snr2,measured); y=lssb.*c; y1=s1.*c; y2=s2.*c; M,m,dfl=fft_seq(m,ts,df); M=M/fs; LSSB,lssb,dfl=fft_seq(lssb,ts,df);频谱 LSSB=LSSB/fs; S1,s1,dfl=fft_seq(s1,ts,df);频谱 S1=S1/fs; S2,s2,dfl=fft_seq(s2,ts,

29、df); % 缩放求下边带信号加小噪声的缩放求下边带信号加大噪声的缩放傅里叶变 换求下边带信号下边带的混频加大噪声后混频加小噪声后混频傅里叶变换求 原始信号频谱时间矢量原始信号载波信号载波信号对应的正弦信号下边带 信号信号保持时间采样时间间隔此为载波频率此为抽样频率滤波器的截止 频率频率分辨率频谱 S2=S2/fs; Y,y,dfl=fft_seq(y,ts,df); Y=Y/fs; Y1,y1,dfl=fft_seq(y1,ts,df); Y1=Y1/fs; Y2,y2,dfl=fft_seq(y2,ts,df); Y2=Y2/fs; n_cutoff=floor(fl/dfl); f=O:

30、dfl:dfl*(length(m)-1)-fs/2; H=zeros(size(f);H(1:n_cutoff)=2*o nes(1,n_cutoff);H(le ngth(f)-n_cutoff+ 1:le ngth(f)=2*o nes(1,n_cutoff);DEM=H.*Y; DEM1=H.*Y1;DEM2=H.*Y2;dem=real(ifft(DEM)*fs;dem1=real(ifft(DEM1)*fs;dem2=real(ifft(DEM2)*fs; figure(1); subplot(3,1,1);plot(t,m(1:length(t); title(调制信号);sub

31、plot(3,1,2);plot(t,lssb(1:le ngth(t);title(下边带信号)subplot(3,1,3);plot(t,dem(1:le ngth(t);title(解调信号）；figure(2)subplot(3,1,1);plot(f,abs(fftshift(M);title(调制信号 频谱)subplot（3,1,2）;plot（f,abs（fftshift（LSSB）; title（下边带信号频谱）;%显示下边带信号频谱 % 显示调制信号频谱显示解调信号显示下边带信号显示调制信号滤波器输出（即 解调后信号）频谱滤波器输出（即加大噪声解调后信号）频谱滤波器输出（即

32、加小噪声解 调后信号）频谱滤波器输出（即解调后）信号滤波器输出（即大噪声解调后）信号滤波 器输出（即小噪声解调后）信号缩放设计滤波器频率矢量缩放傅里叶变换得 出混频频谱缩放傅里叶变换得出混频频谱缩放傅里叶变换得出混频频谱subplot(3,1,3);plot(f,abs(fftshift(DEM); title(解调信号频谱);figure(3); subplot(3,2,1);plot(t,lssb(1:le ngth(t);title( 下 边 带 信 号 );subplot(3,2,2);plot(f,abs(fftshift(LSSB); title(下边带信号频谱); subplot

33、(3,2,3);plot(t,s1(1:length(t); title(加大噪声后下边带信号)； subplot(3,2,4);plot(f,abs(fftshift(S1); title(加大噪声后下边带信号频谱); subplot(3,2,5);plot(t,s2(1:length(t); title(力口 小噪声后下边带信号)； subplot(3,2,6);plot(f,abs(fftshift(S2); title(加小噪声后下边带信号频谱);figure(4); subplot(3,1,1);plot(t,dem(1:le ngth(t); title( 解 调 信 号 ); s

34、ubplot(3,1,2);plot(t,dem1(1:le ngth(t); title(加大噪声解调信号)； subplot(3,1,3);plot(t,dem2(1:length(t); title(力口 小噪声解调信号);figure(5); subplot(3,1,1);plot(f,abs(fftshift(DEM); title(解 调 信 号 频 谱); subplot(3,1,2);plot(f,abs(fftshift(DEM1); title(加大噪声解调信号频谱); subplot(3,1,3);plot(f,abs(fftshift(DEM2); title(加小噪声

35、解调信号频谱)；显示解调信号频 谱显示下边带信号显示下边带信号频谱显示加大噪声后下边带信号显示加大 噪声后下边带信号频谱显示加小噪声后下边带信号显示加小噪声后下边带信号频谱% 显示解调信号显示加大噪声解调信号显示加小噪声解调信号显示解调信号频谱% 显示加大噪声解调信号频谱显示加小噪声解调信号频谱附录2:傅立叶变换函数清单%程序名称：fft_seq.m %程序功能：傅立叶变换函数= fun cti onM,m,df=fft_seq(m,ts,df) fs=1/ts; if nargin=2, else, n1 =fs/df; end n2=length(m); n=2A(max(nextpow2

36、(n 1),nextpow2(n2); M=fft(m,n); m=m,zeros(1,n-n2); df=fs/n; %对 m补零最终频率分辨率n 1=0; %定义函数采样频率设计函数功能附录3： SSB 解调程序方法2清单%程序名称：ssb2.m %程序功能：调用函数 ssbdemod实现 SSB 解 调 %= = t0=0.05; ts=1/2500; fc=500; fs=1/ts; df=0.1; t=0:ts:t0; m=cos(2*pi*50.*t); u=ssbmod(m,fc,fs); y=ssbdemod(u,fc,fs); M,m,dfl=fft_seq(m,ts,df)

37、; M=M/fs; U,u,dfl=fft_seq(u,ts,df); 谱 U=U/fs; Y,y,dfl=fft_seq(y,ts,df); Y=Y/fs; f=0:df:df*(le ngth(m)-1)-fs/2; figure(1); subplot(3,1,1);plot(t,m(1:length(t); title(调制信号);subplot(3,1,2);plot(t,u(1:length(t); title(单边带调制信号);subplot(3,1,3);plot(t,y(1:length(t); title(解调信号);figure(2); subplot(3,1,1);plot(f,abs(fftshift(M); title(调 制 信 号 频 谱); subplot(3,1,2);plot(f,abs(fftshift(U); %显示单边带信号频谱 显示调制信号频谱 显示 解调信号显示单边带信号显示调制信号缩放频率矢量缩放傅里叶变换求 解调信号频谱缩放傅里叶变换求单边带信号频信号保持时间采样时间间隔% 此为载波频率此为抽样频率频率分辨率时间矢量原始信号调制解调% 傅里叶变换求原始信号频谱title(单边带信号频谱);subplot(3,1,3);plot(f,abs(fftshift(Y); title(解调信号频谱);% 显示解 调信号频谱