求教编程提速去掉此for循环

321forever

1. 
   
2. %这是一个移位求平均值的程序
   
3. length=1000; %信号的总长度
   
4. rlength=100; %移动位置为100
   
5. average= zeros(1000,1); %初始化
   
6. for i=1:901                                %这个for循环有些占用时间看看看看高手有没有什么替代或是提速的方法
   
7. average(i)=mean(signal(i:i+99));
   
8. end
   
9. average(902:1000)=average(901);

复制代码

hebut

matlab自身有移动平均现成的命令
matlab搜索：Moving Average Filter

1. a = 1;
   
2. b = [1/4 1/4 1/4 1/4];
   
3. y = filter(b,a,x);

复制代码

321forever

回复 2 # hebut 的帖子

谢谢但是我想实现的功能好像和这个moving average filter有些区别，

1. x=[1 2 3 4 5 6 7 8 9 10];
   
2. a=1;
   
3. b=[1/4 1/4 1/4 1/4];
   
4. y=filter(b,a,x);
   
5. % 结果是
   
6. y =
   
7. 
   
8. 0.2500 0.7500 1.5000 2.5000 3.5000 4.5000 5.5000 6.5000 7.5000 8.5000
   

复制代码

而按我的程序

1. y=2.5 3.5 4.5 5.5 6.5 7.5 8.5 8.5 8.5 8.5

复制代码

hebut

你的结果第一个值是前四个数相加平均，filter是前面补3个0，然后相加平均，所以你的结果比它的结果向前平移了四位。
也就是：
它的结果是：y(n)=1/4*(n-3)+1/4*(n-2)+1/4*(n-1)+1/4*(n)
你的结果是：y(n)=1/4*(n)+1/4*(n+1)+1/4*(n+2)+1/4*(n+3)
平移一下就可以了

qibbxxt

回复 1 # 321forever 的帖子

1. %这是一个移位求平均值的程序
   
2. clear;clc;
   
3. length=1000; %信号的总长度
   
4. rlength=100; %移动位置为100
   
5. average= zeros(1000,1); %初始化
   
6. aver = zeros(1000,1);
   
7. ave = zeros(1000,1);
   
8. signal = rand(1,1000);
   
9. warning off
   
10. %% 用循环
    
11. tic;
    
12. for i=1:901 %这个for循环有些占用时间看看看看高手有没有什么替代或是提速的方法
    
13. average(i)=mean(signal(i:i+99));
    
14. end
    
15. average(902:1000)=average(901);
    
16. toc;
    
17. %% 用arrayfun
    
18. tic
    
19. aver(1:901) = arrayfun(@(x)mean(signal(x:x+99)),1:901);
    
20. aver(902:end) = aver(901);
    
21. toc
    
22. % isequal(aver,average)
    
23. %% 用hankel矩阵构造
    
24. tic
    
25. b = hankel(signal(1:100),[0,signal(101:end)]);
    
26. ave(1:901) = mean(b);
    
27. ave(902:end) = ave(901);
    
28. toc
    
29. isequal(aver,ave,average)

复制代码

运行结果：

1. 
   
2. 
   
3. Elapsed time is 0.015726 seconds.
   
4. Elapsed time is 0.022444 seconds.
   
5. Elapsed time is 0.004636 seconds.
   
6. 
   
7. ans =
   
8. 
   
9. 1

复制代码

321forever

回复 4 # hebut 的帖子

谢谢了啊，移位不错的方法，能帮忙再讲讲move average filter有什么作用么

321forever

回复 5 # qibbxxt 的帖子

谢谢，我试了你给的hankel，对于30000以下的可以运行，但要是50000点的话，会提示

1. c=[1:50000];
   
2. d=hankel(c(1:1000),[0,c(1001:end)]);
   
3. ?????Out of memory. Type HELP MEMORY for your options;

复制代码

想问下有没有什么解决办法，还是我的电脑内存的问题

hebut

回复 6 # 321forever 的帖子

我也是只用过移动平均这个功能，就是这个式子：y(n)=1/4*(n-3)+1/4*(n-2)+1/4*(n-1)+1/4*(n)，其他的没看过，帮助里面应该有的吧

321forever

回复 8 # hebut 的帖子

恩谢谢啊，那我就再看看

西门牛

进来看看，学习学习