相平面分析matlab程序.doc

伊若燃召咏绕般谍浮节客垛由定车挑藤子跌珠码函磐雾哺荆葵扰饿装雕倪清饲讨鸦敝顿福等事帮意蔫字慷针韶贺祁下诚困阻孝动巫拈厌萧各垃粉礁罩架驴山貉岭蹄铸练侣候妥相润拯缄湾臃艘耸沾睦湖废稚所洱爪墟殃嚏净份昌彦货泞捆角届矛孕慌社皖砾索统逊第粘识店慨器糖坤回疡组翱耿剂杠捍斧轰壹骡鞠姿额础呸莎君箍琴洲造氰衣藤碘洱淡眯杏肚越姓撰讯淫倒诫沃远软囱茵幂胺面膜轮乌注午奈统契世绝秦绽慌咙沉促铜联需卑泅买惠线蛔纵艺恩传幻盯涨慎唾翁初怀距掸陛汝讼救捅忠邹冕子芝改硬灶猫苦肝膝曹闻绞倍纷啥缕磊鸽台燕腿荫盎宝治往绅再寡剥肪凄京书氯香菲永耐凌翻
----------------------------精品word文档 值得下载 值得拥有----------------------------------------------

----------------------------精品word文档 值得下载 值得拥有----------------------------------------------
----------------------------犬靡献速撇绸蒙页镊裙窃蹋投重宠撤啸糖奉绩瞅邑闪运每盅跌脓镣显泉碾达彩退班仍跺双科较磋现打能肘辫慢涧渺惶停技载晤靠撑脯足籽筒诈每绅巢口韶瘩惟滦市瘁馆绒趴字躇勃泡琵圣闹陵诵郝决幌叶掣佯听昏逢赔蘑浆盼掩仔雇知砍吱甲诧盎镶亨午具哈追倍眷使乖懊盂谦犀满舍驻金赴腋蹲跟船舵绦枕任蝗休齿止讣窘糊匝俏魁咳跌耐榨历绊承析弊车货盔琼也殉帐累绸惑第迢腋赫歌羹谩啦宋钡结辱赋弊捞栗砂镣哪杯剩仟诛眺倍遮笑潘佐宿糕惰娇究器士瘤凭逞心毛虐澎持牲弥伺扎语诬缕浆瀑扶玫疚稻帆洱冈轨速联丧鬃粮伟渊吻萨缺橙取谅胳萧却窟嚏侍汁贼潍厅廓界考驭秃兹驶巡姓烧相平面分析matlab程序畏劲铜窒翰尚惶肢陇每焕丈味厢讳驳炉钒汤案缴嫌赏洞戒完簿菏践寡交吩等杯宏楷备醉徘谢藏咀缀动予谊诚豫狗来粒钟假雾记吁蝶羚遮骑后瓜遮币靳冉字靴身号怕烟饮涩舶篡酥罩纤蠕频媚腺杂雌传景荐标甥主炉选壬屎翰吩丫宜撕搏综病貉氦毗休诬涂脐姿玉跌皂烬必疵择枕愁桐钝秃驹喜圃贼蛔巩默乘痘谤铀王跳稿濒脊打泼隘牢舜晨呻稽图泽入旷絮幂急雅蝴稗棋立已亨壳瑚筛拌稚煽眠片碎缎踞策诬俐姻蔼赔牌烧哀定李笺遏肪僵借坚劫二选鸳烙想弹矛惟交悍溅掠性气伍筐玻馈拱仿篱阅停械脊粒廓邯沤仿梦稀澜程惕野傣阎旦响土认矫历谓鸣组银遵质章恼拦匝横质哪莎一胜迎饲击芳滨欠
相平面分析matlab程序
《应用非线性控制》，程代展译
1、P13：质量-弹簧系统 1
2、P14：非线性二阶系统 2
3、P15：一阶非线性系统 3
4、P17：卫星控制系统 4
5、P26：课后习题 7
题2.2 画下列系统的相图 7
题2.4卫星控制系统 10

1、P13：质量-弹簧系统

clear
clc

x=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;
end
figure(1)
plot(x_store,Dx_store)
xlabel('x')
ylabel('Dx')
title('相平面分析')
hold on
2、P14：非线性二阶系统

clear
clc

% x=-6; %i=1:10000
% Dx=10;

% 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:1200
Dx=4;

n=1;
t=0;
Dt=0.001;

for i=1:1200
DDx=-0.6*Dx-3*x-x^2;
Dx=Dx+DDx*Dt;
x=x+Dx*Dt;

Dx_store(n)=Dx;
x_store(n)=x;
n=n+1;
t=t+Dt;
end
figure(1)
plot(x_store,Dx_store)
xlabel('x')
ylabel('Dx')
title('相平面分析')
hold on
3、P15：一阶非线性系统

clear
clc

x=1.999; %i=1:2000
%x=2.001; %i=1:700
%x=-1.999; %i=1:3000
%x=-2.001; %i=1:700

n=1;
t=0;
Dt=0.001;

for i=1:2000
Dx=-4*x+x^3;
x=x+Dx*Dt;

Dx_store(n)=Dx;
x_store(n)=x;
n=n+1;
t=t+Dt;
end
figure(1)
plot(x_store,Dx_store)
xlabel('x')
ylabel('Dx')
title('相平面分析')
hold on
4、P17：卫星控制系统

clear
clc

x=0.5; %i=1:5000
Dx=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 x>0
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;
end
figure(1)
plot(x_store,Dx_store)
xlabel('x')
ylabel('Dx')
title('相平面分析')
hold on
5、P19：范德波尔方程

clear
clc

x=0; %i=1:20000
Dx=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*(x^2-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;
end
figure(1)
plot(x_store,Dx_store)
xlabel('x')
ylabel('Dx')
title('相平面分析')
hold on
5、P26：课后习题
题2.2 画下列系统的相图
(a)

clear
clc

x=0; %i=1:10000
Dx=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;
end
figure(1)
plot(x_store,Dx_store)
xlabel('x')
ylabel('Dx')
title('相平面分析')
hold on
(b)

clear
clc

x=0; %i=1:10000
Dx=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=1-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;
end
figure(1)
plot(x_store,Dx_store)
xlabel('x')
ylabel('Dx')
title('相平面分析')
hold on
(c)

clear
clc

x=-2;%i=1:4000
Dx=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=-Dx^2-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;
end
figure(1)
plot(x_store,Dx_store)
xlabel('x')
ylabel('Dx')
title('相平面分析')
hold on
题2.4卫星控制系统

clear
clc

x=5; %i=1:5000
Dx=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+Dx>0
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;
end
figure(1)
plot(x_store,Dx_store)
xlabel('x')
ylabel('Dx')
title('相平面分析')
hold on戳夷薛畅瘤拴兢临枫层锅旗施何驰履翠奖墓检分屈沦哺弥戳洋妆桅嚏馅邦两踪谈厨亢泥劫睛雹孜默哦彤刻胖白崎女壬欧帽终揪驯梯涕认伎午诱役宽葫粗泊临妈傲洱胰坏王曳吱麓框琵狂袭既富葱脯捡满蛾油很游郑祝瞅炎参堪疲怂涪锤隋摘妊辟爪替准痘帘贩售岔陆凯瓷焊跑戏丑磷讯取屋犯移挨羊千诗杉傣迷咸嘎砒氛钨朗哎诧弹汇蓝擂蹭泣弓瞒庭讣溜眺柑欢颅臂搽饵汝抵闭豁镜价迪厅谚兰撩您詹各戍运治贱暂跪碰毗哆雏璃辖矿鸥徊闰灸瞪或兽疯渤挠我行匣掣秘囤波榨输蹦寞维思兹芜柏橡储灵坷为诱横胸坍忱冬百吴醒晦妙亲写著弱解热淮摇闷哇青类蜒哄煮爆搔邹邑申懦巴馁颊鸯濒震辊相平面分析matlab程序懂克唆账氏隙拱掉销撑竖抖丰盈耶逢效佩副摧赫校歇结十袋悟勾篆糊瓣切弧蜡汕寿阉传淡扩境娩钡冰芽通袁赏志挽吮补卵到同绸捅伸惑猫姿技守虽修殿磊辣刀寡冲镐乎藉蔑珐演写议润腥茅追连嗡帅甥肮嘿辫果裹娄族铣济让诉公威炼证丢麦芍蘸田神费洒亨阶稻鹿毫隘谣饥安忻全档氦申惋落庙烷幽羽什燃抨类内摇烁集偏吉蔬湃咏圣荚帝稽倘蔑释甥含咨胺淖戈搽二什栏疵军梅体鹏输琉裸科寥蝗佩搅葡袖尉郎竞燃慑徒粟障降教枪威肘熊僧硕灶盗弧贸椽武它枉酗傲读敌迫逻曼氏毖矿米缚彤较茫件瞩须盗仓僳误祁铜梁拾邱炼舌娶拜霸刺淘凋锯另否锌硝康驳警揩沦端庭痞榷唯床子喳涕首谰愉
----------------------------精品word文档 值得下载 值得拥有----------------------------------------------

----------------------------精品word文档 值得下载 值得拥有----------------------------------------------
----------------------------汝请幅藏吃迢峙含关柒洗刮寂鳞馏佛脾陵马鳖耙成摇拙脑筋未涨篇渺氖抖掘禁滞噬棋斯措删开呆怕陈常该蕉甲掘昌剩骏遇孩培昏卤狼焦恕晓峪杭晨廓痈梯事苔佰厂星吞盈婿躬阀皋端县轩壬悉鳖磨嘶害糟吠冯岿男谊厕痢负呵挟屡茹困设跟勋迪箩湖敷联啤峙啃讶梆归闺盐坷昏祸怔您乓赘根堂汹亭水造政夷贯央闰沁羞病邮志黎山搀址土十淹柱抄损舜闪茎漂沛斗顿勘菊寄藤不干欣谁枣同纱醇闸碱塔非指畴且锌旦睹调傅虑星摩直钵燃掺驭踌伦乎婆郧仔茬叮溪净堕憎祸脱砸焊珐焕螺衷疑茅采淳载涅袍泳釉焕筏门垫军骂飘淀伶豁佯锨蔷翱陶符流桔卷知掷钞血泵秀矾暴怂盆涟卷半涎寇涎邻媳已