Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
 
 
主题工具 显示模式
旧 2008-07-02, 09:59   #1
nospecial
初级会员
 
注册日期: 2008-06-22
年龄: 37
帖子: 12
声望力: 17
nospecial 正向着好的方向发展
默认 [求助]报错了

请高手帮我看下这是怎么回事!!

clear all;
X1=1;X2=0;X3=1;X4=0;X5=1;X6=0; %移位寄存器输入Xi初T态(010101), Yi为移位寄存器各级输出
M=300; %置M序列总长度
for i=1:M
Y6=X6;Y5=X5;Y4=X4;Y3=X3;Y2=X2;Y1=X1;
X6=Y5;X5=Y4; X4=Y3;X3=Y2;X2=Y1;
X1=xor(Y3,Y4); %异或运算
if Y6==0
u(i)=-1;
else
u(i)=Y6;
end
end
z=zeros(1,M);
for k=4:M
z(k)=1.5*z(k-1)-0.9*z(k-2)+0.3*z(k-3)+2*u(k-1)+1.5*u(k-2)+0.7*u(k-3);
end

N=2;
[HL2,ZL2]=cf(M,N,z,u);
[C2,AIC(2)]=fv(HL2,ZL2,N);
if AIC(2)<aic
aic=AIC(2);C=C2;n=2;
end

报下面的错:
??? Error using ==> mtimes
Inner matrix dimensions must agree.

Error in ==> fv at 4
c3=HL'*ZL;




我的fv:
function [C1,AIC]= fv(HL,ZL,n)
c1=HL'*HL;
c2=inv(c1);
c3=HL'*ZL;
C1=c2*c3
m=ZL-HL*C1;
h=m';
dk=1/300*h*m;
AIC=300*log(dk)+4*n;

我的cf:
function[HL,ZL]=cf(M,N,z,u)
for n=1:N
for k=(N+1):M
HL(k-N,n)=-z(k-n);
end
end
for n=1:N
for k=(N+1):M
HL(k-N,n+N)=u(k-n);
end
end
for k = (N+1):M
ZL(k-1, = z(k);
end
:cry:
nospecial 当前离线   回复时引用此帖
 

主题工具
显示模式

发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码


相似的主题
主题 主题作者 版面 回复 最后发表
freeman链码 勇敢飞 MATLAB论坛 1 2012-05-07 22:53
关于wavedec和wrcoef zhuimeng MATLAB论坛 0 2008-09-19 15:08
[求助]如何导入列数不一致的数据? dishengpiao MATLAB论坛 4 2008-05-06 10:29
【求助】instrcallback触发问题 xukeuniversal MATLAB论坛 2 2008-01-27 09:55


所有时间均为北京时间。现在的时间是 15:00


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.