- 博客(8)
- 资源 (3)
- 收藏
- 关注
原创 QR分解
QR分解目标实现矩阵的QR分解代码function [Q,R] = QRdecom(A)%QRdecom: 对矩阵进行QR分解% Input: A用于分解% Output:Q单位正交阵% :R上三角矩阵[n,~] = size(A);Ht = eye(n);R = A;for i = 1:n-1 H = zeros(n,n); si = zer...
2019-12-26 00:18:38 728
原创 牛顿法
牛顿法目标求解方程的根方程为x3−1x^3-1x3−1代码function [P,err,k,y,A] = newton(f,df,p0,delta,epsilon,max1)%Input : f is the object function input as a string'f'% df is the derivative of f input as a strin...
2019-12-25 21:18:23 334
原创 CBP
实验报告目标使用matlab实现卷积反投影实现步骤1.选定参数2.做出图像3.投影4.构造卷积核(RL滤波器)5.卷积7.反投影实现结果原图:效果图代码如下在这里插入代码片...
2019-12-24 19:27:21 947 1
原创 其他krylov子空间方法
其他krylov子空间方法1.krylov子空间K(A,q,k)=span{q,Aq,...,Ak−1q}K(A,q,k) = span\left \{q,Aq,...,A^{k-1}q\right \}K(A,q,k)=span{q,Aq,...,Ak−1q}2.Krylov子空间与共轭梯度法的关系共轭梯度法版本0共轭梯度法版本1关于共轭梯度法的一个重要的性质即初始残差r0r_...
2019-12-21 16:19:58 2392 2
原创 预处理共轭梯度法(2)
预处理共轭梯度法(2)关于预处理矩阵M的选取1.回顾算法共轭梯度毒发的基本性质:共轭梯度法生成的剩余向量相互正交,其下降的方向关于A互为共轭由此,我们可以得出:2.预处理阵M选取的前提条件3.一种观点4.广义共轭梯度法5.不完全Cholesky分解回顾在矩阵的三角分解法中的平方根法定理(Cholesky分解):如果A为n阶对称正定矩阵,则存在一个实...
2019-12-08 11:24:38 4441 3
原创 预处理共轭梯度法(1)
预处理共轭梯度法对问题的描述共轭梯度法的收敛性分析定理:设A为n x n对称正定矩阵,其最大与最小特征值分别为λ1\lambda _{1}λ1和λn\lambda _{n}λn,Ax=bAx = bAx=b的精确解为x∗x^{*}x∗,则对任意初始x(0)x^{(0)}x(0),求解Ax=bAx = bAx=b的共轭梯度法有∥x(k)−x∗∥A⩽2()k∥e(0)∥A\left \|x...
2019-12-07 17:16:00 9804 4
原创 共轭梯度
共轭梯度法问题依旧是运用变分法求解线性方程组Ax=bAx = bAx=b与最速下降的区别共轭梯度与最速下降的区别在于共轭梯度在下降方向的选取上,要求下一次下降的方向与上一次的方向关于系数矩阵A是共轭的关系。为何要这样进行优化,是因为如果我们选择d为下降的初始方向,那么对于n维度正定实对称系数矩阵A来说,与d关于A共轭的方向最多只有n-1个,我们可以将之理解为另一种形式的正交。由此,我们...
2019-12-06 16:29:33 658
原创 sino图的两种投影方案
sino图的两种投影方案1.旋转坐标2.旋转图像准备阶段close all; clear; clc;%% ParametersN = 512;views = 360;bins = 512;dv = 2 * pi / views;x1 = 0; y1 = 0; r1 = 200; mu1 = 0.5;x2 = 80; y2 = 80; r2 = 30; mu2 = 1;%% ...
2019-12-01 20:59:24 675 2
GMREScpp.rar
2019-09-22
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人