牛拉潮流程序(編輯修改稿)
2025-01-09 22:00 本頁面
【文章內(nèi)容簡介】 endendfor i=1:nm1 Delta_Q(i)=Delta_Q(i)(Q(i)QG(i))。endDelta_QDelta_V=zeros(1,m)。for i=1:m for j=1:n if Type(j)==2 Delta_V(i)=U0(j)^2(e(j)^2+f(i)^2)。 end endendDelta_Vnum=0。disp([39。第39。,num2str(num),39。次時的Delta總的失配量為：39。]) %進(jìn)入循環(huán)體判斷是否滿足條件%%先算出最大值，作為判斷是否收斂的依據(jù)%DEL=[Delta_P Delta_Q]。 %Delta_P Delta_Q______________%%%%%%MAX =max(abs(DEL))。 MAXTheta_first=zeros(1,n)。U_f=U_first39。Delta_F_E1=[Theta_first(1:n1) U_f(1:nm1)]。Delta_F=Delta_F_E139。Delta_Cor=Delta_F_E1。 %_______________Delta_the Delta_u________________%%disp([39。第一次最大失配量誤差：39。,num2str(MAX)])%循環(huán)判斷% if MAX1e004 % 判斷依據(jù) disp(39。下面開始下一次迭代過程！39。)endwhile MAX1e004 num=num+1。%%%%%%%%%%%%%%%%%%%形成雅克比矩陣%%%%%%%%%%%%%%%%%%%%%%%%%%%%先求非對角元素(H)%Hik=zeros(n1,n1)。for i=1:n1 for k=1:n1 if i~=k theik=Theta_first(i)Theta_first(k)。 Hik(i,k)=U_first(i)*U_first(k)*(G(i,k)*sin(theik)B(i,k)*cos(theik))。 end endend%再求對角元素（H）%for i=1:n1 for k=1:n if i~=k theik=Theta_first(i)Theta_first(k)。 Hik(i,i)=Hik(i,i)+U_first(k)*(G(i,k)*sin(theik)B(i,k)*cos(theik))。 end end Hik(i,i)=U_first(i)*Hik(i,i)。endHik%先求非對角元素N%Nik=zeros(n1,nm1)。for i=1:n1 for k=1:nm1 if i~=k theik=Theta_first(i)Theta_first(k)。 Nik(i,k)=U_first(i)*U_first(k)*(G(i,k)*cos(theik)+B(i,k)*sin(theik))。 end endend%再求對角元素%for i=1:nm1 for k=1:n if i~=k theik=Theta_first(i)Theta_first(k)。 Nik(i,i)=Nik(i,i)+U_first(k)*(G(i,k)*cos(theik)+B(i,k)*sin(theik))。 end end Nik(i,i)=U_first(i)*Nik(i,i)2*U_first(i)*U_first(i)*G(i,i)。endNik%先求非對角元素（M）%Mik=zeros(nm1,n1)。for i=1:nm1 for k=1:n1 if i~=k theik=Theta_first(i)Theta_first(k)。
公司管理相關(guān)推薦
鄂ICP備17016276號-1