应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 带平方根的（LLT）Cholesky算法分解对称正定矩阵 c语言实现
51/1返回列表
查看: 2747  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

舒马诺

银虫 (初入文坛)

[求助] 带平方根的（LLT）Cholesky算法分解对称正定矩阵 c语言实现

新人，不懂规矩，见谅

大神们好，要求
1先判断任意矩阵A是否为正定对称矩阵，否则，返回输入错误
2若输入为正定对称矩阵，则将其进行带平方根的（LLT）Cholesky算法分解，即实现A=LL^T，其中L为下三角形矩阵。

大致就这意思，求助

定理
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-03 20:37:32
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xiuyouxu

铁杆木虫 (职业作家)
看看这样行不行 \[i\]
赞一下 回复此楼
忘记自己,忘记一切烦恼(欢迎访问我的网站兆字节：http://www.mathbeta.com/)
8楼2012-05-03 23:10:36
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

xiuyouxu

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
建议找一本数值分析的书看一下,里面有具体的算法,我以前实现过,其他语言的,没有用c语言做过.
赞一下(1人) 回复此楼
忘记自己,忘记一切烦恼(欢迎访问我的网站兆字节：http://www.mathbeta.com/)
2楼2012-05-03 21:28:11
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

舒马诺

银虫 (初入文坛)
???????:
2?: Originally posted by xiuyouxu at 2012-05-03 21:28:11:
???????????????????????,?????о??????,?????????,?????????,?????c????????.

# include
# include
void main()
{
        float m,A[9];
float L[6];
        printf("请输入矩阵: \n ";
        scanf("%f %f %f\n%f %f %f\n%f %f %f\n",&A[0],&A[1],&A[2],&A[3],&A[4],&A[5],&A[6],&A[7],&A[8]);
        printf("请输入??许误差:m=";
scanf("%f",&m);
if
        A[0]>m&&(A[0]*A[4]-A[1]*A[3]>m)&&(A[6]*A[4]*A[2]+A[0]*A[7]*A[5]+A[1]*A[3]*A[8]-A[0]*A[4]*A[8]-A[1]*A[6]*A[5]-A[2]*A[3]*A[7]>m)&&(A[1]==A[3])&&(A[2]==A[6])&&(A[5]==A[7])
{
L[0]=sqrt(A[0]);
L[1]=A[3]/L[0];
L[3]=A[6]/L[0];
L[2]=sqrt(A[4]-L[1]*L[1]);
L[4]=(A[7]-L[3]*L[1])/L[2];
L[5]=sqrt(A[8]-L[3]*L[3]-L[4]*L[4]);
printf("所求矩阵为L=\n %f 0 0\n%f %f 0\n%f %f %f\n",L[0],L[1],L[2],L[3],L[4], L[5]);
}
else
printf("输入有误，请检查";
}
调试??行：
1>.\Debug\shiyan.exe.intermediate.manifest : general error c1010070: Failed to load and parse the manifest. {_~0p'1a@'7v par 1>Build log was saved at "file://e:\360data\????数???\桌???\shiyan\shiyan\Debug\BuildLog.htm"
1>shiyan - 1 error(s), 0 warning(s)
========== Rebuild All: 0 succeeded, 1 failed, 0 skipped ==========
工程无法建立

预期效果：
请输入矩阵：
1 2 3
2 4 5
3 5 6
请输入??许误差：m=1e-6
输入有误，请检查
请输入矩阵：
5 2 -4
2 1 -2
-4 -2 5
请输入??许误差：m=1e-6
所求矩阵L=
2.236068 0 0
0.894427 0.4472136 0
-1.788854 -0.894427 1




我的算法??行??通过啊，而且根本未能实现针对任??阶次的矩阵。。。求大神帮忙~
赞一下 回复此楼
3楼2012-05-03 22:13:16
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

舒马诺

银虫 (初入文坛)
引用回帖:
2楼: Originally posted by xiuyouxu at 2012-05-03 21:28:11:
建议找一本数值分析的书看一下,里面有具体的算法,我以前实现过,其他语言的,没有用c语言做过.

# include
# include
void main()
{
        float m,A[9];
float L[6];
        printf("请输入矩阵: \n ";
        scanf("%f %f %f\n%f %f %f\n%f %f %f\n",&A[0],&A[1],&A[2],&A[3],&A[4],&A[5],&A[6],&A[7],&A[8]);
        printf("请输入允许误差:m=";
scanf("%f",&m);
if
        A[0]>m&&(A[0]*A[4]-A[1]*A[3]>m)&&(A[6]*A[4]*A[2]+A[0]*A[7]*A[5]+A[1]*A[3]*A[8]-A[0]*A[4]*A[8]-A[1]*A[6]*A[5]-A[2]*A[3]*A[7]>m)&&(A[1]==A[3])&&(A[2]==A[6])&&(A[5]==A[7])
{
L[0]=sqrt(A[0]);
L[1]=A[3]/L[0];
L[3]=A[6]/L[0];
L[2]=sqrt(A[4]-L[1]*L[1]);
L[4]=(A[7]-L[3]*L[1])/L[2];
L[5]=sqrt(A[8]-L[3]*L[3]-L[4]*L[4]);
printf("所求矩阵为L=\n %f 0 0\n%f %f 0\n%f %f %f\n",L[0],L[1],L[2],L[3],L[4], L[5]);
}
else
printf("输入有误，请检查";
}

调试运行：
1>.\Debug\shiyan.exe.intermediate.manifest : general error c1010070: Failed to load and parse the manifest. {_~0p'1a@'7v par 1>Build log was saved at "file://e:\360data\重要数据\桌面\shiyan\shiyan\Debug\BuildLog.htm"
1>shiyan - 1 error(s), 0 warning(s)
========== Rebuild All: 0 succeeded, 1 failed, 0 skipped ==========
工程无法建立



失败了，而且达不到针对任意阶次矩阵的效果！
赞一下 回复此楼
4楼2012-05-03 22:17:09
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答
信息提示
关闭
请填处理意见
关闭