基于Matlab的数字信号传输系统实验文档格式.docx

基于Matlab的数字信号传输系统实验文档格式.docx

《基于Matlab的数字信号传输系统实验文档格式.docx》由会员分享，可在线阅读，更多相关《基于Matlab的数字信号传输系统实验文档格式.docx（3页珍藏版）》请在冰点文库上搜索。

程序源代码

Fs=1e4;

%采样频率

len=20;

%码元长度

in=randint（1,len,4）;

%产生初始码元序列

sig=[];

out=[];

fort=1:

2000 %产生基带信号

n=fix（t/100）;

ifn==0

in_a（t）=0;

else

in_a（t）=in（n）;

end

end

subplot（2,1,1）;

%基带信号

plot（in_a,'

LineWidth'

3）;

title（'

基带信号'

'

FontWeight'

bold'

FontSize'

20）;

xlabel（'

t/s'

18）;

axis（[100,2100,-0.5,3.5]）;

set（gca,'

XTick'

0:

100:

2000）;

gridon;

cxn=xcorr（in_a,'

unbiased'

）;

%%计算序列的自相关函数

nfft=1024;

CXk=fft（cxn,nfft）;

Pxx=abs（CXk）;

index=0:

round（nfft/2-1）;

k=index*Fs/nfft;

subplot（2,1,2）;

plot_Pxx=10*log10（Pxx（index+1））;

plot（k,plot_Pxx,'

2）;

基带信号功率谱'

axis（[0,5000,-10,40]）;

Hz'

18,'

fori=1:

len%产生单极性归零码信号

ifin（i）==0

ins=[0,0];

elseifin（i）==1

ins=[1,0];

elseifin（i）==2

ins=[2,0];

ins=[3,0];

end

sig=[sig,ins];

4000

s（t）=0;

s（t）=sig（n）;

figure;

%单极性归零码

plot（s,'

单极性归零码'

axis（[100,4100,-0.5,3.5]）;

200:

4100）;

cxn=xcorr（s,'

nfft=1024;

单极性归零码功率谱'

s1=awgn（s,20）;

%添加噪声

plot（s1）;

添加噪声后的信号'

500:

cxn=xcorr（s1,'

添加噪声后的信号功率谱'

%滤波器设计

fp=500;

%通带截止

fs=550;

%阻带截止

ws=fs*2/Fs;

wp=fp*2/Fs;

[N,Wp]=ellipord（wp,ws,1,40）;

[b,a]=ellip（N,1,40,Wp）;

sf0=filter（b,a,s1）;

%滤掉部分噪声后的信号

plot（sf0）;

滤掉部分噪声后的信号'

cxn=xcorr（sf0,'

滤掉部分噪声后的信号功率谱'

form=1:

20%抽样判决

p=round（sf0（200*m-20））;

out=[out,p];

%原始码元序列

stairs（in,'

原始码元序列'

axis（[1,21,-0.5,3.5]）;

1:

20）;

%抽样判决后恢复的信号序列

stairs（out,'

抽样判决后恢复的信号序列'