|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我想把高斯白噪声导成文本,,存成二进制
不知道前后两个噪声的频谱都正确与否??
我感觉频谱应给差不多啊~~
附代码
- t=0:1:1023;%点数
- %%--------生成信号--------%%
- y=wgn(1,length(t),1); %波形表达式
- figure;
- plot(t,y);
- title('生成的高斯白噪声');
- %%--------画频谱图--------%%
- y2=xcorr(y); %自相关
- y1=fft(y2); %傅立叶变换
- y1=fftshift(y1); %取绝对值
- N1=length(y1); %N是y1的长度
- f1=-(N/2-0.5):1:(N/2-0.5);
- figure;
- plot(f1,abs(y1));
- title('生成的高斯白噪声频谱图');
- %%--------取样测试--------%%
- y=abs(y);
- d=fix(y*10000);
- %%--------画波形图--------%%
- figure;
- plot(t,d);
- title('高斯白噪声取整后');
- %%--------转换数制--------%%
- for i=1:1024
- b=dec2bin(d(i),16);
- bout(i,:)=b;
- end;
- %%--------画频谱图--------%%
- y3=xcorr(d); %自相关
- y4=fft(y3); %傅立叶变换
- y4=fftshift(y4); %取绝对值
- N2=length(y4); %N是y1的长度
- f2=-(N/2-0.5):1:(N/2-0.5);
- figure;
- plot(f2,abs(y4));
- title('高斯白噪声取整后频谱图');
- %%--------写入文本--------%%
- fid = fopen('input.txt','wt');
- for j=1:1024
- fprintf(fid,'%s\n',bout(j,:));
- end;
- fclose(fid);
复制代码 |
|