%SSB信号调制解调
clear;clc;
f0 = 1; %信源信号频率(Hz)
E0 = 1; %信源信号振幅(V)
E = 1; %载波分量振幅(V)
fc = 10; %载波分量频率（Hz）
t0 = 1; %信号时长
snr = 15; %解调器输入信噪比dB
dt = 0.003; %系统时域采样间隔
fs = 1/dt; %系统采样频率
df = 0.001; %所需的频率分辨率
t = 0:dt:t0;
Lt = length(t); %仿真过程中，信号长度
snr_lin = 10^(snr/10);%解调器输入信噪比
%-------------画出调制信号波形及频谱
%产生模拟调制信号
m = E*cos(2*pi*f0*t);
L = min(abs(m));%包络最低点
R = max(abs(m));%包络最高点
%画出调制信号波形和频谱
clf;
figure(1);
%%
%画出调制信号波形
subplot(411);
plot(t,m(1:length(t)));
axis([0,t0,-R-0.3,R+0.3]);%设置坐标范围
xlabel('t');title('调制信号');
set(gca,'YTick',-R:1:R);
subplot(412);
[M,m,df1,f] = T2F_new(m,dt,df,fs); %求出调制信号频谱
[Bw_eq] = signalband(M,df,t0); %求出信号等效带宽
f_start_low = fc - Bw_eq; %求出产生下边带信号的带通滤波器的起始频率f_cutoff_low = fc; %求出产生下边带信号的带通滤波器的截止频率f_start_high = fc; %求出产生上边带信号的带通滤波器的起始频率f_cutoff_high = fc + Bw_eq; %求出产生上边带信号的带通滤波器的截止频率
plot(f,fftshift(abs(M))); %画出调制信号频谱%M：傅里叶变换后的频谱序列
xlabel('f');title('调制信号频谱');
axis([-fc-5*f0,fc+5*f0,0,max(M)+0.3]);
set(gca,'XTick', -10:10:10);
set(gca,'XGrid','on');
%%
%载波及其频谱
subplot(413);
c = cos(2*pi*fc*t); %载波
plot(t,c);
axis([0,t0,-E-0.2,E+0.2]);
xlabel('t');title('载波');
subplot(414); %载波频谱
[C,c,df1,f] = T2F_new(c,dt,df,fs);
plot(f,fftshift(abs(C))); %画出载波频谱
xlabel('f');title('载波频谱');
axis([-fc-5*f0,fc+5*f0,0,max(C)+0.3]);
set(gca,'XTick', -10:10:10);
set(gca,'XGrid','on');
%%
%已调信号及其频谱
figure(2);
subplot(321); %画已调信号
u = m(1:Lt).*c(1:Lt);
plot(t,u);
axis([0,t0,-max(u)-0.5,max(u)+0.5]);
xlabel('t');title('DSB信号');
set(gca,'YTick', -max(u):1:max(u));
subplot(322);
[U,u,df1,f] = T2F_new(u,dt,df,fs);
plot(f,fftshift(abs(U))); %画出已调信号频谱
xlabel('f');title('DSB信号频谱');
axis([-fc-5*f0,fc+5*f0,0,max(U)+0.3]);
set(gca,'XTick', -10:10:10);
set(gca,'XGrid','on');
%%
%滤波法产生SSB信号
[H_low,f_low] = bp_f(length(u),f_start_low,f_cutoff_low,df1,fs,1);%求滤波法产生下边带需要的带通滤波器
[H_high,f_high] = bp_f(length(u),f_start_high,f_cutoff_high,df1,fs,1);
%[H,f] = bp_f(length(sam),f_start,f_cutoff,df1,fs,1);
subplot(323);
plot(f_low,fftshift(abs(H_low))); %画出带通滤波器xlabel('f');title('下边带带通滤波器');
axis([-f_cutoff_low-1,f_cutoff_low+1,-0.05,1.05]);
subplot(324);
plot(f_high,fftshift(abs(H_high))); %画出带通滤波器xlabel('f');title('上边带带通滤波器');
axis([-f_cutoff_high-1,f_cutoff_high+1,-0.05,1.05]);
subplot(325);
plot(f_low,fftshift(abs(H_low)));hold on; %画出带通滤波器
plot(f,fftshift(abs(U))); %画出已调信号频谱axis([-fc-5,fc+5,-0.05,1.05]);
xlabel('f');title('下边带信号');
subplot(326);
plot(f_high,fftshift(abs(H_high)));hold on; %画出带通滤波器
plot(f,fftshift(abs(U))); %画出已调信号频谱axis([-fc-5,fc+5,-0.05,1.05]);
xlabel('f');title('上边带信号');
%%
%----------------经过带通滤波器，产生单边带信号（以上边带信号为例）
samuf = H_high.*U; %滤波器输出信号的频谱
[samu] = F2T_new(samuf,fs); %滤波器输出信号的波形
figure(3);
subplot(321);
plot(t,samu(1:Lt));
axis([0,t0,-max(samu)-0.3,max(samu)+0.3]);
xlabel('t');title('上边带信号');
%%
[samuf,samu,df1,f] = T2F_new(samu(1:Lt),dt,df,fs);%上边带信号频谱
subplot(322);
plot(f,fftshift(abs(samuf))); %画出经过理想带通滤波器后信号频谱xlabel('f');title('上边带信号频谱');
axis([-fc-5*f0,fc+5*f0,0,max(samuf)+0.1]);
set(gca,'XTick', -10:10:10);
set(gca,'XGrid','on');
%%
%将已调信号送入信道
%先根据所给信噪比产生高斯白噪声
signal_power = power_x(samu(1:Lt)); %已调信号的平均功率
noise_power = (signal_power * fs)/(snr_lin*4*Bw_eq); %求出噪声方差（噪声均值为0）noise_std = sqrt(noise_power); %噪声标准差
noise = noise_std * randn(1,Lt); %产生噪声
%画出信道高斯白噪声波形及频谱，此时，噪声已实现，为确知信号，可求其频谱
subplot(323);
plot(t,noise);
axis([0,t0,-max(noise),max(noise)]);
xlabel('t');title('噪声信号');
subplot(324);
[noisef,noise,df1,f] = T2F_new(noise,dt,df,fs); %噪声频谱
plot(f,fftshift(abs(noisef))); %画出噪声频谱
xlabel('f');title('噪声频谱');
%%
%信道中的信号%叠加了噪声的已调信号频谱
sam = samu(1:Lt) + noise(1:Lt);
subplot(325);
plot(t,sam);
axis([0,t0,-max(sam),max(sam)]);
xlabel('t');title('信道中的信号');
subplot(326);
[samf,sam,df1,f] = T2F_new(sam,dt,df,fs); %求出叠加了噪声的已调信号频谱plot(f,fftshift(abs(samf))); %画出叠加了噪声的已调信号频谱xlabel('f');title('信道中信号的频谱');
axis([-fc-5*f0,fc+5*f0,0,max(samf)+0.1]);
set(gca,'XTick', -10:10:10);
set(gca,'XGrid','on');
%%
%----------------经过带通滤波器
%经过理想滤波器后的信号及其频谱
DEM = H_high.*samuf; %滤波器输出信号的频谱
[dem] = F2T_new(DEM,fs);%滤波器输出信号的波形
figure(4);
subplot(321); %经过理想带通滤波器后的信号波形
plot(t,dem(1:Lt)); %画出经过理想带通滤波器后的信号波形
axis([0,t0,-max(dem)-0.3,max(dem)+0.3]);
xlabel('t');title('理想BPF输出信号');
%%
[demf,dem,df1,f] = T2F_new(dem(1:Lt),dt,df,fs);%求经过理想带通滤波器后的信号频谱subplot(322);
plot(f,fftshift(abs(demf))); %画出经过理想带通滤波器后信号频谱xlabel('f');title('理想BPF输出信号频谱');
axis([-fc-5*f0,fc+5*f0,0,max(demf)+0.1]);
set(gca,'XTick', [-10:10:10]);
set(gca,'XGrid','on');
%%
%--------------和本地载波相乘，即混频
subplot(323);
plot(t,c(1:Lt));
axis([0,t0,-E-0.2,E+0.2]);
xlabel('t');title('本地载波');
subplot(324); %频谱载波
[C,c,df1,f] = T2F_new(c(1:Lt),dt,df,fs);
plot(f,fftshift(abs(C))); %画出载波频谱
xlabel('f');title('本地载波频谱');
axis([-fc-5*f0,fc+5*f0,0,max(C)+0.3]);
set(gca,'XTick', [-10:10:10]);
set(gca,'XGrid','on');
%再画出混频后信号及其频谱
der = dem(1:Lt).*c(1:Lt); %混频
%%
subplot(325); %画出混频后的信号
plot(t,der);
axis([0,t0,-R,R]);
xlabel('t');title('混频后的信号');
subplot(326);
[derf,der,df1,f] = T2F_new(der,dt,df,fs); %求出混频后的信号频谱
plot(f,fftshift(abs(derf))); %画出混频后的信号频谱
xlabel('f');title('混频后的信号频谱');
axis([-2*fc-5*f0,2*fc+5*f0,0,max(derf)+0.3]);
set(gca,'XTick', [-10:10:10]);
set(gca,'XGrid','on');
%%
%-----------------经过低通滤波器
%画出理想低通滤波器
figure(5);
[LPF,f] = lp_f(length(der),Bw_eq,df1,fs,1); %求出低通滤波器
subplot(411);
plot(f,fftshift(abs(LPF))); %画出理想低通滤波器
xlabel('f');title('理想LPF');
axis([-f0-Bw_eq,f0+Bw_eq,-0.05,1.05]);
%%
%混频信号经过理想低通滤波器后的频谱及波形
DM = LPF.*derf; %理想低通滤波器输出的频谱[dm] = F2T_new(DM,fs); %滤波器的输出波形
subplot(412);
plot(t,dm(1:Lt)); %画出经过低通滤波器后的解调波形axis([0,t0,-max(dm)-0.2,max(dm)+0.2]);
xlabel('t');title('恢复信号');
set(gca,'YTick', [-1:0.5:1]);
set(gca,'YGrid','on');
subplot(413);
[dmf,dm,df1,f] = T2F_new(dm(1:Lt),dt,df,fs); %求LPF输出信号的频谱
plot(f,fftshift(abs(dmf))); %画出LPF输出信号的频谱xlabel('f');title('恢复信号频谱');
axis([-fc,fc,0,max(abs(dmf))+0.1]);
set(gca,'XTick', [-10:10:10]);
set(gca,'XGrid','on');
%%
subplot(414);
plot(t,m(1:Lt)); %画出调制信号波形
xlabel('t');title('调制信号');
set(gca,'YTick', [-R:1:R]);
axis([0,t0,-R-0.3,R+0.3])
xlabel('t');title('调制信号');
子函数
%序列的傅里叶变换
%各参数含义与子函数T2F中的完全相同，完成序列的傅里叶变换
function [M,m,df] = fftseq(m,ts,df)
fs = 1/ts;
if nargin == 2
n1 = 0;
else
n1 = fs / df;
end
n2 = length(m);
n = 2^(max(nextpow2(n1),nextpow2(n2)));
M = fft(m,n);
m = [m,zeros(1,n-n2)];
df = fs / n;
end
%计算信号功率
function p = power_x(x)
%x:输入信号
%p:返回信号的x功率
p = (norm(x).^2)./length(x);
end
%将信号从频域转换到时域
function [m] = F2T(M,fs)
%----------------输入参数
%M：信号的频谱
%fs:系统采样频率
%----------------输出（返回）函数
%m:傅里叶逆变换后的信号，注意其长度为2的整数次幂，利用其画波形时，要注意选取m 的一部分，选取长度和所给时间序列t的长度要一致，plot(t,m(1:length(t))),否则会出错
m = real(ifft(M))*fs;
end
%将信号从时域转换到频域
function [M,m,df1,f] = T2F(m,ts,df,fs)
%----------------输入参数
%m：信号
%ts:系统时域采样间隔、
%df:所需的采样频率
%fs:系统采样频率
%----------------输出（返回）函数
%M：傅里叶变换后的频谱序列
%m:输入信号参与傅里叶变换后对应序列，需要注意的是，该序列与输入信号m的区别，其长度是不一样的，输入的m长度不一定是2的整数次幂，而傅里叶变换要求输入信号长度为2的整数次幂，
%故傅里叶变换前需要对m信号进行补零操作，其长度有所增加，估输出参数中的m为补零后的输入信号，其长度与输入参数m不一样，但与M,f的长度是一样的，
%并且，其与时间序列t所对应的序列m(1:length(t))与输入参数中的m是一致的。
%df1:返回的频率分辨率
%f:与M相对应的频率序列
[M,m,df1] = fftseq(m,ts,df);
f = [0:df1:df1*(length(m)-1)]-fs/2;
M = M / fs;
End
%低通滤波器
function [H,f] = lp_f(n,f_cutoff,df1,fs,p)
%低通滤波器函数输入设计的滤波器参数，产生低通滤波器频率特性函数H和频率向量f %------------输入参数
%n 低通滤波器的输入信号长度
%f_cutoff 低通滤波器截止频率
%df1 频率分辨率
%fs 抽样频率
%p 滤波器振幅
%------------输出（返回）函数
%H 低通滤波器频率响应
%f 频率向量
% 滤波器设计
n_cutoff = floor(f_cutoff/df1);
f = [0:df1:df1*(n-1)]-fs/2;
H = zeros(size(f));
H(1:n_cutoff) = p * ones(1,n_cutoff);
H(length(f)-n_cutoff+1:length(f)) = p * ones(1,n_cutoff);
end
%带通滤波器
function [H,f] = bp_f(n,f_start,f_cutoff,df1,fs,p)
%带通滤波器函数，输入设计的滤波器参数，产生带通滤波器频率特性函数H和频率向量f %------------输入参数
%n 带通滤波器的输入信号长度
%f_start 通带起始频率
%f_cutoff 带通滤波器截止频率
%df1 频率分辨率
%fs 抽样频率
%p 滤波器振幅
%------------输出（返回）函数
%H 带通滤波器频率响应
%f 频率向量
% 滤波器设计
n_cutoff = floor(f_cutoff/df1);
n_start = floor(f_start/df1);
f = [0:df1:df1*(n-1)]-fs/2;
H = zeros(size(f));
H(n_start + 1:n_cutoff) = p * ones(1,n_cutoff - n_start);
H(length(f) - n_cutoff + 1:length(f) - n_start) = p * ones(1,n_cutoff - n_start);
end
function [Bw_eq] = signalband(sf,df,T)
%计算信号等效带宽
%sf:信号频谱
%df:频率分辨率
%T:信号持续时间
sf_max = max(abs(sf));
Bw_eq = sum(abs(sf).^2)*df/T/sf_max.^2;
中南民族大学 软件课程设计报告 电信学院级通信工程专业 题目2PSK数字信号的调制与解调学生学号 42 指导教师 2012年4月21日
基于MATLAB数字信号2PSK的调制与解调 摘要：为了使数字信号在信道中有效地传播，必须使用数字基带信号的调制与解调，以使得信号与信道的特性相匹配。基于matlab实验平台实现对数字信号的2psk的调制与解调的模拟。本文详细的介绍了PSK波形的产生和仿真过程加深了我们对数字信号调制与解调的认知程度。 关键字：2PSK;调制与解调；MATLAB 引言 当今社会已经步入信息时代，在各种信息技术中，信息的传输及通信起着支撑作用。而对于信息的传输，数字通信已经成为重要的手段。因此，数字信号的调制就显得非常重要。 调制分为基带调制和带通调制。不过一般狭义的理解调制为带通调制。带通调制通常需要一个正弦波作为载波，把基带信号调制到这个载波上，使这个载波的一个或者几个参量上载有基带数字信号的信息，并且还要使已调信号的频谱倒置适合在给定的带通信道中传输。特别是在无线电通信中，调制是必不可少的，因为要使信号能以电磁波的方式发送出去，信号所占用的频带位置必须足够高，并且信号所占用的频带宽度不能超过天线的的通频带，所以基带信号的频谱必须用一个频率很高的载波调制，使期带信号搬移到足够高的频率上，才能够通过天线发送出去。 主要通过对它们的三个参数进行调制，振幅，角频率，和相位。使这三个参量都按时间变化。所以基带的数字信号调制主要有三种方式：FSK，PSK，ASK。在这三种调制的基础上为了得到更高的效果也出现了很多其它的调制方式，如：DPSK，MASK，MFSK，MPSK，APK。它们其中有的一些是将基本的调制方式用在多进制上或者引入了一些新的方式来解决基本调制的一些问题如相位模糊和无法提取位定时信号，另外一些由是组合多种基本的调制方式来达到更好的效果。 基带信号的调制主要分为线性调制和非线性调制，线性调制是指已调信号的频谱结构与原基带信号的频谱结构基本相同，只是占用的频率位置搬移了。而非线性调制则是指它们的结构完全不同不仅仅是频谱搬移，在接收方会出现很多新的频谱分量。在三种基本的调制中，ASK 属于线性调制，而FSK和PSK属于非线性调制。已调信号会在接收方通过各种方式通过解调得到，但是由于噪声和码间串扰，总会有一定的失真。所以人们总是在寻找不同的接收方式来降低误码率，其中的接收方式主要有相干接收和非相干接收。在接收方通过载波的相位信号去检测信号的方法称为相干检测，反之若不利用就称为非相干检测，而对于一些特别的调制有特别的解调方式，如过零检测法。 系统的性能好坏取决于传输信号的误码率，而误码率不仅仅与信道、接收方法有关还和发送端采用的调制方式有很大的关系。我们研究的ASK，FSK，PSK等就主要是发送方的调制方式。
******************* 实践教学 ******************* 兰州理工大学 计算机与通信学院 2013年春季学期 信号处理课程设计 题目：基于MATLAB的FIR滤波器语音信号去噪专业班级：通信工程（1）班 姓名：王兴栋 学号：10250114 指导教师：陈海燕 成绩：
摘要 语音信号在数字信号处理中占有极其重要的地位，因此选择通过对语音信号的研究来巩固和掌握数字信号处理的基本能力十分具有代表性。对数字信号处理离不开滤波器，因此滤波器的设计在信号处理中占有极其重要的地位。而MATLAB 软件工具箱提供了对各种数字滤波器的设计。本论文“在MATLAB平台上实现对语音信号的去噪研究与仿真”综合运用了数字信号处理的各种基本知识，进而对不带噪语音信号进行谱分析以及带噪语音信号进行谱分析和滤波处理。通过理论推导得出相应的结论，再通过利用MATLAB作为编程工具来进行计算机实现比价已验证推导出来的结论。在设计过程中，通过设计FIR数字滤波器和IIR数字滤波器来完成滤波处理。在设计过程中，运用了MATLAB对整个设计中的图形的绘制和一些数据的计算以及仿真。 关键字滤波器；MATLAB；仿真；滤波
前言 语音是语言的声学表现，是人类交流信息最自然、最有效、最方便的手段。随着社会文化的进步和科学技术的发展，人类开始进入了信息化时代，用现代手段研究语音处理技术，使人们能更加有效地产生、传输、存储、和获取语音信息，这对于促进社会的发展具有十分重要的意义，因此，语音信号处理正越来越受到人们的关注和广泛的研究。 语音信号是信息技术处理中最重要的一门科学，是人类社会几步的标志。那么什么是语音？语音是人类特有的功能，也是人类获取外界信息的重要工具，也是人与人交流必不可少的重要手段。那么什么又是信号？那信号是什么呢？信号是传递信息的函数。离散时间信号——序列——可以用图形来表示。 语音信号处理是一门用研究数字信号处理研究信号的科学。它是一新兴的信息科学，同时又是综合多个学科领域的一门交叉科学。语音在我们的日常生活中随时可见，也随处可见，语音很大程度上可以影响我们的生活。所以研究语音信号无论是在科学领域上还是日常生活中都有其广泛而重要的意义。 本论文主要介绍的是的语音信号的简单处理。本论文针对以上问题，运用数字信号学基本原理实现语音信号的处理，在matlab7.0环境下综合运用信号提取，幅频变换以及傅里叶变换、滤波等技术来进行语音信号处理。我所做的工作就是在matlab7.0软件上编写一个处理语音信号的程序，能对语音信号进行采集，并对其进行各种处理，达到简单语音信号处理的目的。 对语音信号的研究，本论文采用了设计两种滤波器的基本研究方法来达到研究语音信号去噪的目的，最终结合图像以及对语音信号的回放，通过对比，得出结论。
SSB单边带信号调制 由双边带过渡 双边带信号虽然抑制了载波，提高了调制效率，但调制后的频带宽度仍是基带信号带宽的2倍，而且上、下边带是完全对称的，它们所携带的信息完全相同。因此，从信息传输的角度来看，只用一个边带传输就可以了。我们把这种只传输一个边带的调制方式称为单边带抑制载波调制，简称为单边带调制(SSB)。 原理部分 采用单边带调制，除了节省载波功率，还可以节省一半传输频带，仅传输双边带信号的一个边带(上边带或下边带)。因此产生单边带信号的最简单方法，就是先产生双边带。然后让它通过一个边带滤波器，只传送双边带信号中的一个边带，这种产生单边带信号的方法称为滤波法。由于理想的滤波器特性是不可能作到的，实际的边带滤波器从带通到带阻总是有一个过渡带，随着载波频率的增加，采用一级载波调制的滤波法将无法实现。这时可采用多级调制滤波的办法产生单边带信号。即采用多级频率搬移的方法实现：先在低频处产生单边带信号，然后通过变频将频谱搬移到更高的载频处。产生SSB 信号的方法还有：相移形成法，混合形成法。 SSB移相法原理图
SSB移相法的形成的SystemView仿真 SSB移相法的形成上边带下边带 数学表达式 为简便起见，设调制信号为单频信号f(t)=Amcosωmt，载波为c(t)=cosωct，则调制后的双边带时域波形为：SDSB(t)＝Amcosωmtcost＝[Amcos(ωc+ωm)t+Amcos(ωc-ωm)t]/2 保留上边带，波形为：SUSB(t)＝[Amcos(ωc+ωm)t]/2＝Am(cosωctcosωmt-sinωctsinωmt)/2 保留下边带，波形为：SLSB(t)＝[Amcos(ωc-ωm)t]/2＝Am(cosωctcosωmt+sinωctsinωmt)/2 上两式中的第一项与调制信号和载波信号的乘积成正比，称为同相分量；而第二项的乘积则是调制信号与载波信号分别移相90°后相乘的结果，称为正交分量。由此可以
数字信号处理实验报告 实验名称：基于MATLAB对语音信号 进行分析及滤波处理 院系：物联网工程学院 班级：电子信息工程1101 姓名：
一、实验目的 综合计算运用数字信号处理的理论知识进行频谱分析和滤波器设计，通过理论推导得出相应的结论，培养发现问题、分析问题和解决问题的能力。并利用MATLAB作为工具进行实现，从而复习巩固课堂所学的理论知识，提高对所学知识的综合应用能力，并从实践上初步实现对数字信号的处理。此外，还系统的学习和实现对语音信号处理的整体过程，从语音信号的采集到分析、处理、频谱分析、显示和储存。 二、实验要求 1.分析原始语音信号的时域特性和频谱特性。 2.设计一个IIR滤波器，用该滤波器对语音信号进行滤波处理，分析滤波后信号的时域特性和频谱特性。 3.设计一个FIR滤波器，用该滤波器对语音信号进行滤波处理，分析滤波后信号的时域特性和频谱特性。 三、实验原理 1.采样定理：在进行模拟/数字信号的转换过程中，当采样频率fs.max大于信号中最高频率fmax的2倍时(fs.max>=2fmax)，采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍；采样定理又称奈奎斯特定理。 采样频率越高，即采样的时间间隔越短，对声音波形的表示越精确。 2.时域信号的FFT分析：信号的频谱分析就是计算机信号的傅里叶变换。连续信号与系统的傅里叶分析显然不便于用计算机进行计算，使其应用受到限制。而FFT是一种时域和频域均离散化的变换，适合数值运算，成为用计算机分析离散信号和系统的有力工具。对连续信号和系统，可以通过时域采样，应用DFT进行近似谱分析。 3.IIR数字滤波器设计原理利用双线性变换设计IIR滤波器（巴特沃斯数字低通滤波器的设计），首先要设计出满足指标要求的模拟滤波器的传递函数Ha(s),然后由Ha(s)通过双线性变换可得所要设计的IIR滤波器的系统函数H(z)。如果给定的指标为数字滤波器的指标，则首先要转换成模拟滤波器的技术指标，这里主要是边界频率Wp和Ws的转换，对ap和as指标不作变化。边界频率的转换关系为∩=2/T tan(w/2)。接着，按照模拟低通滤波器的技术指标根据相应设计公式
信息对抗大作业
一、实验目的。 使用 MATLAB构成一个加性高斯白噪声情况下的2psk 调制解系统，仿真分析使用信道编 码纠错和不使用信道编码时，不同信道噪声比情况下的系统误码率。 二、实验原理。 数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数信道具有带通特性 而不能直接传输基带信号。为了使数字信号在带通信道中传输，必须使用数字基带信号对载波 进行调制，以使信号与信道的特性相匹配。这种用数字基带信号控制载波，把数字基带信号变 换为数字带通信号的过程称为数字调制。 数字调制技术的两种方法：①利用模拟调制的方法去实现数字式调制，即把数字调制看成 是模拟调制的一个特例，把数字基带信号当做模拟信号的特殊情况处理；②利用数字信号的离 散取值特点通过开关键控载波，从而实现数字调制。这种方法通常称为键控法，比如对载波的 相位进行键控，便可获得相移键控（PSK）基本的调制方式。 图 1相应的信号波形的示例 101 数字调相：如果两个频率相同的载波同时开始振荡，这两个频率同时达到正最大值，同时达 到零值，同时达到负最大值，它们应处于" 同相 " 状态；如果其中一个开始得迟了一点，就可能不 相同了。如果一个达到正最大值时，另一个达到负最大值，则称为" 反相 " 。一般把信号振荡一次（一周）作为360 度。如果一个波比另一个波相差半个周期，我们说两个波的相位差180 度，也就是反相。当传输数字信号时， "1" 码控制发 0 度相位， "0" 码控制发 180 度相位。载波的初始相位就 有了移动，也就带上了信息。 相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。在2PSK 中，通常用初始相位0 和π分别表示二进制“1”和“ 0”。因此， 2PSK信号的时域表达式为 (t)=Acos t+) 其中，表示第 n 个符号的绝对相位： = 因此，上式可以改写为
MATLAB中用M文件实现SSB解调 一、课程设计目的 本次课程设计是对通信原理课程理论教学和实验教学的综合和总结。通过这次课程设计，使同学认识和理解通信系统，掌握信号是怎样经过发端处理、被送入信道、然后在接收端还原。 要求学生掌握通信原理的基本知识，运用所学的通信仿真的方法实现某种传输系统。能够根据设计任务的具体要求，掌握软件设计、调试的具体方法、步骤和技巧。对一个实际课题的软件设计有基本了解，拓展知识面，激发在此领域中继续学习和研究的兴趣，为学习后续课程做准备。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 二、课程设计内容 （1）熟悉MATLAB中M文件的使用方法，掌握SSB信号的解调原理，以此为基础用M文件编程实现SSB信号的解调。 （2）绘制出SSB信号解调前后在时域和频域中的波形，观察两者在解调前后的变化，通过对分析结果来加强对SSB信号解调原理的理解。 （3）对信号分别叠加大小不同的噪声后再进行解调，绘制出解调前后信号的时域和频域波形，比较未叠加噪声时和分别叠加大小噪声时解调信号的波形有何区别，借由所得结果来分析噪声对信号解调造成的影响。 （4）在老师的指导下，独立完成课程设计的全部内容，并按要求编写课程设计论文，文中能正确阐述和分析设计和实验结果。 三、设计原理 1、 SSB解调原理 在单边带信号的解调中，只需要对上、下边带的其中一个边带信号进行解调，就能够恢复原始信号。这是因为双边带调制中上、下两个边带是完全对称的，它们所携带的信息相同，完全可以用一个边带来传输全部消息。 单边带解调通常采用相干解调的方式，它使用一个同步解调器，即由相乘器和低通
计算机与信息工程学院实验报告 一、实验目的 1.掌握绝对码、相对码概念及它们之间的变换关系。 2.掌握用键控法产生2FSK信号的方法。 3.掌握2FSK过零检测解调原理。 4.了解2FSK信号的频谱与数字基带信号频谱之间的关系。 二、实验仪器或设备 1.通信原理教学实验系统 TX-6（武汉华科胜达电子有限公司 2011.10） 2.LDS20410示波器（江苏绿扬电子仪器集团有限公司 2011.4.1） 三、总体设计 3.1数字调制 3.1.1实验内容： 1、用示波器观察绝对码波形、相对码波形。 2、用示波器观察2FSK信号波形。 3、用频谱仪观察数字基带信号频谱及2FSK信号的频谱。 3.1.2基本原理： 本实验用到数字信源模块和数字调制模块。信源模块向调制模块提供数字基带信号（NRZ码）和位同步信号BS（已在实验电路板上连通，不必手工接线）。调制模块将输入的绝对码AK（NRZ码）变为相对码BK、用键控法产生2FSK信号。调制模块内部只用+5V电压。 数字调制单元的原理方框图如图1-1所示。 图1-1 数字调制方框图 本单元有以下测试点及输入输出点：
? CAR 2DPSK 信号载波测试点 ? BK 相对码测试点 ? 2FSK 2FSK 信号测试点/输出点，V P-P >0.5V 用1-1中晶体振荡器与信源共用，位于信源单元，其它各部分与电路板上主要元器件对 应关系如下： ? ÷2（A ） U8：双D 触发器74LS74 ? ÷2（B ） U9：双D 触发器74LS74 ? 滤波器A V6：三极管9013，调谐回路 ? 滤波器B V1：三极管9013，调谐回路 ? 码变换 U18：双D 触发器74LS74；U19：异或门74LS86 ? 2FSK 调制 U22：三路二选一模拟开关4053 ? 放大器 V5：三极管9013 ? 射随器 V3：三极管9013 2FSK 信号的两个载波频率分别为晶振频率的1/2和1/4，通过分频和滤波得到。 2FSK 信号（相位不连续2FSK ）可看成是AK 与AK 调制不同载频信号形成的两个2ASK 信号相加。时域表达式为 t t m t t m t S c c 21cos )(cos )()(ωω+= 式中m(t)为NRZ 码。 2FSK 信号功率谱 设码元宽度为T S ，f S =1／T S 在数值上等于码速率， 2FSK 的功率谱密度如图所示。多进制的MFSK 信号的功率谱与二进制信号功率谱类似。 本实验系统中m(t)是一个周期信号，故m(t)有离散谱，因而2FSK 也具有离散谱。 3.2 数字解调 3.2.1 实验内容 1、 用示波器观察2FSK 过零检测解调器各点波形。 3.2.2 基本原理 2FSK 信号的解调方法有：包络括检波法、相干解调法、鉴频法、过零检测法等。
通信与信息工程学院 信息处理综合实验报告 班级：电子信息工程1502班 指导教师： 设计时间：2018/10/22-2018/11/23 评语： 通信与信息工程学院 二〇一八年 实验题目：语音信号分析与处理 一、实验内容 1. 设计内容 利用MATLAB对采集的原始语音信号及加入人为干扰后的信号进行频谱分析，使用窗函数法设计滤波器滤除噪声、并恢复信号。 2．设计任务与要求 1. 基本部分
（1）录制语音信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。 （2）对所录制的语音信号加入干扰噪声，并对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。 （3）分别利用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman 窗几种函数设计数字滤波器滤除噪声，并画出各种函数所设计的滤波器的频率响应。 （4）画出使用几种滤波器滤波后信号时域波形和频谱，对滤波前后的信号、几种滤波器滤波后的信号进行对比，分析信号处理前后及使用不同滤波器的变化；回放语音信号。 2. 提高部分 （5）录制一段音乐信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。 （6）利用MATLAB产生一个不同于以上频段的信号；画出信号频谱图。 （7）将上述两段信号叠加，并加入干扰噪声，尝试多次逐渐加大噪声功率，对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。 （8）选用一种合适的窗函数设计数字滤波器，画出滤波后音乐信号时域波形和频谱，对滤波前后的信号进行对比，回放音乐信号。 二、实验原理 1.设计原理分析 本设计主要是对语音信号的时频进行分析，并对语音信号加噪后设计滤波器对其进行滤波处理，对语音信号加噪声前后的频谱进行比较分析，对合成语音信号滤波前后进行频谱的分析比较。 首先用PC机WINDOWS下的录音机录制一段语音信号，并保存入MATLAB软件的根目录下，再运行MATLAB仿真软件把录制好的语音信号用audioread函数加载入MATLAB仿真软件的工作环境中，输入命令对语音信号进行时域，频谱变换。 对该段合成的语音信号，分别用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman窗几种函数在MATLAB中设计滤波器对其进行滤波处理，滤波后用命令可以绘制出其频谱图，回放语音信号。对原始语音信号、合成的语音信号和经过滤波器处理的语音信号进行频谱的比较分析。 2.语音信号的时域频域分析 在Matlab软件平台下可以利用函数audioread对语音信号进行采样,得到了声音数据变量y,同时把y的采样频率Fs=44100Hz放进了MATALB的工作空间。
广西工学院 数字信号处理课程设计设计题目：语音信号的处理与滤波 系别： 学号： 姓名： 班级： 指导教师： 完成日期：
目录 1. 摘要 (3) 2.MATLAB简介 (3) 3.设计目的 (4) 4.设计内容 (4) 5.设计原理 (4) 6.设计步骤 (5) 7.总结与分析 (10) 8.参考资料 (10)
摘要 本课题分析了数字信号处理课程的重要性及特点，可以帮助理解与掌握课程中的基本概念、基本原理、基本分析方法；并利用MATLAB对语音信号进行分析和处理，要求采集语音信号后，在MATLAB软件平台进行频谱分析。用设计的数字滤波器对语音信号进行了滤波 MATLAB简介 MATLAB 是一种对技术计算高性能的语言。它集成了计算，可视化和编程于一个易用的环境中，在此环境下，问题和解答都表达为我们熟悉的数学符号。典型的应用有： ?数学和计算 ?算法开发 ?建模，模拟和原形化 ?数据分析，探索和可视化 ?科学与工程制图 ?应用开发，包括图形用户界面的建立 MATLAB是一个交互式的系统，其基本数据元素是无须定义维数的数组。这让你能解决很多技术计算的问题，尤其是那些要用到矩阵和向量表达式的问题。而要花的时间则只是用一种标量非交互语言(例如C或Fortran)写一个程序的时间的一小部分。 . 名称“MATLAB”代表matrix laboratory（矩阵实验室）。MATLAB最初是编写来提供给对由LINPACK和EINPACK工程开发的矩阵软件简易访问的。今天，MATLA B使用由LAPACK和ARPACK工程开发的软件，这些工程共同表现了矩阵计算的软件中的技术发展。 MATLAB已经与许多用户输入一同发展了多年。在大学环境中，它是很多数学类、工程和科学类的初等和高等课程的标准指导工具。在工业上，MATLAB是高产研究、开发和分析所选择的工具。 MATLAB以一系列称为工具箱的应用指定解答为特征。对多数用户十分重要的是，工具箱使你能学习和应用专门的技术。工具箱是是MATLAB函数（M-文件）的全面的综合，这些文件把MATLAB的环境扩展到解决特殊类型问题上。具有可用工具箱的领域有：信号处理，控制系统神经网络，模糊逻辑，小波分析，模拟等等。
二进制数字调制与解调系统的设计 MATLAB 及SIMULINK 建模环境简介 MATLAB 是美国MathWorks 公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB 和SIMULINK 两大部分。 Simulink 是MATLAB 最重要的组件之一，它提供一个动态系统建模、仿真和综合分析的集成环境。在该环境中，无需大量书写程序，而只需要通过简单直观的鼠标操作，就可构造出复杂的系统。Simulink 具有适应面广、结构和流程清晰及仿真精细、贴近实际、效率高、灵活等优点，并基于以上优点Simulink 已被广泛应用于控制理论和数字信号处理的复杂仿真和设计。同时有大量的第三方软件和硬件可应用于或被要求应用于Simulink 。 Simulink 是MATLAB 中的一种可视化仿真工具， 是一种基于MATLAB 的框图设计环境，是实现动态系统建模、仿真和分析的一个软件包，被广泛应用于线性系统、非线性系统、数字控制及数字信号处理的建模和仿真中。Simulink 可以用连续采样时间、离散采样时间或两种混合的采样时间进行建模，它也支持多速率系统，也就是系统中的不同部分具有不同的采样速率。为了创建动态系统模型，Simulink 提供了一个建立模型方块图的图形用户接口(GUI) ，这个创建过程只需单击和拖动鼠标操作就能完成，它提供了一种更快捷、直接明了的方式，而且用户可以立即看到系统的仿真结果。 数字通信系统的基本模型 从消息传输角度看,该系统包括了两个重要交换,即消息与数字基带信号之间的交换,数字基带信号与信道信号之间的交换.通常前一种交换由发收端设备完成.而后一种交换则由调制和解调完成. 数字通信系统模型 一、2ASK 调制解调 基本原理 2ASK 是利用载波的幅度变化来传递数字信息，而其频率和初始相位保持不变。 其信号表达式为： ，S (t)为单极性数字基带信号。 t t S t e c ωcos )()(0 ?=
《计算机仿真技术》 基于MATLAB的加噪语音信号的滤波学生姓名： 专业：电子信息工程 班级： 学号： 指导教师： 完成时间：2017年12月
一．滤波器的简述 在MATLAB环境下IIR数字滤波器和FIR数字滤波器的设计方法即实现方法，并进行图形用户界面设计，以显示所介绍迷你滤波器的设计特性。 在无线脉冲响应（IIR）数字滤波器设计中，先进行模拟滤波器的设计，然后进行模拟-数字滤波器转换，即采用脉冲响应不变法及双线性Z变化法设计数字滤波器，最后进行滤波器的频带转换。在有限脉冲响应（FIR）数字滤波器设计中，讨论了FIR线性相位滤波的特点和用窗口函数设计FIR数字滤波器两个问题。两类滤波器整个过程都是按照理论分析、编程设计、集体实现的步骤进行的。为方便分析直观者直观、形象、方便的分析滤波器的特性，创新的设计出图形用户界面---滤波器分析系统。整个系统分为两个界面，其内容主要包括四个部分：System（系统）、Analysis（分析）、Tool（工具）、Help（帮助）。 数字滤波在DSP中占有重要地位。数字滤波器按实现的网络结构或者从单位脉冲响应，分为IIR（无限脉冲响应）和FIR（有限脉冲响应）滤波器。如果IRR 滤波器和FIR滤波器具有相同的性能，那么通常IIR滤波器可以用较低的阶数获得高的选择性，执行速度更快，所有的存储单元更少，所以既经济又高效。二．设计要求 1.在matlab平台上录制一段语音信号； 2.完成语音信号的谱分析； 3.对语音信号进行加噪以及加噪后信号的谱分析； 4.选择合适的滤波器进行滤波，确定相关指标； 5.实现滤波过程，显示滤波后的结果，并进行谱分析。 三．实验内容与步骤 1、语音信号的录入
MATLAB仿真信号的相位调制与解调 专业：通信与信息系统 姓名：赵* 学号:********* 指导老师：****教授
摘要 Psk调制是通信系统中最为重要的环节之一，Psk调制技术的改进也是通信系统性能提高的重要途径。本文首先分析了数字调制系统的基本调制解调方法，然后，运用Matlab及附带的图形仿真工具——Simulink设计了这几种数字调制方法的仿真模型。通过仿真，观察了调制解调过程中各环节时域和频域的波形，并结合这几种调制方法的调制原理，跟踪分析了各个环节对调制性能的影响及仿真模型的可靠性。最后，在仿真的基础上分析比较了各种调制方法的性能，并通过比较仿真模型与理论计算的性能，证明了仿真模型的可行性。另外，本文还利用Matlab的图形用户界面（GUI）功能为仿真系统设计了一个便于操作的人机交互界面，使仿真系统更加完整，操作更加方便。 关键词：数字调制；分析与仿真；Matlab；Simulink；PSK；QPSK；
1.数字调制技术 (2) 2.PSK调制系统 (3) 2.1 QPSK调制部分，原理框图如图七所示 (6) 2.2 QPSK解调部分，原理框图如图八所示： (8) 3.用Simulink实现PSK调制 (9) 3.1 2PSK仿真 (9) 3.1.1调制 (9) 3.1.2 解调仿真 (12) 3.2 QPSK仿真 (13) 3.2.1 QPSK调制框图 (13) 参考文献 (18)
1.数字调制技术 通信按照传统的理解就是信息的传输与交换。在当今信息社会，通信则与遥感，计算技术紧密结合，成为整个社会的高级“神经中枢”。没有通信，人类社会是不可想象的。一般来说，社会生产力水平要求社会通信水平与之相适应。若通信水平跟不上，社会成员之间的合作程度就受到限制。可见，通信是十分重要的。 通信传输的消息是多种多样的，可以是符号的，文字的，数据和图像的等等。各种不同的消息可以分为两类：一类称为离散消息；另一类称为连续消息。离散消息的状态是可数的或离散的，比如符号，文字或数据等。离散消息也称数字消息。而连续消息则是其状态连续变化的消息，例如，连续变化的语音，图像等。连续消息也称模拟消息。因此按照信道中传输的是模拟信号还是数字信号可以将通信系统分为模拟通信系统和数字通信系统。 数字通信有以下突出的特点：第一，数字信号传输时，信道噪声或干扰所造成的差错，原则上是可以控制的。第二，当需要保密的时候，可以有效的对基带信号进行人为的“扰乱”，即加上密码。 数字通信系统可以用下图表示： →→→→→→→→信数信信数信 信源 道 字受道源字信 息编编调 解译译信 源 码码调码码者 制 道 器 器 器 器 器 器 图一 数字通信在近20年来得到了迅速的发展，其原因是： （1） 抗干扰能力强 （2） 便于进行各种数字信号处理 （3） 易于实现集成化 （4） 经济效益正赶上或超过模拟通信 （5） 传输与交换可结合起来，传输电话与传输数据也可结合起来，成为一个 统一整体，有利于实现综合业务通信网。
单边带通信的特点 单边带通信的特点在业余无线电短波通信中，单边带 制(SSB)占据着极其重要的位置。与调幅和调频制相比较， 其主要特点如下：一、节约频带。信道容量可增加一倍，从而部分的解决了短波通信中电台拥挤的问题。二、节约功率。与调幅制相比，在一定条件下，要达到相同的通信效果，单边带通讯用到的发射功率仅是调幅通讯用到的发射机功率 的1/16.三、没有门槛效应。即使是在微弱的信号下，仍可 勉强维持通信。四、各信道间相互干扰小。保密性也稍好。 五、网络通信的可能性。不会象调幅制那样出现啸叫声。其缺点则主要是技术难度大。对频率稳定度、滤波器的选择性和放大器的线性要求较高。各种调制式都有它的特点和适用范围，单边带通信和双边带通信相比较也有它突出的特点。下面我们从调制信号占据的带宽，发射机功率的有效利用程度和抵抗传输失真等主要方面，来分析单边带制的特点。(一)、节约频带由于单边带通信只是利用调幅信号中的一个边带 进行通信，由此能节省频带。设被传偷的声频信号的最高频率分量为Fmax，并用Bam、Bssb、Bfm分别表示调幅、单边带和调频时信号的带宽，则它们分别等于： Bam=2FmaxBssb~=FmaxBfm~=2(1+mf)FmaxBfm式中的mf是调频指数，mf=Aω/Ω等于最大频移与最高音频角频率
之比，其值一般取1.6-5视不同用途而异，以Fmax=3KHz 为例，可得：Bam=6KHzBssb=3KHzBfm=18KHz(取mf=2)从比较看出，单边带信号频谱占用的频带宽度最小，因而对高频频谱利用得最经济，在同样的有限高频频段内，就可以使无线电波道容量比用调幅制时增加一倍，从而部分地解决了短波波段空中频谱拥挤的问题。(二)、功率节约在双边带通信中，由于调幅波是由三个分量合成的，因此调幅波的功率就分配在载频和两个边带上。载波成分电流振幅最大，而边频成份电流振幅最小，因此，一个幅调波的总功率的大部分就消耗在不代表信号意义的载频上，而真正含有信号意义的每个边频的功率则是很少的。设载波功率Po=1/2I2cR，式中R为负载，每一边带功率： PΩ=1/2(MIc/2)2R=1/8M2I2cR，调幅波的总功率： P=Po+2PΩ=1/2I2cR+2(1/8M2I2cR) =1/2I2cR(1+M2/2)-Po(1+M2/2).当M=1时，调幅波的总功率是载波功率的150%，其中载波功率占全部功率的2/3，二个边惜共占全部功率的1/3，每个边带只占全部功率的1/6，也就是说，在调幅波中，不代表信号意义的载波分量，白白地消耗了2/3的总功率，而包含有信号意义的每一个边带只分配到1/6的总功率，上面我们假设调幅度M=1，实际上语言调制时M有大有小，最大等于1，不然就产生过调制，所以语言调制时M只能取其平均值，通常M平均=0.3，这样
基于MATLAB的语音信号滤波处理 题目：基于MATLAB的语音信号滤波处理 课程：数字信号处理 学院：电气工程学院 班级： 学生： 指导教师： 二O一三年十二月
目录CONTENTS 摘要 一、引言 二、正文 1.设计要求 2.设计步骤 3.设计内容 4.简易GUI设计 三、结论 四、收获与心得 五、附录
一、引言 随着Matlab仿真技术的推广，我们可以在计算机上对声音信号进行处理，甚至是模拟。通过计算机作图，采样，我们可以更加直观的了解语音信号的性质，通过matlab编程，调用相关的函数，我们可以非常方便的对信号进行运算和处理。 二、正文 2.1 设计要求 在有噪音的环境中录制语音，并设计滤波器去除噪声。 2.2 设计步骤 1.分析原始信号，画出原始信号频谱图及时频图，确定滤波器类型及相关指标； 2.按照类型及指标要求设计出滤波器，画出滤波器幅度和相位响应，分析该滤波器是否符合要求； 3.用所设计的滤波器对原始信号进行滤波处理，画出滤波后信号的频谱图及时频图； 4.对滤波前的信号进行分析比对，评估所设计滤波器性能。 2.3 设计内容 1.原始信号分析
分析信号的谱图可知，噪音在1650HZ和3300HZ附近的能量较高，而人声的能量基本位于1000HZ以下。因此，可以设计低通滤波器对信号进行去噪处理。 2.IIR滤波器设计 用双线性变换法分别设计了巴特沃斯低通滤波器和椭圆低通滤波器和带阻滤波器： ①巴特沃斯滤波器 fp=800;fs=1300;rs=35;rp=0.5; 程序代码如下： fp=800;fs=1300;rs=35;rp=0.5;Fs=44100; wp=2*Fs*tan(2*pi*fp/(2*Fs));ws=2*Fs*tan(2*pi*fs/(2*Fs)); [n,wn]=buttord(wp,ws,rp,rs,'s'); [b,a]=butter(n,wn,'s'); [num,den]=bilinear(b,a,Fs); [h,w]=freqz(num,den,512,Fs);
数字信号载波调制实验指导书 数字信号载波调制实验 一、实验目的 1、运用MATLAB 软件工具仿真数字信号的载波传输．研究数字信号载波调制ASK 、FSK 、PSK 在不同调制参数下的信号变化及频谱。 2，研究频移键控的两种解调方式；相干解调与非相干解调。 3、了解高斯白噪声方差对系统的影响。 4、了解伪随机序列的产生，扰码及解扰工作原理。 二、实验原理 数字信号载波调制有三种基本的调制方式：幅度键控（ASK ），频移键控（FSK ）和相移键控（PSK ）。它们分别是用数字基带信号控制高频载波的参数如振幅、频率和相位，得到数字带通信号。在接收端运用相干或非相干解调方式，进行解调，还原为原数字基带信号。 在幅度键控中，载波幅度是随着调制信号而变化的。最简单的形式是载波在 二进制调制信号1或0的控制下通或断，这种二进制幅度键控方式称为通—断键控（00K ）。二进制幅度键控信号的频谱宽度是二进制基带信号的两倍。 在二进制频移键控中，载波频率随着调制信号1或0而变，1对应于载波频率f 1，0对应于载波频率f 2，二进制频移键控己调信号可以看作是两个不同载频的幅度键控已调信号之和。它的频带宽度是两倍基带信号带宽（B ）与21||f f -之和。 在二进制相移键控中，载波的相位随调制信号1或0而改变，通常用相位0°和180°来分别表示1或0，二进制相移键控的功率谱与通一断键控的相同，只是少了一个离散的载频分量。 m 序列是最常用的一种伪随机序列，是由带线性反馈的移位寄存器所产生的序列。它具有最长周期。由n 级移位寄存器产生的m 序列，其周期为21,n m -序列有很强的规律性及其伪随机性。因此，在通信工程上得到广泛应用，在本实验中用于扰码和解扰。 扰码原理是以线性反馈移位寄存器理论作为基础的。在数字基带信号传输中，将二进制数字信息先作“随机化”处理，变为伪随机序列，从而限制连“0”
山东建筑大学信电学院课程设计说明书 语音信号分析与处理 摘要 用MATLAB对语音信号进行分析与处理，采集语音信号后，在MATLAB软件平台进行频谱分析;并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，设计合适的滤波器滤除噪声，恢复原信号。 数字滤波器是数字信号处理的基础，用来对信号进行过滤、检测和参数估计等处理。IIR数字滤波器最大的优点是给定一组指标时，它的阶数要比相同组的FIR 滤波器的低的多。信号处理中和频谱分析最为密切的理论基础是傅立叶变换（FT）。离散傅立叶变换（DFT）和数字滤波是数字信号处理的最基本内容。 关键词：MATLAB;语音信号；加入噪声；滤波器；滤波 1. 设计目的与要求 （1）待处理的语音信号是一个在20Hz~20kHz频段的低频信号。 （2）要求MATLAB对语音信号进行分析和处理，采集语音信号后，在MATLAB平台进行频谱分析；并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，设计合适的滤波器进行滤除噪声，恢复原信号。 1 山东建筑大学信电学院课程设计说明书
2. 设计步骤 （1）选择一个语音信号或者自己录制一段语音文件作为分析对象； （2）对语音信号进行采样，并对语音信号进行FFT频谱分析，画出信号的时域波形图和频谱图； （3）利用MATLAB自带的随机函数产生噪声加入到语音信号中，对语音信号进行回放，对其进行FFT频谱分析； （4）设计合适滤波器，对带有噪声的语音信号进行滤波，画出滤波前后的时域波形图和频谱图，比较加噪前后的语音信号，分析发生的变化； （5）对语音信号进行回放，感觉声音变化。 3. 设计原理及内容 3.1 理论依据 （1）采样频率：采样频率（也称采样速度或者采样率）定义了每秒从连续信号中提取并组成离散信号的采样个数，它用赫兹（Hz）来表示。采样频率只能用 于周期性采样的采样器，对于非周期采样的采样器没有规则限制。通俗的讲，采样频率是指计算机每秒钟采集多少个声音样本，是描述声音文件的音质、音调，衡量声卡、声音文件的质量标准。采样频率越高，即采样的间隔时间越短，则在单位之间内计算机得到的声音样本数据就越多，对声音波形的表示也越精确。（2）采样位数：即采样值或取样值，用来衡量声音波动变化的参数。 （3）采样定理：在进行模拟/数字信号的的转换过程中，当采样频率f大于信s.max 号中，最高频率f的2倍时，即：f>=2f，则采样之后的数字信号完整的maxmaxs.max 保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的 5~10倍；采样频率又称乃奎斯特定理。 （4）时域信号的FFT分析：信号的频谱分析就是计算信号的傅立叶变换。连续信号与系统的傅立叶分析显然不便于直接用计算机进行计算，使其应用受到限制。而FFT是一种时域和频域均离散化的变换，适合数值计算，成为用计算机分析 离2 山东建筑大学信电学院课程设计说明书 散信号和系统的的有力工具。对连续信号和系统，可以通过时域采样，应用DFT 进行近似谱分析。
一matlab常用函数 1、特殊变量与常数 ans 计算结果的变量名computer 确定运行的计算机eps 浮点相对精 度Inf 无穷大I 虚数单位inputname 输入参数名NaN 非 数nargin 输入参数个数nargout 输出参数的数目pi 圆周 率nargoutchk 有效的输出参数数目realmax 最大正浮点数realmin 最小正浮点数varargin 实际输入的参量varargout 实际返回的参量操作符与特殊字符+ 加- 减* 矩阵乘法 .* 数组乘（对应元素相乘）^ 矩阵幂 .^ 数组幂（各个元素求幂）\ 左除或反斜杠/ 右除或斜面杠 ./ 数组除（对应元素除）kron Kronecker张量积: 冒号() 圆括[] 方括 . 小数点 .. 父目录 ... 继续, 逗号（分割多条命令）; 分号（禁止结果显示）% 注释! 感叹号' 转置或引用= 赋值== 相等<> 不等 于& 逻辑与| 逻辑或~ 逻辑非xor 逻辑异或 2、基本数学函数 abs 绝对值和复数模长acos,acodh 反余弦，反双曲余弦acot,acoth 反余切，反双曲余切acsc,acsch 反余割，反双曲余割angle 相角asec,asech 反正割，反双曲正割secant 正切asin,asinh 反正弦，反双曲正 弦atan,atanh 反正切，双曲正切tangent 正切atan2 四象限反正 切ceil 向着无穷大舍入complex 建立一个复数conj 复数配 对cos,cosh 余弦，双曲余弦csc,csch 余切，双曲余切cot,coth 余切，双曲余切exp 指数fix 朝0方向取整floor 朝负无穷取整*** 最大公因数imag 复数值的虚部lcm 最小公倍数log 自然对数log2 以2为底的对数log10 常用对数mod 有符号的求余nchoosek 二项式系数和全部组合数real 复数的实部rem 相除后求余round 取整为最近的整数sec,sech 正割，双曲正割sign 符号数sin,sinh 正弦，双曲正弦sqrt 平方根tan,tanh 正切，双曲正切 3、基本矩阵和矩阵操作 blkding 从输入参量建立块对角矩阵eye 单位矩阵linespace 产生线性间隔的向量logspace 产生对数间隔的向量numel 元素个数ones 产生全为1的数组rand 均匀颁随机数和数组randn 正态分布随机数和数组zeros 建立一个全0矩阵colon) 等间隔向量cat 连接数组diag 对角矩阵和矩阵对角线fliplr 从左自右翻转矩阵flipud 从上到下翻转矩阵repmat 复制一个数组reshape 改造矩阵roy90 矩阵翻转90度tril 矩阵的下三角triu 矩阵的上三角dot 向量点集cross 向量叉 集ismember 检测一个集合的元素intersect 向量的交 集setxor 向量异或集setdiff 向是的差集union 向量的并集数值分析和傅立叶变换cumprod 累积cumsum 累 加cumtrapz 累计梯形法计算数值微分factor 质因子inpolygon 删除多边形区域内的点max 最大值mean 数组的均 值mediam 中值min 最小值perms 所有可能的转 换polyarea 多边形区域primes 生成质数列表prod 数组元素的乘积rectint 矩形交集区域sort 按升序排列矩阵元 素sortrows 按升序排列行std 标准偏差sum 求
青海师范大学毕业论文 论文题目：数字调制解调的MATLAB仿真 系别：物理系 专业：电子信息工程 班级：05 B 学生姓名：梁俊花 学号：20050811217 指导教师姓名：李文全 职称：教授 最后完成时间：2009-5-10
【内容摘要】 设计了二进制振幅键控(2ASK)、二进制移频键控(2FSK) 、二进制移相键控(2PSK)调制解调系统的工作流程图,并得用了MATLAB软件对该系统的动态进行 了模拟仿真,得用仿真的结果,从而衡量数字信号的传输质量. 【关键词】 调制解调、2ASK、2FSK、2PSK、2DPSK、MATLAB 【Abstract】 The work stream diagrams of 2ASK、2FSK、2PSK are designed .MATLAB softwave is used to simulate the modem system by the scatter diagrams and wave diagrams, then the transmit quality of digital signal can be measured. 【Keys】 Amodulate and ademodulate 、2ASK、2FSK、2PSK、2DPSK、MATLAB 一、数字调制解调的概述 在通信系统中,信道的频段往往是很有限的,而原始的通信信号 的频段与信道要求的频段是不匹配的,这就要求将原始信号进行调制 再进行发送.相应的在接收端对调制的信号进行解调,恢复原始的信号,而且调制解调还可以在一定程度上抑制噪声对通信信号的干扰. 调制解调技术按照通信信号是模拟的还是数字的可分为模拟调 制解调和数字调制解调。数字调制的基本方式可以归结为3类：振幅 键控（ASK）、频率键控（FSK）和移相键控（PSK）。此外还有这3 类的混合方式。 对于数字调制信号，为了提高系统的抗噪声性能，衡量系统性