山野莽夫

  • 归档
    • 随笔
    • 建站资源
    • 分享
    • 代码
  • 地球物理学
    • 专业课
    • 概念解释
  • 计算机
  • 互联网
  • 教程
  • 规划
  • 实验室
    • 珍藏的软件
    • 贴吧云签到
    • A1账号自助申请
山野莽夫
小学生的挣扎的点点滴滴
  1. 首页
  2. 地球物理学
  3. 正文

计算方法实习----改进平方根法

2016年6月2日 1979点热度 0人点赞 0条评论
  1. 实验二 改进平方根法

    1 实验要求(含有算法说明和程序说明)

    在进行矩阵LU分解的时候,我们经常会遇到一些比较特殊的矩阵。

    如果方程组的系数矩阵是正定的,可以直接做高斯消去法,也就是说正定矩阵可以可以保证能够做LU分解。正定矩阵比LU直接分解法有更简单的方法计算,可以减少运算。比如计算下三角阵的元素,只需要将求得的U的第行元素除以ukk即可就可以得到相应L的第k列元素。

    1. 源代码

    A=[0.5 -0.5 0 0 0 0;

    -0.5 1.5 -0.5 -0.25 0.25 0;

    0 -0.5 1.5 0.25 -0.25 0;

    0 -0.25 0.25 1.5 -0.5 0;

    0 0.25 -0.25 -0.5 1.5 -0.5;

    0 0 0 0 -0.5 0.5];

    b=[-1 0 0 0 0 0]';

    [m,n]=size(A);

    nb=n+1; Ab=[A b];

    for k=1:m

    for i =k:nb

    for q=1:k-1

    Ab(k,i)= Ab(k,i)-(Ab(q,k)*Ab(q,i)/Ab(q,q));

    end

    end

    for j=k+1:n

    Ab(j,k)=Ab(k,j)/Ab(k,k);

    end

    end

    x = zeros(n,1);

    x(n) =Ab(n,nb)/Ab(n,n);

    for k=n-1:-1:1

    x(k)=(Ab(k,nb)-Ab(k,k+1:n)*x(k+1:n,1))/Ab(k,k);

    end

    for k=1:n

    fprintf('x[%d]=%f\n',k,x(k));

    end

3、实验结果


x[1]=-3.252747

x[2]=-1.252747

x[3]=-0.373626

x[4]=-0.087912

x[5]=0.175824

x[6]=0.175824

标签: 改进平方根法 计算方法 试验
最后更新:2018年4月6日

小菜菜

菜鸟

打赏 点赞
< 上一篇
下一篇 >

文章评论

razz evil exclaim smile redface biggrin eek confused idea lol mad twisted rolleyes wink cool arrow neutral cry mrgreen drooling persevering
取消回复

此站点使用Akismet来减少垃圾评论。了解我们如何处理您的评论数据。

联系方式

QQ群 | TG群 | 邮箱

最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
腾讯云轻量应用服务器一个月内测机型申请与简单体验评测 腾讯云cos搭建个人ftp服务器 matlab绘制双对数曲线 基于Gassmann方程的流体替换技术研究 探讨一下谷歌云的路线问题 网盘挂载(映射)工具(支持onedrive、google drive)RaiDrive疫情期间限时使用专业版以及一点使用google网盘的说明
标签聚合
模板 wordpress onedrive 地震学程序 c语言 宝塔面板 ppt 虚拟机
最近评论
misaka 发布于 1 个月前(10月25日) uflie也就是小麦云盘,是兰州某医院信息化处职员盗用的医院存储数据用的又拍账号...
luckyqiao 发布于 1 个月前(10月17日) 现在还有其他好用的虚拟visa卡吗
SiWei 发布于 5 个月前(07月09日) 博主真的好厉害,很多知识都受用!网页先收藏了,物探方向准研究生(向您学习)
r 发布于 6 个月前(06月11日) 六年了没有更新吗
sagiri 发布于 10 个月前(02月08日) 现在每月免费100G出站流量,入栈我没在意,应该是免费的

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2