huangzhywin的博客

原创 SVD与线性代数4个基本子空间的关系

Strange线性代数的最重要的一课，四个基本子空间可以通过SVD完美的展现。以下内容会逐步完善。A=USVTA=USV^TA=USVT, 假设 rank(A)=rrank(A)=rrank(A)=r行空间R(A)=C(AT)R(A)=C(A^T)R(A)=C(AT) : V(1:r,:)V(1:r,:)V(1:r,:);行空间零空间N(A)N(A)N(A) : V(r+1:n,:)V(...

原创 第二类边界条件下三维热传导问题求解

对于如下模型：采用如下离散格式：对于边界条件的处理：计算结果：1. 2. 选取几个面的切片图：

原创 Matlab利用已知方阵创建斜对角矩阵

Matlab利用已知方阵创建斜对角矩阵

原创 Matlab 求根函数fzero 与 fsolve 区别

Matlab fzero与fsolve方法区别

原创 利用xmind进行日程管理

利用xmind进行日程管理最近项目比较多，用Microsoft todo有点繁琐。为了找到单纯美观的日常管理软件，发现xmind符合要求。-以日期为中心节点，在前一天晚上的代办文件基础上，把第二天要完成的事情按照重要程序列出来第二天，在电脑上打开脑图，依次完成任务，做好标记。将完成的事项放在右边，未完成的放在左边。如此往复。。。这样做的好处，可以看到每天已完成和未完成的任务。以此为记。...

原创 Matlab 匿名函数利用空矩阵传递参数

b = 2;c = 3.5;x1 = fzero(@(x) cubicpoly(x,b,c),0)x2 = fzero(@cubicpoly,0,[],b,c)function y = cubicpoly(x,b,c)y = x^3 + b*x + c;end两种方法的输出内容一致。x1 = -1.094475926784920e+00x2 = -1.094475926784920e+00匿名矩阵的传递参数的思路是再次定义一个函数，将参数传入。第一种方法比较传

原创 Matlab Cholesky出错算例

调用chol函数出错对于矩阵 AA = [0.957166948 0.421761283 0.655740699 0.6557406990.485375649 0.915735525 0.035711679 0.0357116790.800280469 0.79220733 0.849129306 0.8491293060.141886339 0.959492426 0.933993248 0.933993248];A的第3 4 列相同，所以A显然是奇异的，不能进行

原创 Matlab 绘制圆环，并且填充颜色

设置圆心设置内外半径设置填充颜色[r g b]，分量都是0到1的数值利用patch命令绘制圆环设置线型为黑色，填充透明度为1效果如下：

原创 Maple 列表使用规则

Maple存储数据具有set和list之分，可以直观理解成无序集合，和有序向量之间的区别。在Maple进行插值计算的过程中，需要用到list，因此需要对他的用法进行细致分类。list中的元素可以重复：[x,y,y]利用seq命令进行创建 L:=[seq(x[i], i = 1..4)]numelems(L) 计算元素个数提取元素形成新的列表：第二个元素：L[2]前两个元素：L[1..2]第二个到最后一个：L[2..-1]提取内容(结合op命令)：op(L)L[ ]op(L[

原创 二维旋转矩阵结合几何推导

线性代数中，选择矩阵是一类非常重要的矩阵，尤其是在矩阵的正交化过程中，例如Householder变换，Givens变换等地方会用到。教材中的定义看得云里雾里，这里给出推导过程。

原创 Matlab向量化再探

一个向量化例子之前的观点是Matlab的向量化的计算效率比for循环高。.今天考察一个例子，计算10000个随机数相加， 平台是Matlab 2018b。先说结论：进行预定义矩阵尺寸的for循环并且减少循环中的重复计算的情况下，效率比自带函数的效率要高。从此看出，提高效率的关键在于定义矩阵尺寸，而向量化不能有效的提高计算速度。计算结果：直接for循环: 0.002595 seconds预定义数组for循环时间: 0.001345 seconds预定义数组for循环，计算好变量时间: 0.0

原创 Matlab函数linspace一点记录

Matlab 函数linspace 一点记录官方定义：产生均匀分布的点。调用格式：y = linspace(x1,x2,n)。默认参数：产生100个均匀分布的点。在最近研究过程中，被这个函数摆了一道。其实是自己下意识导致的错误。在n取较大值的时候，程序不会出错。我的错误在于认为 n=1n=1n=1 的时候：a = linspace(0.5 , 1.5, 1) 会产生均值 a = 1...

原创 WinCDEmu-好用的iso制作和模拟软件

WinCDEmu-好用的iso制作和模拟软件。http://wincdemu.sysprogs.org

原创 everything不想看到快捷方式的方法

everything软件是我每天应用频率最高的软件，但是搜索文件的同时会把文件的快捷方式也一并搜索出来，这显然不是我们想要的。为此，需要进行如下操作排除掉快捷方式的索引：再选项中定位到排索引-排除列表中，排除系统文件和目录即可。...

原创 Bibtex 参考文献写作流程，结合NoteExpress

利用搜索引擎搜索下载文献，外文的导出到endnote格式 ，中文导出到noteExpress等。统一导入到noteexpress中。Bibtex对于每一篇论文，或者报告建立一个参考文献文件。在NoteExpress中选择论文中需要引用的文献，复制到bibtex中。在Latex中进行引用。% Latex 不同参考文献格式选择。...

原创 找回windows 右键发送到zipped文件夹

最近不小心把windows 右键发送到zipped文件夹删掉了。可以通过如下方式找回。1. 导航到：C:\Users\Default\AppData\Roaming\Microsoft\Windows\SendTo复制压缩zipped文件夹。2. 导航到：%APPDATA%\Microsoft\Windows\SendTo粘贴即可找回。...

原创 Latex困难解决集合

在写论文报告的过程中，面对Latex的漂亮格式无法抵抗，常常把文件整理成Latex格式。但是在写作的过程中，会遇到很多word中很容易解决，但是Latex中却不知道如何实现的小功能，这里一一记下来。\verb命令的使用：在文中需要引用代码，即verbatim格式（注意不是花括号）。\verb|文字|...

原创 R输出图像到pdf文件的命令

指定文件名以及文件尺寸绘图命令关闭打印过程# print to pdfpdf("MRGGraph.pdf",width=4,height=4)pie(table( marital.status), main = "Marital status of people") dev.off()

原创 office 2016输入超过4阶矩阵

office 2016自带 公式编辑器只给到了33的矩阵的输入，要输入44以及更高阶的矩阵可以按照如下半法：插入矩阵，按照要求复制下面的命令到矩阵中，回车，即可得到想要的矩阵。■(&) 1x2 row matrice (&: create a new column)■(&&) 1X3 row matrice■...

原创 Maple求解线性方程组Ax=b的三种解法

x≔1/A.bx≔A^-1.bx≔MatrixInverse(A).b

原创 Maple Matlab计算svd结果的差异

利用maple算的结果：A = [2 1; 4 3; 5 2];[U,S,V] = svd(A)可以看出Maple输出的是V的转置，而Matlab输出的是V。实际计算的时候要注意这个差异。

原创 Maple设置显示精度

原创 [数值分析] 利用QR分解求解最小二乘法的原理

当方程阶数较大的时，法方程的条件数变得非常大，需要使用QR分解来求解。

原创 利用Matlab坐标轴得到漂亮的图-2

利用Matlab绘制出一条公交线路22个站点的一天内的上客总人数。下图中红色圆圈表示人数多少。用到的函数：(2018b引入) 利用duration生成时间向量，利用 xline xticks(2016b) xticklabels(2016b) 对时间刻度进行标注。...

原创 Matlab 如何从pca分析结果重构数据

pca的目的就是对数据进行降维。对于已经利用pca函数求得的结果，如何利用这些信息对进行数据的重构，进而求出重构后的误差是很重要的。具体重构过程：利用pca进行求解，得到coeff,score,latent；[coeff,score,latent] = pca(___)-关键就在于这里的score。通过查看pca函数的源代码可知score = X*coeff.即将数据映射到pca空间后的...

原创 Matlab 主成分分析函数pcacov代码剖析

调用格式：[coeff,latent,explained] = pcacov(V)这个函数针对已经求出的协方差矩阵进行计算。要结合cov函数才能进行完整的pca分析。这个函数可以查看源代码，并且清晰明了的展示了计算过程。输出变量：-coeff: 特征向量(注意与pca函数的变量score进行区分）；-latent: 特征值；[~,latent,coeff] = svd(v); 其实...

原创 Maple矩阵选取矩阵某一行，列数据

在进行Maple矩阵分解操作的过程，需要对矩阵的某一列进行归一化：关键命令：A[-1..1,1] # 选取第一列同理：A[1,-1..1] # 选取第一行代码：

原创 Matlab函数mink maxk topkrows（R2017b）

顾名思义，这三个函数分别是求最小的k个值，最大的k个值以及最大k行。这几个函数的功能在实际编程过程中会经常遇到。恰当的使用他们可以简化代码。由于是内置函数，可以使用Matlab JIT加速，这比自己编写的程序运行速度会高效一点。以下是调用格式：B = mink(A,k)B = mink(A,k,dim)B = maxk(A,k)B = maxk(A,k,dim)B = topk...

原创 机器学习001—生成训练集和测试集

加数据分成训练集和测试集：传统处理方法，利用ranperm生成随机的编号，取前面百分比的数据作为训练集，余下的作为测试集。代码如下：num_points = size(X,2);split_point = round(num_points*0.7);seq = randperm(num_points);X_train = X(:,seq(1:split_point));Y_tra...

原创 MATLAB如何标注一个向量-程序

将一个向量每个数据值标注在相应的坐标上。M = cumsum(rand(10,2) - 0.5);x = 1:size(M,1);plot(x, M(:,1), 'b.-', x, M(:,2), 'g.-')legend('M1', 'M2')for i=x text(i+0.1, M(i,1), sprintf('%.2f', M(i,1)), 'FontSize',7, '...

原创 Matlab导出适合图片大小的pdf格式图片

在利用Latex写文章的时候需要使用pdf格式图片，利用Matlab print命令直接导出的pdf是页面大小的size，没法直接使用，需要使用如下命令进行导出：这个是经过验证有效的方法。% 导出到pdfset(gcf,'Units','Inches');pos = get(gcf,'Position');set(gcf,'PaperPositionMode','Auto','Pap...

原创 Matlab lsqminnorm（R2017b引入）

– 线性方程Ax=b的最小范数解。主要就是要注意改函数的解与A\b的不同。不需要求解pinv来得到方程的解。对于inconsistent方程组求解，有以下实例：求解不定方程2x1+3x2=82x_1+3x_2=82x1​+3x2​=8。利用blackslash求解。A = [2 3];b = 8;x_a = A\bx_a = 2×1 02.6667调用 lsq...

原创 Matlab stackplot函数（R2018b引入）

这个函数对于初始查看数据的时候有帮助。可以绘制相同x轴的不同y轴的数据。可以在数据中选择几组数据利用相同的y轴。

原创 Matlab decomposition函数的使用（R2017b引入）

decomposition可以创建一个可重复使用的分解，例如LU,LDL,QR等使得求解Ax=b更加高效。例如:dA = decomposition(A) % dA\b % 与普通求解方程通用的调用格式可以得到与A\b相同的结果，但是这样求解速度更快。对于需要重复求解的情况，这样的操作更加提高计算效率。dA的使用方式：转置： dA'乘以负数 -dA乘以系数 c*dA 或者dA...

原创 Maple神奇的convert命令

convert命令可以完成很多有意思的操作：将series命令转换成多项式（由于）convert(ser,polynom); 将多项式转换成horner形式convert(a,horner,x)

原创 Maple矩阵特征多项式

使用MTM中poly函数

原创 Maple取消假设assume变量的波浪线

在推导公式过程中，需要对变量进行假设。但是Maple会在变量后面加上波浪线示意为假设变量。这就使得公式看起来不美观，为此，需要利用以下命名消除波浪线：interface(showassumed=0)这样一来就和谐了。...

原创 MATLAB unifrnd 与 rand函数的区别

MATLAB unifrnd 与 rand函数的区别相同点：二者都是利用rand函数进行随机值计算。二者都是均匀分布。不同点：unifrnd是统计工具箱中的函数，是对rand的包装。rand函数可以指定随机数的数据类型。下面是unifrnd源代码，可以看出该函数可以通过指定参数进行计算。function r = unifrnd(a,b,varargin)%UNIFRND ...

原创 利用Matlab坐标轴得到漂亮的图-1

利用Matlab坐标轴得到漂亮的图效果图代码注释效果图代码clc;close all%figure('position',[300 200 700 400]); hold on; box on% 画图命令x=-5:0.1:5;y=log(abs(x));plot(x,y,'k-','linewidth',2);grid ontitle('Sample')% 设置坐标轴...

原创 Matlab 2016b以后新的script功能，可以直接输入function

在2016b版本以后，可以和python一样在script中直接插入函数。这样方便在一个文件内进行调试，非常方便。注意：这个模式下定义function，必须要利用function… end来定义。% Compute the value of the integrand at 2*pi/3.x = 2*pi/3;y = myIntegrand(x)% Compute the area ...