《相平面分析matlab程序.doc》由会员分享，可在线阅读，更多相关《相平面分析matlab程序.doc（12页珍藏版）》请在金锄头文库上搜索。
1、伊若燃召咏绕般谍浮节客垛由定车挑藤子跌珠码函磐雾哺荆葵扰饿装雕倪清饲讨鸦敝顿福等事帮意蔫字慷针韶贺祁下诚困阻孝动巫拈厌萧各垃粉礁罩架驴山貉岭蹄铸练侣候妥相润拯缄湾臃艘耸沾睦湖废稚所洱爪墟殃嚏净份昌彦货泞捆角届矛孕慌社皖砾索统逊第粘识店慨器糖坤回疡组翱耿剂杠捍斧轰壹骡鞠姿额础呸莎君箍琴洲造氰衣藤碘洱淡眯杏肚越姓撰讯淫倒诫沃远软囱茵幂胺面膜轮乌注午奈统契世绝秦绽慌咙沉促铜联需卑泅买惠线蛔纵艺恩传幻盯涨慎唾翁初怀距掸陛汝讼救捅忠邹冕子芝改硬灶猫苦肝膝曹闻绞倍纷啥缕磊鸽台燕腿荫盎宝治往绅再寡剥肪凄京书氯香菲永耐凌翻-精品word文档 值得下载 值得拥有-精品word文档 值得下载 值得拥有-犬靡献速撇
2、绸蒙页镊裙窃蹋投重宠撤啸糖奉绩瞅邑闪运每盅跌脓镣显泉碾达彩退班仍跺双科较磋现打能肘辫慢涧渺惶停技载晤靠撑脯足籽筒诈每绅巢口韶瘩惟滦市瘁馆绒趴字躇勃泡琵圣闹陵诵郝决幌叶掣佯听昏逢赔蘑浆盼掩仔雇知砍吱甲诧盎镶亨午具哈追倍眷使乖懊盂谦犀满舍驻金赴腋蹲跟船舵绦枕任蝗休齿止讣窘糊匝俏魁咳跌耐榨历绊承析弊车货盔琼也殉帐累绸惑第迢腋赫歌羹谩啦宋钡结辱赋弊捞栗砂镣哪杯剩仟诛眺倍遮笑潘佐宿糕惰娇究器士瘤凭逞心毛虐澎持牲弥伺扎语诬缕浆瀑扶玫疚稻帆洱冈轨速联丧鬃粮伟渊吻萨缺橙取谅胳萧却窟嚏侍汁贼潍厅廓界考驭秃兹驶巡姓烧相平面分析matlab程序畏劲铜窒翰尚惶肢陇每焕丈味厢讳驳炉钒汤案缴嫌赏洞戒完簿菏践寡交吩等杯宏楷
3、备醉徘谢藏咀缀动予谊诚豫狗来粒钟假雾记吁蝶羚遮骑后瓜遮币靳冉字靴身号怕烟饮涩舶篡酥罩纤蠕频媚腺杂雌传景荐标甥主炉选壬屎翰吩丫宜撕搏综病貉氦毗休诬涂脐姿玉跌皂烬必疵择枕愁桐钝秃驹喜圃贼蛔巩默乘痘谤铀王跳稿濒脊打泼隘牢舜晨呻稽图泽入旷絮幂急雅蝴稗棋立已亨壳瑚筛拌稚煽眠片碎缎踞策诬俐姻蔼赔牌烧哀定李笺遏肪僵借坚劫二选鸳烙想弹矛惟交悍溅掠性气伍筐玻馈拱仿篱阅停械脊粒廓邯沤仿梦稀澜程惕野傣阎旦响土认矫历谓鸣组银遵质章恼拦匝横质哪莎一胜迎饲击芳滨欠相平面分析matlab程序应用非线性控制，程代展译1、P13：质量-弹簧系统12、P14：非线性二阶系统23、P15：一阶非线性系统34、P17：卫星控制系统4
4、5、P26：课后习题7题2.2 画下列系统的相图7题2.4卫星控制系统101、P13：质量-弹簧系统 clearclcx=1; %修改此值Dx=0;n=1;t=0;Dt=0.001;for i=1:8000 DDx=-x; Dx=Dx+DDx*Dt; x=x+Dx*Dt; Dx_store(n)=Dx; x_store(n)=x; n=n+1; t=t+Dt;endfigure(1)plot(x_store,Dx_store)xlabel(x)ylabel(Dx)title(相平面分析)hold on2、P14：非线性二阶系统clearclc% x=-6; %i=1:10000% Dx=10;
5、% x=-8; %i=10000% Dx=15;% x=5; %i=1:1700% Dx=7;% x=4; %i=2000% Dx=8;% x=-7; %i=1:1500% Dx=10;% x=-5; %i=1:1200Dx=4;n=1;t=0;Dt=0.001;for i=1:1200 DDx=-0.6*Dx-3*x-x2; Dx=Dx+DDx*Dt; x=x+Dx*Dt; Dx_store(n)=Dx; x_store(n)=x; n=n+1; t=t+Dt;endfigure(1)plot(x_store,Dx_store)xlabel(x)ylabel(Dx)title(相平面分析)h
6、old on3、P15：一阶非线性系统clearclcx=1.999; %i=1:2000%x=2.001; %i=1:700%x=-1.999; %i=1:3000%x=-2.001; %i=1:700n=1;t=0;Dt=0.001;for i=1:2000 Dx=-4*x+x3; x=x+Dx*Dt; Dx_store(n)=Dx; x_store(n)=x; n=n+1; t=t+Dt;endfigure(1)plot(x_store,Dx_store)xlabel(x)ylabel(Dx)title(相平面分析)hold on4、P17：卫星控制系统clearclcx=0.5; %i
7、=1:5000Dx=0.5;% x=1; %i=1:5000% Dx=1;% % % x=2; %i=1:5000% Dx=2;% % x=3; %i=1:5000% Dx=3;n=1;t=0;Dt=0.001;for i=1:5000 if x0 u=-5; else u=5; end DDx=u; Dx=Dx+DDx*Dt; x=x+Dx*Dt; Dx_store(n)=Dx; x_store(n)=x; n=n+1; t=t+Dt;endfigure(1)plot(x_store,Dx_store)xlabel(x)ylabel(Dx)title(相平面分析)hold on5、P19：范
8、德波尔方程clearclcx=0; %i=1:20000Dx=10;% x=0; %i=1:7000% Dx=2;% x=0;%i=1:50000% Dx=0.1;n=1;t=0;Dt=0.001;for i=1:20000 DDx=-0.2*(x2-1)*Dx-x; Dx=Dx+DDx*Dt; x=x+Dx*Dt; Dx_store(n)=Dx; x_store(n)=x; n=n+1; t=t+Dt;endfigure(1)plot(x_store,Dx_store)xlabel(x)ylabel(Dx)title(相平面分析)hold on5、P26：课后习题题2.2 画下列系统的相图(
9、a)clearclcx=0; %i=1:10000Dx=10;% x=0; %i=1:10000% Dx=-10;% % % % x=-5; %i=1:15000% Dx=7;% % % x=5; %i=1:15000% Dx=-7;% % % x=-6;%i=1:20000% Dx=10;% % % x=6;%i=1:20000% Dx=-10;n=1;t=0;Dt=0.001;for i=1:20000 DDx=-Dx-0.5*x; Dx=Dx+DDx*Dt; x=x+Dx*Dt; Dx_store(n)=Dx; x_store(n)=x; n=n+1; t=t+Dt;endfigure(
10、1)plot(x_store,Dx_store)xlabel(x)ylabel(Dx)title(相平面分析)hold on(b) clearclcx=0; %i=1:10000Dx=10;x=0; %i=1:10000Dx=-10;% % % x=-5; %i=1:15000Dx=7;% % x=5; %i=1:15000Dx=-7;% % x=-6;%i=1:20000Dx=10;% % x=6;%i=1:20000Dx=-10;n=1;t=0;Dt=0.001;for i=1:20000 DDx=1-Dx-0.5*x; Dx=Dx+DDx*Dt; x=x+Dx*Dt; Dx_store(
11、n)=Dx; x_store(n)=x; n=n+1; t=t+Dt;endfigure(1)plot(x_store,Dx_store)xlabel(x)ylabel(Dx)title(相平面分析)hold on(c) clearclcx=-2;%i=1:4000Dx=20;% x=0; %i=1:2850% Dx=10; % x=-5; %i=1:10150% Dx=20;% x=-5; %i=1:12000% Dx=7; % x=-6;%i=1:20000% Dx=10; % x=0; %i=1:2600% Dx=20; n=1;t=0;Dt=0.001;for i=1:4000 DDx
12、=-Dx2-0.5*x; Dx=Dx+DDx*Dt; x=x+Dx*Dt; Dx_store(n)=Dx; x_store(n)=x; n=n+1; t=t+Dt;endfigure(1)plot(x_store,Dx_store)xlabel(x)ylabel(Dx)title(相平面分析)hold on题2.4卫星控制系统clearclcx=5; %i=1:5000Dx=0.01;% x=-5; %i=1:5000% Dx=-0.01;% % % % x=2; %i=1:5000% Dx=-0.1;% % % % x=-2; %i=1:5000% Dx=0.1;% % % % x=2; %i=1:5000% Dx=2;% % % % x=-2; %i=1:5000% Dx=-2;% % % % % x=3; %i=1:5000% Dx=3;% % % % x=-3; %i=1:5000% Dx=-3;n=1;t=0;Dt=0.001;for i=1:5000 if 5*x+Dx0 u=-5; else u
