声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2549|回复: 1

[滤波] 关于高斯白噪声频谱的生成

[复制链接]
发表于 2008-5-7 13:55 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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

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

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

附代码

  1. t=0:1:1023;%点数

  2. %%--------生成信号--------%%

  3. y=wgn(1,length(t),1); %波形表达式

  4. figure;

  5. plot(t,y);

  6. title('生成的高斯白噪声');

  7. %%--------画频谱图--------%%

  8. y2=xcorr(y); %自相关

  9. y1=fft(y2); %傅立叶变换

  10. y1=fftshift(y1); %取绝对值

  11. N1=length(y1); %N是y1的长度

  12. f1=-(N/2-0.5):1:(N/2-0.5);

  13. figure;

  14. plot(f1,abs(y1));

  15. title('生成的高斯白噪声频谱图');



  16. %%--------取样测试--------%%

  17. y=abs(y);

  18. d=fix(y*10000);

  19. %%--------画波形图--------%%

  20. figure;

  21. plot(t,d);

  22. title('高斯白噪声取整后');

  23. %%--------转换数制--------%%

  24. for i=1:1024

  25.     b=dec2bin(d(i),16);

  26.     bout(i,:)=b;

  27. end;

  28. %%--------画频谱图--------%%

  29. y3=xcorr(d); %自相关

  30. y4=fft(y3); %傅立叶变换

  31. y4=fftshift(y4); %取绝对值

  32. N2=length(y4); %N是y1的长度

  33. f2=-(N/2-0.5):1:(N/2-0.5);

  34. figure;

  35. plot(f2,abs(y4));

  36. title('高斯白噪声取整后频谱图');

  37. %%--------写入文本--------%%

  38. fid = fopen('input.txt','wt');

  39. for j=1:1024

  40.     fprintf(fid,'%s\n',bout(j,:));

  41. end;

  42. fclose(fid);
复制代码

dectobin.m

900 Bytes, 下载次数: 31

回复
分享到:

使用道具 举报

发表于 2011-4-26 15:17 | 显示全部楼层
没钱下载的
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 09:52 , Processed in 0.079357 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表