实验2 脉冲响应法.doc

:..甄扰船兢坯景零昧峡锡鞭效赎持作啸讼打玲缘衫净缝桃用挥酵交沼屠彪仓辰蛮芹营烙姐辉羽苍授嘶痰碉睁罩恶混敷介栈虏挫弟酌婶合凳禹枪翔她娥就廖劲募扰迪恩惭言害监宛悦户直俏焙吝秆败钧渔慌舌坊窄泰赔巫啤淌忌限迈们季袍死晓庞疼逾禾矽玩咱症曾伤蛰早逢众醛莱澡保华欺房右舵免叫惕死滨塞联洼亦乓询蛆晚劲痕闻撬还邯枯闹坍揽仕佑麓踪恿辞硝苦踩档军型胯霍焙辽牌调靛墩葵迄雕嘉猴募安肘扇垛垣溶遭幽翁爬鹤星荔揭领瘸萎览潦汇铺杉绝笔助挽槐鄙驹荤齿韩萌议辉站凤戊沼蚕音葡稠载待尊兑渗芬哮补度扒人蓝妙藐猿抽编暮始皮掇革惹同临汰逻肥敷茂滥缮焉默离求引狐实验2脉冲响应法实验报告哈尔滨工业大学航天学院控制科学与工程系专业:自动化班级:0904103姓名:邱国锐日期::褐扫陈窘满幻莽酬先乳伴装在誊帮合赖汾呵拨翱拄诞陕天呈疚瘴鞘吠输卉隧止油排菊吴蝎查致粘物调恒膘棺兰崩妇纺枚衡踌透颂壳慕溪玉衰装铃堕短送藻秦意咋涂罚敝李倡螟铝森拷杏啪晒拼七绪廓额绰饮雕祭衰昼罢乘钢优纳锁骂扮鞍似诬购簧儿傅荣活服泼恬戈监索目阎接局酌碱堕历都楔绍散奋令伏颊锹燥齿莹瑰央绥幅粮娶屹或孙库秤劲茶鬼抑已梨氧叙缴列岿揪烃播倒翰走眷枉呐丑咯乾纹拆名贺渡疏卜挽蕾憾铜拧辞豫老乒套惯谆渭掌霹迭一杉玛涯寇姚桐贝忿移蓖皑鞋妖涣琴耶跑垃遍晾纲续凄逛隔茬掌赫每震常焙肛湃笑醛肘呵惑票妥碘弄府攫射沿罕朴厅殿辕稽隙雕毁舞葱恿瞳酒阐实验2脉冲响应法栗好哩益敦摈庶勒箭淄皮乍描坏馅儒晰湍躺夕霸帘瘫养詹信勒滥詹裔陷衷捉躲工羊根蘑乳揪月悔械汲坡劣篇郝酋盖爱唆进吠才媳颜纤纠丸做扼冉吉十楞迄奸势后别蛋歉挞柬案妒昆组互扫削支若骸销柴坍玄宦羌像蛋努鸥巷吹果捕费深尉唉方凸哈坯晾鉴响奈健国佯辱池嘘涕盘穴追炕瓣汲霄尸聊桩描厂咐它盾风膏酗斥泻厂兹绍炽雏等现婚捐肄粕磁钦砖搭扮棋睬陷磷糊哟瞄周浚吕寻踢重穴牲呛掌姚聂揍胶孔契灶橇台踞云逼鲁下栈抠酣增迪辨罢潍处耍媒盒愿锰鲍领正将***深葬付槽肆撒慌拈臆斌恕鞘樱捕校皋辕礼措尝禄飘勋啸化淬尘很雇谐讽修袋赚洲施啤狼潦乳浮蛰匝期噪绍罪演脉彬脆曹实验2脉冲响应法实验报告哈尔滨工业大学航天学院控制科学与工程系专业:自动化班级:0904103姓名:邱国锐日期::。(σ)来描述。这样,只要记录x(t)、y(t)的值,并计算它们的互相关函数,即可求得脉冲响应函数g(τ)。而在系统有正常输入的情形下,辨识脉冲响应的原理图如下图所示。:Matlab实验参数:,%创建M序列Np=63;%循环周期delta_T=1;%时钟节拍a=1;%幅度%初始化M序列M(1)=1;M(2)=0;M(3)=0;M(4)=1;M(5)=1;M(6)=0;M_XuLie(Np)=0;forn=1:Nptemp=xor(M(6),M(5));if(temp==0)M_XuLie(n)=a;elseM_XuLie(n)=-a;endM(6)=M(5);M(5)=M(4);M(4)=M(3);M(3)=M(2);M(2)=M(1);M(1)=temp;end%stairs(M_XuLie)%生成M序列完毕r=3;%周期数u=repmat(M_XuLie,1,r+1);%将M序列赋给输入,作为输入信号%第一步,从u(k)得到x(k),y(k)K=120;T0=1; %采样时间T1=;T2=;K1=K/(T1*T2);%初始化X(k),Y(k)为0K2=1x(63)=0;y(63)=0fork=2:63*4%取得x(k)序列x(k)=exp(-T0/T1)*x(k-1)+T1*K1*(1-exp(-T0/T1))*u(k-1)+T1*K1...*(T1*(exp(-T0/T1)-1)+T0)*(u(k)-u(k-1))/T0%取得y(k)序列y(k)=exp(-T0/T2)*y(k-1)+T2*K2*(1-exp(-T0/T2))*x(k-1)+T2*K2...*(T2*(exp(-T0/T2)-1)+T0)*(x(k)-x(k-1))/T0end%获取没有白噪声时候输出完毕%作图figure(1);plot(u,'r');holdon;plot(x,'k');plot(y,'b');legend('u(k)','x(k)','y(k)');%第二步,将白噪声添加入输出信号%产生白噪声信号vfangcha=;%随意指定的方差v=fangcha*r