灰色预测法程序语言

这个是电气系的编程演算预测程序。用于规划型使用。独家编写。

#include<iostream>//预测第100个数据
#include <math.h>
#include<string>
#include <iomanip>
using namespace std;
int v,d,k1;
void main()
{
cout<<"********************"<<endl;
cout<<" 灰色预测法 "<<endl;
cout<<"********************"<<endl;
for(int pop=1;pop>0;pop++)//是程序可以循环使用
{
double a[100],c[100],J[100],rr[100],uu[100];
string answer;
int dd2=0,nian;
double matrix_k[100][100]={0},matrix_h[100][100]={0},matrix_U[2][2],matrix_Y[2][100]={0},matrix_C[100][1]={0},matrix_X[2][2]={0},matrix_B[100][2]={0};
cout<<"是否使用本软件?【yes/no】"<<endl;
cin>>answer;
if(answer=="yes")
{
cout<<"已知用电量数据对应的末年份"<<endl;
cin>>nian;
cout<<"输入用电量（最多能输入100个数据且结束时【不够100个数据】输入123.123）："<<'\n';
for(int kk=0;kk<10;kk++)
{
cin>>a[kk];
if(a[kk]==123.123)
break;
}
int h=kk;
for(int k1=0;k1<h;k1++)
{
J[k1]=a[k1];
}

doubleb=0;
for(int dd1=0;dd1<h;dd1++)
{
while(dd2<=dd1)
{
b=b+J[dd2];
dd2++;
}
rr[dd1]=b;
}
int H=h-1;
for(int GH=0;GH<H;GH++)
{
c[GH]=0;
}
for(int oo=0;oo<H;oo++)//c[]
{
c[oo]=-0.5*(rr[oo]+rr[oo+1]);
}
cout<<setiosflags(ios::fixed)<<setprecision(3);



for(int j=0;j<H;j++)
{
matrix_B[j][0]=c[j];
}
for(int u=0;u<100;u++)
uu[u]=1;
for(int jp=0;jp<H;jp++)
{
matrix_B[jp][1]=uu[jp];
} //ok
for(int mm=0;mm<2;mm++)
{
for(int mdd=0;mdd<2;mdd++)
{
matrix_k[mm][mdd]=matrix_B[mdd][mm];
}
}
for(int mm0=2;mm0<H;mm0++)
{
for(int mdd0=0;mdd0<2;mdd0++)
{
matrix_k[mdd0][mm0]=matrix_B[mm0][mdd0];
}
}


for(int f=0;f<2;f++)
{
for(int fd=0;fd<2;fd++)
{
for(int r=0;r<H;r++)
{
matrix_h[f][fd]+=matrix_k[f][r]*matrix_B[r][fd];
}
}

}

for(int mt=0;mt<2;mt++)
{
for(int mt1=0;mt1<2;mt1++)
{
matrix_U[mt][mt1]=matrix_h[mt][mt1];
}
}
double t;
for(int yu=0;yu<1;yu++)
{
for(int yu1=0;yu1<1;yu1++)
{
t=matrix_U[yu][yu1]*matrix_U[yu+1][yu+1]-matrix_U[yu][yu1+1]*matrix_U[yu+1][yu1];
}
}
int Q=0;
for(int e=0;e<1;e++)
{
for(int e1=0;e1<1;e1++)
{
Q=matrix_U[e+1][e1+1];
matrix_U[e+1][e1+1]=matrix_U[e][e1];
matrix_U[e][e1]=Q;
}
}
for(int ee=0;ee<1;ee++)
{
for(int ee1=0;ee1<1;ee1++)
{
matrix_U[ee][ee1+1]=-matrix_U[ee][ee1+1];
matrix_U[ee+1][ee1]=-matri
x_U[ee+1][ee1];
}
}
for(int t2=0;t2<2;t2++)
{
for(int t3=0;t3<2;t3++)
{
matrix_U[t2][t3]=matrix_U[t2][t3]/t;
}
}

这个是电气系的编程演算预测程序。用于规划型使用。独家编写。

for(int f22=0;f22<2;f22++)
{
for(int fdd=0;fdd<H;fdd++)
{
for(int rd=0;rd<H;rd++)
{
matrix_Y[f22][fdd]+=matrix_U[f22][rd]*matrix_k[rd][fdd];
}
}
}
cout<<endl;
for(int gg=0;gg<10;gg++)
{
matrix_C[gg][0]=a[gg+1];
if(matrix_C[gg][0]==123.123)
matrix_C[gg][0]=0;
}
for(int z=0;z<2;z++)
{
for(int z1=0;z1<1;z1++)
{
for(int z2=0;z2<H;z2++)
{
matrix_X[z][z1]+=matrix_Y[z][z2]*matrix_C[z2][z1];
}
}
}
double n,nn,tr,GG;
for(int hu=0;hu<2;hu++)
{
if(hu==0)
n=matrix_X[hu][0];
else if(hu==1)
nn=matrix_X[hu][0];
}
int ko=h+5;
for(;h<ko;h++)
{
H=h-1;
for(int jj=H;jj<h;jj++)
{
a[jj+1]=(a[0]-nn/n)*(exp(-n*h))+nn/n;
tr=a[jj+1];
}
for(int je=H;je<h;je++)
{
nian++;
GG=tr-rr[je];
rr[je+1]=tr;
cout<<"预测"<<nian<<"的用电量为："<<GG<<endl;

}
}
}
else if(answer=="
no")
{
cout<<"谢谢使用。下次见!"<<endl;
break;
}
else
{
cout<<"错误!!"<<endl;
break;
}

}
}