c語(yǔ)言計(jì)算閏年程序（推薦）-文庫(kù)吧
2024-11-15 12:03 本頁(yè)面
【正文】 }第三篇：C語(yǔ)言判斷閏年和平年C語(yǔ)言判斷閏年和平年include void fun(int year){ if((year%4==0amp。amp。year%100!=0)||(year%100==0amp。amp。year%400==0)){printf(“%d年是閏年n”,year)。} else {printf(“%d年不是閏年n”,year)。} } void main(){ int Year。scanf(“%d”,amp。Year)。fun(Year)。}第四篇：電力系統(tǒng)通用潮流計(jì)算C語(yǔ)言程序include include include include using namespace std。//節(jié)點(diǎn)號(hào)類型負(fù)荷有功負(fù)荷無功母線數(shù)據(jù)（類型 1＝PV節(jié)點(diǎn)，2＝PQ節(jié)點(diǎn)，3＝平衡節(jié)點(diǎn)）struct BUS { int busno。int type。float Pd。float Qd。}。//發(fā)電機(jī)數(shù)據(jù) 節(jié)點(diǎn)號(hào)有功發(fā)電電壓幅值 struct Generator { int busno。float Pg。float Vg。}。//支路信息 節(jié)點(diǎn)I 節(jié)點(diǎn)J R X B/2 k struct Line { int busi。int busj。float R。float X。float B。float k。}。//deltaP deltaQ deltaV^2// void fun1(double YG[][50],double YB[][50],double e[],double f[],int type[],int N,double W[],double P[],double Q[],double V[]){ double dP=0,dQ=0,dV=0。int i,j。for(i=0。idouble A=0,B=0。for(j=0。jA+=YG[i][j]*e[j]YB[i][j]*f[j]。B+=YG[i][j]*f[j]+YB[i][j]*e[j]。}dV=V[i]*V[i]e[i]*e[i]f[i]*f[i]。dP=P[i]e[i]*Af[i]*B。W[2*i]=dP。dQ=Q[i]f[i]*A+e[i]*B。if(type[i]==1)W[2*i+1]=dQ。else W[2*i+1]=dV。} } //Jacobi矩陣// void Jacobi(double YG[][50],double YB[][50],double e[50],double f[50],int type[50],int N ,double Ja[100][101]){int i,j。for(i=0。ifor(j=0。jif(i!=j){if(type[i]==1){Ja[2*i][2*j]=(YG[i][j]*e[i]+YB[i][j]*f[i])。Ja[2*i][2*j+1]=YB[i][j]*e[i]YG[i][j]*f[i]。Ja[2*i+1][2*j]=Ja[2*i][2*j+1]。Ja[2*i+1][2*j+1]=Ja[2*i][2*j]。}else {Ja[2*i][2*j]=YG[i][j]*e[i]+YB[i][j]*f[i]。Ja[2*i][2*j+1]=YB[i][j]*e[i]YG[i][j]*f[i]。Ja[2*i+1][2*j+1]=Ja[2*i+1][2*j]=0。}}else {double a[50]={0},b[50]={0}。for(int k=0。ka[i]+=(YG[i][k]*e[k]YB[i][k]*f[k])。b[i]+=(YG[i][k]*f[k]+YB[i][k]*e[k])。Ja[2*i][2*j]=a[i]YG[i][i]*e[i]YB[i][i]*f[i]。Ja[2*i][2*j+1]=b[i]+YB[i][i]*e[i]YG[i][i]*f[i]。if(type[i]==1){Ja[2*i+1][2*j]=b[i]+YB[i][i]*e[i]YG[i][i]*f[i]。Ja[2*i+1][2*j+1]=a[i]+YG[i][i]*e[i]+YB[i][i]*f[i]。}else {Ja[2*i+1][2*j]=2*e[i]。Ja[2*i+1][2*j+1]=2*f[i]。}}}}} }//高斯消元法解方程組函數(shù)// void gauss(double a[][101],int n){int i,j,k。double c。for(k=0。kc=a[k][k]。for(j=k。jfor(i=k+1。ic=a[i][k]。for(j=k。j}}a[n1][n]/=a[n1][n1]。for(k=n2。k=0。k)for(j=k+1。jvoid main(){ ifstream f
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1