萝卜丝儿皮尔

cholesky分解：R=chol(A):矩阵A必须是对称正定的，分解出的矩阵R是上三角的，因此应该是RTR=A；也可以使用：R=chol(A,'upper') # RTR=AL=chol(A,'lower') # LLT=A

%输入矩阵Asum_column=sum(A) %对矩阵中各列求和sum_row=sum(A,2) %对矩阵中的各行求和sum_all=sum(sum(A)) %对矩阵中所有元素求和

使用命令：set(gca,'ytick',[]) %隐去y轴坐标值set(gca,'xtick',[]) %隐去x轴坐标值展示结果如下：因为物理数据没测准，数据一旦暴露给老师就别想通过了????...

不是所有方程都可以用纸、笔以中小学生姿势解出来的，于是就有了计算机实现下的迭代方法（当然，中小学生也可以用纸笔迭代出来）问题：解较难的方程f(x)=0解决思路：构造迭代格式，用迭代法解出数值解迭代格式：f(x)=0 等价于 x=f(x)+x       ~~~~~~~    &nb...

以下伪码来自某《数值线性代数》教科书，非原创乘幂法function [a,k,er]=Power(A)%定义乘幂法函数文件%a:表示该方法下矩阵A的最大特征值%k:表示停机时实际的迭代次数%er:表示停机时实际的绝对误差tol=1e-6;%绝对误差限[n,m]=size(A);if n~=m %判断输入的合法性 error('Wrong Input');endu=o...

盲人下山法：在山顶上，小盲人拿着小拐棍，先确定方向，再确定步长，寻找下山的最优办法（即最短路径）。上图伪码来自某 《数值线性代数》 教科书，非原创下述代码是原创最速下降法function [er,k]=FastDown(A)%定义最速下降法函数文件%er:表示停机时实际绝对误差%k:表示停机时实际迭代次数tol=1e-6;[n,m]=size(A);if n~=m %判断输...

选主元LU分解实验内容：列选主元LU分解和全选主元LU分解求解线性方程组计算方法：全选主元消元法1.1 初始化根据参数A、b,记录下矩阵、右端项的尺寸n；以得到的尺寸n初始化解向量x；同时，以尺寸(n-1)初始化一个记录列交换顺序的向量tcolum;1.2 进入循环for k=1:n-1选主元:|A(p1,p2)|=max{|A(i,j)|:i=k:n,j=k:n}记录列交...

1.semilogx()命名理解：semi- ：一半log :取对数x :表示仅x轴采用对数坐标用法：和plot()相同实际结果：对量级较高的函数（结果相差较大），semilogx()函数更能反映图线走势。例如：>> x=10.^(0.1:0.1:4);>> y=1./(x+1000);>> figure>> subplo...

双坐标轴绘制[AX,C1,C2]=plotyy(X1,Y1,X2,Y2)说明：输入参数：[X1,Y1]:曲线1（X、Y为等大数组）在左边[X2,Y2]:曲线2（X、Y为等大数组）在右边返回值：AX:1*2的坐标轴（一左一右）C1:曲线1句柄C2:曲线2句柄参数设置set():set(get(AX(1),‘Ylabel’),‘String’,‘left_name’)%返回曲线...

1.怎么获取绘图框里的n个点的坐标？[X,Y]=ginput(n)其中，n：获取n个数据点坐标X、Y: 返回数据点的坐标2.怎么将文本插入到绘图框的任意位置中？gtext(arg)命令一旦执行，将鼠标移到绘图框中任意位置单击，即可填入其中3.放缩已绘制图形命令zoom on:允许当前图形放缩zoom off:不允许当前图形放缩zoom out:返回原始状态zoom xon:只...

input1.提示用户输入 并将输入参数赋值给变量（数值）a=input(‘Enter a number:’)输入一个数值。如输入发生错误，会接着弹出Enter a numbermatlab中的双引号不同于单引号！2.需要输入参数为字符型时a=input(‘Enter a string:’,‘s’)第二参数 表明接收参数为字符型创建矩阵除了常规创建矩阵，还应注意：与python...

ctrl+c --终止正在执行的命令；在工作区，双击变量 可以根据需要绘制图像；gca --作为句柄，用于获取之前的图像并加以修正，例如：plot(x,y);%绘制图像set(gca,‘xticklabel’);%重新设置图像legend(‘para1’,‘para2’,‘para3’);%为图像设置图例说明在命令行窗口，查看当前文件目录：pwd进入新目录：cd <‘地址’&g...

当矩阵的行数等于列数时，计算矩阵的逆，可直接使用 inv(A);当矩阵的行数不等于列数时，可以考虑计算矩阵的Moore-Penrose逆，有两种方法：第一，直接使用Moore-Penrose逆的而定义B=inv(A’*A)*A’第二，使用命令B=pinv(A)即可求得矩阵A的Moore-Penrose逆。注：非方形矩阵的广义逆不止Moore-Penrose逆这一种情况。...

在命令行窗口输入矩阵A，>> a=[0.780 0.563;0.913 0.659]返回结果输出，a =0.7800 0.56300.9130 0.6590求该矩阵的逆，>>b=inv(a)返回结果输出，b =1.0e+05 *6.5900 -5.6300-9.1300 7.8000注，返回矩阵前的为科学记数法求矩阵的无穷范数，...