Matla系统辨识仿真CH2,CH5,CH7程序注释与剖析.docVIP

第二部分程序注释及运行结果 读者须知: 为了便于读者理解，现将光盘上第一部分可直接在MATLAB6.I下运行的MATLAB程序的编号和书本上的内容对应如下，每个程序题目括号内的file.m是对应书本上的内容在光盘上第一部分的程序编号。 第二章的随机序列产生程序 例2.1 用乘同余法产生随机数 (见光盘 FLch2sjxleg1.m) ①编程如下： A=6; N=100; %初始化； x0=1; M=255; for k=1:N %乘同余法递推100次； x2=A*x0; %x2和x0分别表示xi和xi-1； x1=mod (x2,M); %将x2存储器的数除以M，取余数放x1（xi）中； v1=x1/256; %将x1存储器的数除以256得到小于1的随机数放v1中； v(:,k)=v1; % 将v1中的数（）存放在矩阵存储器v的第k列中，v(:,k) %表示行不变、列随递推循环次数变化； x0=x1; %xi-1= xi； v0=v1; end %递推100次结束； v2=v %该语句末无‘；’，实现矩阵存储器v中随机数放在v2中，%且可直接显示在MATLAB的window