matlab和labview编程的区别

realhappy

这是一个用matlab编写的程序，很多是用循环处理的，用labview处理难度就大了，嘿嘿，不服的写写看！

1. %ZoomFFT谱
   
2. %x-信号序列
   
3. %fs-采样频率
   
4. %N-做谱点数
   
5. %fe-分析中心频率
   
6. %D-细化倍数
   
7. %L-平均段数
   
8. %M-滤波器半阶数
   
9. %f-返回频率向量
   
10. %xz-返回幅值谱
    
11. 
    
12. function [f xz]=ZoomFFT(x,fs,N,fe,D,L,M)
    
13. 
    
14. k=1:M;                          
    
15. w=0.5+0.5*cos(pi*k/M);          %Hanning窗
    
16. 
    
17. fl=max(fe-fs/(4*D),-fs/2.2);
    
18. fh=min(fe+fs/(4*D),fs/2.2);
    
19. 
    
20. yf=D*fl;                       %移频量
    
21. df=fs/D/N;
    
22. f=fl:df:fl+(N/2-1)*df;
    
23. xz=zeros(1,N/2);
    
24. wl=2*pi*fl/fs;
    
25. wh=2*pi*fh/fs;
    
26. hr(1)=(wl-wh)/pi;
    
27. hr(2:M+1)=(sin(wl*k)-sin(wh*k))./(pi*k).*w;
    
28. hi(1)=0;
    
29. hi(2:M+1)=(cos(wl*k)-cos(wh*k))./(pi*k).*w;
    
30. 
    
31. k=0:N-1;
    
32. w=0.5-0.5*cos(2*pi*k/N);
    
33. 
    
34. for i=1:L
    
35.     for k=1:N
    
36.         kk=(k-1)*D+M+(i-1)*N;
    
37.         xrz(k)=x(kk+1)*hr(1)+sum(hr(2:M+1).*(x(kk+2:kk+M+1)+x(kk:-1:kk-M+1)));
    
38.         xiz(k)=x(kk+1)*hi(1)+sum(hi(2:M+1).*(x(kk+2:kk+M+1)-x(kk:-1:kk-M+1)));
    
39.     end
    
40.     xzt=(xrz+j*xiz).*exp(-j*2*pi*(0:N-1)*yf/fs);
    
41.     xzt=xzt.*w;
    
42.     xzt=xzt-sum(xzt)/N;
    
43.     xzt=fft(xzt);
    
44.     xz=xz+(abs(xzt(1:N/2))/N*2).^2;
    
45. end
    
46. xz=(xz/L).^0.5;

复制代码

yangzj

呵呵,这两天在弄labview
正在编一个细化解调
对于labview里面的数据组织还不是很明白,因此做数据处理的时候有点掌控不了
正研究

这个很有挑战性，呵呵！

[ 本帖最后由 realhappy 于 2006-11-7 19:41 编辑 ]

yangzj

呵呵,搞定

realhappy

原帖由 yangzj 于 2006-11-7 21:28 发表
呵呵,搞定

恭喜，呵呵，不容易啊！两种程序能传上来，大家学习学习！

marere_re

REALHAPPY斑竹,其实要讲求方便,其实NI的软件做的很好的,他的那个高级信号处理包里,什么细化谱,倍频谱的函数都有的,还有时频分析,只不过这方面的资料太少,研究起来有点困难而已

FtpAdmin

原帖由 yangzj 于 2006-11-7 21:28 发表
呵呵,搞定

希望能够共享一下

realhappy

原帖由 marere_re 于 2006-11-8 09:18 发表
REALHAPPY斑竹,其实要讲求方便,其实NI的软件做的很好的,他的那个高级信号处理包里,什么细化谱,倍频谱的函数都有的,还有时频分析,只不过这方面的资料太少,研究起来有点困难而已

我有高级信号处理包5.0的，里面的细化谱等都是很复杂，运行效率比较底，小波分析也不是很好用，你给的小波降噪还运行不了，如果自己有时间的话，最好还是自己写程序，一来可以按照自己的要求写，二来可以很快熟悉算法和使用要求。

pigmin

入门中多赐教

syg9696

目前还没有学到那么深，努力！