目录
一．问题描述-----------------------------------------3 二．实验原理-----------------------------------------4 三．源程序-------------------------------------------6 四．数据测试----------------------------------------16 五．调试分析----------------------------------------22 六．用户使用手册------------------------------------23 七．心得体会----------------------------------------24
一、问题描述
1.使用matlab编程完成HDB3的编码与解码。
2.课程设计需要运用MATLAB编程实现2ASK,2FSK,2PSK，
2DPSK调制解调过程，并且输出其源码，调制后码元以及解
调后码元的波形。
二、实验原理
1.HDB3编码解码原理
HDB3码：三阶高密度双极性码。
HDB3码与二进制序列的关系：
(1)二进制信号序列中的“0”码在HDB3码中仍编为“0”码，二进制信号中“1”码，在HDB3码中应交替地成+1和-1码，但序列中出现四个连“0”码时应按特殊规律编码（引入传号交替反转码的“破坏点”V码）；
(2)二进制序列中四个连“0”按以下规则编码：信码中出现四个连“0”码时，要将这四个连“0”码用000V 或B00V 取代节来代替(B 和V 也是“1”码，可正、可负)。这两个取代节选取原则是，使任意两个相邻v 脉冲间的传号数为奇数时选用000V 取代节，偶数时则选用B00V 取代节。 2.二进制数字调制技术原理
数字信号的传输方式分为基带传输和带通传输，在实际应用中，大多数
信道具有带通特性而不能直接传输基带信号。为了使数字信号在带通信道中传输，必须使用数字基带信号对载波进行调制，以使信号与信道的特性相匹配。这种用数字基带信号控制载波，把数字基带信号变换为数字带通信号的过程称为数字调制。通常使用键控法来实现数字调制，比如对载波的振幅、频率和相位进行键控。
(1)2ASK 信号的产生方法通常有两种：模拟调制和键控法。解调有相干解调和非相干解调。P=1时f(t)=Acoswt;p=0时f(t)=0;其功率谱密度是基带信号功率谱的线性搬移
(2) 一个2FSK 信号可以看成是两个不同载波的2ASK 信号的叠加。其解调和解调方法和ASK 差不多。2FSK 信号的频谱可以看成是f1和f2的两个2ASK 频谱的组合。
(3) 2PSK 以载波的相位变化作为参考基准的，当基带信号为0时相位相对于初始相位为0， 当基带信号为1时相对于初始相位为180°。 (4) 2DPSK 调制原理方框图如下图。
间接法信号调制器原理方框图
2DPSK 信号的解调，主要有两种方法，即相位比较法和相干解调法。相干解调法原理方框图如下图：
带通滤
相乘
低通滤波
抽样判决
逆码变换
本地载波
提取
载波
移相
码变换
A(t)
S(t)
相干解调法原理方框图
三、源程序
clear
clc
a=20; %a表示元素个数，可以修改码元为任意个数如1024
r=unidrnd(2,1,a-1); %随机产生离散均匀二值分布（1，2），共a-1个
rn=r-1; %每个元素减1变为（0，1）二值随机数列rn=[1,rn]; %跟书上一致，从1开始随机产生离散均
匀二值分布（1，2），共a个
ori=rn; %将原序列保存起来，便于后面解码后比较
ori
l=1;
for k=1:a-3
if abs(rn(k))==1
m=rn(k); %m记录可能存在的连续4个0前面非零码元的符号
end
if rn(k)==0
if rn(k+1)==0
if rn(k+2)==0
if rn(k+3)==0 %检测是否连续4个码元都是0
rn(k+3)=2*m; %用2表示书上的V
end
end
end
end
end
rn ; %插入V(2)后 checked l=1; %l=1表示两个V之间有偶数个非零码，0是偶数，所以初始为1，l=-1表示两个V之间有奇数个非零码
for s=1:a
if abs(rn(s))==2
for d=s+1:a
if abs(rn(d))==1 %V之间的非零码只有+1，-1遇到一个一使l变一次号，表示偶数或奇数
l=-l;
end
if abs(rn(d))==2
if l~=1
l=1; %不需要插入B时若检测到V 也应把计数器清零
%最初版本的升级处1
else %检测到下一个V时，若为偶数个，则插入B
rn(d-3)=3*(-(rn(d-4))/(abs(rn(d-4)))); %3代表B,插入不带符号的B
rn(d:a)=-rn(d:a);
end %V后的符号再交替
end
end
end
end
rn %到此处完成了插入不带符号的B %for s=1:a
% if abs(rn(s))==3 %找到B
%rn(s)=rn(s)*(-(rn(s-1))/(abs(rn(s-1)))); %B（3）与前一个位置的带符号的归一值相乘再取反，实现符号B(3)的极性与前一非符号的相反
%B后面第三个就是接下来的V,从它开始非零码正负号交替变化
%end
%end
HDB=rn; %给B(3)添加了符号，并且实现了V后的符号再交替
HDB
%以上便实现了HDB3的编码
%下面进行解码
for k=1:a
if abs(rn(k))==2
rn(k-3)=0;rn(k)=0; %每个V（2）都变成0，V前面第三个有可能是B(3)有可能是0也恢复为0
end
if abs(rn(k))==1
rn(k)=1;
end
end
rn %解码后的恢复序列
rn-ori %解码与原码比较全为0则解码正确
s=menu('通信原理','2ASK','2PSK','2FSK','2DPSK')
switch s
case 1,scolor='2ASK';
n=8;N=100;K=4;
a=randint(1,n);
bita=[];sl=[];
bitRate=1e3;fc=1e3;%载频1KHZ
t=linspace(0,1/bitRate,N);
for i=1:length(a)
if a(i)==0
bit1=zeros(1,N);
else
bit1=ones(1,N);
end
bita=[bita,bit1];
c=sin(2*pi*t*fc);
sl=[sl c];
end
figure(1);
subplot(K,1,1);plot(bita,'LineWidth',1.5),title('基带信号'),grid on;axis([0,N*length(a),-2.5,2.5]);
tz=bita*6.*sl;
subplot(K,1,2);plot(tz,'LineWidth',1.5);title('ASK调制后信号');grid on;
signal=awgn(tz,80,'measured');
subplot(K,1,3);plot(signal,'LineWidth',1.5),grid on;title('信号+噪声')
Fs=3e3;
[b,a]=ellip(4,0.1,40,[999.9,1000.1]*2/Fs);%设计IIR带通滤波器，阶数为4，通带纹波0.1，阻带衰减40DB
sf=filter(b,a,signal);%信号通过该滤波器
figure(2);
K1=4;
subplot(K1,1,1);plot(sf,'LineWidth',1.5),grid on;title('BPF')
signal2=abs(sf); %乘同频同相sin
subplot(K1,1,2);plot(signal2,'LineWidth',1.5),grid on;title('全波整流器');
Fs=3e3;%抽样频率400HZ
[b,a]=ellip(4,0.1,40,[50]*2/Fs);%设计IIR低通滤波器
sf1=filter(b,a,signal2);%信号通过该滤波器，输出信号sf
subplot(K1,1,3);plot(sf1,'LineWidth',1.5),grid on;title('LPF');
sf2=[];
LL=fc/bitRate*N;
i=LL/2;
bitb=[];
while (i<=length(sf1)) %判决
sf2=[sf2,sf1(i)>=0.001];
i=i+LL;
end
for i=1:length(sf2)
if sf2(i)==0
bit1=zeros(1,N);
else
bit1=ones(1,N);
end
bitb=[bitb,bit1];
end
figure(1);
subplot(K,1,4); plot(bitb,'LineWidth',1.5),grid on;title('解调后信号'); axis([0,N*length(sf2),-2.5,2.5]);
case 2,scolor='2PSK';
l=linspace(0,pi,50);% 数据初始化
t=linspace(0,9*pi,450);
b=1:1:9;
out=1:1:450;
f=1:1:450;
g=1:1:450;
w1=2 %正弦波f1的频率,可以根据自己想要的频率在此改写
%正弦波f2的频率,可以根据自己想要的频率在此改写
f1=sin(w1*l);
figure(1);
f2=sin(w1*l+pi);
figure(1);
subplot(2,1,1),plot(l,f1),axis([0 pi -1.2 1.2]),xlabel('t'),ylabel('f1');%画出f1信号波形
subplot(2,1,2),plot(l,f2),axis([0 pi -1.2 1.2]),xlabel('t'),ylabel('f2');%画出f2信号波形
a=[0 1 0 0 0 1 1 0 1]
for i=1:9 %2pSK编码
if a(i)==0
for k=1:50 %如果二进制原码为0则输出f1波形
out(k+50*(i-1))=f1(k);
end
else
for j=1:50
out(j+50*(i-1))=f2(j); %r如果二进制原码为1则输出f2波形
end
end
end
for i=1:9 %2PSK解码
n=0;m=0;
for j=1:50
if out(j+50*(i-1))-f1(j)==0
n=n+1;
else if out(j+50*(i-1))-f2(j)==0
m=m+1;
end
end
end
if n>m
b(i)=0;
else b(i)=1;
end
end
b
for i=1:9 %画出解码后的波形 ,包括原码和解码出的码,进行对比
for j=1+50*(i-1):50*i
f(j)=a(i);
g(j)=b(i);
end
end
figure(2);
subplot(3,1,1),plot(t,f),axis([0 9*pi -0.2 1.2]),xlabel('t'),ylabel('数字基带调制原码');
subplot(3,1,2),plot(t,out),axis([0 9*pi -1.2 1.2]),xlabel('t'),ylabel('调制好的波形');
subplot(3,1,3),plot(t,g),axis([0 9*pi -0.2 1.2]),xlabel('t'),ylabel('解码得到的码');
[Num,Rat]=biterr(f,g)
case 3,scolor='2FSK';
%==生成随机码元、基带信号、调制================%
n=8;%随机码元个数
N=100;%模拟一个码元的点数
K=4;%figure1画四个小图
a=randint(1,n)%码元生成
bita=[];%定义空数组，存放基带信号
for i=1:length(a)
if a(i)==0
bit1=zeros(1,N);
else
bit1=ones(1,N);
end
bita=[bita,bit1];%基带信号
end
figure(1);
subplot(K,1,1);plot(bita,'LineWidth',1.5),title('基带信号'),grid on;axis([0,N*length(a),-2.5,2.5]);
bitRate=1e3;%每一个码元中采样点的间隔宽度0.001s
fc=1e3;%载频1KHZ
t=linspace(0,1/bitRate,N);
tz=[];
c1=sin(2*pi*t*fc);%载波
c2=sin(2*pi*t*fc*2);%载波
for i=1:length(a)
if a(i)==1
tz=[tz,c1];
else
tz=[tz,c2];
end
end
subplot(K,1,2);plot(tz,'LineWidth',1.5);title('2FSK已调信号');grid on;
signal=awgn(tz,20,'measured');%加噪
subplot(K,1,3);plot(signal,'LineWidth',1.5),grid on;title('信号+噪声')
%===========解调============%
Fs=5e3;%采样频率
[b1,a1]=ellip(4,0.1,40,[999.9,1000.1]*2/Fs);%设计IIR带通滤波器，阶数为4，通带纹波0.1，阻带衰减40DB
[b2,a2]=ellip(4,0.1,40,[1999.9,2000.1]*2/Fs);
sa=filter(b1,a1,signal);%信号通过该滤波器
sb=filter(b2,a2,signal);
figure(2);
K1=3;%figure2画（3*2)幅图表示解调过程
subplot(K1,2,1);plot(sa,'LineWidth',1.5),grid on;title('BPF')
subplot(K1,2,2);plot(sb,'LineWidth',1.5),grid on;title('BPF')
%======相乘器=============%
t=linspace(0,1/bitRate,N);
c1=sin(2*pi*t*fc);
c2=sin(2*pi*t*fc*2);
sia=[];
sib=[];
for i=1:n
sia=[sia,c1];
sib=[sib,c2];
end
siga=sa.*sia;%乘同频同相sia
sigb=sb.*sib;%乘同频同相sib
subplot(K1,2,3);plot(siga,'LineWidth',1.5),grid on;title('相乘器'); subplot(K1,2,4);plot(sigb,'LineWidth',1.5),grid on;title('相乘器'); %========LPF================%
Fs=5e3;%抽样频率400HZ
[b1,a1]=ellip(4,0.1,40,[50]*2/Fs);%设计IIR低通滤波器
[b2,a2]=ellip(4,0.1,40,[50]*2/Fs);%设计IIR低通滤波器
sfa=filter(b1,a1,siga);%信号通过该滤波器，输出信号sfa
sfb=filter(b2,a2,sigb);%信号通过该滤波器，输出信号sfb
subplot(K1,2,5);plot(sfa,'LineWidth',1.5),grid on;title('LPF'); subplot(K1,2,6);plot(sfb,'LineWidth',1.5),grid on;title('LPF');
%===========抽样判决==================%
s2a=[];s2b=[];
LL=fc/bitRate*N;
i1=LL/2;i2=LL;
bitb=[];
while (i1<=length(sfa)) %判决过程
s2a=[s2a,sfa(i1)>=0];
i1=i1+LL;
end
while (i2<=length(sfb))
s2b=[s2b,sfb(i2)<=0];
i2=i2+LL;
end
for i=1:n
if s2a(i)>s2b(i)
bit1=zeros(1,N);
else
bit1=ones(1,N);
end
bitb=[bitb,bit1];
end
figure(1);
subplot(K,1,4); plot(bitb,'LineWidth',1.5),grid on;title('解调后信号'); axis([0,length(bitb),-2.5,2.5]);
case 4,scolor='2DPSK';
%- 2DPSK 调制与解调
%--------------------------------------------------- %>>>>>>>>>>>>>>>>>>Initial_Part>>>>>>>>>>>>>>>>>>>>>
%---------------------------------------------------
fs = 30000;
Time_Hold_On = 0.1;
Num_Unit = fs * Time_Hold_On;
High_Level = ones ( 1, Num_Unit );
Low_Level = zeros ( 1, Num_Unit );
w = 300;
A = 1;
%--------------------------------------------------- %>>>>>>>>>>>>>>>>>>Initial_The_Signal>>>>>>>>>>>>>>>
%---------------------------------------------------
Sign_Set = [0,1,1,0,1,0,0,1]
Lenth_Of_Sign = length ( Sign_Set );
st = zeros ( 1, Num_Unit * Lenth_Of_Sign );
sign_orign = zeros ( 1, Num_Unit * Lenth_Of_Sign );
sign_result = zeros ( 1, Num_Unit * Lenth_Of_Sign );
t = 0 : 1/fs : Time_Hold_On * Lenth_Of_Sign - 1/fs;
%--------------------------------------------------- %>>>>>>>>>>>产生基带信号>>>>>>>>>>>>
%---------------------------------------------------
for I = 1 : Lenth_Of_Sign
if Sign_Set(I) == 1
sign_orign( (I-1)*Num_Unit + 1 : I*Num_Unit) = High_Level;
else
sign_orign( (I-1)*Num_Unit + 1 : I*Num_Unit) = Low_Level;
end
end
%--------------------------------------------------- %>>>>>>>>>>>>>>>>>>调制部分>>>>>>>>>>>>>>>>>>
%---------------------------------------------------
for I = 1 : Lenth_Of_Sign
if Sign_Set(I) == 1
st( (I-1)*Num_Unit + 1 : I*Num_Unit) = A * cos ( 2 * pi * w * t( (I-1)*Num_Unit + 1 : I*Num_Unit ) + ( pi / 2 ) );
else
st( (I-1)*Num_Unit + 1 : I*Num_Unit) = A * cos ( 2 * pi * w * t( (I-1)*Num_Unit + 1 : I*Num_Unit ) );
end
end
subplot ( 2,1,1 )
plot(t, sign_orign);
axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), - (A / 2), A + (A / 2) ] ); title ( '原始信号' );
grid
subplot ( 2, 1, 2 );
plot ( t, st );
axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), - 3*(A / 2), 3*(A / 2) ] ); title ( '调制后的信号' );
grid
%--------------------------------------------------- %>>>>>>>>>>>>>>>>>>相乘>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
%---------------------------------------------------
dt = st .* cos ( 2 * pi * w * t );
figure
plot ( t, dt );
axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), - 3*(A / 2), 3*(A / 2) ] ); title ( '相乘后的波形' );
grid
%--------------------------------------------------- %>>>>>>>>>>>>>>>>>>>低通滤波部分>>>>>>>>>>>>>>>>>>>>
%---------------------------------------------------
[N,Wn] = buttord( 2*pi*50, 2*pi*150,3,25,'s'); %临界频率采用角频率表示[b,a]=butter(N,Wn,'s');
[bz,az]=impinvar(b,a,fs); %映射为数字的
dt = filter(bz,az,dt);
figure
plot ( t, dt );
axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), - 3*(A / 2), 3*(A / 2) ] ); title ( '低通滤波后的波形' );
grid
%--------------------------------------------------- %>>>>>>>>>>>>>抽样判决 & 逆码变换部分>>>>>>>>>>>>>>>
%---------------------------------------------------
for I = 1 : Lenth_Of_Sign
if dt((2*I-1)*Num_Unit/2) < 0.25
sign_result( (I-1)*Num_Unit + 1 : I*Num_Unit) = High_Level;
sign_result( (I-1)*Num_Unit + 1 : I*Num_Unit) = Low_Level;
end
end
figure
plot ( t, sign_result );
axis( [ 0 , Time_Hold_On *( Lenth_Of_Sign + 1), - 3*(A / 2), 3*(A / 2) ] ); title ( '解调后信号' );
grid
end
四、测试数据
运行程序
命令窗口会显示随进产生的二进制码，对此二进制码进行HDB3编码之后得码元，解码之后得码元以及解码之后得码元与原基带码进行比较得结果如下图：
还会出现调制解调得选择窗口，来进行选择调制解调方式如下
在选择窗口中可以看到有四种调制解调方式
选择2ASK进行调制解调
可以看到随机得基带信号，2ASK调制后的信号以及2ASK解调后的信号进入菜单从新选择2PSK
可以看到随机产生得2PSK基带信号，2PSK调制后的信号以及2PSK解调后的信号
进入菜单栏重新选择2FSK调制解调
可以看到随机产生得基带信号，使用2FSK调制后的信号以及使用2FSK解调后的信号
进入菜单栏从新选择2DPSK调制解调
以上为随机产生得二进制原始信号以及2DPSK调制之后得信号
使用2DPSK解调之后的信号
2DPSK调制与解调系统的仿真 设计原理 (1) 2DPSK信号原理 1.1 2DPSK信号原理 2DPSK方式即是利用前后相邻码元的相对相位值去表示数字信息的一种方式。现假设用Φ表示本码元初相与前一码元初相之差，并规定：Φ＝0表示0码，Φ＝π表示1码。则数字信息序列与2DPSK信号的码元相位关系可举例表示如2PSK信号是用载波的不同相位直接去表示相应的数字信号而得出的，在接收端只能采用相干解调，它的时域波形图如图2.1所示。 图1.1 2DPSK信号 在这种绝对移相方式中，发送端是采用某一个相位作为基准，所以在系统接收端也必须采用相同的基准相位。如果基准相位发生变化，则在接收端回复的信号将与发送的数字信息完全相反。所以在实际过程中一般不采用绝对移相方式，而采用相对移相方式。定义为本码元初相与前一码元初相之差，假设： →数字信息“0”； →数字信息“1”。 则数字信息序列与2DPSK信号的码元相位关系可举例表示如下： 数字信息： 1 0 1 1 0 1 1 1 0 1 DPSK信号相位：0
或 ： 1.2 2DPSK 信号的调制原理 一般来说，2DPSK 信号有两种调试方法，即模拟调制法和键控法。2DPSK 信号的的模拟调制法框图如图1.2.1所示，其中码变换的过程为将输入的单极性不归零码转换为双极性不归零码。 图1.2.1 模拟调制法 2DPSK 信号的的键控调制法框图如图1.2.2所示，其中码变换的过程为将输入的基带信号差分，即变为它的相对码。选相开关作用为当输入为数字信息“0” 时接相位0，当输入数字信息为“1”时接pi 。 图1.2.2 键控法调制原理图 1.3 2DPSK 信号的解调原理 2DPSK 信号最常用的解调方法有两种，一种是极性比较和码变换法，另一种是差分相干解调法。 码变换 相乘 载波 s(t) e o (t)
通信原理课程设计 题目：脉冲编码调制（PCM）系统设计与仿真 院（系）：电气与信息工程学院 班级：电信04-6班 姓名：朱明录 学号： 0402020608 指导教师：赵金宪 教师职称：教授
摘要 : SystemView 仿真软件可以实现多层次的通信系统仿真。脉冲编码调制（PCM ）是现 代语音通信中数字化的重要编码方式。利用SystemView 实现脉冲编码调制(PCM)仿真，可以为硬件电路实现提供理论依据。通过仿真展示了PCM 编码实现的设计思路及具体过程，并加以进行分析。 关键词: PCM 编译码 1、引言 随着电子技术和计算机技术的发展，仿真技术得到了广泛的应用。基于信号的用于通信系统的动态仿真软件SystemView 具有强大的功能，可以满足从底层到高层不同层次的设计、分析使用，并且提供了嵌入式的模块分析方法，形成多层系统，使系统设计更加简洁明了，便于完成复杂系统的设计。 SystemView 具有良好的交互界面，通过分析窗口和示波器模拟等方法，提供了一个可视的仿真过程，不仅在工程上得到应用，在教学领域也得到认可，尤其在信号分析、通信系统等领域。其可以实现复杂的模拟、数字及数模混合电路及各种速率系统，并提供了内容丰富的基本库和专业库。 本文主要阐述了如何利用SystemView 实现脉冲编码调制（PCM ）。系统的实现通过模块分层实现，模块主要由PCM 编码模块、PCM 译码模块、及逻辑时钟控制信号构成。通过仿真设计电路，分析电路仿真结果，为最终硬件实现提供理论依据。 2、系统介绍 PCM 即脉冲编码调制，在通信系统中完成将语音信号数字化功能。PCM 的实现主要包括三个步骤完成：抽样、量化、编码。分别完成时间上离散、幅度上离散、及量化信号的二进制表示。根据CCITT 的建议，为改善小信号量化性能，采用压扩非均匀量化，有两种建议方式，分别为A 律和μ律方式，我国采用了A 律方式，由于A 律压缩实现复杂，常使用 13 折线法编码,采用非均匀量化PCM 编码示意图见图1。 图1 PCM 原理框图 下面将介绍PCM 编码中抽样、量化及编码的原理： (a) 抽样 所谓抽样，就是对模拟信号进行周期性扫描，把时间上连续的信号变成时间上离散的信号。该模拟信号经过抽样后还应当包含原信号中所有信息，也就是说能无失真的恢复原模拟信号。它的抽样速率的下限是由抽样定理确定的。 (b) 量化 从数学上来看，量化就是把一个连续幅度值的无限数集合映射成一个离散幅度值的有限数集合。如图2所示，量化器Q 输出L 个量化值k y ，k=1，2，3，…，L 。k y 常称为重建电
华东交通大学MATLAB程序设计报告书 课题名称：基于MATLAB的粒子群优化算法的实现 姓名： 学号：20160280800014 专业：控制科学与工程 2016年 11月 20日
基于MATLAB的粒子群优化算法的实现 一、课程选题目的 本次课程设计的课题为《基于MATLAB的粒子群优化算法的实现》，主要为学会运用MATLAB对实际算法编程，加深对粒子群优化算法的理解，并为今后熟练使用MA TLAB进行系统的分析仿真和设计奠定基础。数值计算分析可以帮助更深入地理解理论知识，并为将来使用MA TLAB进行各领域数值分析分析和实际应用打下基础。 此次课程主要是为了进一步熟悉对MATLAB软件的使用，以及学会利用MA TLAB对数值运算这种实际问题进行处理，将理论应用于实际，加深对它的理解。 二、粒子群优化算法原理 优化是科学研究、工程技术和经济管理等领域的重要研究工具。它所研究的问题是讨论在众多的方案中寻找最优方案。例如，工程设计中怎样选择设计参数，使设计方案既满足设计要求又能降低成本；资源分配中，怎样分配有限资源，使分配方案既能满足各方面的基本要求，又能获得好的经济效益。在人类活动的各个领域中，诸如此类，不胜枚举。优化这一技术，正是为这些问题的解决，提供理论基础和求解方法，它是一门应用广泛、实用性很强的科学。近十余年来，粒子群优化算法作为群体智能算法的一个重要分支得到了广泛深入的研究，在路径规划等许多领域都有应用。 2.1 粒子群优化算法的起源 粒子群优化（PSO）算法是由Kennedy和Eberhart于1995年用计算机模拟鸟群觅食这一简单的社会行为时，受到启发，简化之后而提出的。 设想这样一个场景：一群鸟随机的分布在一个区域中，在这个区域里只有一块食物。所有的鸟都不知道食物在哪里。但是他们知道当前的位置离食物还有多远。那么找到食物的最优策略是什么呢。最简单有效的方法就是追寻自己视野中目前离食物最近的鸟。如果把食物当作最优点，而把鸟离食物的距离当作函数的适应度，那么鸟寻觅食物的过程就可以当作一个函数寻优的过程。鱼群和鸟群的社会行为一直引起科学家的兴趣。他们以特殊的方式移动、同步，不会相互碰撞，整体行为看上去非常优美。生物学家CargiReynolds提出了一个非常有影响的鸟群聚集模型。在他的模拟模型boids中，每一个个体遵循：避免与邻域个体相冲撞、匹配邻域个体的速度、试图飞向感知到的鸟群中心这三条规则形成简单的非集中控制算法驱动鸟群的聚集，在一系列模拟实验中突现出了非常接近现实鸟群聚集行为的现象。该结果显示了在空中回旋的鸟组成轮廓清晰的群体，以及遇到障碍物时鸟群的分裂和再度汇合过程。由此受到启发，经过简化提出了粒子群优化算法。 2.2粒子群优化算法的原理 在粒子群优化算法中，每个优化问题的潜在解都是搜索空间中的一只鸟，称之为“粒子”。所有的粒子都有一个由被优化的函数决定的适应值，每个粒子还有一个速度决定他们飞翔的方向和距离。然后粒子们就追随当前的最优粒子在解空间中搜索。优化开始时先初始化为一群随机粒子（随机解）。然后通过迭代找到最优解。在每一次迭代中，粒子通过跟踪两个极值来更新自己。第一个极值就是整个种群目前找到的最优解。这个极值是全局极值。另外也可以不用整个种群而只是用其中一部分作为粒子的邻居，那么在所有邻居中的极值就是局部极值。第二个极值是粒子本身所找到的最优解，称为个体极值。这是因为粒子仅仅通过跟踪全局极值或者局部极值来更新位置，不可能总是获得较好的解。这样在优化过程中，粒子在追随全局极值或局部极值的同时追随个体极值则圆满的解决了这个问题。这就是粒子群优化
通信原理课程设计 ＿＿＿＿＿＿ 学号＿＿＿＿___ 班级＿＿＿＿＿
目录 一、目录 (2) 二、任务书 (3) 三、具体容及要求 (4) 3.1 题目一 (4) 3.1.1题目容 (4) 3.1.2设计思想或方法 (4) 3.1.3实现的功能或方法 (4) 3.1.4程序流程图 (4) 3.1.5程序代码 (5) 3.1.6仿真框图 (5) 3.1.7模块描述及参数设置 (5) 3.1.8结果运行…………………………………………………………… 10 3.1.9结果分析…………………………………………………………… 11 3.2 题目二………………………………………………………………… 11 3.2.1题目容 (11) 3.2.2设计思想或方法…………………………………………………… 11 3.2.2程序流程图 (12) 3.2.4程序代码…………………………………………………………… 13 3.2.5仿真框图…………………………………………………………… 13 3.2.6模块描述及参数设置………………………………………………… 14 3.2.7结果运行…………………………………………………………… 20 3.2.8结果分析…………………………………………………………… 20
3.3 题目三………………………………………………………………… 20 3.3.1题目容 (20) 3.3.2设计思想或方法 (20) 3.2.3程序流程图 (21) 3.2.4程序代码 (21) 3.2.5结果运行 (23) 3.2.6结果分析 (23) 四、心得与体会 (23) 五、参考文献 (23) 《通信原理课程设计》任务书 一、目的和要求： 要求学生在熟练掌握MATLAB和simulink仿真使用的基础上，学会通信仿真系统的基本设计与调试。并结合通信原理的知识，对通信仿真系统进行性能分析。 二、实验环境 PC机、Matlab/Simulink 三、具体容及要求 (1)试用Matlab/Simulink研究BPSK在加性高斯白噪声信道下的误码率性能与信 噪比之间的关系； (2)试用Matlab/Simulink研究BPSK+信道编码（取汉明码）在加性高斯白噪声信 道下的误码率性能与信噪比之间的关系；分析不同码率对误码率性能的影响。 (3)试用Matlab编程实现HDB3码的编解码过程，并画出1 1 0 1 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 1 0的原始、编码和解码图形。 四、提交设计报告
通信原理课程设计 --基于FPGA的时分多路数字基带传输系统的设计与开发 指导老师：戴慧洁武卫华 班级：通信111班 组长：徐震震 组员：胡彬、韦景山、谢留香、 徐勇、周晶晶、张秋红 日期：
一、课程设计目的 通信系统课程设计是一门综合设计性实践课程。使大家在综合已学现代通信系统理论知识的基础上，借助可编程逻辑器件及EDA技术的灵活性和可编程性，充分发挥自主创新意识，在规定时间内完成符合实际需求的通信系统电路设计与调试任务。 它不仅能够提高大家对所学理论知识的理解能力，更重要的是能够提高和挖掘大家对所学知识的实际运用能力，为将来进入社会从事相关工作奠定较好的“能力”基础。 二、课程设计内容 时分多路数字电话基带传输系统的设计与开发 三、课程设计要求任务 1、64Kb/S的A律PCM数字话音编译码器的开发设计 2、PCM 30/32一次群时分复接与分接器的开发设计 3、数字基带编码HDB3编译码器的开发设计 4、同步（帧、位、载波同步（可选））电路的开发设计
四、小组分工 小组成员负责项目 徐震震同步（帧同步、位同步） 谢留香PCM 30/32一次群时分复接 韦景山64Kb/S的A律PCM数字话音编码 胡彬PCM 30/32一次群时分分接 徐勇64Kb/S的A律PCM数字话音译码 周晶晶数字基带编码HDB3译码 张秋红数字基带编码HDB3编码 五、时分多路数字电话基带传输系统框图
PCM编码设计 一、设计要求 1、PCM编码器输入信号为: 一个13位逻辑矢量的均匀量化值:D0,D1…D12 其中:D0为极性位,取值范围在-4096～+4096之间； 一个占空比为1/32的8K/S的取样时钟信号； 一个占空比为50%的2.048Mb/S的合路时钟信号； 2、PCM编码器输出信号为: 一个8位逻辑矢量的13折线非均匀量化值:C0,C1…C7 其中:C0为极性位.C0=1为正,C0=0为负； 一个占空比为1/32的8K/S的取样时钟信号； 一个占空比为50%的2.048Mb/S的合路时钟信号； 二、PCM编码分析 脉冲编码调制(PCM)在通信系统中完成将语音信号数字化功能。是一种对模拟信号数字化的取样技术，将模拟信号变换为数字信号的编码方式，特别是对于音频信号。PCM 对信号每秒钟取样8000 次；每次取样为8个位，总共64kbps。PCM的实现主要包括三个步骤完成：抽样、量化、编码。分别完成时间上离散、幅度上离散、及量化信号的二进制表示。根据CCITT的建议，为改善小信号量化性能，采用压扩非均匀量化，有两种建议方式，分别为A 律和μ律方式，本设计采用了A律方式。 在13折线法中，无论输入信号是正是负，均按8段折线（8个段落）进行编码。若用8位折叠二进制码来表示输入信号的抽样量化值，其中用第一位表示量化值的极性，其余七位（第二位至第八位）则表示抽样量化值的绝对大小。具体的做法是：用第二至第四位表示
课程设计任务书 学生姓名：董航专业班级：电信1006班 指导教师：阙大顺，李景松工作单位：信息工程学院 课程设计名称：Matlab应用课程设计 课程设计题目：Matlab运算与应用设计5 初始条件： 1.Matlab6.5以上版本软件； 2.课程设计辅导资料：“Matlab语言基础及使用入门”、“Matlab及在电子信息课程中的应 用”、线性代数及相关书籍等； 3.先修课程：高等数学、线性代数、电路、Matlab应用实践及信号处理类相关课程等。 要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求） 1.课程设计内容：根据指导老师给定的7套题目，按规定选择其中1套完成； 2.本课程设计统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析， 针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结。具体设计要求包括： ①初步了解Matlab、熟悉Matlab界面、进行简单操作； ②MATLAB的数值计算：创建矩阵矩阵运算、多项式运算、线性方程组、数值统计； ③基本绘图函数：plot, plot3, mesh, surf等，要求掌握以上绘图函数的用法、简单图形 标注、简单颜色设定等； ④使用文本编辑器编辑m文件，函数调用； ⑤能进行简单的信号处理Matlab编程； ⑥按要求参加课程设计实验演示和答辩等。 3.课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括： ①目录； ②与设计题目相关的理论分析、归纳和总结； ③与设计内容相关的原理分析、建模、推导、可行性分析； ④程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结； ⑤课程设计的心得体会（至少500字）； ⑥参考文献（不少于5篇）； ⑦其它必要内容等。 时间安排：1.5周（分散进行） 参考文献： [1](美)穆尔,高会生,刘童娜,李聪聪．MA TLAB实用教程(第二版) . 电子工业出版社,2010. [2]王正林,刘明．精通MATLAB(升级版) .电子工业出版社,2011. [3]陈杰. MA TLAB宝典(第3版) . 电子工业出版社,2011. [4]刘保柱,苏彦华,张宏林. MATLAB 7.0从入门到精通(修订版) . 人民邮电出版社,2010. 指导教师签名：年月日 系主任（或责任教师）签名：年月日
通信原理课程设计报告 题目：基于MATLAB 的M-QAM调 制及相干解调的设计与仿真班级：通信工程1411 姓名：杨仕浩(2014111347） 解博文(2014111321） 介子豪(2014111322） 指导老师：罗倩倩 成绩： 日期：2016 年12 月21 日
基于MATLAB的M-QAM调制及相干解调的设计与仿真 摘要：正交幅度调制技术（QAM）是一种功率和带宽相对高效的信道调制技术，因此在自适应信道调制技术中得到了较多应用。本次课程设计主要运用MATLAB软件对M =16 进制正交幅度调制系统进行了仿真，从理论上验证16进制正交幅度调制系统工作原理，为实际应用和科学合理地设计正交幅度调制系统,提供了便捷、高效、直观的重要方法。实验及仿真的结果证明，多进制正交幅度调制解调易于实现，且性能良好，是未来通信技术的主要研究方向之一，并有广阔的应用前景。 关键词：正交幅度调制系统；MATLAB；仿真
目录 1引言 (1) 1.1课程设计的目的 (1) 1.2课程设计的基本任务和要求 (1) 1.3仿真平台Matlab (1) 2 QAM系统的介绍 (2) 2.1正交幅度调制技术 (2) 2.2QAM调制解调原理 (5) 2.3QAM的误码率性能 (7) 3 多进制正交幅度（M-QAM）调制及相干解调原理框图 (9) 4 基于MATLAB的多进制正交幅度（M-QAM）调制及相干解调设计与仿真 (10) 4.1系统设计 (10) 4.2随机信号的生成 (10) 4.3星座图映射 (11) 4.4波形成形（平方根升余弦滤波器） (13) 4.5调制 (14) 4.6加入高斯白噪声之后解调 (15) 5 仿真结果及分析 (20) 6 总结与体会 (23) 6.1总结 (23) 6.2心得体会 (24) 【参考文献】 (25) 附录 (26)
通信原理实验总结 ——电信0906 经过三次的通信原理课程相关实验的学习，我们班同学都受益菲浅。在此仅代表我们班，我对这三次的实验进行了总结。 首先对三次实验中所学到的知识点做个总结。我们一共做了10个实验：数字基带信号，数字调制，模拟锁相环与载波同步，数字解调与眼图，数字锁相环与位同步，帧同步，时分复用数字基带信号，时分复用2DPSK、2FSK通信系统，PCM编译码，时分复用通话与抽样定理。通过这10个实验，我们在理论和实际应用方面都有了一定的提高，比如在理论方面，通过数字基带信号的实验，我们了解了单极性码、双极性码、归零码、不归零码等基带信号波形特点并掌握AMI、HDB3码的编码规则；通过数字调制的实验，我们掌握绝对码、相对码概念及它们之间的变换关系，掌握了相对波形与2PSK信号波形之间的关系、绝对码波形与2DPSK信号波形之间的关系，并对2ASK、2FSK、2DPSK信号的频谱与数字基带信号频谱之间的关系有了进一步了解。又比如在实际应用方面，通过模拟锁相环与载波同步的实验，我们掌握了用平方环法从2DPSK信号中提取相干载波的原理及模拟锁相环的设计方法；通过PCM编译码，我们掌握了PCM基带信号的形成过程及分接过程，并学会了语音信号PCM 编译码系统的动态范围和频率特性的定义及测量方法。 其次在通信原理实验的学习中，我们班同学在学习方法上也受益颇多。总结如下：一、强化了课前预习的好习惯。大多同学觉得自己一直以来就没能养成课前预习的好习惯（虽然一直认为课前预习是很重要的），但通过这三次实验懂得了课前预习的重要。因为只有在课前充分了解了实验原理，才能在课上更好的学习，收获的更多、掌握的更多。二、培养了勤于动手能力。“实验就是为了让你动手做，去探索一些你未知的或是你尚不是深刻理解的东西。”由于10个实验每个步骤大家都必须亲自去做，亲自去调试，同学们都认为动手能力得到了提高。三、懂得了在探索中求得真知。那些伟大的科学家之所以伟大就是他们利用实验证明了他们的伟大。实验是检验理论正确与否的试金石。为了要使你的理论被人接受，你必须用事实（实验）来证明，让那些怀疑的人哑口无言。虽说我们的通信原理实验基本上都是验证性实验，只是对前人的经典实验的重复，但是对于一个知识尚浅、探索能力还不够的人来说，这些探索也非一件易事。对于这些实验，同学们在探索中学习、在模仿中理解、在实践中掌握。通信原理实验让我慢慢开始“摸着石头过河”。学习就是为了能自我学习，这正是实验课的核心，它让我在探索、自我学习中获得知识。 另外，个人方面，经过这一周的通信原理实验课的学习，让我收获多多。但在这中间，我也发现了自身存在的很多不足。我的动手能力还不够强，当有些实验需要很强的动手能力时我还不能从容应对，往往需要组员的帮助；我的探索方式还有待改善，当面对一些复杂的实验时我还不能很快很好的完成，同时面对每个实验后相关的思考题，我有时也会显得非常棘手。 最后，我们班同学在做完实验后也普遍认为，10个实验均为验证性试验，虽然对知识的理解会有很大帮助，但缺乏创新性。因此建议学校能在通信原理课上适当增加一点创新性实验或让同学在完成一定数量的实验且拥有一定的知识积淀后能有所创新。另外，不少同学也表示，他们是在实验结束后，通过写实验报告才对实验原理有了比较清晰的了解，所以大家也建议学校可以先让学生写预习报告，在对实验原理及过程充分了解后，再完成实验。
《计算机仿真及应用》课程设计报告书 学号：08057102，08057127 班级：自动化081 姓名陈婷，万嘉
目录 一、设计思想 二、设计步骤 三、调试过程 四、结果分析 五、心得体会 六、参考文献
选题一、 考虑如下图所示的电机拖动控制系统模型，该系统有双输入，给定输入)(t R 和负载输入)(t M 。 1、 编制MATLAB 程序推导出该系统的传递函数矩阵。 2、 若常系数增益为：C 1＝Ka ＝Km ＝1，Kr ＝3，C2＝0.8，Kb ＝1.5，时间常数T 1＝5， T 2＝0.5，绘制该系统的根轨迹、求出闭环零极点，分析系统的稳定性。若)(t R 和)(t M 分别为单位阶跃输入，绘制出该系统的阶跃响应图。（要求C 1，Ka ，Km ，Kr ，C2，Kb ， T 1，T 2所有参数都是可调的） 一．设计思想 题目分析： 系统为双输入单输出系统，采用分开计算，再叠加。 要求参数均为可调，而matlb 中不能计算未赋值的函数，那么我们可以把参数设置为可输入变量，运行期间根据要求赋值。 设计思路： 使用append 命令连接系统框图。 选择‘参数=input('inputanumber:')’实现参数可调。 采用的方案： 将结构框图每条支路稍作简化，建立各条支路连接关系构造函数，运行得出相应的传递函数。 在得出传递函数的基础上，使用相应的指令求出系统闭环零极点、画出其根轨迹。 通过判断极点是否在左半平面来编程判断其系统是否稳定。 二．设计步骤 （1）将各模块的通路排序编号
（2）使用append命令实现各模块未连接的系统矩阵 （3）指定连接关系 （4）使用connect命令构造整个系统的模型 三．调试过程 出现问题分析及解决办法： 在调试过程出现很多平时不注意且不易寻找的问题，例如输入的逗号和分号在系统运行时不支持中文格式，这时需要将其全部换成英文格式，此类的程序错误需要细心。 在实现参数可调时初始是将其设为常量，再将其赋值进行系统运行，这样参数可调性差，后用‘参数=input('inputanumber:')’实现。 最后是在建立通路连接关系时需要细心。 四．结果分析 源代码： Syms C1 C2 Ka Kr Km Kb T1 T2 C1=input('inputanumber:') C2=input('inputanumber:') Ka=input('inputanumber:') Kr=input('inputanumber:') Km=input('inputanumber:') Kb=input('inputanumber:') T1=input('inputanumber:') T2=input('inputanumber:') G1=tf(C1,[0 1]); G2=tf(Ka*Kr,[0 1]); G3=tf(Km,[T1 1]); G4=tf(1,[T2 1]); G5=tf(1,[1 0]); G6=tf(-C2,1); G7=tf(-Kb,1); G8=tf(-1,1); Sys=append(G1,G2,G3,G4,G5,G6,G7,G8) Q=[1 0 0;2 1 6;3 2 7;4 3 8;5 4 0;6 5 0;7 4 0;8 0 0;]; INPUTS1=1; OUTPUTS=5; Ga=connect(Sys,Q,INPUTS1,OUTPUTS) INPUTS2=8; OUTPUTS=5; Gb=connect(Sys,Q,INPUTS2,OUTPUTS) rlocus(Ga)
目录 第1章绪论 (1) 第2章 SystemView的基本介绍 (2) 第3章二进制振幅键控 2ASK (4) 3.1 2ASK调制系统 (4) 3.2 2ASK调制解调系统 (6) 3.3 2ASK系统仿真结果分析 (9) 第四章二进制频移键控 2FSK (10) 4.1 2FSK调制系统 (10) 4.2 2FSK调制解调系统 (12) 4.3 2FSK仿真结果分析 (17) 第5章二进制移相键控 2PSK (18) 5.1 2PSK调制系统 (18) 5.2 2PSK调制解调系统 (19) 5.3 2PSK仿真结果分析 (23) 第6章二进制差分移相键控 2DPSK (24) 6.1 2DPSK实验原理 (24) 6.2 2DPSK仿真结果分析 (29) 第7章实验总结 (30) 第8章参考文献 (30) 第9章谢辞 (32)
第1章绪论 通信按照传统的理解就是信息的传输，信息的传输离不开它的传输工具，通信系统应运而生，我们此次课题的目的就是要对调制解调的通信系统进行仿真研究。 数字信号的传输方式可以分为基带传输和带通传输。为了使信号在带通信道中传输，必须用数字基带信号对载波进行调制，以使信号与信道特性相匹配。在这个过程中就要用到数字调制。 在通信系统中，利用数字信号的离散取值特点通过开关键控载波，来实现数字调制，这种方法通常称为键控法，主要对载波的振幅，频率，和相位进行键控。键控主要分为：振幅键控，频移键控，相移键控三种基本的数字调制方式。 本次课程设计的目的是在学习以上三种调制的基础上，通过Systemview仿真软件，实现对2ASK,2FSK,2PSK,2DPSK等数字调制系统的仿真，同时对以上系统有深入的了解。 Systemview是美国ELANIX公司于1995年开始推出的软件工具，它为用户提供了一个完整的动态系统设计、仿真与分析的可视化软件环境，能进行模拟、数字、数模混合系统、线性和非线性系统的分析设计，可对线性系统进行拉氏变换和Z变换分析。 SystemView基本属于一个系统级工具平台，可进行包括数字信号处理（DSP）系统、模拟与数字通信系统、信号处理系统和控制系统的仿真分析，并配置了大量图符块（Token）库，用户很容易构造出所需要的仿真系统，只要调出有关图符块并设置好参数，完成图符块间的连线后运行仿真操作，最终以时域波形、眼图、功率谱、星座图和各类曲线形式给出系统的仿真分析结果。 在此次课程设计之前，先学会熟练掌握Systemview的用法，在该软件的配合下完成各个系统的结构图，还有调试结果图。 Systemview对系统的分析主要分为两大块，调制系统的分析和解调系统的分析。由于调制是解调的基础，没有调制就不可能有解调，为了表现解调系统往往需要很高的采样频率来减少滤波带来的解调失真，所以调制的已调信号通过波形模块观察起来不是很清楚，为了更好的弄清楚调制是怎么样的一个过程，在这里，我们把调制单独列出来，用较低的频率实现它，就能从单个周期上观察调制系统的运作模式，更深刻地表现调制系统的调制过程。
通信原理课程设计心得体会 、时分解复用原理 为了提高信道利用率，使多路已抽样的信号组合起来沿同一信道传输而互相不干扰，称时分多路复用。时分复用的解调过程称为时分解复用。目前采用较多的是频分多路解复用和时分多路解复用。频分多路解复用用于模拟通信，而时分多路解复用用于数字通信。为了实现TDM传输，要把传输时间分成若干个时隙，在每个时隙内传输一路信号，将若干个原始的脉冲调制信号在时间上进行交错排列，从而形成一个复合脉冲串，该脉冲串扰码后经信道传输到达接收端。时分解复用通信，是把各路信号在同一信道上占有不同时间间隙进行通信分离出原来的模拟信号。由抽样定理可知，将时间上离散的信号变成时间上连续的信号，其在信道上占用时间的有限性，为多路信号沿同一信道传输提供了条件。时分解复用是建立在抽样定理的基础上的,因为抽样定理连续的基带信号由可能被在时间上离散出现的抽样脉冲所代替.具体说，就是把时间分成一些均匀的时间间隙，将各路信号的传输时间分配在不同的时间间隙，以达到互相分开，互不干扰的目的。抽样脉冲占据时间一般较短,在抽样脉冲之间就留出间隙.利用这些空隙便可以传输其他信号的抽样,因此,就可能用一条信道同时传送若干个基带信号,并且每一个抽
样值占用的时间越短,能够传输的数据也就越多.时分解复用信号在接收端只要在时间上恰当地进行分离,各个信号就能分别互相分开，互不干扰并不失真地还原出原来的模拟信号。 在通信系统中，同步具有相当重要的地位。通信系统能否具有有效、可靠地工作，在很大程度上依赖有无良好的同步系统。同步可分为载波同步、位同步、帧同步和网同步几大类型。他们在通信系统中都具有相当重要的作用。时分解复用通信中的同步技术包括位同步和帧同步，这是数字通信的又一个重要特点。时分解复用的电路原理就是先通过帧同步信号和位同步信号把各路信号数据分开,然后通过移位寄存器构成的并/串转换电路输出串行的数据，把时分复用的调制信号不失真的分离出来。 位同步 位同步的目的是确定数字通信中的个码元的抽样时刻，即把每个码元加以区分，使接受端得到一连串的码元序列，这一连串的码元列代表一定的信息。位同步是最基本的同步，是实现帧同步的前提。位同步的基本含义是收、发两端机的时钟频率必须同频、同相，这样接收端才能正确接收和判决发送端送来的每一个码元。因此,接收端必须提供一个确定抽样判决时刻的定时脉冲序列.
matlab音频降噪课程设计报告
燕山大学 医学软件课程设计说明书 题目：基于MATLAB巴特沃斯滤波器的音频去噪的GUI设计 学院（系）：电气工程学院 年级专业： 13级生物医学工程 2 班 学号： 130103040041 学生姓名：魏鑫 指导教师：许全盛 1
院（系）：电气工程学院基层教学单位：生物医学工程系 学号130103040041 学生 姓名 魏鑫 专业（班 级） 13级生 物医学 工程2 班 设计 题目 基于MATLAB音频去噪的GUI设计设 计 技术参数通带截止频率fp=2700;阻带截止频率fs=3000;采样频率FS=48000; 通带衰减不大于1dB;阻带衰减不小于10dB; 设计要求1.实现用MATLAB导入音频； 2.对音频进行频谱分析； 3.设计滤波器去噪并对含噪信号进行滤 2
波并进行功率谱分析； 4.设计能实现上述功能的GUI； 工作量1.完成音频录入及频谱分析相关程序的编写与调试； 2.设计滤波器去噪； 3.用MATLAB软件做GUI界面的设计； 工作计划11.21-11.24 MATLAB软件中GUIDE 工具箱的使用 11.25-11.29 各处理算法模块的编程实现 11.30-12.1 整体程序联调 12.2 撰写课程设计说明书,答辩 参考资料 1. 陈怀琛吴大正 MATLAB及在电子信息课程中的应用[M] 北京电子工业出版社 2006. 章节2.4； 2. 陈亚勇 MATLAB信号处理详解[M] 北京：人民邮电出版社 2000. 第十 3
章； 3.张康刘雅基于Matlab的巴特沃斯 数字低通滤波器的设计[J] 计算机与现代化 2007年 12期 98-100页 指导 教师签字许全盛 基层教学单 位主任签字 彭勇 目录 一、设计目的意义 (1) 1.1绪论 (1) 1.2设计目的 (1) 1.3意义 (1) 二、设计内容 (2) 2.1 设计原理 (2) 2.2 设计内容 (2) 三、设计过程及结果分析 (3) 3.1 设计步骤 (3) 4
信息工程学院 2014 / 2015学年第一学期 课程设计报告 课程名称：通信原理课程设计 专业班级：统本电信1201 学生学号：12610304152213 12520527151362 学生姓名：陈钰康 夏涛 指导教师：田亚楠
摘要 8PSK(8 Phase Shift Keying，8移相键控)是八进制相移键控，它是一种相位调制算法。相位调制（调相）是频率调制（调频）的一种演变，载波的相位被调整用于把数字信息的比特编码到每一词相位改变（相移）。 8PSK中的“PSK”表示使用移相键控方式，移相键控是调相的一种形式，用于表达一系列离散的状态，8PSK对应8种状态的PSK。如果是其一半的状态，即4种，则为QPSK，如果是其2倍的状态，则为16PSK。因为8PSK拥有8种状态，所以8PSK每个符号(symbol)可以编码3个比特(bits)。8PSK抗链路恶化的能力（抗噪能力）不如QPSK，但提供了更高的数据吞吐容量。本次课程设计过程中，利用了MATLAB7.1仿真实现了8PSK信号的调制与解调，并仿真8PSK载波调制信号在高斯白噪声信道下的误码率及误比特率性能，并用MATLAB仿真出了调制信号、载波信号及已调信号的波形图和频谱图。并在高斯白噪声下，讨论了8PSK 误码率及误比特率性能。 关键字：8PSK;载波的调制；解调;
目录 一．设计内容及要求（PSK信号的仿真） (1) 二．相关理论知识的论述分析 (1) 2. 1.1、8PSK的概念 (1) 2. 1.2、8PSK的特点 (1) 2.2.1、 PSK的调制 (2) 2.2.2、调制的概念 (2) 2.2.3、调制的种类 (2) 2.2.4、调制的作用 (3) 2.2.5、调制方式 (3) 三．系统原理框图及分析(8PSK的原理) (3) 四．完整的设计仿真过程 (4) 五．仿真结果输出及结论 (6) 六．仿真调试中出现的错误、原因及排除方法 (7) 七．总结本次设计，指出设计的核心及应用价值，提出改进意见和展望 (7) 八．收获、体会 (7) 九．参考文献 (8)
实验一 设计任务： 用MatLib仿真一个BFSK通信系统，基本参数： 1）fc=1000Hz； 2）Rb=100bps； 3）信息序列:“Hello world”的ASCII 实验与报告基本要求： 1）Matlab程序，要点旁注（可打印后手写）； 2）绘出信号波形，绘出信号PSD； 3）给出解调后的信息序列； 4）将信息重复3遍以上,FSK信号保存为WAV文件格式，使用音频播放，聆听；M文件： wave.m function[t,mt]=wave(m,dt,fs) l=length(m); mt=[]; ddt=1/fs; n=floor(dt*fs); m_add=ones(1,n); for i=1:l if(m(i)) mt=[mt,m(i),m_add]; else mt=[mt,m(i),m_add*0]; end t=(1:((n+1)*l))*ddt; end my_filter.m function[num,den]=my_filter(wp,ws,ap,as) if nargin<4 as=15; end if nargin<4 ap=3; end [n,wn]=buttord(wp,ws,ap,as); [num,den]=butter(n,wn); end 代码：
f0=800;%‘0’码载波频率 f1=1200;%‘1’码载波频率 fs=4000;%采样频率 Rb=100;%比特率 dt=1/Rb;%一个比特发送时间 A0=2;%调制幅度 A1=2;%相干解调幅度 miu=0;sigma=0.3;%miu:高斯白噪声均值，sigma：高斯白噪声均方差 str='Hello world';%信号字符串 m_dec=abs(str);%将信号字符串转换成ASCII码(十进制) m_bin=dec2bin(m_dec,8); m_bin=abs(m_bin)-48;%将十进制转换成8比特二进制矩阵 m=[]; for i=1:size(m_bin,1) m=[m,m_bin(i,:)]; end%将二进制转换成行向量 [t,m]=wave(m,dt,fs);%对信号采样 mt_f1=m.*cos(2*pi*f1*t)*A0;%频率f1调制 mt_f0=(~m).*cos(2*pi*f0*t)*A0;%频率f0调制 mt=mt_f1+mt_f0;%发送信号 l=length(mt); subplot(2,1,1);plot(t,mt); grid on;xlabel('t/s');title('m(t)');%发送信号波形subplot(2,1,2);periodogram(mt,[],l,fs);grid on;%发送信号PSD
西南科技大学通信原理设计报告 课程名称：通信原理课程设计 设计名称：奇偶校验编码仿真 姓名：王雷 学号: 班级：通信1004 指导教师：秦明伟 起止日期：2013年7月5日星期五 西南科技大学信息工程学院制
方向设计任务书 学生班级：通信1004 学生姓名：王雷学号： 设计名称：奇偶校验编码仿真 起止日期：2013年7月5日星期五指导教师：秦明伟 方向设计学生日志
奇偶校验编码仿真 一、摘要（150-250字） 奇偶校验是一种校验代码传输正确性的方法。根据被传输的一组二进制代码的数位中“1”的个数是奇数或偶数来进行校验。采用奇数的称为奇校验，反之，称为偶校验。采用何种校验是事先规定好的。通常专门设置一个奇偶校验位，用它使这组代码中“1”的个数为奇数或偶数。若用奇校验，则当接收端收到这组代码时，校验“1”的个数是否为奇数，从而确定传输代码的正确性。 二、设计目的和意义 认识matlab软件，学习掌握matlab的基本操作方法，熟悉M文件和simulink的具体实现方法，了解数据奇偶校验的原理和在matlab中的基本仿真，通过对简单的通信实验设计，提高了动手能力和对matlab操作，巩固了课程知识。 三、设计原理 在数据传输前附加一位奇校验位，用来表示传输的数据中"1"的个数是奇数还是偶数，为奇数时，校验位置为"0"，否则置为"1"，用以保持数据的奇偶性不变。例如，需要传输"11001110"，数据中含5个"1"，所以其奇校验位为"0"，同时把"110011100"传输给接收方，接收方收到数据后再一次计算奇偶性，"110011100"中仍然含有5个"1"，所以接收方计算出的奇校验位还是"0"，与发送方一致，表示在此次传输过程中未发生错误。奇偶校验就是接收方用来验证发送方在传输过程中所传数据是否由于某些原因造成破坏。 奇偶校验原理是基于异或的逻辑功能。奇偶校验的编码方法是在原信号码组后面添加以为监督码元，奇偶校验分为奇校验和偶校验，奇校验是原信息码元加上监督码元后，使整个组成的数码组中，1的个数为奇数个。偶校验的工作原理则正好与奇校验相反。 对于n位二进码a1a2a3a4……a n奇校验有如下表示： a1⊕a2⊕a3⊕a4……⊕a n⊕C=1 偶校验的表达式为： a1⊕a2⊕a3⊕a4……⊕a n⊕C =1 其中，C为监督码元，在本设计中n为8，可以推出C的表达式为： C =a1⊕a2⊕a3⊕a4……⊕a8 在发送端让其监督码和信息码一起发送，在信息接收端，计算校验因子的表达式为： 、 S=a1⊕a2⊕a3⊕a4……⊕a n⊕C
Matlab程序设计任务书 分院（系）信息科学与工程专业 学生姓名学号 设计题目车牌识别系统设计 内容及要求： 车牌定位系统的目的在于正确获取整个图像中车牌的区域，并识别出车牌号。通过设计实现车牌识别系统，能够提高学生 分析问题和解决问题的能力，还能培养一定的科研能力。 1.牌照识别系统应包括车辆检测、图像采集、牌照识别等几 部分。 2.当车辆检测部分检测到车辆到达时,触发图像采集单元，采 集当前的视频图像。 3.牌照识别单元对图像进行处理，定位出牌照位置，再将牌 照中的字符分割出来进行识别，然后组成牌照号码输出。 进度安排： 19周：Matlab环境熟悉与基础知识学习 19周：课程设计选题与题目分析 20周：程序设计编程实现 20周：课程设计验收与答辩 指导教师（签字）： 年月日学院院长（签字）： 年月日 目录
一．课程设计目的 (3) 二．设计原理 (3) 三．详细设计步骤 (3) 四. 设计结果及分析 (18) 五. 总结 (19) 六. 设计体会 (20) 七. 参考文献 (21) 一、课程设计目的 车牌定位系统的目的在于正确获取整个图像中车牌的区域，并识别出车牌号。通过
设计实现车牌识别系统，能够提高学生分析问题和解决问题的能力，还能培养一定的科研能力。 二、设计原理: 牌照自动识别是一项利用车辆的动态视频或静态图像进行牌照号码、牌照颜色自动识别的模式识别技术。其硬件基础一般包括触发设备、摄像设备、照明设备、图像采集设备、识别车牌号码的处理机等，其软件核心包括车牌定位算法、车牌字符分割算法和光学字符识别算法等。某些牌照识别系统还具有通过视频图像判断车辆驶入视野的功能称之为视频车辆检测。一个完整的牌照识别系统应包括车辆检测、图像采集、牌照识别等几部分。当车辆检测部分检测到车辆到达时触发图像采集单元，采集当前的视频图像。牌照识别单元对图像进行处理，定位出牌照位置，再将牌照中的字符分割出来进行识别，然后组成牌照号码输出。 三、详细设计步骤: 1. 提出总体设计方案: 牌照号码、颜色识别 为了进行牌照识别，需要以下几个基本的步骤： a.牌照定位，定位图片中的牌照位置；
通信原理实验报告 一、实验名称 MATLAB验证低通抽样定理 二、实验目的 1、掌握抽样定理的工作原理。 2、通过MATLAB编程实现对抽样定理的验证，加深抽样定理的理解。同时训练应用计算机分析问题的能力。 3、了解MATLAB软件，学习应用MATLAB软件的仿真技术。它主要侧重于某些理论知识的灵活运用，以及一些关键命令的掌握，理解，分析等。 4、计算在临界采样、过采样、欠采样三种不同条件下恢复信号的误差，并由此总结采样频率对信号恢复产生误差的影响，从而验证时域采样定理。 三、实验步骤及原理 1、对连续信号进行等间隔采样形成采样信号，采样信号的频谱是原连续信号的频谱以采样频率为周期进行周期性的延拓形成的。 2、设连续信号的的最高频率为Fmax，如果采样频率Fs>2Fmax，那么采样信号可以唯一的恢复出原连续信号，否则Fs<=2Fmax会造成采样信号中的频谱混叠现象，不可能无失真地恢复原连续信号。 四、实验内容 1、画出连续时间信号的时域波形及其幅频特性曲线，信号为 x=cos(4*pi*t)+1.5*sin(6*pi*t)+0.5*cos(20*pi*t)
2、对信号进行采样，得到采样序列，画出采样频率分别为10Hz，20 Hz，50 Hz时的采样序列波形； 3、对不同采样频率下的采样序列进行频谱分析，绘制其幅频曲线，对比各频率下采样序列和的幅频曲线有无差别。 4、对信号进行谱分析，观察与3中结果有无差别。 5、由采样序列恢复出连续时间信号，画出其时域波形，对比与原连续时间信号的时域波形。 五、实验仿真图 (1) x=cos(4*pi*t)+1.5*sin(6*pi*t)+0.5*cos(20*pi*t)的时域波 形及幅频特性曲线。 clear; close all; dt=0.05; t=-2:dt:2 x=cos(4*pi*t)+1.5*sin(6*pi*t)+0.5*cos(20*pi*t); N=length(t); Y=fft(x)/N*2; fs=1/dt; df=fs/(N-1); f=(0:N-1)*df; subplot(2,1,1) plot(t,x)