发布时间 : 星期日 文章通信原理课程设计(1)更新完毕开始阅读
用合适的滤波器对多元码流进行滤波实现，由于生成的是基带信号，因此这一过程又称“基带成形滤波”。 代码如下： bshape.m
%基带升余弦成形滤波器
function y=bshape(x,fs,fb,N,alfa,delay); %设置默认参数 if nargin<6; delay=8; if nargin<5; alfa=0.5; if nargin<4; N=16;
end; end; end;
b=firrcos(N,fb,2*alfa*fb,fs); y=filter(b,1,x); 4.5调制
已调号进行串并转换分离出I分量、Q分量，然后再分别进行电平映射，经过成形滤波后与载波相乘，最后相加得到已调信号。 调制代码： qam.m
function [y,I,Q]=qam(x,Kbase,fs,fb,fc); T=length(x)/fb; m=fs/fb;
nn=length(x); dt=1/fs; t=0:dt:T-dt;
%串/并变换分离出I分量、Q分量，然后再分别进行电平映射
14
I=x(1:2:nn-1);
[I,In]=two2four(I,4*m); Q=x(2:2:nn);
[Q,Qn]=two2four(Q,4*m);
if Kbase==2; %基带成形滤波 I=bshape(I,fs,fb/4); Q=bshape(Q,fs,fb/4); end;
y=I.*cos(2*pi*fc*t)-Q.*sin(2*pi*fc*t); %调制 4.6加入高斯白噪声之后解调
为了简化程序和得到可靠的误码率，我们在解调时并未从已调信号中恢复载波，而是直接产生与调制时一模一样的载波来进行信号解调。 加入不同强度的高斯白噪声代码：
SNR_in_dB=8:2:24; %AWGN信道信噪比 for j=1:length(SNR_in_dB)
y_add_noise=awgn(y2,SNR_in_dB(j)); %加入不同强度的高斯白噪声 y_output=qamdet(y_add_noise,fs,fb,fc); %对已调信号进行解调
解调时先设计一个巴特沃斯滤波器，然后将I分量、Q分量并/串转换，最终恢复成码元序列xn，然后进行解调。 解调的代码如下： %QAM信号解调
function [xn,x]=qamdet(y,fs,fb,fc); dt=1/fs;
t=0:dt:(length(y)-1)*dt;
15
I=y.*cos(2*pi*fc*t); Q=-y.*sin(2*pi*fc*t);
[b,a]=butter(2,2*fb/fs); %设计巴特沃斯滤波器 I=filtfilt(b,a,I); Q=filtfilt(b,a,Q); m=4*fs/fb; N=length(y)/m; n=(.6:1:N)*m; n=fix(n); In=I(n); Qn=Q(n);
xn=four2two([In Qn]);
%I分量Q分量并/串转换，最终恢复成码元序列xn nn=length(xn); xn=[xn(1:nn/2); xn(nn/2+1:nn)]; xn=xn(:); xn=xn';
四进制转换成二进制代码如下： %四进制转换成二进制 function xn=four2two(yn); y=yn; ymin=min(y); ymax=max(y);
ymax=max([ymax abs(ymin)]);
16
ymin=-abs(ymax); yn=(y-ymin)*3/(ymax-ymin); %设置门限电平，判决 I0=find(yn< 0.5); yn(I0)=zeros(size(I0)); I1=find(yn>=0.5 & yn<1.5); yn(I1)=ones(size(I1)); I2=find(yn>=1.5 & yn<2.5); yn(I2)=ones(size(I2))*2; I3=find(yn>=2.5); yn(I3)=ones(size(I3))*3;
%一位四进制码元转换为两位二进制码元 T=[0 0;0 1;1 1;1 0]; n=length(yn); for i=1:n;
xn(i,:)=T(yn(i)+1,:); end; xn=xn';
xn=xn(:); xn=xn'; 误码率代码如下： numoferr=0;
for i=1:N if (y_output(i)~=info(i)), numoferr=numoferr+1;
17