Matlb_系统辨识_仿真_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中； k1=k; %grapher %以下是绘图程序； k=1:k1; plot(k,v,k,v,r); xlabel(k), ylabel(v);tktle((0-1)均匀分布的随机序列) ② 程序运行结果如图2.5所示。 图2.5 采用MATLAB产生的(0,1)均匀分布的随机序列图 ③ 产生的(0-1)均匀分布的随机序列 在程序运行结束后，产生的(0，1)均匀分布的随机序列，直接从MATLAB的window界面中copy出来如下（v2中每行存6个随机数）： v2 = 0.0234 0.1406 0.8438 0.0820 0.4922 0.9609 0.7852 0.7266 0.3750 0.2578 0.5508 0.3164 0.9023 0.4336 0.6094 0.6680 0.0234 0.1406 0.8438 0.0820 0.4922 0.9609 0.7852 0.7266 0.3750 0.2578 0.5508 0.3164 0.9023 0.4336 0.6094 0.6680 0.0234 0.1406 0.8438 0.0820 0.4922 0.9609 0.7852 0.7266 0.3750 0.2578 0.5508 0.3164 0.9023 0.4336 0.6094 0.6680 0.0234 0.1406 0.8438 0.0820 0.4922 0.9609 0.7852 0.7266 0.3750 0.2578 0.5508 0.3164 0.9023 0.4336 0.6094 0.6680 0.0234 0.1406 0.8438 0.0820 0.4922 0.9609 0.7852 0.7266 0.3750 0.2578 0.5508 0.3164 0.9023 0.4336 0.6094 0.6680 0.0234 0.1406 0.8438 0.0820 0.4922 0.9609 0.7852 0.7266 0.3750 0.2578 0.5508 0.3164 0.9023 0.4336 0.6094 0.6680 0.0234 0.1406 0.8438 0.0820 第二章的白噪声产生程序 例2.2 用乘同余法产生(见光盘 FLch2bzsheg2.m) 编程如下： A=6; x0=1; M=255; f=2; N=100； %初始化； x0=1; M=255; for k=1: N %乘同余法递推100次； x2=A