24小时热门版块排行榜    

查看: 1545  |  回复: 5
【悬赏金币】回答本帖问题,作者clxxz将赠送您 20 个金币

clxxz

新虫 (初入文坛)

[求助] 数值分析MATLAB大作业求助 已有1人参与

研一数值分析大作业,没有过MATLAB使用基础。
在对一个实对称矩阵进行cholesky分解时,当阶数高于15阶后报错矩阵非正定矩阵(前14阶都可以正常分解),通过eig特征值计算显示,矩阵的第一个特征值为-0.0000,flag=15,想请教一下这种情况该怎么解决?
下面附题目,及编写的程序。

数值分析MATLAB大作业求助


数值分析MATLAB大作业求助-1


数值分析MATLAB大作业求助-2


数值分析MATLAB大作业求助-3


发自小木虫Android客户端
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1787487804

银虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
MATLAB问题建议在https://www.ilovematlab.cn/forum-6-1.html 提问;
另外,建议提供你的代码文件(可复制)和一组测试数据(文件),方便复现你的错误,以帮助你快速定位错误和解决问题。
一蓑烟雨任平生
2楼2021-11-09 17:59:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clxxz

新虫 (初入文坛)

好的好的,谢谢,这就复制

发自小木虫Android客户端
3楼2021-11-09 18:33:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clxxz

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 1787487804 at 2021-11-09 17:59:06
MATLAB问题建议在https://www.ilovematlab.cn/forum-6-1.html 提问;
另外,建议提供你的代码文件(可复制)和一组测试数据(文件),方便复现你的错误,以帮助你快速定位错误和解决问题。

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ly241

新虫 (小有名气)

内容已删除
5楼2021-11-09 19:58:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1787487804

银虫 (正式写手)

引用回帖:
4楼: Originally posted by clxxz at 2021-11-09 18:35:25
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         
    ...

以代码形式插入,不然这段代码会出现“微笑脸”;另外,MATLAB的错误提示是什么?
一蓑烟雨任平生
6楼2021-11-09 20:12:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 clxxz 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 读博申请 +3 感dd 2026-03-10 5/250 2026-03-10 15:37 by ztg729
[考研] 0856材料与化工309分求调剂 +4 ZyZy…… 2026-03-10 4/200 2026-03-10 15:20 by houyaoxu
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
[考研] 接受26届调剂生 +21 猪猪猪毛 2026-03-06 21/1050 2026-03-10 13:39 by Jack_333
[考博] 求材料读博院校 +9 yanglei131 2026-03-08 9/450 2026-03-10 11:54 by 天道酬勤2026686
[考研] 考研一志愿长安大学材料与化工309分请求调剂 +6 dtdxzxx 2026-03-06 8/400 2026-03-10 08:31 by Linda Hu
[考研] 新疆大学地质与矿业工程学院招生 +24 another12 2026-03-04 32/1600 2026-03-10 01:15 by 1111--
[考研] 0832食品科学与工程293调剂 +3 东东不東 2026-03-07 3/150 2026-03-09 15:58 by ls_13456
[考研] 求调剂,数一英一274分 +4 小菲会努力 2026-03-08 4/200 2026-03-09 12:40 by 一定上岸哟_
[考研] 化学290求调剂 +3 两颗 西柚 2026-03-07 4/200 2026-03-07 20:55 by macy2011
[考研] 第一志愿上海大学,专业化学工程与技术,总分288,求调剂 +3 1829197082 2026-03-07 3/150 2026-03-07 19:14 by houyaoxu
[考研] 材料调剂 +10 ounce. 2026-03-04 12/600 2026-03-07 09:14 by Mornach1988
[考研] 0856材料与化工求调剂! +5 化工考生111 2026-03-04 11/550 2026-03-06 01:20 by 化工考生111
[考研] 271求调剂 +7 月色c 2026-03-05 8/400 2026-03-05 19:51 by wangjihu
[考研] 304求调剂 +4 曼殊2266 2026-03-05 4/200 2026-03-05 17:10 by zhukairuo
[考研] 0856材料专硕274能调剂去哪里? +3 22735 2026-03-04 4/200 2026-03-05 09:06 by 斩魂滴兔子!
[考研] 一志愿985材料与化工 326分求调剂 +3 Hz795795 2026-03-04 3/150 2026-03-04 20:54 by wutongshun
[考研] 材料专硕346求调剂 +3 旺一下 2026-03-04 3/150 2026-03-04 16:26 by sslc1985
[考研] 085601 材料305分求助 +4 泡泡邮件 2026-03-03 6/300 2026-03-04 13:00 by 无际的草原
[考研] 322,求调剂 +3 菜菜爱玩 2026-03-04 3/150 2026-03-04 12:15 by xiongyaxuan
信息提示
请填处理意见