清水小鱼 发表于 2007-6-14 09:50

下标出现了0怎么处理?

具体是这样的
clear all
close all
sample_N=16;
n=1:sample_N;
x(n)=sin(2*pi*n/sample_N);
r(n)=randn(1,sample_N)
u(n)=x(n)+r(n)
figure(1);
plot(u(n));
pr=0.01;
d(n)=2*cos(2*pi*n/sample_N);
h1=mean(u(n)*u(n-1));
g1=mean(d(n)*u(n-1));

错误提示是这个:??? Subscript indices must either be real positive integers or logicals.
Error in ==> lms at 12
h1=mean(u(n)*u(n-1))

求均方的公式出现了u(0)这种不存在的情况

如果改成
u(n+1)u(n)的话矩阵又溢出了,
谁帮下忙,谢谢

[ 本帖最后由 eight 于 2007-6-14 10:39 编辑 ]

pengweicai 发表于 2007-6-14 10:13

clear all
close all
sample_N=16;
n=1:sample_N;
x=sin(2*pi*n/sample_N);
r=randn(1,sample_N);
u=x(n)+r(n);

figure(1);
plot(u);
pr=0.01;
d=2*cos(2*pi*n/sample_N);
for n=2:sample_N
    h1=mean(u(n)*u(n-1));
    g1=mean(d(n)*u(n-1));
end

花如月 发表于 2007-6-14 10:20

clear all
sample_N=16;
n=1:sample_N;
x(n)=sin(2*pi*n/sample_N);
r(n)=randn(1,sample_N);
u(n)=x(n)+r(n);
figure(1);
plot(u(n));
pr=0.01;
d(n)=2*cos(2*pi*n/sample_N);
for i=2:length(u)
    h1(i-1)=u(i)*u(i-1);
    g1(i-1)=d(i)*u(i-1);
end
h=mean(h1)
g=mean(g1)
慢了,还不简练,功力不够啊

[ 本帖最后由 花如月 于 2007-6-14 10:23 编辑 ]

eight 发表于 2007-6-14 10:44

原帖由 清水小鱼 于 2007-6-14 09:50 发表 http://www.chinavib.com/forum/images/common/back.gif
具体是这样的
clear all
close all
sample_N=16;
n=1:sample_N;
x(n)=sin(2*pi*n/sample_N);
r(n)=randn(1,sample_N)
u(n)=x(n)+r(n)
figure(1);
plot(u(n));
pr=0.01;
d(n)=2*cos(2*pi*n/sample_N) ...

看这个帖子:常见的程序出错问题整理

ps:看在有两个回帖的份上这个帖子就不删除了,否则一定不会手软。下次发贴前先自己动手
页: [1]
查看完整版本: 下标出现了0怎么处理?