电力系统稳态实验报告(7)

34. deta_p(i)=pis(i)-p(i); 35. end 36. for i=1:3 37. for j=1:5

38. q(i)=q(i)+f(i)*(G(i,j)*e(j)-B(i,j)*f(j))-e(i)*(G(i,j)*f(j)+B(i,j)*e(j)); 39. end

40. deta_q(i)=qis(i)-q(i); 41. end

42. deta_UU=1.05*1.05-(e(4)*e(4)+f(4)*f(4)); 43. jacobi=zeros(8,8); 44. for i=1:4 45. for j=1:4

46. jacobi(2*i-1,2*j-1)=-G(i,j)*e(i)-B(i,j)*f(i); 47. jacobi(2*i-1,2*j)=B(i,j)*e(i)-G(i,j)*f(i); 48. end 49. end 50. for i=1:3 51. for j=1:4

52. jacobi(2*i,2*j-1)=B(i,j)*e(i)-G(i,j)*f(i); 53. jacobi(2*i,2*j)=G(i,j)*e(i)+B(i,j)*f(i); 54. end 55. end

56. for i=1:2:7

57. jacobi(i,i)=0; 58. jacobi(i,i+1)=0; 59. for j=1:5

60. jacobi(i,i)=jacobi(i,i)-(G((i+1)/2,j)*e(j)-B((i+1)/2,j)*f(j)); 61. jacobi(i,i+1)=jacobi(i,i+1)-(G((i+1)/2,j)*f(j)+B((i+1)/2,j)*e(j)); 62. end 63.

jacobi(i,i)=jacobi(i,i)-G((i+1)/2,(i+1)/2)*e((i+1)/2)-B((i+1)/2,(i+1)/2)*f((i+1)/2); 64.

jacobi(i,i+1)=jacobi(i,i+1)+B((i+1)/2,(i+1)/2)*e((i+1)/2)-G((i+1)/2,(i+1)/2)*f((i+1)/2); 65. end

66. for i=2:2:6

67. jacobi(i,i-1)=0; 68. jacobi(i,i)=0; 69. for j=1:5

70. jacobi(i,i-1)=jacobi(i,i-1)-(G(i/2,j)*f(j)-B(i/2,j)*e(j)); 71. jacobi(i,i)=jacobi(i,i)-(G(i/2,j)*e(j)-B(i/2,j)*f(j)); 72. end

73. jacobi(i,i-1)=jacobi(i,i-1)+B(i/2,i/2)*e(i/2)-G(i/2,i/2)*f(i/2); 74. jacobi(i,i)=jacobi(i,i)+G(i/2,i/2)*e(i/2)+B(i/2,i/2)*f(i/2);