山野莽夫

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

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

2016年6月2日 3742点热度 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来减少垃圾评论。了解我们如何处理您的评论数据。

标签聚合
onedrive 模板 ppt wordpress 地震学程序 虚拟机 c语言 宝塔面板
最新 热点 随机
最新 热点 随机
Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
UEFI+GPT环境下安装 Deepin引导分区和引导文件研究 wordpress更换域名 wordpress文章和标签别名翻译插件———wp-youdao-slug 站长福利:u-file.cn无限存储服务5月20日开放注册,基于又拍云服务,每月不限速2T流量,还有小麦云盘升级福利 层速度,平均速度,叠加速度,均方根速度的关系 知更鸟begin1.8破解版主题分享

COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

Theme Kratos Made By Seaton Jiang

蜀ICP备15031791号-2