matlab为什么找不到常用函数,这个说函数找不到，怎么办呢，求助

SNR_dB=0:1:30;  %信噪比范围

LOOP=500;

%------------------------------------------

%--------------信源产生---------------------

%------------------------------------------

for i=1:length(SNR_dB)

snr=10^(SNR_dB(i)/10);     %信噪比转换

sgma=sqrt(8/snr/3);        %噪声方差

error0=0;

error1=0;

error2=0;

error3=0;

for a=1:LOOP

N=2560;         %仿真序列长度

S=randsrc(1,N,[0,1]);  %二进制序列信源      产生n个随机数

%-------------------------------------------

%------------有信道编码----------------------

%-------------------------------------------

S00=convolutionalcode(S);

%--------------QPSK调制---------------------

s=QPSK(S00);

s0=zeros(1,3*N/4);  %将调制后的信息奇偶拆分成两个信息以进行MIMO发射

s1=zeros(1,3*N/4);

for b1=1:2:(3*N/2)

s0(ceil(b1/2))=s(b1);

s1(ceil(b1/2))=s(b1+1);

end

N1=size(s0,2);

%--------有信道编码(2发2收)-----------------

h0=sqrt(0.5).*(randn(1,N1)+1i.*randn(1,N1)); %设置信道冲激响应

h1=sqrt(0.5).*(randn(1,N1)+1i.*randn(1,N1));

h2=sqrt(0.5).*(randn(1,N1)+1i.*randn(1,N1));

h3=sqrt(0.5).*(randn(1,N1)+1i.*randn(1,N1));

%------接收到的信号(冲激响应与发送序列卷积)---

R0=h0.*s0+h1.*s1;                      %接收信号矢量(不含高斯噪声)

R1=-conj(s1).*h0+h1.*conj(s0);         %移动通信P166

R2=h2.*s0+h3.*s1;

R3=-conj(s1).*h2+h3.*conj(s0);

%通过AWGN信道加入噪声后的冲击序列

r0=(R0+sgma.*(randn(1,N1)+1i*randn(1,N1)));%接收信号矢量(有噪声)

r1=(R1+sgma.*(randn(1,N1)+1i*randn(1,N1)));

r2=(R2+sgma.*(randn(1,N1)+1i*randn(1,N1)));

r3=(R3+sgma.*(randn(1,N1)+1i*randn(1,N1)));

s_0=conj(h0).*r0+h1.*conj(r1)+conj(h2).*r2+h3.*conj(r3);%STBC空时译码

s_1=conj(h1).*r0-h0.*conj(r1)+conj(h3).*r2-h2.*conj(r3);

for l0=1:3*N/4     %将接收到的信号合并

Y0(2*l0-1)=s_0(l0);

Y0(2*l0)=s_1(l0);

end

y0=DEQPSK(Y0);  %QPSK解调

m0=viterbi(y0); %维特比译码

[n0,e0]=biterr(m0,S);%计算误比特率

error0=error0+n0;

%有信道编码2发1收

hr0=sqrt(0.5).*(randn(1,N1)+1i.*randn(1,N1)); %瑞利信道函数

hr1=sqrt(0.5).*(randn(1,N1)+1i.*randn(1,N1));

Rr0=hr0.*s0+hr1.*s1;                      %接收信号矢量(不含高斯噪声)

Rr1=-conj(s1).*hr0+hr1.*conj(s0);

rN0=(Rr0+sgma.*(randn(1,N1)+1i*randn(1,N1)));%接收信号矢量(含复高斯噪声)

rN1=(Rr1+sgma.*(randn(1,N1)+1i*randn(1,N1)));

sr_0=conj(hr0).*rN0+hr1.*conj(rN1);%STBC空时译码

sr_1=conj(hr1).*rN0-hr0.*conj(rN1);

for l1=1:3*N/4   %将接收到的信号合并

Y1(2*l1-1)=sr_0(l1);

Y1(2*l1)=sr_1(l1);

end

y1=DEQPSK(Y1);   %QPSK解调

m1=viterbi(y1);  %维特比译码

[n1,e1]=biterr(m1,S);%计算误比特率

error1=error1+n1;

%---------------------------

%无信道编码

%---------------------------

sn=QPSK(S); %QPSK调制

for b2=1:2:(N/2)  %将调制后的信息奇偶拆分成两个信息以进行MIMO发射

sn0(ceil(b2/2))=sn(b2);

sn1(ceil(b2/2))=sn(b2+1);

end

N2=size(sn0,2);

%无信道编码(2发2收)

hn0=sqrt(0.5).*(randn(1,N2)+1i.*randn(1,N2));  %瑞利信道函数

hn1=sqrt(0.5).*(randn(1,N2)+1i.*randn(1,N2));

hn2=sqrt(0.5).*(randn(1,N2)+1i.*randn(1,N2));

hn3=sqrt(0.5).*(randn(1,N2)+1i.*randn(1,N2));

Rn0=hn0.*sn0+hn1.*sn1;                          %接收信号矢量(无噪声)

Rn1=-conj(sn1).*hn0+hn1.*conj(sn0);

Rn2=hn2.*sn0+hn3.*sn1;

Rn3=-conj(sn1).*hn2+hn3.*conj(sn0);

rn0=(Rn0+sgma.*(randn(1,N2)+1i*randn(1,N2))); %接收信号矢量(有噪声)

rn1=(Rn1+sgma.*(randn(1,N2)+1i*randn(1,N2)));

rn2=(Rn2+sgma.*(randn(1,N2)+1i*randn(1,N2)));

rn3=(Rn3+sgma.*(randn(1,N2)+1i*randn(1,N2)));

sn_0=conj(hn0).*rn0+hn1.*conj(rn1)+conj(hn2).*rn2+hn3.*conj(rn3);

sn_1=conj(hn1).*rn0-hn0.*conj(rn1)+conj(hn3).*rn2-hn2.*conj(rn3);

for l2=1:N/4   %将接收到的信号合并

Y2(2*l2-1)=sn_0(l2);

Y2(2*l2)=sn_1(l2);

end

y2=DEQPSK(Y2); %QPSK解调

y2=y2';

[n2,e2]=biterr(y2,S);%计算误比特率

error2=error2+n2;

%无信道编码(2发1收)

hN0=sqrt(0.5).*(randn(1,N2)+1i.*randn(1,N2));   %瑞利信道函数

hN1=sqrt(0.5).*(randn(1,N2)+1i.*randn(1,N2));

RN0=hN0.*sn0+hN1.*sn1;                         %发射序列与接收序列的卷积

RN1=-conj(sn1).*hN0+hN1.*conj(sn0);

rN0=(RN0+sgma.*(randn(1,N2)+1i*randn(1,N2)));  %接收信号矢量(有噪声)

rN1=(RN1+sgma.*(randn(1,N2)+1i*randn(1,N2)));

sN_0=conj(hN0).*rN0+hN1.*conj(rN1);

sN_1=conj(hN1).*rN0-hN0.*conj(rN1);

for l3=1:N/4                                  %将接收到的信号合并

Y3(2*l3-1)=sN_0(l3);

Y3(2*l3)=sN_1(l3);

end

y3=DEQPSK(Y3);    %QPSK解调

y3=y3';

[n3,e3]=biterr(y3,S);%计算误比特率

error3=error3+n3;

end

BER0(i)=error0/(LOOP*N);

BER1(i)=error1/(LOOP*N);

BER2(i)=error2/(LOOP*N);

BER3(i)=error3/(LOOP*N);

end

semilogy(SNR_dB,BER0,'-b+');

hold on;

semilogy(SNR_dB,BER1, 'r');

grid on;

semilogy(SNR_dB,BER2,'--g');

grid on;

semilogy(SNR_dB,BER3,'-.k');

grid on;

legend('2Tx2R:卷积信道编码','2Tx1R有信道编码','2Tx2R无信道编码','2Tx1R无信道编码');

xlabel('SNR/dB');

ylabel('BER');