牛頓-拉夫遜法潮流計(jì)算課程設(shè)計(jì)說(shuō)明書(shū)(參考版)
2025-08-05 21:06本頁(yè)面
【正文】 33 。printf(P3+jQ3=%f+j%f,P3,Q3)。}P3=e[3]*a1+f[3]*b1。j++){a1+=G[3][j]*e[j]B[3][j]*f[j]。}for(j=1。k1=k1+1。 printf( ,)。i++){printf(f[%d]=,i)。}for(i=1。 printf(%f,e[i])。i3。 f[i]+=xf[i]。i3。 printf(節(jié)點(diǎn)電壓:\n)。 xf[1]=x[2]。}xe[1]=x[1]。n++)sigma2+=U[i][n]*x[n]。for(n=i+1。i=1。y[i]=b[i]sigma1。n=i1。i++){sigma1=0。for(i=1。b[3]=P[2]。}}b[1]=P[1]。s++)sigma2+=L[i][s]*U[s][n]。for(s=0。iM。U[n][j]=J[n][j]sigma1。s=n1。j++){sigma1=0。n++){for(j=n。}for(n=2。i++){U[1][i]=J[1][i]。}for(i=1。iM。}printf(\n)。j++){printf(%f,J[i][j])。i++){for(j=1。for(i=1。J[2*N][j*N1]=0。J[2*N][j*N]=0。J[2*N][j*N]=2*f[2]。J[2*N1][j*N]=d+B[2][j]*e[2]G[2][j]*f[2]。 d+=G[2][m]*f[m]+B[2][m]*e[m]。m4。int m。j3。 J[1*N][j*N1]=B[1][j]*e[1]G[1][j]*f[1]。 J[1*N][j*N]=G[1][j]*e[1]+B[1][j]*f[1]。J[1*N][j*N]=c+G[1][j]*e[1]+B[1][j]*f[1]。J[1*N1][j*N]=d+B[1][j]*e[1]G[1][j]*f[1]。 d+=G[1][m]*f[m]+B[1][m]*e[m]。m4。int m。j3。printf(max=%f\n,max)。max=maxfabs(Q[1])?max:fabs(Q[1])。printf(\n)。printf(\n)。printf(無(wú)功功率增量Q[1]=%f,Q[1])。 printf( ,)。printf(\n)。}printf(有功功率增量P[1]=%f,P[1])。//計(jì)算有功功率的增量Q[i]=Qs[i](f[i]*ae[i]*b)。 b+=G[i][j]*f[j]+B[i][j]*e[j]。j4。i++){float a=0,b=0。for (i=1。//形成節(jié)點(diǎn)導(dǎo)納矩陣//*******************************************}printf(\n)。printf( )。j4。i4。B[2][3]=B[3][2]=4。B[1][3]=B[3][1]=。B[1][2]=B[2][1]=3。B[3][3]=。B[2][2]=7。B[1][1]=。float L[M][M]={0},U[M][M]={0},sigma1,sigma2,b[M],y[M],x[M]。float a1=0,b1=0。float e[4]={0,1,1,1},f[4]={0},P[4],Q[4],Ps[4]={0,2,},xe[3],xf[3]。做好每一個(gè)細(xì)節(jié)，不斷完善自我，提高自身的學(xué)習(xí)的水平。經(jīng)過(guò)了自己的改造將程序完整的編輯了出來(lái)，并實(shí)現(xiàn)的預(yù)期的功能!通過(guò)本次課程設(shè)計(jì)，我知道了自己學(xué)習(xí)的知識(shí)還不夠扎實(shí)！很多方面只是應(yīng)付考試，到了讓你做東西的時(shí)候確實(shí)還是相當(dāng)困難的!尤其是在編程方面的缺陷！現(xiàn)如今是一個(gè)軟硬件相結(jié)合的時(shí)代，其中軟件更具有競(jìng)爭(zhēng)性。最后，遇到了最難難的地方：用程序來(lái)實(shí)現(xiàn)牛頓拉夫遜算法！起初是想用matlab程序來(lái)進(jìn)行程序的編寫(xiě)并進(jìn)行仿真，但是查找了一些資料以后，由于自己對(duì)于matlab的了解并不深而且自己學(xué)習(xí)的知識(shí)也不夠扎實(shí)，給自己帶來(lái)了不小的麻煩最終無(wú)法完成牛頓拉夫遜的算法。大體上知道了運(yùn)用牛頓拉夫遜法的各個(gè)步驟！其次，讀題的過(guò)程中也遇到了一些麻煩:“”和“1”分別指代的是什么？？通過(guò)自己上網(wǎng)查資料，在網(wǎng)上看到了一些人有關(guān)的說(shuō)法以及相似的題目的解答，這給我不小的啟發(fā)！我認(rèn)識(shí)到了查找資料的重要性！經(jīng)過(guò)了自己的努力，我知道了以上的答案：“”是有功功率，“1”代表的是電壓幅值。首先，在拿到本次課設(shè)的題目時(shí)，就看不懂題目！這就給了自己一個(gè)不小的打擊。因此，一般認(rèn)為算法不是導(dǎo)致不收斂的原因，潮流不收斂產(chǎn)生的主要原因是計(jì)算的初始條件給得不合理，導(dǎo)致潮流方程無(wú)解。給定節(jié)點(diǎn)電壓初值,經(jīng)過(guò)五次迭代過(guò)程后，得到程序的顯示結(jié)果如下（?。赫砜傻霉?jié)點(diǎn)電壓和不平衡功率的變化情況，分別于表1和表2所示：迭代計(jì)數(shù)k節(jié)點(diǎn)電壓1121314151表1 迭代過(guò)程中節(jié)點(diǎn)電壓變化情況迭代計(jì)數(shù)k節(jié)點(diǎn)不平衡量012345表2 迭代過(guò)程中節(jié)點(diǎn)不平衡量變化情況進(jìn)行了五次迭代，結(jié)果仍然沒(méi)有收斂。 printf( ,)。i++){printf(f[%d]=,i)。}for(i=1。 printf(%f,e[i])。i3。 f[i]+=xf[i]。i3。 printf(節(jié)點(diǎn)電壓:\n)。 xf[1]=x
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1