数字基带信号
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
单极性不归零
function y=snrz(x)
t0=200;
t=0:1/t0:length(x);
for i=1:length(x)
if x(i)==1
for j=1:t0
y((i-1)*t0+j)=1;
end
else
for j=1:t0
y((i-1)*t0+j)=0;
end
end
end
y=[y,x(i)];
plot(t,y);
title('1 0 1 1 0 0 1 0');
grid on;
axis([0,i,-0.1,1.1]);
单极性归零
function y=srz(x)
t0=200;
t=0:1/t0:length(x);
for i=1:length(x)
if x(i)==1
for j=1:t0/2
y((2*i-2)*to/2+j)=1;
y((2*i-1)*t0/2+j)=0;
end
else
for j=1:t0
y((i-1)*t0+j)=0;
end
end
end
y=[y,x(i)];
plot(t,y);
grid on;
axis([0,i,-0.1,1.1]);
x=[10101010]
srz(x)
蒙特卡罗
EbN0dB=0:0.5:10;
N0=10.^(-EbN0dB/10);
sigma=sqrt(N0/2);
%理论计算误码率
pb=0.5*erfc(sqrt(1./N0));
%仿真误码率
for n=1:length(EbN0dB)
%产生等概率信源
a=sign(randn(1,100000));
%离散等效接收模型
rk=a+sigma(n)*randn(1,100000);
dec_a=sign(rk);%判决
%计算误码率
ber(n)=sum(abs(a-dec_a)/2)/length(a); end
semilogy(EbN0dB,pb);
hold;
semilogy(EbN0dB,ber,'rd-');
legend('理论值','仿真结果');
xlabel('Eb/N0(dB)');
ylabel('pb');
双极性不归零
function y=shuangsnrz(x)
t0=200;
t=0:1/t0:length(x);
for i=1:length(x)
if x(i)==1
for j=1:t0
y((i-1)*t0+j)=1;
end
else
for j=1:t0
y((i-1)*t0+j)=-1;
end
end
end
y=[y,x(i)];
plot(t,y);
title('1 0 1 1 0 0 1 0');
grid on;
axis([0,i,-1.1,1.1]);
双极性归零
function y=shuangsrz(x)
t0=200;
t=0:1/t0:length(x);
for i=1:length(x)
if x(i)==1
for j=1:t0/2
y((2*i-2)*t0/2+j)=1;
y((2*i-1)*t0/2+j)=0;
end
else
for j=1:t0/2
y((2*i-2)*t0/2+j)=-1;
y((2*i-1)*t0/2+j)=0;
end
end
end
y=[y,x(i)];
plot(t,y);
grid on;
axis([0,i,-1.1,1.1]);