窄带随机信号性能分析

窄带随机信号性能分析

一．摘要

窄带信号在通信系统中有着重要的意义，信号处理技术及通信网络系统与计算机分析技术的相互融合，都要求我们对研究分析随机信号经过系统的响应有一个深入的了解。本实验包括四部分：窄带信号及包络和相位检波分析，窄带随机信号的仿真与分析，希尔伯特变换在单边带系统中的应用，随机信号的DSB分析。主要涉及窄带滤波器的设计，高斯窄带信号包络的均值，均方值和方差的测定，相位概率密度函数的测定等。通过实验了解窄带信号在信号处理领域的应用。

复杂的实际通信系统可以通过抽象与仿真来研究它的特性。本实验通过MATLAB中的仿真出理想信号，并对其进行分析与测量。

二．实验特点与原理

1.窄带信号及包络和相位检波分析

一般无线电接收机中，通常都有高频或中频放大器，它们的通频带往往远小于中心频率f0，既有

?f??1 f0这种线性系统通称为窄带线性系统。

在通信、雷达等许多电子系统中，都常常用一个宽带平稳随机过程来激励一个窄带滤波器，这是在滤波器输出端得到的便是一个

窄带随机过程。若用示波器观测此波形，则可看到，它接近一个正弦波，但此正弦波的幅度和相位都在缓慢的随机变化。我们可以证明，任何一个实窄带随机过程X(t)都可以表示为： X(t)?A(t)cos(?0t??(t))

式中，?0 是固定值，对于窄带随机过程来说，?0一般取窄带滤波器的中心频率或载波频率。

在实际应用中，常常需要检测出包络A(t)和?(t)的信息。若将窄带随机过程X(t)送入包络检波器，则在检波器的输出端可得到包络

A(t)；若将窄带随机过程X(t)送入一个相位检波器，便可检测出相位

信息?(t)。如图10所示：

包络检波器 x(t) w(t) A(t)窄带滤波限幅器× 2cosωt低通滤波器 φ(t)

图10 窄带信号及包络和相位检波器 图10中，在相位检波器之前加入一个理想限幅器，其作用是消除包络起伏对相位检波器的影响。 2.窄带随机信号的仿真与分析 （1）.窄带随机过程

一个实平稳随机过程X(t)，若它的功率谱密度具有下述性质：

中心频率为ωc，带宽为△ω=2ω0，当△ω<<ωc 时，就可认为满足窄带条件。若随机过程的功率谱满足该条件则称为窄带随机过程。若带通滤波器的传输函数满足该条件则称为窄带滤波器。随机过程通过窄带滤波器传输之后变成窄带随机过程。

图1 为典型窄带随机过程的功率谱密度图。若用一示波器来观测次波形，则可看到，它接近于一个正弦波，但此正弦波的幅度和相位都在缓慢地随机变化，图2所示为窄带随机过程的一个样本函数。

图1 典型窄带随机过程的功率谱密度图

图2 窄带随机过程的一个样本函数

(2).窄带随机过程的数学表示 1).用包络和相位的变化表示

由窄带条件可知,窄带过程是功率谱限制在ωc附近的很窄范围内的一个随机过程,从示波器观察(或由理论上可以推知):这个过程中的一个样本函数(一个实现)的波形是一个频率为?c且幅度和相位都做缓慢变化的余弦波。

写成包络函数和随机相位函数的形式：

X(t)=A(t)*cos[ωct+ Φ(t)]

其中:A(t)称作X(t)的包络函数; Φ(t)称作X(t)的随机相位函数。包络随时间做缓慢变化，看起来比较直观，相位的变化，则看不出来。 2).莱斯（Rice）表示式

任何一个实平稳随机过程X(t)都可以表示为：

X(t)=Ac(t) cosωct-AS(t) sinωct 其中同相分量：

Ac(t)= X(t) cosφt= X(t) cosωct＋ sinωct=LP[X(t) *2cosωct] 正交分量：

AS(t) = X(t)sinφt= cosωct— X(t) sinωct= LP[-X(t) *2sinωct] （LP[A]表示取A的低频部分）。Ac(t)和AS(t)都是实随机过程，均值为0，方差等于X(t)的方差。 3. 希尔伯特变换子单边带系统中的应用

在单边带幅度调制中，可以保留上边带，也可以保留下边带。信号单边带调制可以提高信道的利用率。信号单边调制（SSB）有上边带（USB）和下边带（LSB）两种，一般利用希尔伯特变换来实现。 1).利用希尔伯特实现单边带调制的原理框图如下所示：

Xt(t)

图1 利用希尔伯特变换实现单边带调制框图

其中输入信号x(t)： x(t)=s(t)+n(t)。s(t)为频率为1KHz、幅值为1v的正弦波信号。载波为4 KHz、幅值为1v的正弦波信号。n(t)为高斯噪声。

2).希尔伯特变换器的时域特性h(t)为

h(t)?1 ?t低通滤 波器 a点 希尔伯特变换 单边输出

对上式进行傅里叶变化，可得希尔伯特变换器的频率特性H(jw)为：

由以上可知，希尔伯特变换器的幅度响应为︱H(jw)︱=1,相 位响应为φ（w）=-

?2sgn(w)，因此，希尔伯特变换器是一个全通

系统，称为90度相移器。

3).如果调制信号的频谱为X(jw)，则对ｙUSB(t)及ｙLSB(t)的时域表达式两边进行傅里叶变换可得下式： YUSB(jw)=

1111??（jw）??[?(w?wc)??(w?wc)]???h(jw)?j?[?(w?wc)??(w?wc)]= 2?22?2

1j[?(j(w?wc))??(j(w?wc))]?[?h(j(w?w))??h(j(w?w))] 44利用上面Xh(jw)与X(jw)的关系，将Xh(jw)用X(jw)替换得：

YUSB(j?)?sgn(w?wc)]11??（j(w?wc)）??(j(w?wc))??[?(j(w?wc))sgn(w?w)??(j(w?w))*44

YLSB(j?)?1111??（jw）??[?(w?wc)??(w?wc)]???h(jw)?= 2?22?2j?[?(w?wc)??(w?wc)]

1j[?(j(w?wc))??(j(w?wc))]?[?h(j(w?w))??h(j(w?w))]= 44

11??（j(w?wc)）??(j(w?wc))??[?(j(w?wc))sgn(w?w)??(j(w?w))sgn(w?wc)]44

4.随机信号DSB分析 1）.DSB信号的模型

在AM信号中，载波分量并不携带信息，信息完全由边带传送。如果将载波抑制，只需在将直流A0去掉，即可输出抑制载波双边带信号，简称双边带信号（DSB）。 DSB调制器模型如图1所示。

图1 DSB调制器模型

其中，设正弦载波为

c(t)?Acos(?ct??0)

式中，A为载波幅度；?c为载波角频率；?0为初始相位（假定?0为0）。

调制过程是一个频谱搬移的过程，它是将低频信号的频谱搬移到载频位置。而解调是将位于载频的信号频谱再搬回来，并且不失真地恢复出原始基带信号。

双边带解调通常采用相干解调的方式，它使用一个同步解调器，即由相乘器和低通滤波器组成。在解调过程中，输入信号和噪声可以分别单独解调。相干解调的原理框图如图2所示：

图2 相干解调器的数学模型

信号传输信道为高斯白噪声信道，其功率为?2。 2).DSB信号调制过程分析

假定调制信号m(t)的平均值为0，与载波相乘，即可形成DSB信号，其时域表达式为

sDSB?m(t)cos?ct

式中，m(t)的平均值为0。DSB的频谱为

1sDSB(?)?[M(???c)?M(???c)]

2DSB信号的包络不再与调制信号的变化规律一致，因而不能采用简单的包络检波来恢复调制信号， 需采用相干解调(同步检波)。另外，在调制信号m(t)的过零点处，高频载波相位有180°的突变。

除了不再含有载频分量离散谱外，DSB信号的频谱与AM信号的频谱完全相同，仍由上下对称的两个边带组成。所以DSB信号的带宽与AM信号的带宽相同，也为基带信号带宽的两倍， 即

BDSB?BAM?2fH

式中，fH为调制信号的最高频率。 3).高斯白噪声信道特性分析

在实际信号传输过程中，通信系统不可避免的会遇到噪声，例如自然界中的各种电磁波噪声和设备本身产生的热噪声、散粒噪声等，

它们很难被预测。而且大部分噪声为随机的高斯白噪声，所以在设计时引入噪声，才能够真正模拟实际中信号传输所遇到的问题，进而思考怎样才能在接受端更好地恢复基带信号。信道加性噪声主要取决于起伏噪声，而起伏噪声又可视为高斯白噪声，因此我在此环节将对双边带信号添加高斯白噪声来观察噪声对解调的影响情况。

为了具体而全面地了解噪声的影响问题，我将分别引入大噪声（信噪比为20dB）与小噪声（信噪比为2dB）作用于双边带信号，再分别对它们进行解调，观察解调后的信号受到了怎样的影响。

在此过程中，我用函数randn来添加噪声，此函数功能为向信号中添加噪声功率为其方差的高斯白噪声。

正弦波通过加性高斯白噪声信道后的信号为

r(t)?Acos(?ct??)?n(t)

故其有用信号功率为

A2S?

2噪声功率为

N??2

信噪比N满足公式

B?10log10(SN)

S则可得到公式

??2A22?10B10

我们可以通过这个公式方便的设置高斯白噪声的方差。

4).DSB解调过程分析

所谓相干解调是为了从接收的已调信号中，不失真地恢复原调制信号，要求本地载波和接收信号的载波保证同频同相。相干解调的一般数学模型如图所示。

图5 DSB相干解调模型

设图四的输入为DSB信号

Sm(t)?SDSB(t)?m(t)cos(?ct??0)

乘法器输出为

?(t)?SDSB(t)?m(t)cos(?ct??0)cos(?ct??)?1m(t)[cos(???0)?cos(2?ct??0??)]2

通过低通滤波器后

m0(t)?1m(t)cos(?0??) 2当?0???常数时，解调输出信号为

m0(t)?1m(t) 2三．实验的设计思想与实现

1.窄带信号及包络和相位检波分析 程序实现

clear all

Wp1=2*pi*1*10^3;Ws1=4*pi*1*10^3;rp=1;rs=30;

wp1=1;ws1=Ws1/Wp1;

[N1,wc1]=buttord(wp1,ws1,rp,rs,'s'); [z1,p1,k1]=buttap(N1); [B1,A1]=zp2tf(z1,p1,k1); w1=0:pi/512:2*pi;

[h1,w1]=freqs(B1,A1,w1);

plot(w1,20*log10(abs(h1)),'k');grid; xlabel('\\lambda');

ylabel('A(\\lambda)/dB'); title('低通滤波器幅度谱') grid on;

plot(w1,angle(h1)); %在 0~2π 区间内作相位谱 title('低通滤波器相位谱'); grid on;

fs=20000; %设定采样频率 N=100; %取的样本点数 n=0:N-1; %取样点序列

t=n/fs; %获得以1/16000为时间间隔采样序列

noisy=wgn(1,N,0); %产生高斯白噪声 a=cos(2*pi*1000*t); %获取A(t)的采样点

x=(1+a).*cos(2*pi*4000*t+2)+noisy/10; %获取X(t)的采样点 figure(1); %画图 subplot(2,1,1);

plot(n,x);%以t为横坐标画出x(t)的时域图型 axis([0 100 -2.5 2.5]); xlabel('采样点'); ylabel('X(t)/V');

title('图1-窄带随机信号波形');

%窄带系统检测

xc=filter(B1,A1,x); %对信号进行滤波 Ac=xc.*2.*cos(2*pi*4000*t); As=xc.*2.*sin(2*pi*4000*t);

y=Ac.*cos(2*pi*4000*t)-As.*sin(2*pi*4000*t); %作输出信号时域波形 figure(10)

subplot(2,1,1);

n=0:N-1; plot(n,y); %作信号时域波形 axis([0 100 -3 3]);

%axis tight;

xlabel('采样点' ); ylabel('幅值') title('Y时域谱'); grid on;

%进行 FFT 变换并做频谱图

yy=fft(y,length(y)); %对滤波后信号进行 fft 变换 longy=length(yy); %求傅里叶变换后的序列长度

labelx=(0:longy-1)*16000/longy;

magn3=abs(yy); % 求滤波后信号幅值 subplot(2,1,2);

plot(labelx,magn3); %做频谱图 axis tight;

xlabel('频率(Hz)'); ylabel('幅值' ); title('Y频谱图'); grid on;

%求输出信号的自相关函数

[c3,lags3]=xcorr(y,'coeff'); %求出自相关序列 figure(11);

subplot(2,1,1);

plot(lags3/fs,c3); %在时域内画自相关函数 xlabel('τ'); ylabel('Rx(τ)'); axis tight;

title('Y的自相关函数'); grid on;

%%.作窄带随即过程y(t)的相位图 xiangwei=angle(y); figure(9);

plot(n,xiangwei);%以t为横坐标画出y(t)的时域图型 axis([0 50 0 4]);

title('图9-y(t)的相位'); grid on;

.做相位得一维概率函数 figure(6);

hist(xiangwei);

title('相位的一维概率密度'); grid on;

each1=linspace(min(xiangwei),max(xiangwei),10); %将最大最小区间分成14等份,然后分别计算各个区间的个数

nr1=hist(xiangwei,each1); %计算各个区间的个数

nr1=nr1/length(xiangwei); %计算各个区间的个数归一化 subplot(2,1,2); %%%多项式拟合 bar(each1,nr1); hold on;

plot(each1,nr1,'g') hold on;

p1=polyfit(each1,nr1,5); eachi1=0:0.1:3.5;

nri1=polyval(p1,eachi1); plot(eachi1,nri1,'r') axis tight;

title('相位概率密度分布'); xlabel('xiangwei'); ylabel('P1(X)'); grid on;

%9.窄带随即过程y(t)的包络 A=abs(y); figure(7); plot(n,A);

axis([0 100 0 5]); xlabel('时间t/s'); ylabel('包络A(t)/v'); title('图7-y(t)的包络'); grid on;

%%.作包络的一维概率密度图 figure(8); hist(A);

title('包络的一维概率密度'); grid on;

each=linspace(min(A),max(A),14); %将最大最小区间分成14等份,然后分别计算各个区间的个数

nr=hist(A,each); %计算各个区间的个数

nr=nr/length(A); %计算各个区间的个数归一化 subplot(2,1,2); %%%多项式拟合 bar(each,nr); hold on;

plot(each,nr,'g') hold on;

p=polyfit(each,nr,6); eachi=0:0.1:3.5;

nri=polyval(p,eachi); plot(eachi,nri,'r') axis tight;

title('包络概率密度分布'); xlabel('A');

ylabel('P(X)')

2.窄带随机信号的仿真与分析 1）.程序设计

clc

fs=16000; %设定采样频率 N=1300;

n=0:N-1; %取的样本点数

t=n/fs; %获得以1/16000为时间间隔采样序列 noisy=wgn(1,N,0); %产生高斯白噪声 w1=2*pi*1000;%低频信号角频率 w2=2*pi*4000;%载波角频率

a=cos(2*pi*1000*t); %获取A(t)的采样点

x=(1+a).*cos(2*pi*4000*t+2)+noisy/10; %获取x(t)的采样点

%以t为横坐标画出x(t)的时域图型

figure(1); subplot(2,1,1); plot(n,x);

axis([0 140 -3 3]);xlabel('采样点');ylabel('X(t)/V');title('窄带随机信号波形');grid on;

%求X(t)的统计特性 并画出来

disp('X(t)的均值为'); Ex=mean(x); disp(Ex);%求X(t)均值 hold on; plot(n,Ex,'r.');

disp('X(t)的方差为');Dx=var(x); disp(Dx);%求x(t)方差 hold on; plot(n,Dx,'g.'); %画出X(t)的概率分布函数

each=linspace(min(x),max(x),14); %将最大最小区间分成14等份,然后分别计算各个区间的个数

nr=hist(x,each); %计算各个区间的个数

nr=nr/length(x); %计算各个区间的个数归一化 subplot(2,1,2); p=polyfit(each,nr,9); %画出概率分布直方图 bar(each,nr); %多项式拟合 hold on; plot(each,nr,'g') eachi=-2:0.1:2;

nri=polyval(p,eachi); plot(eachi,nri,'r')

axis tight;title('X(t)概率密度分布');xlabel('X(t)');ylabel('P(x)');grid on;

%对X(t)进行频谱分析

Fx=fft(x,N); %对x(t)进行fft变换，在0~16000区间内得到2N-1个频率值 magn=abs(Fx); %求x(t)幅值 xangle=angle(Fx); %求X(t)相位

labelang=(0:length(x)-1)*16000/length(x); %在0~16000区间内求横坐标刻度 figure(2); plot(labelang,magn*10); %在0~16000区间内做频谱和相位图 axis([0 16000 -0.5 600]); xlabel('频率/Hz');ylabel('幅值');title('X(t)频谱图');grid on;

%求X(t)的自相关函数

[c,lags]=xcorr(x,'coeff'); %求出自相关序列

figure(3); subplot(2,1,1); plot(lags/fs,c); %在时域内画自相关函数 axis tight; xlabel('T');ylabel('Rx(T)');title('X(t)的自相关函数');grid on;

%求X(t)的功率谱密度 long=length(c);

Sx=fft(c,long); labelx=(0:long-1)*2*pi;

plot_magn=10*log10(abs(Sx));

subplot(2,1,2); plot(labelx,plot_magn); %画功率谱密度 axis tight;xlabel('w');ylabel('Sx(w)');title('X(t)的功率谱密度');grid on; %窄带系统检测

z1=2.*cos(2*pi*4000*t); z2=-2.*sin(2*pi*4000*t);

Ac=z1.*x; %滤波后生成Ac(t) As=z2.*x; %滤波后生成As(t)

y=Ac.*cos(2*pi*4000*t)-As.*sin(2*pi*4000*t); %滤波器设计

f_p=1000;f_s=1600;R_p=1;R_s=35; %设定滤波器参数; 通、阻带截止频率,通、阻带衰减

Ws=2*f_s/fs;Wp=2*f_p/fs; %频率归一化

[n,Wn]=buttord(Wp,Ws,R_p,R_s); %采用巴特沃思滤波器 [b,a]=butter(n,Wn); %求得滤波器传输函数的多项式系数 figure(4);

[H,W]=freqz(b,a); %求得滤波器传输函数的幅频特性

subplot(2,1,1); plot(W*fs/(2*pi),abs(H)); %在0～2pi区间内作幅度谱 title('低通滤波器幅度谱'); grid on;

subplot(2,1,2); plot(W*fs/(2*pi),angle(H)); %在0～2pi区间内作相位谱 title('低通滤波器相位谱'); grid on;

%求Ac(t)滤波后的统计特性

mc=filter(b,a,Ac); %上支路通过滤波器 Ac(t) disp('Ac(t)的均值');Eh=mean(mc) %求Ac(t)的均值

disp('Ac(t)的均方值是');E2h=mc*mc'/N %求Ac(t)的均方值 disp('Ac(t)的方差');Dh=var(mc) %求Ac(t)的方差 %画Ac(t)的时域波形

figure(6); subplot(2,1,1); n=0:N-1; plot(n,mc);

axis([0 300 -1 1]);xlabel('采样点');ylabel('幅值');title('Ac(t)的时域波形');grid on; %画Ac(t)的频谱图

yc=fft(mc,length(mc)); %对Ac(t)进行fft变换 longc=length(yc); %求傅里叶变换后的序列长度 labelx=(0:longc-1)*16000/longc;

magnl=abs(yc); %求Ac(t)的幅值

subplot(2,1,2); plot(labelx,magnl); %画Ac(t)的频谱图

axis tight; xlabel('频率(Hz)'); ylabel('幅值'); title('Ac(t)频谱图'); grid on; %求Ac(t)的自相关函数

[c1,lags1]=xcorr(mc,'coeff'); %求出Ac(t)的自相关序列

figure(7); subplot(2,1,1); plot(lags1/fs,c1); %在时域内画Ac(t)的自相关函数 xlabel('T');ylabel('Rx(T)');axis tight; title('Ac(t)的自相关函数'); grid on;

%求Ac(t)的双边功率谱

Sac=fft(c1,length(c1)); %对Ac(t)的自相关函数进行傅里叶变换 magnc=abs(Sac); %求Ac(t)的双边功率谱幅值 long=length(Sac); %求傅里叶变换后的序列长度 labelc=(0:long-1)*16000/long;

subplot(2,1,2); plot(labelc,10*log10(magnc)); %画Ac(t)的自相关函数频谱 即为Ac(t)的双边功率谱

xlabel('频率(Hz)');ylabel('功率谱(dbW)');axis tight;title('Ac(t)的双边功率谱');grid on;

%求得As(t)的统计特性

ms=filter(b,a,As); %对下支路信号进行滤波得As(t) disp('As(t)的均值'); Eh=mean(ms) %求As(t)的均值

disp('As(t)的均方值是'); E2h=ms*ms'/N %求As(t)的均方值 disp('As(t)的方差'); Dh=var(ms) %求As(t)的方差 %作As(t)的时域波形

figure(8);subplot(2,1,1); n=0:N-1;plot(n,ms); %画出As(t)的时域波形

axis([0 300 -0.5 2]); xlabel('采样点');ylabel('幅值');title('As(t)的时域波形');grid on; %对As(t)进行FFT变换并做频谱图

ys=fft(ms,length(ms)); %对As(t)进行fft变换 longs=length(ys); %求傅里叶变换后的序列长度 labelx=(0:longs-1)*16000/longs;

magn2=abs(ys); %求As(t)的幅值

subplot(2,1,2); plot(labelx,magn2); %画出As(t)的频谱图 axis tight;xlabel('频率(Hz)');ylabel('幅值');title('As(t)的频谱图');grid on; %求As(t)的自相关函数

[c2,lags2]=xcorr(ms,'coeff'); %求出As(t)的自相关序列

figure(9);subplot(2,1,1);plot(lags2/fs,c2); %画出As(t)自相关函数的时域波形 xlabel('T');ylabel('Rx(T)');axis tight;title('As(t)的的自相关函数');grid on; %求As(t)的双边功率谱

Sas=fft(c2,length(c2)); %对As(t)的自相关函数进行傅里叶变换 magnc=abs(Sac); %求As(t)的双边功率谱幅值 long=length(Sas); %求傅里叶变换后的序列长度 labels=(0:long-1)*16000/long;

subplot(2,1,2); plot(labelc,10*log10(magnc)); %画As(t)的自相关函数频谱 xlabel('频率(Hz)');ylabel('功率谱(dbW)');axis tight;title('As(t)的双边功率谱');

% 求y(t)的统计特性

disp('输出信号Y(t)的均值');Eh=mean(y) %求输出信号Y(t)的均值 disp('输出信号Y(t)的均方值');E2h=y*y'/N %求输出信号Y(t)的均方值 disp('输出信号Y(t)的方差');Dh=var(y) %求输出信号Y(t)的方差 %作输出信号Y(t)的时域波形

figure(10); subplot(2,1,1);n=0:N-1;plot(n,y);

axis([0 150 -2 2]);xlabel('采样点');ylabel('幅值');title('Y(t)的时域波形');grid on;

6）输出下边带信号的统计特性

均值mean= -0.0023

均方差meansquare = 0.5066 方差var1 = 0.5066

自相关函数图象如下：

7）输出上边带信号的统计特性

均值mean= -8.2231e-004 均方差meansquare = 0.5065 方差var1 = 0.5065

自相关函数图象和功率谱密度图象如下：

4.随机信号的DSB分析

DSB信号调制过程：

clf; %清除窗口中的图形 clc; clear;

ts=0.01;%定义变量区间步长 t0=2;%定义变量区间终止值

t=-t0:ts:t0;%定义变量区间取值情况 fc=10;%给出相干载波的频率 A=1;%定义调制信号幅度 fa=1;%定义调制信号频率

mt=A*cos(2*pi*fa.*t);%输入调制信号表达式 ct=cos(2*pi*fc.*t);%输入载波信号表达式 psnt=mt.*ct;%输出调制信号表达式 subplot(5,1,1);%划分画图区间

plot(t,mt,'g');%画出调制信号波形 title('输入信号波形'); xlabel('Variable t'); ylabel('Variable mt'); subplot(5,1,2);

plot(t,ct,'b');%画出载波信号波形 title('输入载波波形'); xlabel('Variable t'); ylabel('Variable ct'); subplot(5,1,3);

plot(1:length(psnt),psnt,'r');%length用于长度匹配 title('已调信号波形');%画出已调信号波形 xlabel('Variable t'); ylabel('Variable psnt'); mt1=fftshift(fft(mt));

mt2=abs(mt1.^2); df=1/(2*t0); ff=length(mt1);

f=-ff/2*df:df:ff/2*df-df; subplot(5,1,4);

plot(f,mt2);%调制信号功率谱密度 psnt1=fftshift(fft(psnt)); psnt2=abs(psnt1.^2); subplot(5,1,5);

plot(f,psnt2);%已调信号功率谱密度 高斯白噪声信道特性： clf;%清除窗口中的图形 ts=0.01;%定义变量区间步长 t0=2;%定义变量区间终止值

t=-t0+0.0001:ts:t0;%定义变量区间 fc=10;%给出相干载波的频率 A=1;%定义输入信号幅度 fa=1;%定义调制信号频率

mt=A*cos(2*pi*fa.*t);%输入调制信号表达式 xzb=2;%输入小信噪比(dB) snr=10.^(xzb/10);

[h,l]=size(mt);%求调制信号的维数 fangcha=A*A./(2*snr);%由信噪比求方差

nit=sqrt(fangcha).*randn(h,l);%产生小信噪比高斯白躁声 psmt=mt.*cos(2*pi*fc.*t);%输出调制信号表达式 psnt=psmt+nit;%输出叠加小信噪比已调信号波形 xzb=20;%输入大信噪比(dB) snr1=10.^(xzb/10);

[h,l]=size(mt);%求调制信号的维数

fangcha1=A*A./(2*snr1);%由信噪比求方差

nit1=sqrt(fangcha1).*randn(h,l);%产生大信噪比高斯白噪声 psnt1=psmt+nit1;%输出已调信号波形 subplot(2,2,1);%划分画图区间

plot(t,nit,'g');%画出输入信号波形 title('小信噪比高斯白躁声'); xlabel('Variable t'); ylabel('Variable nit'); subplot(2,2,2); plot(t,psnt,'b');

title('叠加小信噪比已调信号波形'); xlabel('Variable t'); ylabel('Variable psnt'); subplot(2,2,3);

plot(t,nit1,'r');%length用于长度匹配

title('大信噪比高斯白躁声');%画出输入信号与噪声叠加波形 xlabel('Variable t'); ylabel('Variable nit'); subplot(2,2,4); plot(t,psnt1,'k');

title('叠加大信噪比已调信号波形');%画出输出信号波形 xlabel('Variable t'); ylabel('Variable psmt');

调制解调仿真过程：

clf;%清除窗口中的图形 ts=0.01;%定义变量区间步长 t0=2;%定义变量区间终止值

t=-t0+0.0001:ts:t0;%定义变量区间 fc=10; %给出相干载波的频率 A=1; %定义输入信号幅度 fa=1; %定义调制信号频率

mt=A*cos(2*pi*fa.*t); %输入调制信号表达式 xzb=20;%输入信噪比(dB) snr=10.^(xzb/10);

[h,l]=size(mt);%求调制信号的维数 fangcha=A*A./(2*snr);%由信噪比求方差

nit=sqrt(fangcha).*randn(h,l); %产生高斯白噪声 snit=mt+nit; %调制信号与噪声叠加

psmt=mt.*cos(2*pi*fc.*t); %输出调制信号表达式 pnit=nit.*cos(2*pi*fc.*t); %输出噪声表达式 psnt=psmt+pnit; %输出已调信号波形

jic=psnt.*cos(2*pi*fc.*t); %调制信号乘以相干载波

ht=(2*pi*fc.*sin(2*pi*fc.*t)./(2*pi*fc.*t))./pi; %低通滤波器的时域表达式

htw=abs(fft(ht)); %低通滤波器的频域表达式 jt=conv(ht,jic); %解调信号的时域表达式 subplot(3,3,1); %划分画图区间

plot(t,mt,'g'); %画出输入信号波形 title('输入信号波形'); xlabel('Variable t'); ylabel('Variable mt'); subplot(3,3,2); plot(t,nit,'b');

title('输入噪声波形'); xlabel('Variable t'); ylabel('Variable nit'); subplot(3,3,3);

plot(1:length(snit),snit,'r'); %length用于长度匹配

title('输入信号与噪声叠加波形'); %画出输入信号与噪声叠加波形 xlabel('Variable t'); ylabel('Variable snit'); subplot(3,3,4); plot(t,psmt,'k');

title('输出信号波形'); %画出输出信号波形 xlabel('Variable t'); ylabel('Variable psmt'); subplot(3,3,5); plot(t,pnit,'k');

title('输出噪声波形'); %画出输出噪声波形 xlabel('Variable t'); ylabel('Variable pnit'); subplot(3,3,6); plot(t,psnt,'k');

title('输出信号与输出噪声叠加波形'); %画出输出信号与输出噪声叠加波形 xlabel('Variable t'); ylabel('Variable psnt'); subplot(3,3,7);

plot(1:length(htw),htw,'k');

title('低通滤波器频域波形'); %画出低通滤波器频域波形 xlabel('Variable w'); ylabel('Variable htw'); axis([0 60 0 150]); subplot(3,3,8);

plot(1:length(ht),ht,'k');

title('低通滤波器时域波形');%画出低通滤波器时域波形 xlabel('Variable t'); ylabel('Variable psnt');

axis([150 250 -20 25]); %给出坐标轴范围 subplot(3,3,9);

plot(1:length(jt),jt,'k');

title('输出信号与输出噪声叠加波形');%画出输出信号与输出噪声叠加波形 xlabel('Variable t'); ylabel('Variable jt'); axis([200 600 -50 50]);

仿真结果及分析 DSB信号调制过程

本文来源：https://www.bwwdw.com/article/fudd.html