MATLAB程序出现 Index exceeds matrix dimensions. 错误 为什么？求高手解答

Question

Numusers=1;
Nc=16; %扩频因子
ISI_Length=1; %每径延时为ISI_Length/2（码间干扰）
Tlen=5000;%数据长度
EbN0db = [0:2:100];
Bit_Error_Number1=0;%误比特率的初始值
a1=0.8;%每径幅度衰减
a2=0.2;
a3=0.1;
power_unitary_factor1=sqrt(a1);%每径功率因子
power_unitary_factor2=sqrt(a2);
power_unitary_factor3=sqrt(a3);
s_initial=randsrc(1,Tlen);%数据源，以1/2的概率随机产生1或-1的1行Tlen列矩阵
%产生Walsh矩阵
Wal2=[1 1;1 -1];
Wal4=[Wal2 Wal2;Wal2 Wal2*(-1)];
Wal8=[Wal4 Wal4;Wal4 Wal4*(-1)];
Wal16=[Wal8 Wal8;Wal8 Wal8*(-1)];
%扩频
s_spread=zeros(Numusers,Tlen*Nc);
ray1=zeros(Numusers,2*Tlen*Nc);
ray2=zeros(Numusers,2*Tlen*Nc);
ray3=zeros(Numusers,2*Tlen*Nc);
ray4=zeros(Numusers,2*Tlen*Nc);
ray5=zeros(Numusers,2*Tlen*Nc);
for i=1:Numusers
x0=s_initial(i,:).'*Wal16(8,:);
x1=x0.';
s_Spread(i,:)=(x1(:)).';
end
%将每个扩频后输出重复为两次，有利于下面的延迟（延迟了半个码元）
ray1(1:2:2*Tlen*Nc-1)=s_Spread(1:Tlen*Nc);
ray1(2:2:2*Tlen*Nc)=ray1(1:2:2*Tlen*Nc-1);

%产生第二径和第三径信号
ray2(ISI_Length+1:2*Tlen*Nc)=ray1(1:2*Tlen*Nc-ISI_Length);
ray3(2*ISI_Length+1:2*Tlen*Nc)=ray1(1:2*Tlen*Nc-2*ISI_Length);
ray4(1:2*Tlen*Nc-ISI_Length)=ray1(2:2*Tlen*Nc);%表示减小1个时延
ray5(1:2*Tlen*Nc-2*ISI_Length)=ray1(3:2*Tlen*Nc);
for nEN = 1:length(EbN0db)
en = 10^(EbN0db(nEN)/10);% 将 Eb/N0 的dB值转换为普通十进制数值
sigma = sqrt((32/(2*en)));
sig=32/(2*en);
%接收到的信号demp
demp1=a1*ray1 + a2*ray2 + a3*ray3 + (randn(1,2*Tlen*Nc) + randn(1,2*Tlen*Nc)*i)*sig;
%接收到的信号进行时间反转技术处理
demp=(a1*a1+a2*a2+a3*a3)*ray1+(a1*a2+a2*a3)*ray2+a1*a3*ray3+(a1*a2+a2*a3)*ray4+a1*a3*ray5+(randn(1,2*Tlen*Nc) + randn(1,2*Tlen*Nc)*i)*sigma;

dt=reshape(demp,32,Tlen)';
%将Walsh码重复为两次
Wal16_d(1:2:31)=Wal16(8,1:16);
Wal16_d(2:2:32)=Wal16(8,1:16);
%解扩后的输出
rdata1=dt*Wal16_d(1,:).';
rd_m=0.2*32;
%判决输出
r_Data1=sign(rd_m)';
%计算误比特率
Bit_Error_Number1=length(find(r_Data1(1:Tlen)~=s_initial(1:Tlen)));
Bit_Error_Rate1(nEN)=Bit_Error_Number1/(Tlen);
end
semilogy(EbN0db,Bit_Error_Rate1,'*-');hold on;
xlabel('信噪比');
ylabel('误比特率');
title('时间反转技术接收性能');
我想知道为什么是在最后出的问题 在at 59出错的{%计算误比特率
Bit_Error_Number1=length(find(r_Data1(1:Tlen)~=s_initial(1:Tlen)));
Bit_Error_Rate1(nEN)=Bit_Error_Number1/(Tlen)};，求高手帮忙改改
谢谢

Answer 1

Bit_Error_Number1=length(find(r_Data1(1:Tlen)~=s_initial(1:Tlen)));
你的r_Data是一个1乘1的数据，没有那么多。
Tlen=5000,明显是不行的。
帮你从逻辑上简单debug了一下：
Numusers=1;
Nc=16; %扩频因子
ISI_Length=1; %每径延时为ISI_Length/2（码间干扰）
Tlen=5000;%数据长度
EbN0db = [0:2:100];
Bit_Error_Number1=0;%误比特率的初始值
a1=0.8;%每径幅度衰减
a2=0.2;
a3=0.1;
power_unitary_factor1=sqrt(a1);%每径功率因子
power_unitary_factor2=sqrt(a2);
power_unitary_factor3=sqrt(a3);
s_initial=randsrc(1,Tlen);%数据源，以1/2的概率随机产生1或-1的1行Tlen列矩阵
%产生Walsh矩阵
Wal2=[1 1;1 -1];
Wal4=[Wal2 Wal2;Wal2 Wal2*(-1)];
Wal8=[Wal4 Wal4;Wal4 Wal4*(-1)];
Wal16=[Wal8 Wal8;Wal8 Wal8*(-1)];
%扩频
s_spread=zeros(Numusers,Tlen*Nc);
ray1=zeros(Numusers,2*Tlen*Nc);
ray2=zeros(Numusers,2*Tlen*Nc);
ray3=zeros(Numusers,2*Tlen*Nc);
ray4=zeros(Numusers,2*Tlen*Nc);
ray5=zeros(Numusers,2*Tlen*Nc);
for i=1:Numusers
x0=s_initial(i,:).'*Wal16(8,:);
x1=x0.';
s_Spread(i,:)=(x1(:)).';
end
%将每个扩频后输出重复为两次，有利于下面的延迟（延迟了半个码元）
ray1(1:2:2*Tlen*Nc-1)=s_Spread(1:Tlen*Nc);
ray1(2:2:2*Tlen*Nc)=ray1(1:2:2*Tlen*Nc-1);

%产生第二径和第三径信号
ray2(ISI_Length+1:2*Tlen*Nc)=ray1(1:2*Tlen*Nc-ISI_Length);
ray3(2*ISI_Length+1:2*Tlen*Nc)=ray1(1:2*Tlen*Nc-2*ISI_Length);
ray4(1:2*Tlen*Nc-ISI_Length)=ray1(2:2*Tlen*Nc);%表示减小1个时延
ray5(1:2*Tlen*Nc-2*ISI_Length)=ray1(3:2*Tlen*Nc);
for nEN = 1:length(EbN0db)
en = 10^(EbN0db(nEN)/10);% 将 Eb/N0 的dB值转换为普通十进制数值
sigma = sqrt((32/(2*en)));
sig=32/(2*en);
%接收到的信号demp
demp1=a1*ray1 + a2*ray2 + a3*ray3 + (randn(1,2*Tlen*Nc) + randn(1,2*Tlen*Nc)*i)*sig;
%接收到的信号进行时间反转技术处理
demp=(a1*a1+a2*a2+a3*a3)*ray1+(a1*a2+a2*a3)*ray2+a1*a3*ray3+(a1*a2+a2*a3)*ray4+a1*a3*ray5+(randn(1,2*Tlen*Nc) + randn(1,2*Tlen*Nc)*i)*sigma;

dt=reshape(demp,32,Tlen)';
%将Walsh码重复为两次
Wal16_d(1:2:31)=Wal16(8,1:16);
Wal16_d(2:2:32)=Wal16(8,1:16);
%解扩后的输出
rdata1=dt*Wal16_d(1,:).';
rd_m=0.2*32;
%判决输出
r_Data1=sign(rd_m)';
%计算误比特率
Bit_Error_Number1=length(find(rdata1(1:Tlen)'~=s_initial(1:Tlen))); %考虑到你可能输错了
%变量，再就是差一个“转置”
Bit_Error_Rate1(nEN)=Bit_Error_Number1/(Tlen);
end
semilogy(EbN0db,Bit_Error_Rate1,'*-');hold on;
xlabel('信噪比');
ylabel('误比特率');
title('时间反转技术接收性能');

追问

问题还是没有解决 不过还是谢谢你了
不过程序里面的find(r_data1(1:Tlen)~=s_initial(1:Tlen))这里好像是不能加转置的，这里的两个相同行列的矩阵才能够找到不相同的地方，要不然这样出来的结果就是这误比特率为0了，那这个程序就一点意义都没有了
还有就是关于这个Tlen=5000应该是没有问题的，我这个程序是参照另外的一个程序编写的，那个程序能够完好的运行的 能加个qq253778040帮忙解答下么 我是个matlab新手

Answer 2

j取值超出A的长度了。
i是从1到length(A)
j始终比i大1，所以从2到length(A)+1了，超出范围了。

循环那里用for i=1:length(A)-1就好了。

追问

我想知道为什么是在最后出的问题 在at 59出错的，求高手帮忙改改

Answer 3

Index exceeds matrix dimensions 问题补充：把从一个矩阵中提取某些元素作为子函数怎么写程序啊？ index exceeds matrix dimensions 。说的是下标超过矩阵