开发工具:
文件大小: 1kb
下载次数: 0
上传时间: 2015-11-28
详细说明： #include #include using namespace std; int main() { double a[] = { 4, 2, -2, 2, 2, -3, -2, -3, 14 }; double *L = NULL; double *U = NULL; double tmp=0; int n; //矩阵内总数据个数 int s; //矩阵的阶数 n = sizeof(a)/sizeof(double); s = sqrtl(n); L = new double[n]; U = new double[n]; for (int i = 0; i < s; i++) { for (int j = 0; j < s; j++) { if (i = = j) L[i*s+j] = 1; if (i < j) L[i*s+j] = 0; if (i > j) U[i*s+j] = 0; U[0*s+j] = a[0*s+j]; L[i*s+0] = a[i*s+0] / U[0*s+0]; } } for (int k = 1; k < s; k++) { for (int j = k; j < s; j++) { tmp = 0; for (int m = 0; m < k; m++) { tmp += L[k*s+m] * U[m*s+j]; } U[k*s+j] = a[k*s+j] - tmp; } for (int i = k+1; i < s; i++) { tmp = 0; for (int m = 0; m < k; m++) { tmp += L[i*s+m] * U[m*s+k]; } L[i*s+k] = ( a[i*s+k] - tmp ) / U[k*s+k]; } } //这里就得到L和U矩阵的值了 delete [] L; delete [] U; system("pause"); return 0; } ...展开收缩
(系统自动生成,下载前可以参看下载内容)
下载文件列表
相关说明
- 本站资源为会员上传分享交流与学习,如有侵犯您的权益,请联系我们删除.
- 本站是交换下载平台，提供交流渠道，下载内容来自于网络，除下载问题外，其它问题请自行百度。
- 本站已设置防盗链，请勿用迅雷、QQ旋风等多线程下载软件下载资源，下载后用WinRAR最新版进行解压.
- 如果您发现内容无法下载，请稍后再次尝试；或者到消费记录里找到下载记录反馈给我们.
- 下载后发现下载的内容跟说明不相乎，请到消费记录里找到下载记录反馈给我们，经确认后退回积分.
- 如下载前有疑问，可以通过点击"提供者"的名字，查看对方的联系方式，联系对方咨询.
