clear;close?all;
gamma=0.75;%設定慣性因子

eps1=0.005;%設定停止訓練的條件參數
m1=8;%設定隸屬函數個數
m2=8;

a=-1;b=1;
w0=a+（b-a）*rand（1m1*m2）;%初始化權值陣

for?i=1:2
????switch?i
????????case?1beta=0.75;%設定學習率
????????otherwisebeta=0.25;
????end
????c=[2/7*（0:m1-1）-1;2/7*（0:m2-1）-1];%初始化Cij
????sigma=0.1213*ones（2m1）;%初始化σij　
????w=w0;
????mu=zeros（2m1）;%初始化隸屬度層
????alpha=zeros（1m1*m2）;%初始化規則層
????alpha_=zeros（1m1*m2）;%初始化歸一化層

????delta2=zeros（2m1）;%初始化δ2
????dw=zeros（1m1*m2）;%初始化權值變化量
????dc=zeros（2m1）;%初始化c變化量
????dsigma=zeros（2m1）;%初始化σij變化量

%----------------訓練過程----------------%
????err=1;%初始化誤差
????er=[];
????counter=0;%統計訓練步數
????while（err>=eps1）
????????Par_E_w=zeros（1m1*m2）;%誤差對權值陣的偏導數
????????Par_E_c=zeros（2m1）;%誤差對c的偏導數
????????Par_E_sigma=zeros（2m1）;%誤差對σ的偏導數
????????E=0;
????????for?x1=-1:2/19:1
????????????for?x2=-1:2/19:1
????????????????yd=sin（pi*x1）*cos（pi*x2）;%期望輸出，導師信號
????????????????%正向傳播
????????????????mu（1:）=exp（-（x1-c（1:））.^2./sigma（1:）.^2）;%計算隸屬度層
????????????????mu（2:）=exp（-（x2-c（2:））.^2./sigma（2:）.^2）;
????????????????s=zeros（2m1m1*m2）;%初始化Sij
????????????????for?m=1:m1
????????????????????for?n=1:m2
????????????????????????alpha（（m-1）*m2+n）=min（mu（1m）mu（2n））;%計算規則層（取小運算）
????????????????????????if?mu（1m）<=mu（2n）
????????????????????????????s（1m（m-1）*m2+n）=1;
????????????????????????end
????????????????????????if?mu（1m）>=mu（2n）
????????????????????????????s（2n（m-1）*m2+n）=1;
????????????????????????end
????????????????????end
????????????????end
????????????????alpha_=alpha/sum（alpha）;%計算歸一化層
????????????
????????????????y=alpha_*w.‘;%計算網絡輸出
????????????????E=E+1/2*（yd-y）^2;%計算誤差
????????????????%反向傳播
????????????????delta5=yd-y;
????????????????delta4=delta5*w;

????????????????for?k=1:m1*m2
????????????????????delta3（k）=delta4（k）*（sum（alpha）-alpha（k））./sum（alpha）^2;
????????????????end
????????????????for?m=1:2
????????????????????for?n=1:m1
????????????????????????delta2（mn）=0;
????????????????????????for?l=1:m1*m2
????????????????????????????delta2（mn）=delta2（mn）+delta3（l）*s（mnl）*mu（mn）;
????????????????????????end
????????????????????end
????????????????end
????????????
????????????????Par_E_w=Par_E_w-delta5*alpha_;%計算偏導數
????????????????Par_E_c（1:）=Par_E_c（1:）-2*delta2（1:）.*（x1-c（1:））./sigma（1:）.^2;
????????????????Par_E_c（2:）=Par_E_c（2:）-2*delta2（2:）.*（x2-c（2:））./sigma（2:）.^2;
????????????????Par_E_sigma（1:）=Par_E_sigma（1:）-2*delta2（1:）.*（x1-c（1:））.^2./sigma（1:）.^3;
????????????????Par_E_sigma（2:）=Par_E_sigma（2:）-2*delta2（2:）.*（x2-c（2:））.^2./sigma（2:）.^3;
????????????end
????????end
??
????????num=20*20;
????????Par_E_w=Par_E_w/num;
????????Par_E_c=Par_E_c/num;
????????Par_E_si