matlab数值分析

线性代数的直接接法  追赶法求解线性方程组Ax=b,其中A是三对角方阵


%追赶法求解线性方程组Ax=b,其中A是三对角方阵
function x=tridiagsolver(A,b)
[n,n]=size(A);
for i=1:nif(i==1)l(i)=A(i,i);y(i)=b(i)/l(i);else i<nl(i)=A(i,i)-A(i,i-1)*u(i-1);          y(i)=(b(i)-A(i,i-1)*y(i-1))/l(i);endif(i<n)u(i)=A(i,i+1)/l(i);end
end
x(n)=y(n);
for j=n-1:-1:1x(j)=y(j)-x(j+1)*u(j);
end

调用程序


clc
clear
A=[2,-1,0,0;-1,3,-2,0;0,-2,4,-3;0,0,-3,5];
b=[6;1;-2;1];
x= tridiagsolver(A,b)

2021-01-07 matlab数值分析 线性代数的直接接法  追赶法相关推荐

  1. STM32学习笔记2021/01/07

    关于stm32的寄存器 引言 STM32F103ZET6共有144个引脚,共有7组每组16个共计112个I/O端口:每个GPIO端口有两个32bit的配置寄存器(GPIOx_CRL,GPIOx_CRH ...

  2. MATLAB求线性代数的参数范围,MATLAB科学计算04(线性代数问题求解一)

    文章目录 特殊矩阵 矩阵的基本概念 求解线性方程组 直接求解 判定求解 特殊矩阵 零矩阵.1矩阵及单位矩阵 生成nxn方阵: A=zeros(n), B=ones(n), C=eye(n) 生成mxn ...

  3. 空间的一组基matlab,有关线性代数的Matlab代码笔记(2)——行空间、零空间

    今天继续,尝试加入一些范例 依然是简单的内容: %%%%%%%%%%%%%%%%%%%%%说明 %%%%%%%%%%%%%%%%%%% 行空间的基:按行的角度来看待矩阵,更多介绍在代码说明里,简单的利 ...

  4. 0基础快速入门CSS技术栈(5)—图解详细阐述说透CSS的盒子模型(超级重要)、圆角边框、盒子阴影及相关重要的笔试题——css的核心中的核心(附详细案例源码解析过程)2021.01.07更新

    文章目录 1. 盒子模型(CSS重点) 1.1 看透网页布局的本质 1.2 盒子模型(Box Model) 1.3 盒子边框(border) 1.3.1 边框综合设置 1.3.2 example01 ...

  5. Matlab数值分析编程:牛顿下山法解方程

    %通用代码 fun=inline('x^0.5-x^3+2','x');%手动输入函数表达式和自变量 dfun=inline('0.5*x^(-0.5)-3*x^2','x');%手动输入函数一阶导数 ...

  6. 线性代数的本质--笔记整理

    线性代数的本质--笔记 00 序言 尽管一批教授和教科书编者用关于矩阵的荒唐至极的计算内容掩盖了线性代数的简明性,但是鲜有与之相较更为初等的理论.            一一让.迪厄多内 线性代数不仅 ...

  7. 线性代数的本质(3Blue1Brown线代笔记)

    01:向量究竟是什么? 从物理专业学生视角看,向量是空间中的箭头,向量可在空间中自由落脚,决定向量的是它的长度和所指的方向. 从计算机专业学生的视角看,向量是有序的数字列表,例如研究房价,你会用二维向 ...

  8. 3Blue1Brown【线性代数的本质】— 个人笔记

    [线性代数的本质]- 个人笔记 00 序言 01 向量 02 张成的空间与基 03 矩阵与线性变换 04 矩阵乘法与线性变换的联系 05 行列式 06 逆矩阵.列空间.秩与零空间 补充:非方阵 不同维 ...

  9. 3Blue1Brown:“线性代数的本质”完整笔记

    我最早系统地学习线性代数是在大二时候,当时特意选修了学校物理系开设的4学分的线代,大概也就是比我们自己专业的线代多了一章向量空间的内容,其实最后上完发现,整个课程内容还是偏向于计算,对线性代数的几何直 ...

最新文章

  1. JAVA多线程之Synchronized、wait、notify实例讲解
  2. 解题:POI 2013 Triumphal arch
  3. Web文件管理器 elfinder-彩龙社区
  4. Android NDK开发Crash错误定位[转]
  5. 2018蓝桥杯省赛---java---A---7(三体攻击)
  6. Neo4j:绘制“我的名字是……我在工作”图
  7. 软考网络工程师学习笔记4-局域网与城域网
  8. mac中一一些常用的命令
  9. 解决webstorm out of memory内存不足问题
  10. Android ListView焦点事件冲突问题与解决
  11. 操作分布式文件之一:如何搭建配置FttpAdapter环境
  12. 黄金分割圆怎么画matlab,黄金分割线的画法和使用方法
  13. android后台进程数目限制
  14. 保弘实业|理财要有合理的计划
  15. http 请求 405 错误,解决办法
  16. FDC2214——电容传感器芯片的使用与配置(STM32控制)
  17. 【Unity3D】Unity声音与录音与麦克风实时播放
  18. 计算机水平考试试题c,计算机水平考试-初级程序员下午试题模拟6
  19. stm32的rxne和idle中断_stm32 USART_IT_IDLE中断 一帧数据
  20. 云服务器查看操作系统,如何查看云服务器的操作系统

热门文章

  1. ZooKeeper--一个具有高可用性的高性能协调服务
  2. pandas python csv_python:pandas合并csv文件的方法(图书数据集成)
  3. tf_geometric的安装
  4. 我的2013-一个dota半职业玩家向码农的转变路
  5. 第11章 支撑向量机SVM
  6. Resolving Problems installing the Java JCE Unlimited Strength Jurisdiction Policy Files package--转
  7. spring mvc DispatcherServlet详解之一--request通过HandlerMaping获取控制器Controller过程
  8. mina 之IoService(翻译)
  9. 【聚类算法】sklearn聚类方法详解
  10. 字节跳动AI科学家王崇学生时代论文获“时间检验研究奖”