%零均值的白噪聲輸入
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');