为什么用matlab画出的高斯白噪声功率谱密度不是直线 而是如下的图，按理论来讲最后一个应该是直线，可是并非如此，请问为什么？是啊，理论归理论，现实归现实。。
谁会matlab编程显示功率谱密度呢，跪求。 t=0：0.0001：0.1；时间间隔为0.0001，说明采样频率为10000Hzx=square(2*pi*1000*t)；产生基e68a8462616964757a686964616f31333239303261频为1000Hz的方波信号n=randn(size(t))；白噪声f=x+n；在信号中加入白噪声figure(1)；subplot(2，1，1)；plot(f)；画出原始信号的波形图ylabel('幅值(V)')；xlabel('时间(s)')；title('原始信号')；y=fft(f，1000)；对原始信号进行离散傅里叶变换，参加DFT采样点的个数为1000subplot(2，1，2)；m=abs(y)；f1=(0：length(y)/2-1)'*10000/length(y)；计算变换后不同点对应的幅值plot(f1，m(1：length(y)/2))；ylabel('幅值的模')；xlabel('时间(s)')；title('原始信号傅里叶变换')；用周期图法估计功率谱密度p=y.*conj(y)/1000；计算功率谱密度ff=10000*(0：499)/1000；计算变换后不同点对应的频率值figure(2)；plot(ff，p(1：500))；ylabel('幅值')；xlabel('频率(Hz)')；title('功率谱密度（周期图法）')；
怎么在matlab中仿真功率谱密度为1的高斯白噪声？ 功率谱2113为1，则在一个周期内的功率为2*pi，那5261么由函数wgn产生4102高斯白噪声即可。例如：y=wgn(100，100，2*pi，‘1653linear’)；产生高斯白噪声，100*100为高斯白噪声的矩阵大小，2*pi为其功率，linear表示其输出的幅度为功率形式。则 imshow(y)
