




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、M序列是工程中常用的输入信号,它的性质类似于白噪声,而白噪声是理论上最好的输入信号,可见M序列的价值。下面介绍M序列的matlab产生方法。 idinput函数产生系统辨识常用的典型信号。格式u = idinput(N,type,band,levels)u,freqs = idinput(N,'sine',band,levels,sinedata)N产生的序列的长度,如果N=N nu,则nu为输入的通道数,如果N=P nu M,则nu指定通道数,P为周期,M*P为信号长度。默认情况下,nu=1,M=1,即一个通道,一个周期。Type指定产生信号的类型,可选类型如下rg
2、s高斯随机信号rbs (默认)二值随机信号prbs二值伪随机信号(M序列)sine正弦信号和Band指定信号的频率成分。对于rgs、rbs、sine,band = wlow, whigh指定通带的范围,如果是白噪声信号,则band=0, 1,这也是默认值。指定非默认值时,相当于有色噪声。对于prbs,band=0, B,B表示信号在一个间隔1/B(时钟周期)内为恒值,默认为0, 1。Levels指定输入的水平。Levels=minu, maxu,在type=rbs、prbs、sine时,表示信号u的值总是在minu和maxu之间。对于type=rgs,minu指定信号的均值减标准差,maxu指
3、定信号的均值加标准差,对于0均值、标准差为1的高斯白噪声信号,则levels=-1, 1,这也是默认值。说明对于PRBS信号,如果M>1,则序列的长度和PRBS周期会做调整,使PRBS的周期为对应一定阶数的最大值(即2n-1,n为阶数);如果M=1,PRBS的周期是大于N的相应阶数的值。在多输入的情形时,信号被最大平移,即P/nu为此信号能被估计的模型阶次的上界。上面的意思可如下理解:对于M=1时,ms = idinput(12, 'prbs', 0 1, 0 1);figurestairs(ms)title('M序列')ylim(-0.5 1.5)结果如
4、下同时,matlab给出如下警告Warning: The PRBS signal delivered is the 12 first values of a full sequence of length 15.即函数的输出为周期为15(大于12的第一个2n-1的值)PRBS信号的前12个值组成的序列。如ms = idinput(15, 'prbs', 0 1, 0 1);figurestairs(ms)title('M序列')ylim(-0.5 1.5)可以看到指定12时的序列为指定15时的序列的前面部分。 对于M>1时,ms = idinpu
5、t(12,1,2, 'prbs', 0 1, 0 1);figurestairs(ms)title('M序列')ylim(-0.5 1.5)结果如下 Matlab给出的响应警告为Warning: The period of the PRBS signal was changed to 7. Accordingly, the length of thegenerated signal will be 14.对于正弦信号和的产生,貌似用的不多,语法还挺复杂,等用的时候再看吧。 方法产生rgs信号的带通信号使用的是一个8阶巴特沃斯滤波器,使用idf
6、ilt做的非因果滤波,这个是可信赖的方法。对于rbs信号,使用的是相同的滤波器,但是是在二值化之前,这意味着频率成分并不保证是精确的。 产生高斯随机信号clcclear allclose all% 高斯随机信号u = idinput(1000, 'rgs');figurestairs(u)title('高斯随机信号')figurehist(u, -4:4)title('高斯随机信号的分布')产生二值随机信号clcclear allclose all% 二值随机信号u = idinput(100, 'rbs');figur
7、estairs(u)title('二值随机信号')ylim(-1.5 1.5)产生二值伪随机信号(M序列)合理的选择输入激励信号,能有效的激励起系统的动态信号。白噪声的平稳谱的性质决定了它是一个很好的输入信号,但它在工程中不易实现,而M序列具有近似白噪声的性质,可保证良好的辨识精度。clcclear allclose all% 二值伪随机信号(M序列)n = 8; % 阶次p = 2n -1; % 循环周期ms = idinput(p, 'prbs');figurestairs(ms)title('M序列')ylim(-1.5 1.5)结果验证M
8、序列的性质如下· -1和1的个数差1sum(ms=1) % 1的个数sum(ms=-1) % -1的个数ans =127ans =128· 存在直流分量mean(ms) % 直流分量ans =-0.0039· 相关函数a = zeros(length(ms)*10, 1); % 采样for i = 1:10a(i:10:end) = ms;endc = xcorr(a, 'coeff'); % 自相关函数figureplot(c)title('相关函数')自相关函数接近于函数。· 谱密度figurepwelch(a) %
9、谱密度说明M序列不含基频的整数倍的频率成分。 产生逆M序列谱分析表明,M序列含有直流成分,将造成对辨识系统的“净扰动”,这通常不是所希望的。而逆M序列将克服这一缺点,是一种比M序列更为理想的伪随机码序列。clcclear allclose all% 二值伪随机信号(M序列)n = 8; % 阶次p = 2n -1; % 循环周期ms = idinput(p, 'prbs', , 0 1);figurestairs(ms)title('M序列')ylim(-0.5 1.5)% 产生逆M序列s = 0;ims = zeros(2*p, 1);mstemp = ms; ms;for i = 1:2*pims(i) = xor(mstemp(i), s);s = not(s);endims(ims=0) = -1;figurestairs(ims)title('逆M序列')ylim(-1.5 1.5)· -1和1的个数差1sum(ims=1) % 1的个数sum(ims=-1) % -1的个数ans =255ans =255· 无直流分量mean(ims) % 直流分
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 兼职工作合同协议
- 消防系统检测合同
- 小数的意义(教学设计)-2023-2024学年四年级下册数学人教版
- 管理软件系统购买合同范文格式7篇
- 吨的认识(教学设计)-2024-2025学年三年级上册数学人教版
- 双手胸前传接球 教学设计-2023-2024学年高二下学期体育与健康人教版必修第一册
- 小学三年级数学几百几十加减几百几十水平练习习题
- 简易家用活动平台施工方案
- Unit 1 Lesson 3 The Sun Is Rising教学设计 -2024-2025学年冀教版八年级英语下册
- 第9课 两宋的政治和军事 教学设计-2023-2024学年高一上学期统编版(2019)必修中外历史纲要上
- 南宋北京大学历史学系课件
- 重庆市房屋建筑与装饰工程计价定额2018-建筑工程
- 三年级数学-解决问题策略(苏教版)
- 不吃路边摊精品课件
- 《网络服务器搭建、配置与管理-Linux(RHEL8、CentOS8)(微课版)(第4版)》全册电子教案
- 心理评估与诊断简介
- 无痛病房管理课件
- 让孩子变成学习的天使——由《第56号教室的奇迹》读书分享
- 球泡检验标准
- 公安笔录模板之询问嫌疑人(书面传唤治安案件)
- 振动分析基础讲义1
评论
0/150
提交评论