- 1、本文档共2页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
functionber=sim_ASK2(snrIndB,N,flag)
%二进制ASK调制解调仿真程序
%输入
%snr:信噪比
%N:仿真比特数
%flag:绘制调制解调波形标志位
%输出
%y:误码率
x=round(rand(1,N));%产生信号比特
fc=200;%载波频率
fs=2000;%采样频率
fb=50;%比特频率
nSamples=fs/fb;%每比特采样数
signal_in_one_bit=cos(2*pi*fc*(0:nSamples-1)/fs);%传号
space_in_one_bit=zeros(1,nSamples);%空号
limitation=sum(signal_in_one_bit.^2)/2;%判决门限
sig_ori=reshape(repmat(x,nSamples,1),1,N*nSamples);%原始信号
%%ASK调制
y=zeros(1,nSamples*N);
fori=1:length(x)
if(x(i)==1)
y((i-1)*nSamples+1:i*nSamples)=signal_in_one_bit;
else
y((i-1)*nSamples+1:i*nSamples)=space_in_one_bit;
end
end
mod_sig=y;%调制信号
%%高斯信号传输
%rec_sig=awgn(mod_sig,snr,0);%通过高斯信道后的接收信号
snr=10^(0.1*snrIndB)/(fb/4);sigma=1/sqrt(snr);
rec_sig=awgn_gauss(mod_sig,sigma);
%%ASK解调
bits_judged=zeros(1,N);
fori=1:N%逐比特处理
rec_seg=rec_sig((i-1)*nSamples+1:(i*nSamples));
sum_rec=rec_seg*signal_in_one_bit;%相关累积
if(sum_rec=limitation)%判决
bits_judged(i)=1;
else
bits_judged(i)=0;
end
end
x_pie=bits_judged;%接收信号判决值
ber=length(find(x~=x_pie))/N;%误码率
%绘制调制解调结果图
if(1==flag)
bits_display=30;%绘图显示的比特数
figure;
subplot(411);stem(x(1:bits_display));title(原始数据比特),xlabel(t)
subplot(412);plot(mod_sig(1:bits_display*nSamples));title(调制信号
),xlabel(t)
subplot(413);plot(rec_sig(1:bits_display*nSamples));title(经过信道后的接
收信号),xlabel(t)
subplot(414);stem(x_pie(1:bits_display));title(解调以后的数据比特
),xlabel(t)
%%计算和绘制信号频谱
[f1,powerIndB1]=getFre(sig_ori,fs);
[f2,powerIndB2]=getFre(mod_sig,fs);
[f3,powerIndB3]=getFre(rec_sig,fs);
figure
subplot(311);plot(f1,powerIndB1);xlim([0fs/2]);ylim([-8010]);title(原始
数据的功率谱密度);xlabel(频率(Hz));ylabel(功率谱密度(w/Hz));
subplot(312);plot(f2,powerIndB2);xlim([0fs/2]);ylim([-8010
文档评论(0)