关于高斯白噪声频谱的生成

HolySaint

我想把高斯白噪声导成文本，，存成二进制

不知道前后两个噪声的频谱都正确与否？？

我感觉频谱应给差不多啊～～

附代码

1. t=0:1:1023;%点数
   
2. 
   
3. %%--------生成信号--------%%
   
4. 
   
5. y=wgn(1,length(t),1); %波形表达式
   
6. 
   
7. figure;
   
8. 
   
9. plot(t,y);
   
10. 
    
11. title('生成的高斯白噪声');
    
12. 
    
13. %%--------画频谱图--------%%
    
14. 
    
15. y2=xcorr(y); %自相关
    
16. 
    
17. y1=fft(y2); %傅立叶变换
    
18. 
    
19. y1=fftshift(y1); %取绝对值
    
20. 
    
21. N1=length(y1); %N是y1的长度
    
22. 
    
23. f1=-(N/2-0.5):1:(N/2-0.5);
    
24. 
    
25. figure;
    
26. 
    
27. plot(f1,abs(y1));
    
28. 
    
29. title('生成的高斯白噪声频谱图');
    
30. 
    
31. 
    
32. 
    
33. %%--------取样测试--------%%
    
34. 
    
35. y=abs(y);
    
36. 
    
37. d=fix(y*10000);
    
38. 
    
39. %%--------画波形图--------%%
    
40. 
    
41. figure;
    
42. 
    
43. plot(t,d);
    
44. 
    
45. title('高斯白噪声取整后');
    
46. 
    
47. %%--------转换数制--------%%
    
48. 
    
49. for i=1:1024
    
50. 
    
51.     b=dec2bin(d(i),16);
    
52. 
    
53.     bout(i,:)=b;
    
54. 
    
55. end;
    
56. 
    
57. %%--------画频谱图--------%%
    
58. 
    
59. y3=xcorr(d); %自相关
    
60. 
    
61. y4=fft(y3); %傅立叶变换
    
62. 
    
63. y4=fftshift(y4); %取绝对值
    
64. 
    
65. N2=length(y4); %N是y1的长度
    
66. 
    
67. f2=-(N/2-0.5):1:(N/2-0.5);
    
68. 
    
69. figure;
    
70. 
    
71. plot(f2,abs(y4));
    
72. 
    
73. title('高斯白噪声取整后频谱图');
    
74. 
    
75. %%--------写入文本--------%%
    
76. 
    
77. fid = fopen('input.txt','wt');
    
78. 
    
79. for j=1:1024
    
80. 
    
81.     fprintf(fid,'%s\n',bout(j,:));
    
82. 
    
83. end;
    
84. 
    
85. fclose(fid);

复制代码

zhaohua5269

没钱下载的