稳定性评价

1. %*****************控制系统ITAE性能计算*****************%  
2. function [Ts,CT,Ess,ITAE]=analysis(y,r,t)  
3. mp=max(y);  
4. tp=spline(y,t,mp);  %峰值时间  
5. cs=length(t);  
6. yss=y(cs);          %稳态值  
7. CT=(mp-yss)/yss;    %超调量  
8. i=length(t);  
9. while(y(i)>0.98*yss)&(y(i)<1.02*yss)  
10.     i=i-1;  
11. end  
12. Ts=t(i);            %调节时间  
13. %稳态误差  
14. Ess=abs(r(length(t))-yss);  
15.   
16. e=0;  
17. dt=0.01;  
18. ITAE=0;  
19. for j=1:length(t)  
20.     tdt=t(j)*dt;   
21.     e=y(j)-r(j);  
22.     ITAE=ITAE+abs(e)*tdt; %ITAE指标计算   
23. end  
24. ITAE=roundn(ITAE,-4);  

1. %************构造评价矩阵**************%  
2. function [R1,R2]=evaluate()  
3. % 参数初始化  
4. inputfile = 'C:\Users\beartree\Desktop\GUI\ITAEsample.xls';   
5. %% 读取数据并进行fcm计算  
6. [num,txt]=xlsread(inputfile);  
7.   
8. for i=1:11  
9.     data(:,1)=num(:,(2*i-1));  
10.     data(:,2)=num(:,(2*i));  
11.     options=[2;100;1e-5;1];  
12.     [center,U,obj_fcn]=fcm(data,3,options);  
13.     %%结果处理  
14.     maxC=min(center(:,2));  
15.     minC=max(center(:,2));  
16.     stab=find(center(:,2)==maxC);  
17.     lose=find(center(:,2)==minC);  
18.     for k=1:3  
19.         if k~=stab && k~=lose  
20.             wave=k;  
21.         end  
22.     end  
23.     R(i,1)=U(stab,1);  
24.     R(i,2)=U(wave,1);  
25.     R(i,3)=U(lose,1);  
26. end  
27. %%输出结果  
28. for j=1:11  
29.     fprintf('评价指标参数C%d对应稳定、波动、失控状态隶属度分别为：',j);  
30.     fprintf('%1.4f,%1.4f,%1.4f \n',R(j,1),R(j,2),R(j,3));  
31. end  
32.  fprintf('评价矩阵为:\n');  
33.  R1=R(1:7,:);  
34.  R2=R(8:11,:);  
35. clc;  

1. %*************模糊综合合成*****************%  
2. function [S1,S2,S]=AHPFCE()  
3. W1=[0.4614 0.0973 0.0801 0.0818 0.1297 0.0901 0.0595];  
4. W2=[0.1915 0.5495 0.0081 0.1408];  
5. W=[0.55 0.45];  
6. [R1,R2]=evaluate();  
7. S1=W1*R1;  
8. S2=W2*R2;  
9. R=[S1;S2];  
10. S=W*R;