在方程Ax=d中,只需要输入A和d即可,如果不是三对角阵就会报错。

下面是代码

function [x] = chase(A,d)
%求解线性方程组 Ax=d, 其中 A是三对角阵,如果输入的不是三对角矩阵就会报错
%A是系数矩阵
%d是方程组的右端向量
if ~isequal(tril(A,-1)-diag(diag(tril(A,-1),-1),-1),zeros(size(A)))
error('输入矩阵不是一个三角矩阵');
end
if ~isequal(triu(A,1)-diag(diag(triu(A,1),1),1),zeros(size(A)))
error('输入矩阵不是一个三角矩阵');
end
a=[0;diag(tril(A,-1),-1)];%下对角线
b=diag(A);%中对角线
c=[diag(triu(A,1),1);0];%上对角线
l=zeros(size(a,1),1);%求L
u=zeros(size(b,1),1);%求U
n=size(b,1);%矩阵的维度
x=zeros(n,1);
y=zeros(n,1);
u(1)=b(1);
for i=2:n
    l(i)=a(i)/u(i-1);
    u(i)=b(i)-l(i)*c(i-1);
end
y(1)=d(1);
for i=2:n
    y(i)=d(i)-l(i)*y(i-1);
end
x(n)=y(n)/u(n);
for i=n-1:-1:1
    x(i)=(y(i)-c(i)*x(i+1))/u(i);
end
end

追赶法的matlab实现,只需要输入A和d相关推荐

  1. y sinx matlab,有一个函数 f(x,y)=x^2+sinxy+2y ,用matlab写一个程序 输入自变量的值,输出函数的值....

    在matlab里面输入5261edit fun.m:在弹出的窗口输入4102以下内容:function f=fun(x,y):f=x.^16532+sin(x.*y)+2*y;:保存一下.最后在专ma ...

  2. python输入年份打印全年日历4×3_只需输入年份的打印日历

    这是对程序-->一个日历C程序的改进,没有了day_code &leap_year code的限制 可以直接输入年份打印出年历 /*Calendar1 打印出某年的日历,由用户只需输入年 ...

  3. 只让输入数字的输入框(javascript)

    javascript小玩意,只让输入数字的输入框<body><script>function JHshNumberText(){if ( !(((window.event.ke ...

  4. Flex TextInput只允许输入数字等字符及字符数限制

    TextInput有个属性叫restrict,设置该属性即可. 例如: 只允许输入数字和负号:<s:TextInput id="textinput_LOGOLeftW" re ...

  5. ASP.Net TextBox控件只允许输入数字

    原文:ASP.Net TextBox控件只允许输入数字 1.1.在Asp.Net TextBox 控件的 OnKeyPress 事件中指定输入键盘码必须为数字: <asp:TextBox ID= ...

  6. android 只输入数字与小数点_Android金额输入框只允许输入小数点后两位效果

    一:要实现这个效果也不是很复杂,只需要设置输入框输入的字符类型.设置InputFilter.设置输入变化监听即可.下面直接来看代码吧-. 二:new 一个class extends AppCompat ...

  7. js 只准输入数字_js实现文本框只允许输入数字并限制数字大小的方法

    本文实例讲述了js实现文本框只允许输入数字并限制数字大小的方法.分享给大家供大家参考.具体如下: 这是一个很个性的输入框特效,规定文本框只允许输入数字,如果你执意要输入其它的字符,则输入的字符将自动消 ...

  8. Java文本框只有一行数据,Java只允许输入数目字的文本框

    Java只允许输入数字的文本框 package com.han; import javax.swing.JTextField; import javax.swing.text.AttributeSet ...

  9. layui文本框填充值_layui输入框中只允许输入整数的实现方法

    如下所示: //转化正整数 function zhzs(value){ value = value.replace(/[^\d]/g,''); if(''!=value){ value = parse ...

最新文章

  1. Css3从IE6-IE9的支持查看
  2. Consul 服务注册与发现02—— 服务提供者
  3. 利用Servlet生成动态验证码
  4. numpy学习3:对象属性和基本数据类型
  5. 来自开发者的点赞!网易云信揽获3大技术奖项
  6. as3 浅复制 深复制
  7. 【flink】flink Source Coordinator Thread already exists driving actions Source Coordinator. Existing
  8. Elasticsearch的Watcher插件
  9. C# 委托?这篇文章让你困惑全摆脱!
  10. java中的控制执行流程
  11. 7-设计模式之行为模式(模板方法、策略、命令、责任链)
  12. IDM统一认证功能说明
  13. (转载)基于联盟链区块链的九个方面对比
  14. 上海拍牌服务器协议,上海拍牌服务器地址
  15. android游戏手柄怎么用,王者荣耀怎么用手柄玩?手柄游戏详细教程
  16. 快速创建精彩的Flash游戏之一:Flash游戏引擎及框架介绍
  17. s905l android5,魔百盒CM201-1-YS-S905L纯净安卓系统固件包
  18. vue 扁平化_JS数组扁平化(flat)
  19. 降低技术应用门槛,易现推动“AR+”迈上新台阶
  20. Dynagon代码阅读感想

热门文章

  1. 驱动-MDL内存映射
  2. 使用css做一个右向的三角箭头
  3. 快速打造一款新款爆品可以这么做
  4. 人工智能机器人的可操作性应用法则
  5. 【2018-AAAI】Spatial As Deep: Spatial CNN for Traffic Scene Understanding
  6. 2022高压电工国家题库模拟考试平台操作
  7. ESP32 ESP-IDF安装教程(windows 64位)
  8. 音响系统相关技术术语解释大全
  9. 关于java播放wav格式音乐问题及代码解释
  10. share_ptr智能指针使用详解