在MATLAB中輸入下述語句，調(diào)用randn函數(shù)產(chǎn)生白噪聲信號，并將噪聲信號序列存儲為一個激勵文件。PCI9054-AC50PI每行語句中“%”及后面的文字起注釋作用。

   openfile=fopen(’mystl.txt‘，’w')；%以寫模式打開并創(chuàng)建激勵文件mystl.txt

   t=[1:1:4000l;%噪聲序列中產(chǎn)生4000個數(shù)據(jù)點

   y=0.1；randn(size(t))；%調(diào)用randn函數(shù)產(chǎn)生白噪聲序列

   for i=l:length(t)

   fprintf(openfle，’+(%de．10%f)＼r＼nI，i，y(i))；%將數(shù)據(jù)寫到mystl.txt文件中

   end

   fciose(openfile)

   (2)在MATLAB命令窗中運行上述語句后，產(chǎn)生描述白噪聲信號的數(shù)據(jù)文件mystl.txt。

   該文件中前10組數(shù)據(jù)如下：

   +(le-10 -0.102491)

   +(2e-10 -0.011376)

   +(3e-10 -0.167544)

   +(4e-10 -0.055026)

   +(5e-10 -0.068639)

   +(6e-10 0.103039)

   +(7e-10 -0.117071)

   +(8e-10 -0.126143)

   +(9e-10 -0.065016)

   +(lOe-10 0.106448)

   步驟四：將MATLAB產(chǎn)生的信號數(shù)據(jù)復制到步驟二生成的激勵信號文件中

   打開noise.stl文件，將步驟三中產(chǎn)生的mystl.txt文件中描述波彤的數(shù)據(jù)添加到noise.stl文件中，作為圖8-29中描述信號波形時間和信號值的坐標數(shù)據(jù)。為了便于比較添加噪聲后的仿真結(jié)果，本例中將默認的(0，0)作為噪聲信號的初始值。添加了MATLAB中產(chǎn)生數(shù)據(jù)之后的noise.stl文件如圖8-30所示，其中只列出了前11個數(shù)據(jù)坐標。