接收机抽取的Matlab仿真代码,接收机的MATLAB仿真程序代码

clear all; close all;

Tc = 1;

N = 128; %扩频增益

gx= '1000010111000101'; %g(x) = x^15+x^13+x^9+x^8+x^7+x^5+1 g = bin2dec(gx);

state = 1;

L = 2^13;

EcN0dB = -21:-14;

for k=1:length(EcN0dB)

error(k)=0; %计数错误比特数

total(k)=0; %计数总的传输比特数

sigma(k) = sqrt( 10.^(-EcN0dB(k)/10) /2 );

while( error(k)<100 )

%多径结构

p1 = sqrt(0.5/2)*( randn(1,L)+j*randn(1,L) );

p2 = sqrt(0.3/2)*( randn(1,L)+j*randn(1,L) );

p3 = sqrt(0.2/2)*( randn(1,L)+j*randn(1,L) );

t1 = 0;

t2 = 1;

t3 = 2;

[pt state]= mgen(g,state,L+t3);

pt = 2*pt-1;

%数据产生

d = sign( randn(1,L/N) ); %一次64个

%扩频,先将数据扩展,然后与pt点积

dd = sigexpand(d,N);

s = conv( dd,ones(1,N) );

st = s(1:L+t3).*pt(1:L+t3); %扩频

%经过多径信道,加入噪声

z = sigma(k)* ( randn(1,L)+j*randn(1,L) );

rt = st(1:L).*p1 + st(t2+1:L+t2).*p2 + st(t3+1:L+t3).*p3 + z;

%rake接收

r1 = rt.*conj(p1).*pt(1:L);

r2 = rt.*conj(p2).*pt(t2+1:L+t2);

r3 = rt.*conj(p3).*pt(t3+1:L+t3);

%积分

  • 0
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值