滤波的程序

dlxzw

1. clc;
   
2. % load L20X.txt;
   
3. % y = L20X;
   
4. load M1;
   
5. load L1;
   
6. figure(1)
   
7. subplot(2,1,1),plot(M1);
   
8. subplot(2,1,2),plot(L1);
   
9. 
   
10. 
    
11. N = size(M1);
    
12. figure(2)
    
13. % n1 = find(y<=min(y));
    
14. % plot(abs(fft(y)));
    
15. % y1  =ones(N(1,1),1);
    
16. % a = y(1);
    
17. % y1 = a.*y1;
    
18. % y = y - y1;
    
19. b = fir1(100,0.2,0.3,blackman(101));
    
20. m1 = filter(b,1,M1);
    
21. l1 = filter(b,1,L1);
    
22. subplot(2,1,1),plot(m1);
    
23. subplot(2,1,2),plot(l1);
    
24. 
    
25. % x = x + y1;
    
26. % n2 = find(x<=min(x));
    
27. % n = n2-n1;
    
28. % figure(3);
    
29. % % plot(abs(fft(x)));
    
30. % plot(x);
    
31. % save y.txt y -ASCII
    
32. 
    
33. p_a = find(l1>=max(l1))
    
34. p_b = find(l1<=min(l1))
    
35. figure(3)
    
36. plot([l1(31:(p_b-p_a+1)*2); l1(31)],[m1(31:(p_b-p_a+1)*2); m1(31)])
    
37. hold on

复制代码

[ 本帖最后由 suffer 于 2006-10-9 20:34 编辑 ]

suffer

能否对这个程序作一个简单的介绍？

zonglei

光秃秃的一个程序肯定运行不起来的
像load M1;
load L1;
M1，L1是什么？？

suffer

原帖由 zonglei 于 2006-11-1 01:15 发表
光秃秃的一个程序肯定运行不起来的
像load M1;
load L1;
M1，L1是什么？？

这是载入数据文件，你没有该文件当然不能运行了

pigmin

介绍还是需要的

叶川

希望大家以后在发代码的时候有一些简单的注释,毕竟这里新人比较多

hbyyj2008

没有说明，让我很受伤！想说爱你不容易！不过还是要谢谢楼主！

wangyongq246

看别人的代码其实是一种考验，也是一种折磨，还是自己先写！

suffer

原帖由 wangyongq246 于 2006-11-16 16:11 发表
看别人的代码其实是一种考验，也是一种折磨，还是自己先写！

呵呵，确实如此，看别人的代码确实比自己写还要困难一些，个人感觉

fenghai222

原帖由 wangyongq246 于 2006-11-16 16:11 发表
看别人的代码其实是一种考验，也是一种折磨，还是自己先写！

是呀，不过别人能说明白的话是最好的了，这样就可以多听取一点经验了。

suffer

原帖由 fenghai222 于 2006-11-21 21:58 发表

是呀，不过别人能说明白的话是最好的了，这样就可以多听取一点经验了。

这个其实要求很高~

yangji

不错，运行试试

lym

很好 好好学习