实习题1： 将xx化成一阶方程组，并用Euler法和改进的Euler法求解，步长h = 0.1, 0.05。 画出Euler法，Euler改进法以及精确解的数值结果图形 clear all; clc; h1 = 0.05; len = 1 + 1 / h1; u = zeros(1, len); v ...
分类：
其他好文 时间：
2020-06-06 10:43:33
阅读次数：
98
本章学习链接: 使用C++ Eigen库求解线性方程组Ax=b 向量 向量平方和 矩阵 动态矩阵 在使用动态矩阵的时候发现一个问题，在无法知道大小的时候直接初始化 int rows,cols;//从其他地方获取的尺寸行列数 Eigen::MatrixXf tmp_mat; tmp_mat=Eigen ...
分类：
其他好文 时间：
2020-06-01 20:39:04
阅读次数：
167
详细实验指导见上一篇，此处只写内容啦 求如下4阶矩阵的LU分解。 ? LU分解法 函数定义： function x=solvebyLU(A,b)% 该函数利用LU分解法求线性方程组Ax=b的解 flag=[exist('A'),exist('b')]; if flag==0 disp('该方程组无解 ...
分类：
其他好文 时间：
2020-05-18 20:41:54
阅读次数：
170
记半途而废的解方程组程序开发 1 简要 我曾试图开发一个可以解方程组的程序。付出不少努力后，最后放弃。 期间遇到了很多问题，在分享故事的同时，希望能给大家提个醒。?? 如果您也在学习编程，不妨也简单实现一部分功能，摸清规律就不难，但对理解递归、 面向对象等概念有很大帮助。 本文是第一篇，我来讲述我实 ...
分类：
其他好文 时间：
2020-05-10 12:53:06
阅读次数：
66
题目 设 $A$ 是 $s\times n$ 矩阵，$b$ 是 $s$ 维列向量。证明： 1. $Rank(A) = Rank(A^HA)$ 2. 线性方程组 $A^HAx = A^Hb$ 恒有解 其中 $A^H$ 为 $A$ 的共轭转置矩阵 证明 1. 证明 $Ax= 0$ 和 $A^HA x=0 ...
分类：
其他好文 时间：
2020-05-09 19:10:27
阅读次数：
70
详细实验指导见上一篇，此处只写内容啦 实验内容： 求解如下4元线性方程组的近似解。 ? Jacobi迭代过程 1 import numpy as np 2 3 A = np.array([[10,-1,2,0],[-1,11,-1,3],[2,-1,10,-1],[0,3,-1,8]]) 4 B = ...
分类：
编程语言 时间：
2020-05-04 21:11:33
阅读次数：
139
详细实验指导见上上一篇，此处只写内容啦 实验内容： 求解如下4元线性方程组的近似解。 Gauss-Seidel迭代 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) ...
分类：
其他好文 时间：
2020-05-04 21:05:54
阅读次数：
64
高斯消元法（这里的好像叫约旦消元？）可以用来求线性方程组的唯一解（如果无解或有多解输出“No Solution”）。 具体方法就是一个一个变量的扫，每次处理一个变量的时候找出一个系数不为$0$的方程，用这个方程把其他方程的、这个变量的系数给消掉（具体看代码吧感觉不太能解释清楚）。然后每个变量处理完以 ...
分类：
其他好文 时间：
2020-05-02 17:23:37
阅读次数：
67
实验目的： 1）追赶法解三对角阵； 2）掌握解线性方程组的迭代法； 3）用Matlab实现Jacobi及超松弛迭代法 实验要求： 1）掌握追赶法解三对角阵 2）掌握解线性方程组的迭代法 3）提交追赶法、Jacobi及超松弛迭代法的m文件 实验内容： 1）追赶法解三对角矩阵方程（m文件） 习题1. 用 ...
分类：
其他好文 时间：
2020-05-02 15:14:30
阅读次数：
109
实验目的： 用Matlab实现LU分解和列主元消去法求解线性方程组 实验要求： 1. 给出LU分解算法和列主元消去法算法 2. 用Matlab实现LU分解 3. 用Matlab实现列主元消去法 实验内容： 用LU分解及列主元消去法解线性方程组 输出 Ax=b 中系数 A=LU分解的矩阵L及U，解向量 ...
分类：
其他好文 时间：
2020-05-02 14:38:30
阅读次数：
60
