乔利斯基三角分解_杜利特尔及乔利斯基三角分解

/**********改进乔利斯基三角分解**********/

void Improved_Cholesky()

{

int i,j,k;

float t;

L=(float

**)malloc(sizeof(float)*(n+1));

for(i=0;i

L[i]=(float

*)malloc(sizeof(float)*(n+1));

T=(float

**)malloc(sizeof(float)*(n+1));

for(i=0;i

T[i]=(float

*)malloc(sizeof(float)*(n+1));

D=(float

*)malloc(sizeof(float)*(n+1));

D[1]=a[1][1];

for(i=2;i

{

for(j=1;j

{

t=0.0;

for(k=1;k

t=t+T[i][k]*L[j][k];

T[i][j]=a[i][j]-t;

L[i][j]=T[i][j]/D[j];

}

t=0.0;

for(k=1;k

t=t+T[i][k]*L[i][k];

D[i]=a[i][i]-t;

}

for(i=1;i

for(j=i;j

if(j==i)

T[i][j]=D[i];

else

T[i][j]=0.0;

for(i=1;i

for(j=i;j

if(j==i)

L[i][j]=1.0;

else

L[i][j]=0.0;

printf("\nArray T is:\n");

for(i=1;i

{

for(j=1;j

printf("%f ",T[i][j]);

printf("\n");

}

printf("\nArray L is:\n");

for(i=1;i

{

for(j=1;j

printf("%f ",L[i][j]);

printf("\n");

}

printf("\nArray D is:\n");

for(i=1;i

{

for(j=1;j

if(j==i)

printf("%f ",D[i]);

else

printf("%f ",D[i]-D[i]);

printf("\n");

}

/**********改进乔利斯基三角分解法解方程组**********/

void Improved_Method_Cholesky()

{

int i,j,k;

float t;

y[1]=b[1];

for(i=2;i

{

t=0.0;

for(k=1;k

t=t+L[i][k]*y[k];

y[i]=b[i]-t;

}

x[n]=y[n]/D[n];

for(i=n-1;i>0;i--)

{

t=0.0;

for(k=i+1;k

t=t+L[k][i]*x[k];

x[i]=y[i]/D[i]-t;

}

乔利斯基三角分解_杜利特尔及乔利斯基三角分解相关推荐

数值分析-杜利特尔分解法C语言代码
直接三角分解法如果方程组的系数矩阵可以分解为一个下三角矩阵和一个上三角矩阵的乘积,即:A=LU. Doolittle(杜利特尔)分解在对系数矩阵进行三角分解时,如果L为单位下三角矩阵(主对角线系数 ...
紧凑存储的杜利特尔分解法Doolittle(LU分解法)_解线性方程组的直接解法
紧凑存储的杜利特尔分解法Doolittle(LU分解法)_解线性方程组的直接解法标签:计算方法实验 /* 紧凑存储的杜利特尔分解法Doolittle:如果初始矩阵不要求保留的话,可以紧凑存储.因为每 ...
杜利特尔 (Doolittle)矩阵分解法求线性方程组的解
简介若方阵 A 可以分解为一个下三角矩阵 L 和一个上三角矩阵 U的乘积,即 A = LU,则这种分解称为 A 的一种三角分解或 LU分解.如果 L 为单位下三角矩阵,则称为杜利特尔 (Doolit ...
2 shell 锂基脂_【壳牌爱万利EP2润滑脂/Shell Alvania EP2锂基脂/黄油】价格_厂家 - 中国供应商...
壳牌爱万利 EP00 EP0 EP1 EP2 EP3 润滑脂(黄油) Shell Alvania EP 包装:16KG 180KG 壳牌爱万利EP润滑脂是同高粘 ...
贝索斯每年信_沙特王储黑杰夫·贝索斯的电话分析表明
贝索斯每年信 By Greg Farrell, Alyza Sebenius, and Matt Day Greg Farrell,Alyza Sebenius和Matt Day Jeff Bez ...
设计桑基图_教你用pyecharts制作交互式桑基图，赶快学起来吧！
由于工作的需要,经常需要进行可视化展示,除了一些常用的BI工具,我也会使用python对数据进行可视化. python的第三方可视化库有很多,比如matplotlib.seaborn.plotly.b ...
[oeasy]python0095_乔布斯求职_雅达利_atari_breakout_打砖块_布什内尔_游戏机_Jobs
编码进化回忆上次内容上次我们回顾了电子游戏的历史从电子游戏鼻祖双人网球到视频游戏 PingPong 再到街机游戏 Pong 雅达利公司来了嬉皮士捣乱?
图灵在布莱切利园的女朋友_极客参观布莱切利公园
图灵在布莱切利园的女朋友 Jeffrey Osier-Mixon和Dawn Foster共同撰写的文章当一些极客发现自己在都柏林参加LinuxCon Europe的途中会做什么? 他们当然去了布莱奇 ...
2 shell 锂基脂_壬二酸和癸二酸制备的复合锂基脂到底有那些差别！
2019年12月21-22日深圳<润滑脂生产技术培训班>烃基脂,钙基脂,锂基脂,复合锂基脂.铝基脂,复合铝,磺酸钙脂等润滑脂的生产配方设计及生产工艺和流程,专家现场答疑,报名咨询:189 ...

乔利斯基三角分解_杜利特尔及乔利斯基三角分解

乔利斯基三角分解_杜利特尔及乔利斯基三角分解相关推荐

最新文章

热门文章