[求助]帮忙看看这个matlab程序

coldstar88 · 发表于 2005-11-4 09:25

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

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

x

<DIV>clc clear tic; fd=1;fs=3;flag1=0;flag2=0;suma=0;sumb=0; k=10000;er=1;M=4;snr_db=0:12;snr1_db=0:0.1:12; for i=1:length(snr_db) flag1=0;flag2=0; for n=1:k suma=0;sumb=0; matrix=randint(1071,1,2); syn=[1 1 1 0 0 1 0]; matrix_syn=[syn,matrix']'; modata=dmodce(matrix_syn,fd,fs,'psk',M); modata_nois=awgn(modata,snr_db(i)); y=ddemodce(modata_nois,fd,fs,'psk',M); suma=y(1)+y(2)+y(3)+~y(4)+~y(5)+y(6)+~y(7); if suma>=7; flag1=flag1+1; else flag2=flag2+1; end end pm(i)=flag2/k; %实际漏同步概率 end for j=1:length(snr1_db) eb_n0=10^(snr1_db(j)/10); pm_theory(j)=2*qfunct(sqrt(2*log2(M)*eb_n0)*sin(pi/M));%理论误符号率 per(j)=pep(pm_theory(j),13,er); %理论漏同步概率 end semilogy(snr_db,pm,'*',snr1_db,per); xlabel('信噪比/dB'); ylabel('漏同步概率'); legend('实际漏同步概率','理论漏同步概率'); toc</DIV>
<DIV>我有几个地方不明白，请教你一下，1.取1071位是什么意思呢，就算加上7位巴克码，是1078位，说是5个语音帧插入一个同步帧，可是语音帧不是137bit吗？那同步帧又是多少位呢？2.这个程序的信噪比就是设计成等差数列吗？3.if suma>=7，这里为什么是7呢，下面理论漏同步概率中er取的是1，既是允许码组中错误码元的最大数是1，那这里就应该是6啊，因为允许有一位出错阿。4.per(j)=pep(pm_theory(j),13,er)中为什么取13呢？ %漏同步概率 function [y]=pep(pe,n,k) % pe 误码率 % n 码长 y=0; for i=k:n y=y+factorial(n)/(factorial(i)*factorial(n-i))*pe^i*(1-pe)^(n-i); end</DIV>
<DIV>n是同步码组的码元数，这里取13说明同步帧是13位吗？还是别的什么原因？</DIV>

suffer · 发表于 2005-11-4 15:26

这个问题好像不是matlab问题，好像是语音信号处理问题吧

账号		自动登录	找回密码
密码			我要加入

[求助]帮忙看看这个matlab程序

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

回复：（coldstar88）[求助]帮忙看看这个matlab程序