



全文预览已结束
下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
递推的极大似然法辨识程序例5.2系统模型如图5.5所示。试用递推的极大似然法对系统辨识的参数集 图5.5 例5.2系统模型u(k)+ e(k)y(k)+ z(k)v(k) u v(k)随机信号,输入信号为幅值为的M序列或随机信号,要求画出程序流程图,打印出程序(程序中带有注释)和辨识中的参数、误差曲线。解: 首先解释编程所用的部分字母:由于在MATLAB语言中无法用希腊字母描述、无法用上标及下标,故用 o和o1表示;令;产生M序列时,a(i),b(i),c(i),d(i)表示四级移位寄存器的第1,2,3,4级寄存器的输出; 编程如下(光盘上该程序:FLch5RMLeg2.m,可在MATLAB6.I下直接运行):编程如下:clear %清零a(1)=1;b(1)=0;c(1)=1;d(1)=0;u(1)=d(1);z(1)=0;z(2)=0; %初始化for i=2:1200 %产生m序列u(i) a(i)=xor(c(i-1),d(i-1); b(i)=a(i-1); c(i)=b(i-1); d(i)=c(i-1); u(i)=d(i);endu; %若取去;可以在程序运行中观测到m序列v=randn(1200,1); %产生正态分布随机数V=0; %计算噪声方差for i=1:1200 V=V+v(i)*v(i);endV1=V/1200; for k=3:1200 %根据v和u计算z z(k)=1.2*z(k-1)-0.6*z(k-2)+u(k-1)+0.5*u(k-2)+v(k)-v(k-1)+0.2*v(k-2);endo1=0.001*ones(6,1);p0=eye(6,6); %赋初值 zf(1)=0.1;zf(2)=0.1;vf(2)=0.1;vf(1)=0.1;uf(2)=0.1;uf(1)=0.1; %迭代计算参数值和误差值for k=3:1200 h=-z(k-1);-z(k-2);u(k-1);u(k-2);v(k-1);v(k-2); hf=h; K=p0*hf*inv(hf*p0*hf+1); p=eye(6,6)-K*hf*p0; v(k)=z(k)-h*o1; o=o1+K*v(k) ; p0=p; o1=o; a1(k)=o(1); a2(k)=o(2); b1(k)=o(3); b2(k)=o(4); d1(k)=o(5); d2(k)=o(6); e1(k)=abs(a1(k)+1.2); e2(k)=abs(a2(k)-0.6); e3(k)=abs(b1(k)-1.0); e4(k)=abs(b2(k)-0.5); e5(k)=abs(d1(k)+1.0); e6(k)=abs(d2(k)-0.2); zf(k)=z(k)-d1(k)*zf(k-1)-d2(k)*zf(k-2); uf(k)=u(k)-d1(k)*uf(k-1)-d2(k)*uf(k-2); vf(k)=v(k)-d1(k)*vf(k-1)-d2(k)*vf(k-2); hf=-zf(k-1);-zf(k-2);uf(k-1);uf(k-2);vf(k-1);vf(k-2); end o1 %若取去;可以在程序运行中观测到参数 V1 %绘图 subplot(4,1,1)k=1:1200;plot(k,a1,k:,k,a2,b,k,b1,r,k,b2,m:,k,d1,g,k,d2,k);xlabel(k)ylabel(parameter)legend(a1=-1.2,a2=0.6,b1=1.0,b2=0.5,d1=-1.0,d2=0.2); %图标炷title(The parameter idendification of the RML);endsubplot(4,1,2)k=1:1200;plot(k,e1,k,k,e2,b,k,e3,r,k,e4,m,k,e5,g,k,e6,k);xlabel(k)ylabel(error)%title(误差曲线)endsubplot(4,1,3)k=1:1200;plot(k,u);xlabel(k)ylabel(input)%title(系统输入信号)endsubplot(4,1,4)k=1:1
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 桂花蜜饯铺创新创业项目商业计划书
- 读书笔记演讲稿模版模板
- 夏威夷果果汁制造创新创业项目商业计划书
- 果蔬罐头厂创新创业项目商业计划书
- 酱油酱类制作工特殊工艺考核试卷及答案
- 时尚科技融合加速器创新创业项目商业计划书
- 智能垃圾分类系统孵化器创新创业项目商业计划书
- 丁辛醇装置操作工晋升考核试卷及答案
- 乡镇医院安全知识培训课件
- 2025-2030年新能源汽车电池回收与再生利用产业链投资机会报告
- Unit 2 Home Sweet Home 重点短语和句式-人教版英语八年级上册
- 黄体破裂护理常规课件
- 防治大气污染课件
- 环境监测质量管理课件
- 国际音标教学课件
- 2025-2030中国可变磁阻旋转变压器行业产销状况与应用趋势预测报告
- (高清版)DB31∕T 310001-2020 船舶水污染物内河接收设施配置规范
- 诊所治疗室管理制度
- 2025年高考英语全国二卷听力试题真题及答案(含MP3+原文)
- 常用概率分布题目及答案
- 2025年互联网营销师考试题及答案
评论
0/150
提交评论