声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1676|回复: 2

[编程技巧] 以下QPSK相干解调的代码是什么意思?

[复制链接]
发表于 2010-4-6 11:16 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
t=0:0.001:0.999;%输入比特流
b=randint(1,10,2);
s=b(ceil(10*t+0.01));%ceil值为>=10*t+0.01 的最小整数,s为矩阵的第几列

%QPSK调制
m=cos(2*pi*20*t+s*pi/2);%调制信号
r=awgn(m,0);
T=0:0.001:0.099;
Icorrelator=cos(2*pi*20*T);
Qcorrelator=-sin(2*pi*20*T);
for n=0:9
result(n+1,:)=[trapz(T,r(n*100+1:n*100+100).*Icorrelator),trapz(T,r(n*100+1:n*100+100).*Qcorrelator)]; %这一行是什么意思?
end
s0=[trapz(T,cos(2*pi*20*T).*Icorrelator),trapz(T,cos(2*pi*20*T).*Qcorrelator)]; %这一行是什么意思?
s1=[trapz(T,cos(2*pi*20*T+pi/2).*Icorrelator),trapz(T,cos(2*pi*20*T+pi/2).*Qcorrelator)];
for n=0:9
[C,I(n+1)]=min([norm(result(n+1,:)-s0),norm(result(n+1,:)-s1)]);
end
%从for开始,这段代码实现的是相干解调,但是具体的是什么意思,我不大明白
input=b;
output=I-1;

figure
subplot(2,1,1);
plot(t,input(ceil(10*t+0.01)));
axis([0,1,-0.2,1.2]);
xlabel('Time')
ylabel('Amplitude')
title('Input Waveform')
subplot(2,1,2);
plot(t,output(ceil(10*t+0.01)));
axis([0,1,-0.2,1.2]);
xlabel('Time')
ylabel('Amplitude')
title('Output Waveform')
回复
分享到:

使用道具 举报

发表于 2010-4-6 12:16 | 显示全部楼层
个人水平专业有限, 并不懂LZ的东西!
若是函数使用不清楚, 可善用help功能!
 楼主| 发表于 2010-4-6 12:18 | 显示全部楼层

回复 沙发 ChaChing 的帖子

谢谢!我也还在继续研究中。。。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-13 18:58 , Processed in 0.055334 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表