关于高斯白噪声频谱的生成
我想把高斯白噪声导成文本，，存成二进制不知道前后两个噪声的频谱都正确与否？？
我感觉频谱应给差不多啊～～
附代码
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); 没钱下载的
页:
[1]