用MATLAB实现的非线性频率响应分析.docx

广布石油1 ：高等专杓学校学抱
JOURNAL CW' 皿LEGE
i JF P］：T^l X'HFCMICAl 「"IN〕】」
用MATLAB实现的非线性频率响应分析
司注3
:'：1.■:.J：i_ . &可缺少的旧力11, MATLAB时以在—触的机 器类Y I使周.，：除「，逐琨供「更降迄算、13妣始制■.数据姓原:、 Wmcl-wvs^S等便孩〕具*它提IU -个开放式的墟成环境羽山町世运行系筮提供的大景血令戒用起 塞极其方便二［
2非线匹系统的频率恂应分析及MATLAB实现
严格世说，艳仄部弁控制系统都是样戮性系州旧我率分析法是工程中虚用…泛的-特分析方凯所
审3期 讶徒莹，的MATLAR实现的非线宅巍率吟电分析 37
以IE绒性羸统的频率响应务折是十分重罢的。#蚊性系统鼠率-响直分航的结果不但取决「辎人法肆的晦 率，」.；• 求取，这神方法这但能院指茗七绑折车线性环节系统的频率响应，还可以求出带有分授税性壮II浅性环节 菸琉的厥率响应。
2 1仿真算由，与思踌分析
设某罪螭性环节的摘入函数为『(〔)= MxWi),其中M fit w可以在MATLAB坏境中用变量祖和 ww来表示D ft St过程中雨璀用定£长4/3阶Range-Kurta方法，其中输入住寸rC t j为土批信号,r(x)~ .=l2…=1,，••「*)：.如果股定壬弦周明T=h~向仿真n ML则怙真步长 幻=】Kn-l);若取丁 =却,则总定的非攻性系统在有限次的仿真之后「 岭F两个周期曲线布覆于-寮「若轴人的僖峙为L.'fM」，便可叩*下分」71集幅值
A,., -〔m敬(跟」-nu-；iyp i|) )/2A,
可见击鲁个周期氏至少有两个过零点,如果有备于两个过零点』J存在求堰响应相世伉的何瓯耳果 只有两&迂零雅，就要职出处T_r升段批豹相技值。暇设这洋的过零焉，茨％在第k和郃*卜I仿耳立迎， 则可以通隽曜性捕『［得出相技值
S| = ：km ,6-._k，：)苒。.•T】1 - L)
拉一舞试虽然前季：四即使对「一个简单的系统栾说，殷率宿底也需要几打个同期的仿世分F「，所以 当血很多撮率值和■值时,十|*量逐是相当大的「为r提高计算牧率，可以在MATLAB中漆用以下措
1｝巡免彳；，子要的计算周期 定义一个收散精则占捕足血敛街帅I,则级止罚芯(代.，略对的访真， 存站苴-得&的蠕值冗电曲值，卅始对上一-个g,w徭吁伯真“收匙准则叫定义为
.1 = mfljcm-y . - ：. J -*丽以|>- f =e
其山：q= I疽，…， I为第］】周期内的所有仿真点一 ，
2)使用在(A”叫一］)时时收散处的状态值作仿真初直在实际仿真中，系统的收敏通度取映看初始 状槛的遂取，所以叫一 |对的收散状茶可能对垠取下-个领率下的响廉有用，因为当频牢11勺变化值不是很 X3>,(Alrw ,｝ T'g过程的初帷

根拇理上思路，