信号处理：频域转时域

%*******时域反演*******%
rng default;  % For reproducibility
Phi = 2*pi*rand(NN,1);  %生成NN个在（0，2*pi）上均匀分布的相互独立的随机变量
fv = zeros(1,NN);
for k = 1:length(LS)
    %采用三角级数以余弦波进行时域样本的模拟
    for i = 1:NN
        fv(k) = fv(k) + av(i)*cos(Omega(i)*LS(k)+Phi(i)); 
    end
end
%*******时域反演结束*******%
subplot(212)  
plot(LS,fv)
xlabel('时间/s')  
ylabel('加速度 m/s')
[f,y0]=change_fft(frecy,fv); %#codegen
figure(3)
plot(f,y0,'r',ff1,F,'b')
xlim([2.5 75])
legend('反演的时域然后傅里叶变换得到的频域','原始频域')
title('频域谱')  
xlabel('频率/hz')  
ylabel('加速度功率密度g^2/hz')