王朝网络
分享
 
 
 

[数值算法]线性方程组求解算法---基于LU分解法的追赶法

王朝other·作者佚名  2006-02-01
宽屏版  字体: |||超大  

[数值算法]线性方程组求解算法---基于LU分解法的追赶法

By EmilMatthew

05/9/08

基于LU分解法的追赶法主要用于求解一类形式上特殊的线性方程组。

Ax=b,其中,当|i-j|>1时,有aij=0.求解这类形式的线性方程组,在样条函数,常微分方程边值问题的数值解中都会有用.

原理上仍是考虑将这个矩阵分成L和U两个矩阵,不过,由于这个矩阵存大着大量的0值元素,所以,给最终的递推求解带来了相当大的便利.

更加详细的内容请去参考相应的数值算法书籍.

/*This function :chasingMethod is coded by EmilMatthew,you can use and modify it as you wish,

But I do not pomise it can fit all your need .*/

void chasingMethod(Type** inMatrixArr,Type* bList,Type* xAnsList,int size)

{

/*Core think of this algorithm,this algorithm is just a special version of the LUPationMethod:

matrix_L*yAnsList=bList;

matrix_U*xAnsList=yAnsList;

*/

Type* listA,* listB,* listC;/*to crack the in matrix in three one demision list*/

Type* listAlpha,* listBeta;/*The tween list for resolve the problem*/

Type* yAnsList;

int i,j;

/*pointer data assertion*/

assertF(inMatrixArr!=NULL,"in LUPationMethod,matrixArr is NULL\n");

assertF(bList!=NULL,"in LUPationMethod,bList is NULL\n");

assertF(xAnsList!=NULL,"in LUPationMethod,xAnsList is NULL\n");

/*The Assertion of the pass in matrix's format.*/

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

for(j=i;j<size;j++)

if(j-i>1)

assertF(inMatrixArr[i][j]==0,"In chasingMethod,the pass in matrix format is not correct\n");

/*Mem Apply*/

listArrMemApply(&listA,size);

listArrMemApply(&listB,size);

listArrMemApply(&listC,size);

listArrMemApply(&listAlpha,size);

listArrMemApply(&listBeta,size);

listArrMemApply(&yAnsList,size);

/*locate the inMatrixArr data into correct position*/

listB[0]=inMatrixArr[0][0];

listC[0]=inMatrixArr[0][1];

for(i=1;i<size-1;i++)

{

listA[i]=inMatrixArr[i][i-1];

listB[i]=inMatrixArr[i][i];

listC[i]=inMatrixArr[i][i+1];

}

listA[size-1]=inMatrixArr[size-1][size-2];

listB[size-1]=inMatrixArr[size-1][size-1];

/*main program of the chasingMethod*/

listAlpha[0]=listB[0];

listBeta[0]=listC[0]/listAlpha[0];

for(i=1;i<size-1;i++)

{

listAlpha[i]=listB[i]-listA[i]*listBeta[i-1];

listBeta[i]=listC[i]/listAlpha[i];

}

listAlpha[i]=listB[i]-listA[i]*listBeta[i-1];

/*To get the yAnsList*/

yAnsList[0]=bList[0]/listAlpha[0];

for(i=1;i<size;i++)

yAnsList[i]=(bList[i]-listA[i]*yAnsList[i-1])/listAlpha[i];

showArrListFloat(yAnsList,0,size);

/*To get the xAnsList*/

xAnsList[size-1]=yAnsList[size-1];

for(i=size-2;i>=0;i--)

xAnsList[i]=yAnsList[i]-listBeta[i]*xAnsList[i+1];

/*Mem Free*/

free(listA);

free(listB);

free(listC);

free(listAlpha);

free(listBeta);

free(yAnsList);

}

/*算法简单测试*/

//inputData

3;

2,1,0,3;

1,3,1,5;

0,1,2,3;

//ouputData:

after the chasingMethod:the ans x rows is:(from x0 to xn-1)

1.00000 1.00000 1.00000

 
 
 
免责声明:本文为网络用户发布,其观点仅代表作者个人观点,与本站无关,本站仅提供信息存储服务。文中陈述内容未经本站证实,其真实性、完整性、及时性本站不作任何保证或承诺,请读者仅作参考,并请自行核实相关内容。
2023年上半年GDP全球前十五强
 百态   2023-10-24
美众议院议长启动对拜登的弹劾调查
 百态   2023-09-13
上海、济南、武汉等多地出现不明坠落物
 探索   2023-09-06
印度或要将国名改为“巴拉特”
 百态   2023-09-06
男子为女友送行,买票不登机被捕
 百态   2023-08-20
手机地震预警功能怎么开?
 干货   2023-08-06
女子4年卖2套房花700多万做美容:不但没变美脸,面部还出现变形
 百态   2023-08-04
住户一楼被水淹 还冲来8头猪
 百态   2023-07-31
女子体内爬出大量瓜子状活虫
 百态   2023-07-25
地球连续35年收到神秘规律性信号,网友:不要回答!
 探索   2023-07-21
全球镓价格本周大涨27%
 探索   2023-07-09
钱都流向了那些不缺钱的人,苦都留给了能吃苦的人
 探索   2023-07-02
倩女手游刀客魅者强控制(强混乱强眩晕强睡眠)和对应控制抗性的关系
 百态   2020-08-20
美国5月9日最新疫情:美国确诊人数突破131万
 百态   2020-05-09
荷兰政府宣布将集体辞职
 干货   2020-04-30
倩女幽魂手游师徒任务情义春秋猜成语答案逍遥观:鹏程万里
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案神机营:射石饮羽
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案昆仑山:拔刀相助
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案天工阁:鬼斧神工
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案丝路古道:单枪匹马
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案镇郊荒野:与虎谋皮
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案镇郊荒野:李代桃僵
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案镇郊荒野:指鹿为马
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案金陵:小鸟依人
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案金陵:千金买邻
 干货   2019-11-12
 
>>返回首页<<
推荐阅读
 
 
频道精选
 
静静地坐在废墟上,四周的荒凉一望无际,忽然觉得,凄凉也很美
© 2005- 王朝网络 版权所有