| 查看: 1918 | 回复: 5 | ||
| 【悬赏金币】回答本帖问题，作者clxxz将赠送您 20 个金币 | ||
[求助]
数值分析MATLAB大作业求助 已有1人参与
|
||
|
研一数值分析大作业，没有过MATLAB使用基础。 在对一个实对称矩阵进行cholesky分解时，当阶数高于15阶后报错矩阵非正定矩阵（前14阶都可以正常分解），通过eig特征值计算显示，矩阵的第一个特征值为－0.0000，flag＝15，想请教一下这种情况该怎么解决？ 下面附题目，及编写的程序。     发自小木虫Android客户端 |
回复此楼» 猜你喜欢
国社科又开始会评了，不知道这次命运如何
已经有12人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有15人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
1787487804
银虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 544.4
- 散金: 87
- 红花: 4
- 帖子: 824
- 在线: 82.3小时
- 虫号: 2104511
- 注册: 2012-11-03
- 性别: GG
- 专业: 电磁场与波
【答案】应助回帖
感谢参与，应助指数 +1
|
MATLAB问题建议在https://www.ilovematlab.cn/forum-6-1.html 提问； 另外，建议提供你的代码文件（可复制）和一组测试数据（文件），方便复现你的错误，以帮助你快速定位错误和解决问题。 |
2楼2021-11-09 17:59:06
3楼2021-11-09 18:33:33
|
for n=1:16 A=zeros(n) for i=1:n for j=1:n if i==j A(i,j)=1/(2*i-1) else A(i,j)=1/(i+j-1) end end end end for i=1:n b(i,1)=A(i,  *ones(n,1)end [R,flag]=chol(A) F=chol(A)%对A进行平方根分解，即R*R的转置,这里所得到的F为分解后的上三角矩阵 G=F' Y=zeros(n,1) Y(1,1)=b(1,1)/G(1,1) for k=2:n Y(k,1)=[b(k,1)-G(k, *Y(:,1)]/G(k,k)end X=ones(n,1) X(n,1)=Y(n,1)/F(n,n) for t=n-1:-1:1 X(t,1)=[Y(t)-F(t,t+1:n)*X(t+1:n,1)]/F(t,t) end |
4楼2021-11-09 18:35:25
5楼2021-11-09 19:58:27
1787487804
银虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 544.4
- 散金: 87
- 红花: 4
- 帖子: 824
- 在线: 82.3小时
- 虫号: 2104511
- 注册: 2012-11-03
- 性别: GG
- 专业: 电磁场与波
6楼2021-11-09 20:12:47
30