Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

lai5011 发表于 2006-4-12 12:11

求助，改了好久改不过来

s=; % m序列1初始化值 b=; % m序列2初始化值 %function bm = subsindex(s) for i=1:p al=bm(s(2),s(11),0,0,0,0,0,0); s(2),s(11)可以为s的第2与第11个元素啊 为什么报错 ??? Subscript indices must either be real positive integers or logicals.
Error in ==> squense1 at 16 al=bm(s(2),s(11),0,0,0,0,0,0);

happy 发表于 2006-4-12 15:47

回复：（lai5011）求助，改了好久改不过来

<DIV class=quote>以下是引用lai5011在2006-4-12 12:11:24的发言： 
s=; % m序列1初始化值 b=; % m序列2初始化值 %function bm = subsindex(s) for i=1:p al=bm(s(2),s(11),0,0,0,0,0,0); s(2),s(11)可以为s的第2与第11个元素啊 为什么报错 ??? Subscript indices must either be real positive integers or logicals.
Error in ==> squense1 at 16 al=bm(s(2),s(11),0,0,0,0,0,0);</DIV>
你这里的bm是什么东西啊

cdwxg 发表于 2006-4-12 23:55

你上面的bm=subsindex（s）又有什么用呢？ 那你下面为什么不直接用subsindex代替bm呢，这样不是绕路走么？ 再说，好象subsindex不能有太多的输入吧。 你的错误报的是下标等要么是正整数或者逻辑值吧。 不晓得究竟错在哪里，你要不全部拿来，好让看下全部运行后的错误。

lai5011 发表于 2006-4-13 17:57

1、bm是一函数，另一m文件 2、%function bm = subsindex(s) 这条语句不要，隐掉 % 用matlab构建直接序列扩频通信系统，pn码用gold码，m序列用11级斐波纳契结构（SSRG）移位寄存器实现 % 通过加入加性高斯白噪声来判定系统性能，得出不同信噪比情况下的误码率 % 程序结构：gold码产生-》传送数据产生-》扩频-》加入噪声-》解扩-》样本求和、信息判决-》误码率记数 function=squense(N,s,p) clear all global bm Lc=input('please input the value of Lc:'); % 每一信息比特的码片数目，即处理增益 N=input('please input the value of N:'); % 传送的信息比特数目 g=input('please input the value of g:'); % 是否要加高斯噪声 p=2047; % 扩频序列长度 n=11; % 序列级数 s=; % m序列1初始化值 b=; % m序列2初始化值 %function bm = subsindex(s) for i=1:p al1=bm(s(2),s(11),0,0,0,0,0,0); for j=n:-1:2 s(j)=s(j-1) end s(1)=al1; l1(i)=s(n); % m序列1 end for i=1:p al2=bm(b(2),b(3),b(8),b(11),0,0,0,0); for j=n:-1:2 b(j)=b(j-1) end b(1)=al2; l2(i)=b(n); % m序列2 end 
 l3=xor(l1,l2); % 生成gold码序列l33（i） for i=1:p if (l3(i)>0); l33(i)=-1; else l33(i)=1; end end k=cell(N*Lc/p)-1; for i=0:k for j=1:p l4(j-p*i)=l33(j); % 多个gold码序列 end end sgma=1 for SNRindb=1:30 snr=10^(SNRindb/10); % 单位转换db-》一般 Eb=2*sgma^2*snr; % 信号能量 E_chip=Eb/Lc; % 每码片的信号能量 num_of_err=0; % 误码率初始值 for i=0:N-1 % 产生均匀分布的数据信号 temp=rand; if (temp<0.5) data=-1; else data=1; end for j=1:Lc % 数据重复lc次 Tb_data(j)=data; PN(j)=l4(j-i*Lc); end Tx_data=sqrt(E_chip)*Tb_data.*PN; % 扩频 noise=g*sgma*randn(1,Lc); % 产生加性高斯白噪声 %h=(i*Lc-1):(i-1)*Lc; Rx_data=Tx_data+noise; % 加入高斯白噪声 plot(Rx_data) keyboard temp=Rx_data.*PN; % 解扩 decision_variable=sum(temp); % 样本求和 if(decision_variable<0) % 输出判决 decision=-1; else decision=1; end if(decision~=data) num_of_err= num_of_err+1; % 误码率记数 else num_of_err=num_of_err; end end Err(SNRindb)=num_of_err/N; % 单位数据误码率 double(err); keyboard; end subplot(2,1,2); Xlim(); Ylim(); Err plot(err);
%存成另一文件bm.m function=bm(x,y,z,h,r,t,p,q) % m序列生成算法 add=x+y+z+h+r+t+p+q; if(add==0|add==2|add==4|add==6|add==8) add=0; else add=1; end 
[此贴子已经被作者于2006-4-14 10:41:39编辑过]

AaronSpark 发表于 2006-4-13 23:06

回复：（lai5011）1、bm是一函数，另一m文件2、%funct...

<DIV class=quote>以下是引用lai5011在2006-4-13 17:57:56的发言： 
1、bm是一函数，另一m文件 2、%function bm = subsindex(s) 这条语句不要，隐掉 % 用matlab构建直接序列扩频通信系统，pn码用gold码，m序列用11级斐波纳契结构（SSRG）移位寄存器实现 % 通过加入加性高斯白噪声来判定系统性能，得出不同信噪比情况下的误码率 % 程序结构：gold码产生-》传送数据产生-》扩频-》加入噪声-》解扩-》样本求和、信息判决-》误码率记数 function=squense(N,S,s,p,bm) clear all global bm </DIV>
既然bm是函数，你在红色标记部分又把它作为变量使用，明显就错了

lai5011 发表于 2006-4-14 10:55

我改了下，问题还没解决，同样的错误，各位大哥，看看还错在哪里

其中： l1是变化的量，会随al1变化，l2随al2， al1、al2随s（i）变化

happy 发表于 2006-4-14 10:57

回复：（lai5011）求助，改了好久改不过来

见5楼

cdwxg 发表于 2006-4-14 11:00

他可能想知道错了怎么改吧，偶不会..

lai5011 发表于 2006-4-14 11:10

问题还没解决，

谢谢各位，不过问题还没解决 各路好手能不能全看看， 看有没有其他定义上的错误

happy 发表于 2006-4-14 11:54

回复：（lai5011）问题还没解决，

<DIV class=quote>以下是引用lai5011在2006-4-14 11:10:10的发言： 谢谢各位，不过问题还没解决 各路好手能不能全看看， 看有没有其他定义上的错误</DIV>
 自己调吧，别人没那么多功夫，有问题再问

lai5011 发表于 2006-4-17 21:19

有点急了，无奈，只好求助

改了好多地方，没有效果 还是知识太浅了 不知道看哪方面的细节知识更有效 希望各位好手帮我指出程序中的错误 并说说看哪方面的知识有助于编程调测 谢谢

happy 发表于 2006-4-17 21:39

回复：（lai5011）求助，改了好久改不过来

把global bm去掉，四楼不是说过了嘛
[此贴子已经被作者于2006-4-17 21:41:35编辑过]

lai5011 发表于 2006-4-18 11:25

谢谢了

是我理解错了 现在结果接近理想了

页: [1]

声振论坛's Archiver

求助，改了好久改不过来

回复：（lai5011）求助，改了好久改不过来

回复：（lai5011）1、bm是一函数，另一m文件2、%funct...

我改了下，问题还没解决，同样的错误，各位大哥，看看还错在哪里

回复：（lai5011）求助，改了好久改不过来

问题还没解决，

回复：（lai5011）问题还没解决，

有点急了，无奈，只好求助

回复：（lai5011）求助，改了好久改不过来

谢谢了