亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频

蟲蟲首頁| 資源下載| 資源專輯| 精品軟件
登錄| 注冊

您現(xiàn)在的位置是：蟲蟲下載站 > 資源下載 > 軟件 > 道理特分解法

道理特分解法

  • 資源大?。?/b>18 K
  • 上傳時(shí)間： 2018-05-20
  • 上傳用戶：Aa123456789
  • 資源積分：2 下載積分
  • 標(biāo)      簽： 道理特分解法

資 源 簡 介

#include "iostream" using namespace std;

class Matrix

{

private:

double** A; //矩陣A

double *b; //向量b

public:

int size;

Matrix(int );

~Matrix();

friend double* Dooli(Matrix& );

void Input();

void Disp();

};

Matrix::Matrix(int x) {

size=x;

//為向量b分配空間并初始化為0

b=new double [x];

for(int j=0;j<x;j++)

b[j]=0;

//為向量A分配空間并初始化為0

A=new double* [x];

for(int i=0;i<x;i++)

A[i]=new double [x];

for(int m=0;m<x;m++)

for(int n=0;n<x;n++)

A[m][n]=0;

}

Matrix::~Matrix() {

cout<<"正在析構(gòu)中~~~~"<<endl;

delete b;

for(int i=0;i<size;i++)

delete A[i];

delete A;

}

void Matrix::Disp()

{

for(int i=0;i<size;i++)

{

for(int j=0;j<size;j++)

cout<<A[i][j]<<" ";

cout<<endl;

}

}

void Matrix::Input()

{

cout<<"請輸入A:"<<endl;

for(int i=0;i<size;i++)

for(int j=0;j<size;j++){

cout<<"第"<<i+1<<"行"<<"第"<<j+1<<"列:"<<endl;

cin>>A[i][j];

}

cout<<"請輸入b:"<<endl;

for(int j=0;j<size;j++){

cout<<"第"<<j+1<<"個(gè):"<<endl;

cin>>b[j];

}

}

double* Dooli(Matrix& A) {

double *Xn=new double [A.size];

Matrix L(A.size),U(A.size);

//分別求得U，L的第一行與第一列

for(int i=0;i<A.size;i++)

U.A[0][i]=A.A[0][i];

for(int j=1;j<A.size;j++)

L.A[j][0]=A.A[j][0]/U.A[0][0];

//分別求得U，L的第r行，第r列

double temp1=0,temp2=0;

for(int r=1;r<A.size;r++){

//U

for(int i=r;i<A.size;i++){

for(int k=0;k<r-1;k++)

temp1=temp1+L.A[r][k]*U.A[k][i];

U.A[r][i]=A.A[r][i]-temp1;

}

//L

for(int i=r+1;i<A.size;i++){

for(int k=0;k<r-1;k++)

temp2=temp2+L.A[i][k]*U.A[k][r];

L.A[i][r]=(A.A[i][r]-temp2)/U.A[r][r];

}

}

cout<<"計(jì)算U得:"<<endl;

U.Disp();

cout<<"計(jì)算L的:"<<endl;

L.Disp();

double *Y=new double [A.size];

Y[0]=A.b[0];

for(int i=1;i<A.size;i++ ){

double temp3=0;

for(int k=0;k<i-1;k++)

temp3=temp3+L.A[i][k]*Y[k];

Y[i]=A.b[i]-temp3;

}

Xn[A.size-1]=Y[A.size-1]/U.A[A.size-1][A.size-1];

for(int i=A.size-1;i>=0;i--){

double temp4=0;

for(int k=i+1;k<A.size;k++)

temp4=temp4+U.A[i][k]*Xn[k];

Xn[i]=(Y[i]-temp4)/U.A[i][i];

}

return Xn;

}

int main()

{

Matrix B(4);

B.Input();

double *X;

X=Dooli(B);

cout<<"~~~~解得:"<<endl;

for(int i=0;i<B.size;i++)

cout<<"X["<<i<<"]:"<<X[i]<<" ";

cout<<endl<<"呵呵呵呵呵";

return 0;

} 

相 關(guān) 資 源

主站蜘蛛池模板： 晋宁县| 梓潼县| 库伦旗| 松阳县| 娄底市| 达拉特旗| 邵阳县| 陆川县| 台江县| 荣成市| 扶沟县| 五大连池市| 汝阳县| 小金县| 饶河县| 吴川市| 神农架林区| 台东县| 丰城市| 庆云县| 东光县| 日喀则市| 赣榆县| 嘉定区| 兰坪| 梁河县| 若尔盖县| 茂名市| 鲜城| 兴宁市| 扶风县| 东阿县| 屏边| 临猗县| 沾化县| 齐齐哈尔市| 光泽县| 百色市| 潢川县| 唐河县| 嵊泗县|