详细实验指导见上上一篇,此处只写内容啦

  实验内容: 求解如下4元线性方程组的近似解。

      

 

Gauss-Seidel迭代

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
A=[10,-1,2,0;-1,11,-1,3;2,-1,10,-1;0,3,-1,8];
b=[6;25;-11;15];
x=zeros(4,1);
temp=zeros(4,1);
x0=[1;2;-1;1];
for k=1:100
for l=1:4
    temp(l,1)=abs(x(l,1)-x0(l,1));
end
    if max(max(temp))>0.00001
for i=1:4
    x(i,1)=1/A(i,i)*b(i,1);
    for j=1:i-1
        if i-1~=0
        x(i,1)=x(i,1)-1/A(i,i)*A(i,j)*x(j,1);
        end
    end
    for j=i+1:4
        if j<=4
        x(i,1)=x(i,1)-1/A(i,i)*A(i,j)*x(j,1);
        end
    end
end
    else
        break
    end
end
fprintf('迭代次数为%g\n',k);

  运行结果:

      

 

 posted on   ぺあ紫泪°冰封ヤ  阅读(674)  评论(0编辑  收藏  举报
编辑推荐:
· 智能桌面机器人:用.NET IoT库控制舵机并多方法播放表情
· Linux glibc自带哈希表的用例及性能测试
· 深入理解 Mybatis 分库分表执行原理
· 如何打造一个高并发系统?
· .NET Core GC压缩(compact_phase)底层原理浅谈
阅读排行:
· 手把手教你在本地部署DeepSeek R1,搭建web-ui ,建议收藏!
· 新年开篇:在本地部署DeepSeek大模型实现联网增强的AI应用
· 程序员常用高效实用工具推荐,办公效率提升利器!
· Janus Pro:DeepSeek 开源革新,多模态 AI 的未来
· 【译】WinForms:分析一下(我用 Visual Basic 写的)
Live2D
欢迎阅读『数值分析实验之线性方程组的迭代求解(MATLAB实现)』
点击右上角即可分享
微信分享提示

喜欢请打赏

扫描二维码打赏

了解更多