exm080652_1.m

同上

%exm080652_1.m 		用于ITAE标准型系数的寻优
clc
disp('按任意键，开始运行。')
disp('(1) 计算 6 阶 ITAE 传递函数分母多项式系数。')
pause
disp(' ')
disp('    正在计算中，请稍等 ！ ')
clear
global a ss
a0=[3.3 6.6 8.6 7.5 3.9];	
Tspan=(0:500)/10;			
options=optimset('LargeScale','off');
a=fminunc(@itae,a0,options);		
coeff=[1,a,1]					
ss							


disp(' ')
disp('按任意键，继续运行。')
disp('(2) 新老 6 阶标准型的单位阶跃响应比较。')
pause
old=tf(1,[1 3.25  6.60  8.60  7.45  3.95  1]);
new=tf(1,coeff);
[yold,told]=step(old,50);
[ynew,tnew]=step(new,50);
plot(told,yold,'b','LineWidth',1)
axis([3,18,0.95,1.05])
hold on,plot(tnew,ynew,'r','LineWidth',3),hold off
legend('Old','New',4),grid on
shg