北邮通原软件实验
北邮通信原理软软件实验报告
实验八：一、实验目的假设基带信号为m(t)=sin(2000*pi*t)+2cos(1000*pi*t)，载波频率为20kHz，请仿真出AM，DSB-SC、SSB信号，观察已调信号的波形和频谱。
二、实验模型基带信号m(t)可以分成两个信号的叠加，分别记为m1(t),m2(t)。
借助公式s DSB-SC=m(t)cos(2*pi*fc*t),S AM=(1+m(t))cos(2*pi*fc*t),s SSB=m(t)cos(2*pi*fc*t)+H[m(t)]sin(2*pi*fc*t)分别仿真出m1(t)和m2(t)的信号波形，然后叠加便可以得到m(t)的波形和频谱三、仿真设计设计程序时先确定采样点、采样频率，然后分别表示出m1(t)和m2(t)的表达式，然后表示出后面仿真SSB信号所需要的两个信号的希尔伯特变换表达式。
其中表示希尔伯特变换时，采用的方法是先表示出频域的形式MH1和MH2，然后再傅里叶反变换得出对应的mh1和mh2。
对应代码如下：m1=sin(2*pi*fm1*t);M1=t2f(m1,fs);MH1=-j*sign(f).*M1;mh1=real(f2t(MH1,fs));m2(t)信号做相同的处理。
处理完信号后，就利用上述的三个公式，表示出AM、DSB-SC和SSB信号s1、s2和s3和其对应傅里叶变换得到其频谱S1 、S2、S3。
为了方便实验结果的观察与对比，将这三组图处理在一张图内，利用的函数是subplot。
四、实验结果五、分析讨论由实验结果可见，AM与DSB-SC相比，频谱多了一个离散的大载波直流分量，而且DSB-SC信号波形会有相位翻转的现象出现；而DSB-SC和SSB相比，SSB信号的频谱是DSB-SC的一个边带，本实验中采用的上边带滤波。
可见实验结果与理论结果是相一致的。
六、思考题1.如何仿真VSB系统？答：将残留边带滤波器用M文件实现，然后当做函数使用，在程序中调用。
北邮通信原理实验报告
北京邮电大学通信原理实验报告学院：信息与通信工程学院班级：姓名：姓名：实验一：双边带抑制载波调幅（DSB-SC AM）一、实验目的1、了解DSB-SC AM 信号的产生以及相干解调的原理和实现方法。
2、了解DSB-SC AM 信号波形以及振幅频谱特点,并掌握其测量方法。
3、了解在发送DSB-SC AM 信号加导频分量的条件下,收端用锁相环提取载波的原理及其实现方法。
4、掌握锁相环的同步带和捕捉带的测量方法,掌握锁相环提取载波的调试方法。
二、实验原理DSB 信号的时域表达式为()()cos DSB c s t m t tω=频域表达式为1()[()()]2DSB c c S M M ωωωωω=-++其波形和频谱如下图所示DSB-SC AM 信号的产生及相干解调原理框图如下图所示将均值为零的模拟基带信号m(t)与正弦载波c(t)相乘得到DSB—SC AM信号，其频谱不包含离散的载波分量。
DSB—SC AM信号的解调只能采用相干解调。
为了能在接收端获取载波，一种方法是在发送端加导频，如上图所示。
收端可用锁相环来提取导频信号作为恢复载波。
此锁相环必须是窄带锁相，仅用来跟踪导频信号。
在锁相环锁定时，VCO输出信号sin(2πf c t+φ)与输入的导频信号cos(2πf c t)的频率相同，但二者的相位差为(φ+90°)，其中很小。
锁相环中乘法器的两个输入信号分别为发来的信号s(t)（已调信号加导频）与锁相环中VCO的输出信号，二者相乘得到[A C m(t)cos(2πf c t)+A p cos(2πf c t)]∙sin(2πf c t+φ)=A c2m(t)[sinφ+sin(4πf c t+φ)]+A p2[sinφ+sin(4πf c t+φ)]在锁相环中的LPF带宽窄，能通过A p2sinφ分量，滤除m(t)的频率分量及四倍频载频分量，因为很小，所以约等于。
LPF的输出以负反馈的方式控制VCO,使其保持在锁相状态。
北邮通信原理软件实验报告
通信原理软件实验报告学院：信息与通信工程学院班级：一、通信原理Matlab仿真实验实验八一、实验内容假设基带信号为m(t)=sin(2000*pi*t)+2cos(1000*pi*t),载波频率为20kHz，请仿真出AM、DSB-SC、SSB信号，观察已调信号的波形和频谱。
二、实验原理1、具有离散大载波的双边带幅度调制信号AM该幅度调制是由DSB-SC AM信号加上离散的大载波分量得到，其表达式及时间波形图为：应当注意的是，m(t)的绝对值必须小于等于1，否则会出现下图的过调制：AM信号的频谱特性如下图所示：由图可以发现，AM信号的频谱是双边带抑制载波调幅信号的频谱加上离散的大载波分量。
2、双边带抑制载波调幅（DSB—SC AM）信号的产生双边带抑制载波调幅信号s(t)是利用均值为0的模拟基带信号m(t)和正弦载波c(t)相乘得到，如图所示：m(t)和正弦载波s(t)的信号波形如图所示：若调制信号m(t)是确定的，其相应的傅立叶频谱为M(f)，载波信号c(t)的傅立叶频谱是C(f),调制信号s(t)的傅立叶频谱S(f)由M(f)和C(f)相卷积得到，因此经过调制之后，基带信号的频谱被搬移到了载频fc处，若模拟基带信号带宽为W，则调制信号带宽为2W，并且频谱中不含有离散的载频分量，只是由于模拟基带信号的频谱成分中不含离散的直流分量。
3、单边带条幅SSB信号双边带抑制载波调幅信号要求信道带宽B=2W, 其中W是模拟基带信号带宽。
从信息论关点开看，此双边带是有剩余度的，因而只要利用双边带中的任一边带来传输，仍能在接收机解调出原基带信号，这样可减少传送已调信号的信道带宽。
单边带条幅SSB AM信号的其表达式：或其频谱图为：三、仿真设计1、流程图：Array2、实验结果&分析讨论实验仿真结果从上至下依次是AM信号、DSB信号、SSB信号。
从仿真结果看，AM调制信号包络清晰，可利用包络检波恢复原信号，接收设备较为简单。
北邮通信原理软件实验报告
北邮通信原理软件实验报告北邮通信原理软件实验报告一、实验目的本次实验旨在加深对通信原理知识的理解，并通过实际操作掌握常用通信原理技术。
二、实验内容及原理本次实验分为三项实验内容：1、FSK解调实验FSK是一种通信调制技术，在数字信号传输领域中比较常用。
FSK解调实验中，我们使用MATLAB软件编写程序，模拟FSK解调过程，了解解调过程中的基本原理。
FSK是通过改变载波频率来传输信息的调制技术。
在数字信号的传输中，我们一般将数字信号分为两种，0和1，然后分别将它们对应到两个频率上，再将这两个频率进行交错发送，接收方通过检测频率的变化来判断发送方的信息。
在FSK解调实验中，我们使用的解调技术是匹配滤波器法。
解调的过程是将接收到的信号经过低通滤波器，合并成一个信号。
2、QAM解调实验QAM是一种把两路模拟信号叠加的数字调制技术，它是组合了ASK和PSK的数字传输技术。
QAM解调实验中，通过MATLAB软件仿真的方法，模拟QAM解调过程，了解解调过程中的基本原理。
QAM技术是将两路数模合成的模拟信号进行数字化处理，将两路模拟信号进行分别调制成两个独立的数字信号，然后将这两个数字信号通过载波同步合成一个数字信号进行传输。
在QAM解调实验中，我们使用的解调技术是相干解调。
解调的过程是将接收到的信号经过相干解调器解调，得到原始的时域信号，然后通过低通滤波器进行滤波。
3、OFDM调试实验OFDM技术是目前广泛应用于高速数据传输的一种技术，它是通过将信号分成多个子载波进行传输，提高频率利用率，并实现抗多径衰落的效果。
OFDM调试实验中，我们通过软件界面和Matlab代码相结合，模拟OFDM调制和解调过程，了解其中的基本原理。
OFDM技术是通过将原始信号分成多个子信道，每个子信道独立传输，最终将其合并成整个信号。
因此，在OFDM模式下，每个子信道的公共频率就成为可利用的带宽，提高了传输率并减少了所需的带宽。
在OFDM调试实验中，我们使用了MATLAB软件进行调制和解调。
北京邮电大学信息与通信工程学院实验报告
第二学期《通信原理软件》实验报告专业班级姓名学号开课系室报告日期目录实验一声音播放和滤波 (6)实验二时域仿真精度分析 (17)实验三频域仿真精度分析 (21)实验四噪声产生 (27)实验五取样和重建 (37)实验八抑制双边带的调制与解调 (46)实验十二ASK调制与解调 (53)综合实验数字基带系统仿真 (61)实验一 声音播放和滤波实验目的掌握声音播放模块，FIR 滤波器的使用方法。
深入理解频率的大小如何影响声音的变化，FIR 滤波器的特性和模块各参数的作用。
主要功能：利用FIR 滤波器实现对声音信号中某些频率的滤除，频带内的信号分量通过，频带外的滤除，来观察声音的变化。
实验原理声音频率的大小影响声音的音调，即频率越高我们听到的声音越尖，频率越低听到的声音越低沉；因此设置声音信号的频率为不同的频率时，我们可以通过声音播放器发出的声音感知声音音调的变化。
FIR 滤波器又名为有限冲激响应滤波器，利用离散时间系统的特性来对输入信号的波形或频谱进行加工处理，本实验中即利用该滤波器实现对声音信号中某些频率的滤除，频带内的信号分量通过，频带外的滤除，来观察声音的变化。
10||2,111()1cos ,||22220,1||2ss s s r s s s s a f T T T T a a a H f f f T T T a f T πα-≤≤⎧⎪⎧⎫⎡⎤⎛⎫--+⎪⎪⎪=+-≤≤⎨⎨⎬⎢⎥ ⎪⎝⎭⎪⎪⎣⎦⎪⎩⎭⎪+⎩> 下图所示是满足上式的理想的FIR 滤波器的时域冲激响应和传递函数。
声音播放和滤波系统框图如图1.1所示：图1.1 声音播放和滤波系统框图实验方案所需元件：●正弦波发生器（sinusoid generator）（两个，①产生频率为200Hz的正弦波，②频率为800Hz的正弦波）●音频播放模块（Play sound，选自Scicom_sinks元件库）●触发时钟（CLOCK_c）●FIR滤波器（FIR Filter，选自Scicom_Filter原件库）●频谱示波器模块（FFT，选自Scicom_sinks元件库）具体步骤如下：将正弦波发生器、音频播放模块、触发时钟、频谱示波器模块按下图连接：打开Diagram菜单栏中的Context输入框，输入下图内容：设置正弦波模块，产生频率为200Hz的信号，设置Play sound模块的参数。
通原软件实验报告
信息与通信工程学院通信原理软件实验报告班级：201121xxxx姓名：xxx学号：序号：目录实验八 (4)一、实验内容 (4)二、实验原理 (4)三、仿真设计 (6)1. 仿真思路 (6)2. 程序框图 (6)3. 源程序 (7)四、实验结果及分析 (10)1．实验仿真结果 (10)2. 结果分析 (12)五、实验总结 (12)实验九 (13)一、实验内容 (13)二、实验原理 (13)三、仿真设计 (13)1. 仿真思路 (13)2. 程序框图 (14)3. 源代码 (14)四、实验结果及分析 (16)1. 实验仿真结果 (16)2. 结果分析 (17)五、实验总结 (17)实验十一 (18)一、实验内容 (18)二、实验原理 (18)1、单极性归零码 (18)2、双极性归零码 (18)3、各种码的比较 (19)三、仿真设计 (19)1. 仿真思路 (19)2. 程序框图 (20)3. 源代码 (20)四、实验结果及分析 (23)1. 实验仿真结果 (23)2. 结果分析 (24)五、实验总结 (24)实验十二 (25)一、实验内容 (25)二、实验原理 (25)三、仿真设计 (26)1. 仿真思路 (26)2. 程序框图 (26)3. 源程序 (26)四、实验结果及分析 (28)1. 实验仿真结果 (28)2. 结果分析 (29)五、实验总结 (29)实验八一、实验内容假设基带信号为()sin(2000)2cos(1000)m t t t ππ=+，载波频率为20kHz ，请仿真出AM 、DSB-SC 、SSB 信号，观察已调信号的波形及频谱。
二、实验原理1. 具有离散大载波的双边带幅度调制信号AM该幅度调制是由DSB-SC AM 信号加上离散的大载波分量得到，其表达式及时间波形图为：应当注意的是，m(t)的绝对值必须小于等于1，否则会出现下图的过调制：AM 信号的频谱特性如下图所示：由图可以发现，AM信号的频谱是双边带抑制载波调幅信号的频谱加上离散的大载波分量。
北邮通信原理软件实验报告
北邮通信原理软件实验报告一、实验目的本次实验的目的是通过使用软件进行通信原理实验，探究数字通信系统的原理和性能。
二、实验内容1.利用软件计算并绘制理想低通滤波器的频率响应曲线。
2.通过软件模拟并比较维纳滤波器与理想低通滤波器的频率响应曲线。
3.仿真带通调制器和解调器在理想信道中的性能。
三、实验步骤1.理想低通滤波器的设计：（1）利用软件，设置滤波器参数，如截止频率和滤波器类型。
（2）计算并绘制理想低通滤波器的频率响应曲线。
2.维纳滤波器与理想低通滤波器的比较：（1）利用软件设置维纳滤波器参数，如截止频率和信噪比。
（2）仿真并比较维纳滤波器与理想低通滤波器的频率响应曲线。
3.带通调制器和解调器的性能仿真：（1）设置带通调制器和解调器的参数，如载波频率和调制系数。
（2）仿真并分析带通调制器和解调器的性能，如频率响应和误码率。
四、实验结果五、实验分析通过本次实验，我们对数字通信系统的原理和性能有了更深入的了解。
首先，理想低通滤波器的频率响应曲线能够更清晰地展现滤波器的特性，帮助我们更好地了解滤波器的设计和应用。
其次，维纳滤波器相对于理想低通滤波器而言，频率响应存在一定的失真，但对于噪声有一定的抑制作用。
在实际应用中，需要根据具体需求选择适合的滤波器类型。
最后，带通调制器和解调器的性能仿真结果能够帮助我们评估系统的性能，如误码率和频率响应。
通过调整调制系数和载波频率，可以使系统在一定范围内具有较好的性能。
六、实验结论通过本次实验的软件仿真，我们探究了数字通信系统的原理和性能。
实验结果显示，理想低通滤波器具有良好的频率响应特性；维纳滤波器能对噪声进行一定的抑制，但频率响应存在一定的失真；带通调制器和解调器在适当的调制系数和载波频率下能够获得较低的误码率和良好的频率响应。
总之，本次实验通过软件仿真，使我们更好地理解了通信原理中的数字通信系统及其性能分析方法，提高了我们的实践能力和理论知识水平。
北邮通原软件实验报告
北邮通原软件实验报告北京邮电大学实验报告题目：班级：专业：姓名：成绩：实验1：抽样定理一．实验目的（1）掌握抽样定理（2）通过时域频域波形分析系统性能二．实验原理抽样定理：设时间连续信号m（t），其最高截止频率为fm ，如果用时间间隔为T抽样过程原理图（时域）重建过程原理图（频域）具体而言：在一个频带限制在（0，f h）内的时间连续信号f（t），如果以小于等于1/(2 f h)的时间间隔对它进行抽样，那么根据这些抽样值就能完全恢复原信号。
或者说，如果一个连续信号f（t）的频谱中最高频率不超过f h，这种信号必定是个周期性的信号，当抽样频率f S≥2 f h时，抽样后的信号就包含原连续信号的全部信息，而不会有信息丢失，当需要时，可以根据这些抽样信号的样本来还原原来的连续信号。
根据这一特性，可以完成信号的模-数转换和数-模转换过程。
三．实验步骤1.将三个基带信号相加后抽样，然后通过低通滤波器恢复出原信号。
实现验证抽样定理的仿真系统，同时在必要的输出端设置观察窗。
如下图所示2.设置各模块参数三个基带信号频率从上至下依次为10hz、20hz、40hz。
抽样信号频率fs设置为80hz，即2*40z。
（由抽样定理知，fs≥2fH）。
低通滤波器频率设置为40hz 。
设置系统时钟，起始时间为0，终止时间设为1s.抽样率为1khz。
3.改变抽样速率观察信号波形的变化。
四．实验结果五．实验建议、意见将抽样率fs设置为小于两倍fh的值，观察是否会产生混叠失真。
实验2：验证奈奎斯特第一准则一．实验目的（1）理解无码间干扰数字基带信号的传输；（2）掌握升余弦滚降滤波器的特性；（3）通过时域、频域波形分析系统性能。
二．实验原理基带传输系统模型奈奎斯特准则提出：只要信号经过整形后能够在抽样点保持不变，即使其波形已经发生了变化，也能够在抽样判决后恢复原始的信号，因为信息完全恢复携带在抽样点幅度上。
无码间干扰基带传输时，系统冲击响应必须满足x(nTs)=1(n=0); x(nTs)=0(n=!0)。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验一实验目的：假设基带信号为m(t)=sin(2000πt)+2cos(1000πt)，载波频率为20kHz,请仿真出AM，DSB-SC，SSB信号，观察已调信号的波形和频谱。
1．AM信号：（1）信号的表达式（3）流程图AM信号s= (1+0.3*m).*cos(2*pi*fc*t);绘制时域波形及频谱傅氏变换S= t2f(s,fs)（2）源代码%AM信号的产生fs= 800; %采样频率KHzT= 200; %截短时间msN= T*fs; %采样点数dt= 1/fs;t= [-T/2:dt:T/2-dt];df= 1/T;f=[-fs/2:df:fs/2-df];fm= 1; % kHzfc= 20; % kHzm= sin(2*pi*fm*t)+2*cos(1*fm*pi*t);s= (1+0.3*m).*cos(2*pi*fc*t); %AM 信号S= t2f(s,fs);figure(1)plot(f,abs(S1))title('AM信号频谱')xlabel('f')ylabel('S(f)')axis([-25,25,0,max(abs(S1))]);%xset('window',2)figure(2)plot(t,s1)title('AM信号波形')xlabel('t')ylabel('s(t)')axis([-3,3,-3,3]);（4）实验结果精选文库-3-2-10123-3-2-1123AM 信号波形t(ms)s (t )-25-20-15-10-505101520250102030405060708090100AM 信号频谱f(kHz)S (f )精选文库2．DSB-SC信号（1）信号的产生和表达式（2）流程图DSC-SB信号s= m.*cos(2*pi*fc*t);傅氏变换S= t2f(s,fs)绘制时域波形及频谱（3）源代码fs= 800; %KHzT= 200; %msN= T*fs;dt= 1/fs;t= [-T/2:dt:T/2-dt];df= 1/T;f=[-fs/2:df:fs/2-df];fm= 1; % kHzfc= 20; % kHzm= sin(2*pi*fm*t)+2*cos(1*fm*pi*t);s= m.*cos(2*pi*fc*t); %DSB-SC 信号S= t2f(s,fs);figure(1)plot(f,abs(S2))title('DSB-SC信号频谱')xlabel('f')ylabel('S(f)')axis([-25,25,0,max(abs(S2))]);figure(2)plot(t,s2)title('DSB-SC信号波形')xlabel('t')ylabel('s(t)')axis([-1,4,-3,3]);（4）实验结果()()()()()cosc c cs t m t c t m t A tωϕ==+精选文库-1-0.500.511.522.533.54-3-2-1123DSB-SC 信号波形t(ms)s (t )-25-20-15-10-50510152025DSB-SC 信号频谱f(kHz)S (f )2 DSC-SB 频谱精选文库3．SSB 信号（1）信号的产生和表达式（2）流程图（3）源代码：%SSB 信号的产生 fs= 800; %KHz T= 200; %ms N= T*fs; dt= 1/fs;t= [-T/2:dt:T/2-dt]; df= 1/T;f=[-fs/2:df:fs/2-df];fm= 1; % kHz fc= 20; % kHz m= sin(2*pi*fm*t)+2*cos(1*fm*pi*t); M= t2f(m,fs);MH=-j*sign(f).*M; %在频域进行希尔伯特变换 mh= real(f2t(MH,fs)); %希尔伯特变换后的信号 s= m.*cos(2*pi*fc*t)-mh.*sin(2*pi*fc*t); %SSB signal S= t2f(s,fs); figure(1)plot(f,abs(S3))title('SSB 信号频谱') xlabel('f') ylabel('S(f)')axis([-25,25,0,max(abs(S3))]) figure(2) plot(t,s3)title('SSB 信号波形') xlabel('t') ylabel('s(t)') axis([0,6,-3,3])()()()()µ()()()µ()cos 2sin 2cos 2sin 2DSB SC AM c c c c SSB c c c cs t s t s t A m t f t A m t f t s t A m t f t A m t f t ππππ--=-=-=m 下上（4）实验结果0123456-3-2-1123SSB 信号波形t(ms)s (t )-25-20-15-10-50510152025020406080100120140160180200SSB 信号频谱f(kHz)S (f )实验二实验目的：假设基带信号为m(t)=sin(2000πt)+2cos(1000πt)+4sin(500πt+π/3),载波频率为40kHz，仿真产生FM信号，观察波形与频谱，并与卡松公式作对照。
FM的频率偏移常数是5kHz/V。
（1）信号表达式（2）流程图FM信号phi= 2*pi*Kf*cumsum(m)*dt;s= cos(2*pi*fc*t+phi);傅氏变换S= t2f(s,fs)绘制时域波形及频谱设定采样频率、采样点数、时间截短等（3）源代码fs= 800; %kHzT= 16; %msN= T*fs;dt= 1/fs;t= [-T/2:dt:T/2-dt];df= 1/T;f= [-fs/2:df:fs/2-df];fm= 1; %kHzKf= 5; %kHz/Vfc= 40; %kHzm= sin(2*pi*fm*t)+2*cos(1*pi*fm*t)+4*sin(0.5*pi*fm*t+pi/3);phi= 2*pi*Kf*cumsum(m)*dt; %求相位s= cos(2*pi*fc*t+phi); % s(t)S= t2f(s,fs);figure(1)plot(f,abs(S).^2)title('FM信号功率谱')xlabel('f')ylabel('S(f)')axis([-80,80,0,max(abs(S).^2)]); %功率谱密度为|S|^2 figure(2) plot(f,abs(S))title('调制信号频谱') xlabel('f') ylabel('S(f)')axis([-80,80,0,max(abs(S))]); figure(3) plot(t,s)title('FM 信号波形') xlabel('t') ylabel('s(t)') axis([0,3,-2,2]);（4）实验结果00.511.522.53-2-1.5-1-0.500.511.52FM 信号波形t(ms)s (t )-80-60-40-20020406080调制信号频谱f(kHz)S (f )试验结论：fm 取1kHz,用卡松公式计算得到FM 信号带宽: Bfm=2*(Kf*max(abs(m))+1) = 66.8325与FM 频谱图比较，基本相等，说明实验FM 信号带宽与理论值基本相符。
时域图也可看到疏密不同的波形，符合FM 信号的特点。
实验三实验目的：通过仿真测量占空比为25%、50%、75%以及100%的单、双极性归零码波形及其功率谱。
clear all ; close all ;L=32; %每个码元间隔内的采样点数 N=2^13; %总采样点数 M=N/L; %总码元数 Rb=2; %码元速率 Ts=1/Rb; %比特间隔 fs=L/Ts; %采样速率 T=N/fs; %截断时间 Bs=fs/2; %系统带宽t=-T/2+[0:N-1]/fs; %时域采样点精选文库f=-Bs+[0:N-1]/T; %频域采样点L0=input('请输入占空比（0~1）：')EP=zeros(1,N);ch=input('请选择要观察的码型：1-单极性；2-双极性：')for loop=1:1000 %1000次样本函数取平均if ch==1a=(rand(1,M)>0.5); %生成单极性序列elsea=sign((rand(1,M)>0.5)-0.5); %生成双极性序列endtmp=zeros(L,M); %一个码元的归零部分取零L1=L*L0; %占空比，求出一个码元不归零部分的采样点数tmp([1:L1],:)=ones(L1,1)*a; %将一个码元不归零部分的取样点置为1s=tmp(:)';S=t2f(s,fs);P=abs(S).^2/T; %样本功率谱密度EP=EP*(1-1/loop)+P/loop; %随机部分的功率谱是各个样本功率谱的期望endfigure(1)plot(t,s)grid ontitle('时域图')xlabel('t')ylabel('S(t)')axis([-3,3,-1.5,1.5]);figure(2)plot(f,abs(EP+eps))grid ontitle('功率谱图形')xlabel('f')ylabel('功率')axis([-35,35,-5,max(EP+eps)]);figure(3)plot(f,10*log10(EP+eps))grid ontitle('功率谱图形(dB)')xlabel('f')ylabel('功率')实验结果：（1）.单极性-3-2-10123-1.5-1-0.50.511.5时域图t(ms)S (t ) (V )-30-20-100102030功率谱图形f(kHz)功率 (V 2/k H z )-40-30-20-10010203040-160-140-120-100-80-60-40-20020功率谱图形(dB)f(kHz)功率 (d B /k H z )修改占空比可得到以下图形-3-2-10123-1.5-1-0.50.511.5时域图t(ms)S (t ) (V )-30-20-100102030-551015功率谱图形f(kHz)功率 (v 2/k H z )-40-30-20-10010203040-160-140-120-100-80-60-40-20020功率谱图形(dB)f(kHz)功率 (d B /k H z )-3-2-10123-1.5-1-0.50.511.5时域图t(ms)S (t ) (V )-30-20-100102030-5051015202530功率谱图形f(kHz)功率 (v 2/k H z )-40-30-20-10010203040-160-140-120-100-80-60-40-20020功率谱图形(dB)f(kHz)功率 (d B /k H z )从上至下依次是占空比为50%、75%、100%的波形图及功率谱密度图。