evluation.cpp

数据结构 数学计算 数 据结构 数学计算 数据结构 数学计算

#include <iostream.h> 
#include<iomanip.h>
struct element
{
	int I,J,K,IT;
	double VALUE;
}e[10];
void main()
{
	int n=2;
	double h=0.01,Tmax=10;
	cout<<setw(14)<<"T(s)"<<setw(13)<<"IL3"<<setw(15)<<"UC5";
	cout<<endl;
	cout<<setw(13)<<0<<setw(15)<<0<<setw(16)<<1<<endl;
	e[1].J=1;e[1].K=0;e[1].IT=2;e[1].VALUE=1;
	e[2].J=1;e[2].K=0;e[2].IT=1;e[2].VALUE=20;
	e[3].J=1;e[3].K=0;e[3].IT=5;e[3].VALUE=2;
	e[3].IT=1;
	e[3].VALUE=2*e[3].VALUE/h;
	e[4].J=e[3].J;
	e[4].K=e[3].K;
	e[4].IT=2;
	e[4].VALUE=-30.00/(11*e[3].VALUE);
	e[5].J=1;e[5].K=2;e[5].IT=1;e[5].VALUE=2;
	e[6].J=2;e[6].K=0;e[6].IT=4;e[6].VALUE=0.2;
	e[6].VALUE=h/(2*e[6].VALUE);
	e[7].J=e[6].J;
	e[7].K=e[6].K;
	e[7].IT=2;
	e[7].VALUE=19.00/22.00+1/e[6].VALUE;
	double x[2];
	x[0]=30.00/11.00;
	x[1]=1;
	double a[2][3];
	for(int i=0;i<2;i++)
		for(int j=0;j<3;j++)
			a[i][j]=0;
	for(double c=0;c<Tmax;c+=h)
	{
		if(c>0)
		{
			e[7].VALUE=(x[0]-x[1])/2+x[1]/e[6].VALUE;
			e[4].VALUE=-(1-(x[0]-x[1])/2-x[0]/20+x[0]/e[3].VALUE);
		}else;
		a[0][0]=1/e[2].VALUE+1/e[3].VALUE+1/e[5].VALUE;
		a[0][1]=-1/e[5].VALUE;
		a[0][2]=e[1].VALUE+e[4].VALUE;
		a[1][0]=a[0][1];
		a[1][1]=1/e[5].VALUE+1/e[6].VALUE;
		a[1][2]=e[7].VALUE;
		for(int j=1;j<=2;j++)
			a[1][j]=a[1][j]-a[0][j]*a[1][0]/a[0][0];			
		x[1]=a[1][2]/a[1][1];
		x[0]=(a[0][2]-a[0][1]*x[1])/a[0][0];
		if(int(c/h+1)%20==0)
		{
			if(c>2)	
				cout<<setw(13)<<c<<setw(15)<<1-(x[0]-x[1])/2-x[0]/20<<setw(16)<<x[1];
			else
				cout<<setw(13)<<c+h<<setw(15)<<1-(x[0]-x[1])/2-x[0]/20<<setw(16)<<x[1];
			cout<<endl;
		}
	}
}