请教下大家,如何用visual c++进行高斯消元法求解呢?谢谢!
2条回答
- threenewbee 2015-07-25 15:10关注
列选主元消元法 /* *Gauss.cpp *功能:列选主元消元法 课本3.5为例 *时间:2013.4.1 */ #include<stdio.h> #include"Gass.h" //高斯消元法(列选主元) void Gauss (double a[][MAXNUM],int n) { int i,j; SelectColE(a,n); //列选主元并消元成上三角 //回代求解 printf("上三角的结果\n"); printM(a,3); for(i=n;i>=1;i--) { for(j=i+1;j<=n;j++) a[i][n+1]-=a[i][j]*a[j][n+1]; a[i][n+1]/=a[i][i]; } return ; } //选择列主元并进行消元 void SelectColE(double a[][MAXNUM],int n) { int i,j,k,maxRowE; double temp; //用于记录消元时的因数 for(j=1;j<=n;j++) { maxRowE=j; for(i=j;i<=n;i++) if(fabs(a[i][j])>fabs(a[maxRowE][j])) maxRowE = i; if(maxRowE!=j) swapRow(a,j,maxRowE,n); //与最大主元所在行交换 //消元 for(i=j+1;i<=n;i++) { temp =a[i][j]/a[j][j]; for(k=j;k<=n+1;k++) a[i][k]-=a[j][k]*temp; } printf("第%d列消元后:\n",j); printM(a,3); } return; } void swapRow(double a[][MAXNUM],int m,int maxRowE,int n) { int k; double temp; for(k=m;k<=n+1;k++) { temp = a[m][k]; a[m][k] = a[maxRowE][k]; a[maxRowE][k] = temp; } return ; } //测试函数 int main() { double a[4][MAXNUM]; int i,n,j; a[1][1] = 0.001; a[1][2]=2.000;a[1][3]=3.000;a[1][4]=1.000; a[2][1] = -1.000;a[2][2]=3.712;a[2][3]=4.623;a[2][4]=2.000; a[3][1] = -2.000;a[3][2]=1.070;a[3][3]=5.643;a[3][4]=3.000; Gauss (a,3); for(i=1;i<=3;i++) printf("X%d=%f\n",i,a[i][4]); return 0; } //输出矩阵 void printM(double a[][MAXNUM], int n) { int i,j; for(i=1;i<=n;i++) { for(j=1;j<=n+1;j++) { printf("%f ,",a[i][j]); } printf("\n"); } }
本回答被题主选为最佳回答 , 对您是否有帮助呢?解决 无用评论 打赏 举报
悬赏问题
- ¥15 关于#matlab#的问题:在模糊控制器中选出线路信息,在simulink中根据线路信息生成速度时间目标曲线(初速度为20m/s,15秒后减为0的速度时间图像)我想问线路信息是什么
- ¥15 banner广告展示设置多少时间不怎么会消耗用户价值
- ¥16 mybatis的代理对象无法通过@Autowired装填
- ¥15 可见光定位matlab仿真
- ¥15 arduino 四自由度机械臂
- ¥15 wordpress 产品图片 GIF 没法显示
- ¥15 求三国群英传pl国战时间的修改方法
- ¥15 matlab代码代写,需写出详细代码,代价私
- ¥15 ROS系统搭建请教(跨境电商用途)
- ¥15 AIC3204的示例代码有吗,想用AIC3204测量血氧,找不到相关的代码。