关于王济、胡晓《matlab在振动信号处理中的应用》书中白噪声激励

左言

关于王济、胡晓《matlab在振动信号处理中的应用》书中白噪声激励程序的问题，我按照书中给出的程序运行，生成的图形与书上的怎么不一样呢，请给位大虾帮忙看一下

左言

这是书中的给出的图形

左言

这是我按照书中给出的程序生成的图形

impulse

你的数据哪来的？

Cena马

每次产生的噪声不一样，分析结果当然也不一样了

左言

回复 Cena马 的帖子

这个我知道，但是两个图中的频率范围是不一样的

Cena马

回复 左言 的帖子

你采样频率多少？

左言

回复 Cena马 的帖子

采样频率是400

左言

1. %生成白噪声信号
   
2. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   
3. clear
   
4. clc
   
5. close all hidden
   
6. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   
7. % fni='out.txt';
   
8. fid=fopen('out.txt','r');
   
9. sf=fscanf(fid,'%f',1);%采样频率
   
10. fi=fscanf(fid,'%f',1);%最小截止频率
    
11. fa=fscanf(fid,'%f',1);%最大截止频率
    
12. t1=fscanf(fid,'%f',1)%时间长度（秒）
    
13. am=fscanf(fid,'%f',1);%最大幅值
    
14. fno=fscanf(fid,'%f',1);%输出数据文件名;
    
15. status=fclose(fid);
    
16. %计算白噪声信号数据长度n
    
17. n=round(sf*t1)+1
    
18. t=0:1/sf:(n-1)/sf;%建立信号离散时间向量t
    
19. nfft=2^nextpow2(n);%大于并接近n的2的幂次方为FFt长度
    
20. ni=round(fi*nfft/sf+1);%四舍五入取整求最小截止频率对应数组元素的下标
    
21. na=round(fa*nfft/sf+1);%四舍五入取整求最大截止频率对应数组元素的下标
    
22. r=zeros(1,nfft/2);%建立长度为nfft/2元素全为0的向量
    
23. r(ni:na)=1;%将r的正频带通内的元素赋值为1，生成幅值谱
    
24. 
    
25. a=2*pi*rand(1,nfft/2);%生成0~2pi的随机数为随机相位
    
26. b=r.*exp(i*a)*nfft/2;%将幅值谱和相位谱转换为实部和虚部
    
27. a=[b,b(nfft/2:-1:1)];%将正负圆频率向量合成一个向量
    
28. x=real(ifft(a,nfft));%IFFT变换，并取实部为生成的白噪声信号
    
29. y=am*x(1:n)/max(abs(x(1:n)));%去指定长度的数据并使最大值为指定幅值
    
30. 
    
31. nft=1024;%定义自功率谱的FFt长度
    
32. y1=psd(y,nft,sf,boxcar(nft),nft/2);%计算生成的白噪声信号的自功率谱
    
33. f=0:sf/nft:(nft/2-1)*sf/nft;%定义自功率谱的离散频率向量
    
34. 
    
35. %绘制白噪声信号随时间变化的图形
    
36. subplot(2,1,1);
    
37. plot(t,y);
    
38. xlabel('时间 （s）');
    
39. ylabel(' 幅值 ');
    
40. grid on
    
41. 
    
42. %绘制自功率谱图形
    
43. subplot(2,1,2);
    
44. plot(f(1:nft/4),y1(1:nft/4));
    
45. xlabel('频率 (HZ)');
    
46. ylabel('幅值');
    
47. grid on;
    
48. save force y
    
49. % %打开文件输出生成的白噪声信号数据
    
50. % fid=fopen(fno,'w');
    
51. % for k=1:n
    
52. %     fprintf(fid,'%f %f/n',t(k),y(k));%每一行输出两个实型数据，t为时间，y为白噪声信号值
    
53. % end
    
54. % status=fclose(fid);
    

复制代码

这是书中给出的程序

VibrationMaster

这个程序需要输入频率参数。你的输入对吗？

左言

回复 VibrationMaster 的帖子

恩，就是按照书中给出的输入的，频率范围为0.5-90Hz，时间长度40s，最大幅值1.1

rogen

我也出现同样的问题，不知道怎么解决，还有其中有一句：“ni=round(fi*nfft/sf+1);%四舍五入取整求最小截止频率对应数组元素的下标",这句话怎么理解啊？