通信原理实验--数字基带传输仿真实验(4)

采用匹配滤波方式下系统的可靠性很高。 2、（不采用匹配滤波器） %输入符号序列,形成发送信号 M=8;%符号数 k=4;

T0=1;Ts=k*T0; Rs=1/Ts;

fs=1/T0;%抽样频率 Bs=fs/2;%折叠频率

x=1-2*(rand(1,M)>0.5); x0=x>0; n=0:M-1;

subplot(4,2,1),stem(x0,'b.'); axis([0 M+1 0 1.5]) title('输入符号序列');

n=0:k*M-1;

x1=zeros(1,k*M); for i=1:M

x1(k*i-(k-1))=x(i); end

subplot(4,2,2);stem(n,x1,'.'); title('发送信号');

%升余弦滚降滤波器 N=31; alpha=1;

[h,n]=syx(alpha,N);

subplot(4,2,3);stem(n,h,'.') title('升余弦滚降滤波器');

%发送滤波器输出 y=conv(x1,h);

n=-(N-1)/2:k*M-1+(N-1)/2; subplot(4,2,4);plot(n,y);

axis([-(N-1)/2,k*M-1+(N-1)/2,1.1*min(y),1.1*max(y)]); title('发送滤波器输出');

%计算平均每比特功率 Eb=0;

for i=1:length(y)

Eb=Eb+abs(y(i))*abs(y(i));

end

Eb=Eb/(M*fs/Rs); %信噪比 SNR=20; %产生噪声

sgma=sqrt(Eb/(10^(SNR/10))/2); P=M*fs/Rs+N-1;

n0=gnguass(0,sgma,P); t=0:P-1;

subplot(425);plot(t,n0);

axis([0 P-1 -12*sgma 12*sgma]) title('噪声'); %形成接收信号 r=y+n0;

subplot(426);plot(n,r);

axis([-(N-1)/2,k*M-1+(N-1)/2,1.1*min(r),1.1*max(r)]); title('接收信号'); %抽样判决

sam=zeros(1,M); for i=0:M-1

c=find(n==i*fs/Rs); sam(i+1)=r(c); end

n1=0:M-1;

subplot(4,2,7),stem(n1,sam,'.');

title('抽样值') b=zeros(1,M); for i=1:M

if sam(i)>0 b(i)=1; end end

subplot(4,2,8),stem(b,'b.'); axis([0 M+1 0 1.5]) title('判决结果'); %眼图

eyediagram(y,k,1,3);

title('发送滤波器输出信号眼图'); eyediagram(r,k,1,3);

title('接收滤波器输出信号眼图'); %星座图

scatterplot(y,k,0,'r+');

title('发送滤波器输出信号星座图');

scatterplot(r,k,0,'r+');

title('接收滤波器输出信号星座图'); 仿真结果；