求助：请高手帮忙看下这个白噪声信号生成程序

lvpin

程序如下

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear
clc
close all hidden
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fni=input('生成白噪声信号-输入数据文件名：','s');
fid=fopen(fni,'r');
sf=fscanf(fid,'%f',1); 采样频率
fi=fscanf(fid,'%f',1);  最小截止频率
fa=fscanf(fid,'%f',1); 最大
tl=fscanf(fid,'%f',1);  时间
am=fscanf(fid,'%f',1); 最大幅值
fno=fscanf(fid,'%s',1);  输出数据文件名
fclose(fid);
%计算白噪声信号数据长度n
n=round(sf*tl)+1;
%建立信号离散时间向量t
t=0:1/sf:(n-1)/sf;
%大于并最接近n的2的幂次方为fft长度
nfft=2^nextpow2(n);
%四舍五入取整求最小截止频率对应数组元素的下标
ni=round(fi*nfft/sf+1);
%四舍五入取整求最大截止频率对应数组元素的下标
na=round(fa*nfft/sf+1);
%建立一个长度为nfft/2元素全为0的向量
r=zeros(1,nfft/2);
%将r的正频率带通内的元素赋值1，生成赋值铺
r(ni:na)=1;
%生成0~2pi的随机数为随机相位
a=2*pi*rand(1,nfft/2);
%将幅值谱和相位谱转换成实部和虚部
b=r.*exp(i*a)*nfft/2;
%将正负圆频率向量组合成一个向量
a=[b,b(nfft/2:-1:1)];
%IFFT变换，并取实部为生成白噪声信号
x=real(ifft(a,nfft));
%取指定长度的数据并使数据的最大值为指定的幅值
y=am*(1:n)/max(abs(x(1:n)));
%定义自功率谱的FFT长度
nft=1024;
%计算生成白噪声信号自功率谱
y1=psd(y,nft,sf,boxcar(nft),nft/2);
%定义自功率谱的离散频率向量
f=0:sf/nft:(nft/2-1)*sf/nft;
%绘制白噪声信号随时间变化的曲线图
subplot(2,1,1);
plot(t,y);
xlabel('时间(s)');
ylabel('幅值');
grid on;
%绘制自功率谱曲线图
subplot(2,1,2);
plot(f(1:nft/4),y1(1:nft/4));
xlabel('频率(Hz)');
ylabel('幅值');
grid on;
fid=fopen(fno,'w');
for k=1:n
    fprintf(fid,'%f%f\n',t(k),y(k));
end
status=fclose(fid);

最后得不到白噪声曲线，看看是不是算法或者程序编写有误，谢谢！