c語言計(jì)算閏年程序推薦(編輯修改稿)
2024-11-15 12:03 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 in。int N=0,GS=0,LD=0,ZLs=0。//節(jié)點(diǎn)數(shù)發(fā)電機(jī)數(shù)負(fù)荷數(shù) BUS*B。Generator*G。Line*L。//從文本中讀入原始數(shù)據(jù)到數(shù)組中//(“C:”)。if(!fin){coutgetchar()。}int m1[50]={0},m2[50]={0}。float m3[50],m4[50],m5[50],m6[50]。int i,j,l。支路數(shù)//for(i=0。i++){finm1[i]。if(m1[i]==0)break。finm2[i]m3[i]m4[i]。N++。} B=new BUS[N]。for(i=0。iB[i].busno=m1[i]。B[i].type=m2[i]。B[i].Pd=m3[i]。B[i].Qd=m4[i]。}for(i=0。i++){finm1[i]。if(m1[i]==0)break。finm4[i]m3[i]。GS++。} G=new Generator[GS]。for(i=0。iG[i].busno=m1[i]。G[i].Pg=m4[i]。G[i].Vg=m3[i]。}for(i=0。i++){finm1[i]。if(m1[i]==0)break。finm2[i]m3[i]m4[i]m5[i]m6[i]。ZLs++。}L=new Line[ZLs]。for(i=0。iL[i].busi=m1[i]。L[i].busj=m2[i]。L[i].R=m3[i]。L[i].X=m4[i]。} L[i].B=m5[i]。L[i].k=m6[i]。LD=NGS。()。//節(jié)點(diǎn)導(dǎo)納矩陣形成// double YB[50][50],YG[50][50],BB[50][50],K[50][50]。for(i=0。iYB[i][j]=0。YG[i][j]=0。BB[i][j]=0。K[i][j]=1。} } for(l=0。l} for(i=0。ifor(j=i。jK[i][j]=K[j][i]。K[j][i]=1。}for(j=0。jif(i!=j){YG[i][i]=YG[i][i]+(YG[i][j]*K[i][j]*K[i][j])。YB[i][i]=YB[i][i]+(YB[i][j]*K[i][j]*K[i][j]+BB[i][j])。}}} //修正后// for(l=0。l// 求A=e+f// double e[50]={0},f[50]={0}。double C[100]={0},D[100]={0}。for(i=0。iC[2*i]=1。} else C[2*i]=V[i]。}double W[100]={0},Ja[100][101]={0}。//調(diào)用Jacobi函數(shù)和高斯函數(shù)//for(int t=1。tfor(i=0。ie[i]=C[2*i]。f[i]=C[2*i+1]。}fun1(YG,YB,e,f,type,N,W,P,Q,V)。double it=fabs(W[0])。for(i=1。iif(itit=fabs(W[i])。j=i。}}//中間迭代過程//coutcoutif(itJacobi(YG,YB,e,f,type,N,Ja)。for(i=0。iJa[i][2*N2]=W[i]。}//高斯消元法解方程//gauss(Ja,2*N2)。for(i=0。iD[i]=Ja[i][2*(N1)]。C[i]+=D[i]。}}//平衡節(jié)點(diǎn)//for(i=0。idouble a=0,b=0。for(int j=0。ja+=(YG[i][j]*e[j]YB[i][j]*f[j])。b+=(YB[i][j]*e[j]+YG[i][j]*f[j])。}P[i]=e[i]*a+f[i]*b。Q[i]=f[i]*ae[i]*b。}//支路//double PZL[100][101]={0},QZL[100][101]={0},pr[100][101]={0},qx[100][101]={0}。double x1=0,x2=0,y1=0,y2=0,I2=0。for(int k=0。ki=L[k].busi1。j=L[k].busj1。x1=e[i]/L[k].ke[j]。y1=f[i]/L[k].kf[j]。x2=e[i]*YG[i][j]f[i]*YB[i][j]。y2=f[i]*YG[i][j]+e[i]*YB[i][j]。QZL[i][j]=(x1*y2x2*y1)。PZL[i][j]=(x1*x2+y1*y2)。I2=(PZL[i][j]*PZL[i][j]+QZL[i][j]*QZL[i][j])/(e[i]*e[i]+f[i]*f[i])。pr[i][j
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1