MATLAB仿真在通信和电子工程中应用第6章信号处理办法.ppt

2.眼图和星座图 工程上，眼图是观察带限信道传输码间串扰程度的方便手段。MATLAB通信工具箱中提供了眼图显示函数“eyediagram”。其常用调用格式是： ? eyediagram(x,n); eyediagram(x,n,period); eyediagram(x,n,period,offset); eyediagram(x,n,period,offset,plotstring); h=eyediagram(...); 其中，x是输入信号序列，n为每一波形轨迹的样值数，period为作图的时间范围刻度，在-period/2到period/2之间，offset为时间偏移量，在0到n-1之间，plotstring是作图线型选项。 MATLAB中观察多元数字调制输出的星座图可使用指令“scatterplot”。其调用格式是：? scatterplot(x); scatterplot(x，n); scatterplot(x，n，offset); scatterplot(x，n，offset，plotstring); scatterplot(x，n，offset，plotstring，h); 其中，x是输入信号序列。在x序列中每隔n个值在星座图上作出一个点。offset为偏移量，plotstring是作图线型选项。scatterplot命令的详细用法参见在线帮助文档。举例说明如下： 程序6-37 M=16;Fd=1;Fs=10;%定义M元数和采样率 Pd=100;%计算的轨迹数 msg-d=randint(Pd,1,M);%产生在［0,M-1］范围的随机整数 %使用矩形星座图的QASK调制 msg-a=modmap(msg-d,Fd,Fd,′qask′,M); %假定信道等价于升余弦滤波器 delay=3;%升余弦滤波器的延时 rcv=rcosflt(msg-a,Fd,Fs,′fir/normal′,0.5,delay);%通过滤波器 %去掉升余弦滤波器的暂态响应部分的数据 propdelay=delay.*Fs/Fd+1;%滤波器的适合的时延 rcv1=rcv(propdelay:end-(propdelay-1),:); N=Fs/Fd; %作出接收信号的眼图和星座图（无噪也无串扰） offset1=0; h1=eyediagram(rcv1,N,1/Fd,offset1); h2=scatterplot(rcv1,N,0,′bx′); rcv2=rcv1+0.1*randn(size(rcv1));%接收QASK信号加高斯噪声 %作出接收信号的眼图和星座图（加入噪声的情况） h3=eyediagram(rcv2,N,1/Fd,offset1); h4=scatterplot(rcv2,N,0,′bx′); 程序作出的QASK的I路和Q路眼图如图6-45所示，星座图如图6-46所示。 图6-45 QASK同相支路和正交支路的眼图 （左：无噪情况，右：有噪情况） 图6-46 QASK星座图（左：无噪情况，右：有噪情况） 6.6 数字滤波器设计与实现 6.6.1 数字滤波器的设计参数 滤波器的4个重要的通带、阻带参数是:  fp：通带截止频率(Hz);  fs：阻带起始频率(Hz);  Rp:通带内波动(dB),即通带内所允许的最大衰减; Rs:阻带内最小衰减(dB)。 设采样率(即奈奎斯特速率)为fN，那么可将以上参数中的频率参数转换为归一化角频率参数:  ωp:通带截止角频率(rad/s), ωp =fp/(fN/2);  ωs:阻带起始角频率(rad/s),ωs =fs/(fN/2)。 通过这些参数就可以进行离散滤波器的设计了。例如： ·低通滤波器情况：在采样频率为8000Hz的条件下设计一个低通滤波器，要求通带截止频率为1500Hz，阻带起始频率为2000Hz，通带内波动3dB，阻带内最小衰减50dB，则  ωp=1500/4000，ωs=2000/4000，Rp=3, Rs=50。 ·高通滤波器情况：在采样频率为8000Hz的条件下设计一个高通滤波器，要求通带截止频率为1500Hz，阻带起始频率为1000Hz，通带内波动3dB，阻带内最小衰减65dB,则  ωp=1500/4000，ωs=1000/4000，Rp=3, Rs=65。 ·带通滤波器情况：在采样频率为8000Hz的条件下设计一个带