clc;clear;close?all;%清除變量、窗體、及工作區間
global?rin?yout?timef

%*********第一步：單純形替換法變量準備及設定***********
x0?=?[1?0.8?0.8];%Kp，Ti，Td初始值
l?=?1e-6;%單純形棱長
r?=?1;%反射系數Gama，通常取1
e?=?2;%延伸系數，通常取2
n?=?3;%n?=?3表示問題為三維空間最優點求解
c?=?0.5;%收縮系數，通常取0.5
Maxstep?=?1000;%迭代最大次數
MarginErr?=?5e-13;%誤差限，一部西路
Bestv?=?zeros（1?3）;%最優解
Bestf?=?0;%最優解對應的函數值
[v?f]?=??Initialize（x0?n?l）;%調用初始化函數

%**第二步：單純形反射，延伸，收縮，減小棱長得到最優點****
Deltarecord?=?[];%誤差記錄矩陣
frecord?=?[];%函數值記錄矩陣
for?i?=?1?:?Maxstep
????%調用FYSJ函數求的下一次迭代所需要的單純形
????[Nextv?Nextf?Delta?Meanf]?=?FYSJ（v?f?r?e?c?n）;
????Deltarecord?=?[Deltarecord?Delta];%記錄誤差
????if?（Delta?????????for?i?=?1?:?3
???????????Bestv（i）?=?sum（Nextv（:?i））?/?（n?+?1）;
????????end
????????Bestf?=?Targetf（Bestv）;
????????frecord?=?[frecord?Bestf];%記錄函數值
????????break;
????else
????????v?=?Nextv;
????????f?=?Nextf;
????????frecord?=?[frecord?Meanf];%記錄函數值
????end
end
%********第三步：做出誤差收斂曲線，函數值變化曲線*******
figure;
%誤差收斂曲線
[msize?nsize]?=?size（Deltarecord）;
t?=?1?:?nsize;
plot（t?Deltarecord?‘b‘）;
xlabel（‘時間‘）;ylabel（‘誤差‘）;
title（‘誤差收斂曲線‘）;
figure;
%函數值變化曲線
[msize?nsize]?=?size（frecord）;
t?=?1?:?nsize;
plot（t?frecord?‘b‘）;
xlabel（‘時間‘）;ylabel（‘函數值‘）;
title（‘二次型性能指標收斂曲線‘）;
%系統響應圖
figure;
hold?on;
plot（timef?yout）;
xlabel（‘時間‘）;ylabel（‘輸出‘）;
title（‘整定后系統階躍響應圖‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1416??2012-11-19?17:14??piddanchn\danchunmain.asv

?????文件???????1614??2012-12-12?09:56??piddanchn\danchunmain.m

?????文件????????586??2012-11-01?11:51??piddanchn\Error.m

?????文件???????1905??2012-11-01?11:04??piddanchn\FYSJ.m

?????文件????????726??2012-11-01?10:41??piddanchn\Initialize.m

?????文件????????908??2012-11-19?17:29??piddanchn\Targetf.asv

?????文件????????910??2012-11-22?14:59??piddanchn\Targetf.m

?????目錄??????????0??2012-11-19?22:10??piddanchn

-----------?---------??----------?-----??----

?????????????????8065????????????????????8