f14_4.m

《Matlab时频分析及其应用》的详细代码

%零均值的白噪声输入
load wnoise.dat -ascii
N=length(wnoise);
w=wnoise;
x=w;
%路面不平度的模拟
for i=8:N
    x(i)=1.0658*x(i-1)-0.4533*x(i-2)+0.3940*x(i-3)-0.1735*x(i-4)+...
        0.1375*x(i-5)-0.0551*x(i-6)+0.0799*x(i-7)+w(i);
end
%路面不平度输出
x=(x-mean(x));
%显示白噪声
subplot(2,1,1);
plot(1:N,w);
xlabel('样本序号 n');
ylabel('幅值 A');
%显示路面不平度
subplot(2,1,2);
plot(1:N,x);
xlabel('样本序号 n');
ylabel('幅值 A');