完善资料让更多小伙伴认识你,还能领取20积分哦, 立即完善>
clear all
close all t=[-1:0.01:7-0.01]; %变量赋值 tt=length(t); %循环的长度 x1=ones(1,800); %定义1*1800全1矩阵x1 for i=1:tt if(t(i)>=-1&t(i)<=1)| (t(i)>=5&t(i)<=7); %如果-1= else x1(i)=-1; %否则,矩阵x1的对应元素=-1 end end t1=[0:0.01: 8-0.01]; t2=0:0.01: 7-0.01; t3=-1:0.01:7-0.01; t4=0:0.01: 8-0.01; %定义4个变量 tt1=length(t1); %循环的长度 x2=ones(1,800); %定义1*1800全1矩阵x2 for i=1:tt1 if(t1(i)>=0&t1(i)<=2)| (t1(i)>=4&t1(i)<=8); %如果0= else x2(i)=-1; %否则,矩阵x2的对应元素=-1 end end f=0:0.1:1; %定义频率变量f xrc=0.5+0.5*cos(pi*f); %给函数xrc赋值 y1_rc=conv(x1,xrc)/5.5; %计算x1和xrc的卷积 y2_rc=conv(x2,xrc)/5.5; %计算x2和xrc的卷积 n0=randn(size(t2)); %产生与t2一样大的随机矩阵n0 f1=1; i=x1.*cos(2*pi*f1*t); %计算i分量 q=x2.*sin(2*pi*f1*t); %计算q分量 I=i(101:800); Q=q(1:700); %截取部分数据作为I、Q QPSK=sqrt(1/2).*I+sqrt(1/2).*Q; %计算正交相移键控 QPSK_n=(sqrt(1/2).*I+sqrt(1/2).*Q)+n0; %叠加前面的随机数n0 n1=randn(size(t2)); %产生与t2一样大的随机矩阵n1 i_rc=y1_rc.*cos(2*pi*f1*t3); %计算i分量 q_rc=y2_rc.*sin(2*pi*f1*t4); %计算q分量 I_rc=i_rc(101:800); Q_rc= q_rc (1:700); QPSK_rc =(sqrt(1/2).*I_rc +sqrt(1/2).*Q_rc); QPSK_rc_n1= QPSK_rc +n1; %叠加前面的随机数n1 figure(1) subplot(4,1,1);plot(t3, i_rc);axis([-1 8 -1 1]);ylabel('a序列'); %图形1有4个子图 %画变量t3的i分量序列 subplot(4,1,2);plot(t4, q_rc);axis([-1 8 -1 1]);ylabel('b序列');%子图2%画变量t4的q分量序列 subplot(4,1,3);plot(t2, QPSK_rc);axis([-1 8 -1 1]);ylabel('合成序列'); %子图3 %画变量t2的QPSK subplot(4,1,4);plot(t2, QPSK_rc_n1);axis([-1 8 -1 1]);ylabel('加入噪声'); %子图4%画变量t2的A加入噪声的QPSK 为什么程序 i_rc=y1.*cos(2*pi*f1*t3); %计算i分量 q_rc=y2.*sin(2*pi*f1*t4); %计算q分量 会有错呢,已经看了好久还是没有发现啊! |
|
相关推荐 |
|
你正在撰写答案
如果你是对答案或其他答案精选点评或询问,请使用“评论”功能。
小黑屋| 手机版| Archiver| 电子发烧友 ( 湘ICP备2023018690号 )
GMT+8, 2024-9-2 04:13 , Processed in 0.553945 second(s), Total 75, Slave 56 queries .
Powered by 电子发烧友网
© 2015 bbs.elecfans.com
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 (电路图) 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证:合字B2-20210191