y=yn; ymin=min(y); ymax=max(y); ymax=max([ymax abs(ymin)]); 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';
constel。m %画出星座图
function c=constel(x,fs,fb,fc); N=length(x); m=2*fs/fb; n=fs/fc;
i1=m-n; i=1; ph0=(i1-1)*2*pi/n;
while i <= N/m; xi=x(i1:i1+n-1); y=2*fft(xi)/n; c(i)=y(2); i=i+1; i1=i1+m; end;
%如果无输出,则作图 if nargout<1;
cmax=max(abs(c)); ph=(0:5:360)*pi/180;
plot(1。414*cos(ph),1。414*sin(ph),'c'); hold on;
for i=1:length(c); ph=ph0-angle(c(i));
44
a=abs(c(i))/cmax*1。414; plot(a*cos(ph),a*sin(ph),'r*'); end;
plot([-1。5 1。5],[0 0],'k:',[0 0],[-1。5 1。5],'k:'); hold off; axis equal; axis([-1。5 1。5 -1。5 1。5]); end;
[5] 沈越泓,高媛媛,魏以民.通信原理[M].北京:机械工业出版社,2004: 141-143.[6] Instructor’s Manual to accompany TIMS[M].Tim Hooper,2001:104-105. [7] 李倜,宁金辉. 数字分量图像信号系统的检测(上)[J].世界广播电视,2003,9. [8] 谢筱.试论基带数字与模拟视频测量的区别[J].广播与电视技术,2004,1.
45
搜索“diyifanwen.net”或“第一范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读,第一范文网,提供最新人文社科毕业设计(论文)-基于MATLAB仿真QAM调制与解调的设计 (5)全文阅读和word下载服务。
相关推荐: