西电 电院 随机信号分析

一、用matlab语言产生一个随机白噪声序列的样本序列X(n)，要求
3.用遍历性估计X(n)的自相关序列R X(m)，画出R X(m)的图像。

二、将一中产生的序列通过一个线性系统，其单位脉冲响应为h(n)=0.9n，n=0，
1，…，100
三、比较X(n)与Y(n)的幅度分布直方图，发生了什么变化。

分析其变化的原
因。

随机信号经过线性系统后，不会增加新的频率分量，但是输出的幅度和相位会发生变化。

白噪声X(n)的幅度基本相同，而Y(n)的幅度基本呈正态分布。

因为均匀白噪声是一种宽带非正态过程，所以通过一有限带宽线性系统后，输出Y(n)近似呈正态分布。

——via 1402011 赵春昊。

[键入公司名称]第一章1.23 上机题：设有随机初相信号X(t)=5cos（t+），其中相位是在区间（0,2）上均匀分布的随机变量。

试用Matlab编程产生其三个样本函数。

MATLAB源代码：clc,clear;o=2*pi*rand(1,3)for n=1:3t=0:.01:10;y=5*cos(t+o(n));figure(1);plot(t,y),grid on;hold on;endtitle('三个来自随机初相信号的样本函数');第二章2.22 上机题：利用MATLAB程序设计一正弦信号加高斯白噪声的复合信号。

（1）.分析复合信号的功率谱密度、幅度分布特性；MATLAB源代码：clccleart=0:0.001:0.1;s=3*sin(1e3*t);%正弦信号（幅度分布）figure(1);subplot(3,1,1),plot(t,s),grid on;title('原正弦信号');xlabel('t/s');ylabel('s');s1=fft(s);subplot(3,1,2),plot(t,abs(s1)),grid on;title('正弦信号幅度谱');xlabel('t/s');ylabel('s_fft');n=100;f=100;window=boxcar(length(s));[p1,f1]=periodogram(s,window,n,f)subplot(3,1,3),plot(f1,10*log10(p1));xlabel('f/Hz');ylabel('Gs');title('正弦信号功率谱');（2）分析复合信号通过RC积分电路后的功率谱密度和相应的幅度分布特性；MATLAB源代码：clccleart=0:0.001:0.1;s=3*sin(1e3*t);%正弦信号（幅度分布）y1=awgn(s,10,'measured');r=(1/(pi*pi))*10e6;c=(1/16)*10e-4;hw=1/(1+1i*2*pi*r*c);h1=abs(hw);f1=fft(y1);f11=fftshift(f1);y2=f11*h1;a2=ifft(y2,length(t));fs=100;n=100;window=boxcar(length(a2));[p2,f2]=periodogram(a2,window,n,fs);figure(1);subplot(3,1,1),plot(t,y2),grid on;xlabel('t/s');ylabel('u/v')title('复合信号通过RC积分电路波形')subplot(3,1,2),plot(t,abs(y2)),grid on;xlabel('t/s');ylabel('s_fft');title('复合信号通过RC积分电路幅度谱');subplot(3,1,3),plot(f2,10*log10(p2));xlabel('f/Hz');ylabel('Ga');title('复合信号通过RC积分电路功率谱');（3）分析复合信号通过理想低通系统后的功率谱密度和相应的幅度分布特性；MATLAB源代码：clccleart=0:0.001:0.1;s=3*sin(1e3*t);%正弦信号（幅度分布）y1=awgn(s,10,'measured');f=0:length(t)-1/200:10;f1=fft(y1);f11=fftshift(f1);g=(heaviside(f+20)-heaviside(f-20));y3=f11*g;a3=ifft(y3,length(t));n=100;fs=100;window=boxcar(length(y3));[p3,f3]=periodogram(a3,window,n,fs);figure(1);subplot(3,1,1),plot(t,abs(y3)),grid on;xlabel('t/s');ylabel('u/v');title('复合信号理想低通系统波形');subplot(3,1,2),plot(t,abs(f11));title('复合信号理想低通系统幅度谱');xlabel('t/s');ylabel('a3_fft');subplot(3,1,3),plot(f3,10*log10(p3)),grid on;title('复合信号通过理想低通系统功率谱');第三章3.11 上机题：利用Matlab程序设计一正弦型信号、高斯白噪声信号。

实验六、FIR数字滤波器设计及其网络结构班级: 学号: 姓名: 成绩:1实验目得(1)熟悉线性相位FIR数字滤波器得时域特点、频域特点与零极点分布;(2)掌握线性相位FIR数字滤波器得窗函数设计法与频率采样设计法;(3)了解IIR数字滤波器与FIR数字滤波器得优缺点及其适用场合。

2 实验内容(1)设计计算机程序,根据滤波器得主要技术指标设计线性相位FIR数字低通、高通、带通与带阻滤波器;(2)绘制滤波器得幅频特性与相频特性曲线,验证滤波器得设计结果就是否达到设计指标要求;(3)画出线性相位FIR数字滤波器得网络结构信号流图。

3实验步骤(1)设计相应得四种滤波器得MATLAB程序;(2)画出幅频相频特性曲线;(3)画出信号流图。

4 程序设计%% FIR低通f=[0、2,0、35];m=[1,0];Rp=1;Rs=40;dat1=(10^(Rp/20)-1)/(10^(Rp/20)+1);dat2=10^(-Rs/20);rip=[dat1,dat2];[M,f0,m0,w]=remezord(f,m,rip);M=M+2;hn=remez(M,f0,m0,w);w=0:0、001:pixn=[0:length(hn)-1];H=hn*exp(-j*xn'*w);figuresubplot(2,1,1)plot(w/pi,20*log10(abs(H)));gridon;xlabel('\omega/\pi'),ylabel('|H(e^j^w)|/dB')subplot(2,1,2)plot(w/pi,angle(H)/pi);xlabel('\omega/\pi'),ylabel('\phi(\omega)/\pi') %% FIR高通f=[0、7,0、9];m=[0,1];Rp=1;Rs=60;dat1=(10^(Rp/20)-1)/(10^(Rp/20)+1);dat2=10^(-Rs/20);rip=[dat2,dat1];[M,f0,m0,w]=remezord(f,m,rip);hn=remez(M,f0,m0,w);w=0:0、001:pixn=[0:length(hn)-1];H=hn*exp(-j*xn'*w);figuresubplot(2,1,1)plot(w/pi,20*log10(abs(H)));gridon;xlabel('\omega/\pi'),ylabel('|H(e^j^w)|/dB')subplot(2,1,2)plot(w/pi,angle(H)/pi);xlabel('\omega/\pi'),ylabel('\phi(\omega)/\pi') %% FIR带通f=[0、2,0、35,0、65,0、8];m=[0,1,0];Rp=1;Rs=60;dat1=(10^(Rp/20)-1)/(10^(Rp/20)+1);dat2=10^(-Rs/20);rip=[dat2,dat1,dat2];[M,f0,m0,w]=remezord(f,m,rip);M=M+3hn=remez(M,f0,m0,w);w=0:0、001:pixn=[0:length(hn)-1];H=hn*exp(-j*xn'*w);figuresubplot(2,1,1)plot(w/pi,20*log10(abs(H)));gridon;xlabel('\omega/\pi'),ylabel('|H(e^j^w)|/dB')subplot(2,1,2)plot(w/pi,angle(H)/pi);xlabel('\omega/\pi'),ylabel('\phi(\omega)/\pi') %% FIR带阻f=[0、2,0、35,0、65,0、8];m=[1,0,1];Rp=1;Rs=60;dat1=(10^(Rp/20)-1)/(10^(Rp/20)+1);dat2=10^(-Rs/20);rip=[dat1,dat2,dat1];[M,f0,m0,w]=remezord(f,m,rip);hn=remez(M,f0,m0,w);w=0:0、001:pixn=[0:length(hn)-1];H=hn*exp(-j*xn'*w);figuresubplot(2,1,1)plot(w/pi,20*log10(abs(H)));gridon;xlabel('\omega/\pi'),ylabel('|H(e^j^w)|/dB')subplot(2,1,2)plot(w/pi,angle(H)/pi);xlabel('\omega/\pi'),ylabel('\phi(\omega)/\pi') 5实验结果及分析(1)FIR低通滤波器自动得到得M值不满足要求,故我们将M加上2 在w=0、2π时,H=-0、5dB;w=0、35π时,H=-41dB。

2ASK调制器的设计与实现一、实验目的（1）掌握2FSK 调制的原理及实现方法；（2）学习与熟悉Quratus II 软件的使用；（3）掌握如何应用仿真器来开发应用系统及仿真调试的过程，了解FPGA 开发的基本流程。

二、实验仪器或软件QuartusII 9.0、FPGA实验板、WD990 微机电源、双踪示波器。

三、实验原理3.1 2ASK调制原理2ASK 是数字调制技术的基础，是一种实用的二进制振幅监控方式。

2ASK调制解调器系统框图如图3.1-1所示：图3.1-1 2ASK调制解调器系统3.1.1 2ASK调制器原理及波形图在幅度键控中载波幅度是随着调制信号而变化的。

最简单的形式是载波在二进制调制信号1或0的控制之下通或断。

典型的波形如图3.1-2所示：图3.1-2 2ASK调制原理与波形3.1.2 实验步骤2ASK调制器的调制器电路如图3.1-3所示：图3.1-3 2ASK调制器电路① m序列信号：设计一个周期为15的M序列作为基带信号（信源）（见例1）。

信源码率5000bit/s。

②时钟信号：由实验板提供20MHz时钟clk，仿真时可设20MHz时钟分频。

③ Start信号：开始信号。

在实验板中需要硬件实现。

Start 信号波形图如下:③分频器：编写一个通用的奇偶通用分频程序。

分频次数与载波频率有关。

本实验用数字电路完成分频器设计。

④计数器：使用QuartusⅡ中的lpm_counter器件实现。

⑤开关电路：由基带信号来控制它的输出。

使用QuartusⅡ中的lpm_latch器件实现。

当M序列输出为“1”时输出“f”，为“0”时输出“0“。

3.1.3 测试．按系统方框图，模块化设计，在顶层文件中调用各模块，最终完成一个2ASK调制器。

实现此系统可分3步完成：①电路设计或程序设计。

②QuartusⅡ软件仿真。

③在FPGA实验板中下载并用示波器观察2ASK调制信号。

四、实验结果及分析第一部分电路图图3.1.1 分频器分频器的设计：采用Verilog语言编写程序，再转换为逻辑符号，接入总电路中。

随机信号分析实验报告目录随机信号分析 (1)实验报告 (1)理想白噪声和带限白噪声的产生与测试 (2)一、摘要 (2)二、实验的背景与目的 (2)背景： (2)实验目的: (2)三、实验原理 (3)四、实验的设计与结果 (4)实验设计: (4)实验结果: (5)五、实验结论 (12)六、参考文献 (13)七、附件 (13)1理想白噪声和带限白噪声的产生与测试一、摘要本文通过利用MATLAB软件仿真来对理想白噪声和带限白噪声进行研究。

理想白噪声通过低通滤波器和带通滤波器分别得到低通带限白噪声和帯通带限白噪声。

在仿真的过程中我们利用MATLAB工具箱中自带的一些函数来对理想白噪声和带限白噪声的均值、均方值、方差、功率谱密度、自相关函数、频谱以及概率密度进行研究，对对它们进行比较分析并讨论其物理意义。

关键词：理想白噪声带限白噪声均值均方值方差功率谱密度自相关函数、频谱以及概率密度二、实验的背景与目的背景：在词典中噪声有两种定义：定义1：干扰人们休息、学习和工作的声音，引起人的心理和生理变化。

定义2：不同频率、不同强度无规则地组合在一起的声音。

如电噪声、机械噪声，可引伸为任何不希望有的干扰。

第一种定义是人们在日常生活中可以感知的，从感性上很容易理解。

而第二种定义则相对抽象一些，大部分应用于机械工程当中。

在这一学期的好几门课程中我们都从不同的方面接触到噪声，如何的利用噪声，把噪声的危害减到最小是一个很热门的话题。

为了加深对噪声的认识与了解，为后面的学习与工作做准备，我们对噪声进行了一些研究与测试。

实验目的:了解理想白噪声和带限白噪声的基本概念并能够区分它们，掌握用MATLAB 或c/c++软件仿真和分析理想白噪声和带限白噪声的方法，掌握理想白噪声和带限白噪声的性质。

三、实验原理所谓白噪声是指它的概率统计特性服从某种分布而它的功率谱密度又是均匀的。

确切的说，白噪声只是一种理想化的模型，因为实际的噪声功率谱密度不可能具有无限宽的带宽，否则它的平均功率将是无限大，是物理上不可实现的。