1、系統(tǒng)辨識與MATLAB仿真附錄例2.1用乘同余法產(chǎn)生隨機數(shù) (見光盤FLch2sjxleg1 .m)編程如下：A=6;N=100;x0=1;M=255;%初始化；for k=1:N x2=A*x0;x仁 mod (x2,M); v仁 X1/256;%乘冋余法遞推100次;%x2和x0分別表示Xi和xi-1 ;%將x2存儲器的數(shù)除以M，取余數(shù)放x1 (xi)中；%將x1存儲器的數(shù)除以256得到小于1的隨機數(shù)放v1中；v(:,k)=v1;%將v1中的數(shù)(i )存放在矩陣存儲器v的第k列中，v(:,k)%表示行不變、列隨遞推循環(huán)次數(shù)變化；x0=x1;v0=v1;%Xi-1= Xi ；end v2=v

2、%遞推100次結(jié)束；%該語句末無；'實現(xiàn)矩陣存儲器v中隨機數(shù)放在v2中， 且可直接顯示在 MATLAB的window中；k仁k;%grapher k=1:k1; plot(k,v,k,v,'r');%以下是繪圖程序；xlabel('k'), ylabel('v');tktle('(0-1)均勻分布的隨機序列') 程序運行結(jié)果如圖2.5所示。(0.1均為-分布的葫枷序列圖2.5采用MATLAB產(chǎn)生的(0,1)均勻分布的隨機序列圖 產(chǎn)生的(0-1)均勻分布的隨機序列在程序運行結(jié)束后，產(chǎn)生的(0 , 1)均勻分布的隨機序列，直接

3、從MATLAB的window界面中copy出來如下(v2中每行存6個隨機數(shù))：v2 =0.02340.14060.84380.08200.49220.96090.78520.72660.37500.25780.55080.31640.90230.43360.60940.66800.02340.14060.84380.08200.49220.96090.78520.72660.37500.25780.55080.31640.90230.43360.60940.66800.02340.14060.84380.08200.49220.96090.78520.72660.37500.25780.550

4、80.31640.90230.43360.60940.66800.02340.14060.84380.08200.49220.96090.78520.72660.37500.25780.55080.31640.90230.43360.60940.66800.02340.14060.84380.08200.49220.96090.78520.72660.37500.25780.55080.31640.90230.43360.60940.66800.02340.14060.84380.08200.49220.96090.78520.72660.37500.25780.55080.31640.902

5、30.43360.60940.66800.02340.14060.84380.0820第二章的白噪聲產(chǎn)生程序例22用乘同余法產(chǎn)生(見光盤FLch2bzsheg2 .m) 編程如下：A=6; x0=1; M=255; f=2;N=100 ；%初始化;x0=1;M=255;for k=1: N x2=A*xO;x仁 mod (x2,M); v仁 x1/256;%乘冋余法遞推100次；%分別用x2和x0表示xi+1和X/ ；%取x2存儲器的數(shù)除以M的余數(shù)放x1 (Xi)中；%將x1存儲器中的數(shù)除以256得到小于1的隨機數(shù)放v1中；v(:,k)=(v1-0.5 )*f;%將v1中的數(shù)(J減去0.5再乘

6、以存儲器f中的系數(shù)，存放 在矩陣存儲器 v的第k列中，v(:,k)表示行不變、列隨遞推 循環(huán)次數(shù)變化；x0=x1;v0=v1;% Xi-1= Xi ；end v2=v%遞推100次結(jié)束；%該語句后無；'實現(xiàn)矩陣存儲器v中隨機數(shù)放在v2中, 且可直接顯示在 MATLAB的window中；k1=k;%grapher k=1:k1; plot(k,v,k,v,'r');%以下是繪圖程序；xlabel('k'), ylabel('v');tktle(' (-1,+1) 均勻分布的白噪聲') 程序運行結(jié)果如圖2.6所示。均勻分布的白

7、毘言圖2.6采用MATLAB產(chǎn)生的(-1,+1)均勻分布的白噪聲序列 產(chǎn)生的(-1，1)均勻分布的白噪聲序列在程序運行結(jié)束后，產(chǎn)生的(-1 ,1)均勻分布的白噪聲序列， 直接從MATLAB的window 界面中copy出來如下(v2中每行存6個隨機數(shù))：v2 =-0.9531-0.71880.6875-0.8359-0.01560.92190.57030.4531-0.2500-0.48440.1016-0.36720.8047-0.13280.21880.3359-0.9531-0.71880.6875-0.8359-0.01560.92190.57030.4531-0.2500-0.484

8、40.1016-0.36720.8047-0.13280.21880.3359-0.9531-0.71880.6875-0.8359-0.01560.92190.57030.4531-0.2500-0.48440.1016-0.36720.8047-0.13280.21880.3359-0.9531-0.71880.6875-0.8359-0.01560.92190.57030.4531-0.2500-0.48440.1016-0.36720.8047-0.13280.21880.3359-0.9531-0.71880.6875-0.8359-0.01560.92190.57030.4531-

9、0.2500-0.48440.1016-0.36720.8047-0.13280.21880.3359-0.9531-0.71880.6875-0.8359-0.01560.92190.57030.4531-0.2500-0.48440.1016-0.36720.8047-0.13280.21880.3359-0.9531-0.71880.6875-0.8359*另外，書中圖2.3白噪聲的產(chǎn)生如下：顯然，只要在例22程序的初始化部分中給 N=300, f=6，運行程序就可以得到如圖 2.3 所示的(-3，3)的白噪聲過程. 編程如下：A=6; xO=1; M=255; f=6;N=300 ；%

10、初始化；xO=1;M=255;for k=1: N x2=A*xO;x仁 mod (x2,M); v仁 x1/256;%乘冋余法遞推100次；%分別用x2和x0表示Xi+1和X/ ；%取x2存儲器的數(shù)除以M的余數(shù)放x1 (Xi)中；%將x1存儲器中的數(shù)除以256得到小于1的隨機數(shù)放v1中；v(:,k)=(v1-0.5 )*f;%將v1中的數(shù)(i)減去0.5再乘以存儲器f中的系數(shù)，存放 在矩陣存儲器 v的第k列中，v(:,k)表示行不變、列隨遞推 循環(huán)次數(shù)變化；x0=x1;v0=v1;% Xi_1= Xi ；end v2=v%遞推100次結(jié)束；%該語句后無；'實現(xiàn)矩陣存儲器v中隨機數(shù)放在

11、v2中， 且可直接顯示在 MATLAB的window中；k仁k;%grapher k=1:k1; plot(k,v,k,v,'r');%以下是繪圖程序；xlabel('k'), ylabel('v');tktle(' (-1,+1) 均勻分布的白噪聲') 程序運行結(jié)果如圖2.3所示。圖2.3白噪聲過程第二章的M序列產(chǎn)生程序例2.3用移位寄存器 產(chǎn)生M序列的MATLAB軟件實現(xiàn)（見光盤FLch2bzsheg3 .m）編程如下：X1=1;X2=0;X3=1;X4=0;%移位寄存器輸入 Xi初T態(tài)（0101）， Yi為移位寄存器各級輸出

12、m=60;%置M序列總長度for i=1:m%1#Y4=X4;Y3=X3;Y2=X2;Y1=X1;X4=Y3;X3=Y2;X2=Y1;X1= xor(Y3,Y4);% 異或運算if Y4=0U(i)=-1;elseU(i)=Y4;endendM=U%繪圖i1=ik=1:1:i1;plot(k,U,k,U,'rx') xlabel('k') ylabel('M 序列')title（'移位寄存器產(chǎn)生的 M序列'）程序運行結(jié)果如圖2.8所示。移拉寄存器產(chǎn)生的M序列110.80.60402國 n-0.2-0.4-06武 E *7*7-0.8*寸0102030405060k圖2.8軟件實現(xiàn)的移位寄存器產(chǎn)生的M序列圖 四級移位寄存器產(chǎn)生的 M序列M =Columns 1 through 106系統(tǒng)辨識與MATLAB仿真附錄#系統(tǒng)辨識與MATLAB仿真附錄-11-1111 1 -1 -1 -1#系統(tǒng)辨識與MATLAB仿真附錄#系統(tǒng)辨識與MAT