• 新记

    • 请matlab高手帮忙指点 以下是程序 望高手解释一下程序的意思及运行方法，这是LU分解 问题补充里面有程序

    function [x,L,U]= Doolittle (A,b)
    N = size(A);
    n = N(1);
    L = eye(n,n); %L的对角元素为1
    U = zeros(n,n);

    U(1,1:n) = A(1,1:n); %U的第一行
    L(1:n,1) = A(1:n,1)/U(1,1); %L的第一列

    for k=2:n
    for i=k:n
    U(k,i) = A(k,i)-L(k,1:(k-1))*U(1:(k-1),i);
    %U的第k行
    end
    for j=(k+1):n
    L(j,k) = (A(j,k)-L(j,1:(k-1))*U(1:(k-1),k))/U(k,k);
    %L的第k列
    end
    end

    y = SolveDownTriangle(L,b);
    x = SolveUpTriangle(U,y); %求解方程

    举报该文章
    其他看法
    第1个回答  2014-10-31
    先做LU分解,A=L*U。
    再解 L*y=b
    U*x=y
    还缺 SolveDownTriangle.m 和SolveUpTriangle.m 2个程序文件。
    否则 [x,L,U]= Doolittle (A,b) 就出结果。追问

    SolveDownTriangle.m 和SolveUpTriangle.m 2个程序文件我没有，还有就是我不知道怎么运行这个程序，怎么调用这程序呢。直接在命令窗口输入Doolittle (A,b)，老提示错误，怎么回事？请高手指点

    追答

    没有 SolveDownTriangle.m 和SolveUpTriangle.m 2个程序文件，程序不全。当然错。

    追问

    SolveDownTriangle.m 和SolveUpTriangle.m 2个程序文件怎么编写，有的话麻烦你给我发过来，谢谢。我也想知道到的是该怎么运行这个程序？

    帮忙解释一下matlab中这些程序的意义
    xlim([mu-3*sigma,mu+3*sigma]) %设置x轴显示范围 a2=axes; %新建一个坐标轴 ezplot(@(x)normpdf(x,mu,sigma),[mu-3*sigma,mu+3*sigma]) %在新坐标轴里画拟合得到的正态分布，显示范围在x取[mu-3*sigma,mu+3*sigma]之间 set(a2,'box','off','yaxislocation','right','color','...

    求直接三角分解法的matlab程序,每一步都需要解释,谢谢!
    你的代码我帮你解释了，果然是个复杂的活，如下function hl=zhjLU(A) %函数名为zhiLU，输入矩阵A,来进行LU分解，返回值hl为A的各阶主子式的行列式[n n] =size(A);%返回矩阵A的维数 RA=rank(A); %返回矩阵A的秩if RA~=n %判断矩阵A的秩RA是否不等于A的维数n，当不等于n时，即小于n...

    matlab里的一段程序 求大神解释一下每句的含义。尤其是程序的前面5行...
    使用fir2构造一个带限输入序列 clf %清除当前图形窗口图形。freq=[0 0.45 0.5 1] %定义一个1*4的矩阵freq。mag=[0 1 0 0] %定义一个1*4的矩阵mag。x=fir2(99,freq,mag) %设置滤波器。fir2(n,f,m)函数fir2用来设计多通带任意响应FIR滤波器，该滤波器的幅频特性由向量对f...

    求高手帮忙解两道MATLAB题SOS!在线等
    方法一:（矩阵左除法）x=A\\b 方法二：（矩阵求逆法）x=inv(A)*b 方法三：（LU分解法）[L,U]=lu(A);y=L\\b;x=U\\y 实验二：是说随便画个曲面和平面吗？ezmesh('x^2+y^2+5')hold on,ezmesh('1')画出来后，在图形上可以用工具按钮进行旋转视角的 ...

    求问下面matlab程序的解释,每一行都是什么意思
    三个for 对数组的每个元素进行循环 if 是如果 下标之和等于9 那么 d(m)就定义为 对应的a b c 的元素的和 并把当前下表存在E的第m行 然后把 m加一, 用来存储下一个 等于9 的位置 MAXnum用来存数组d的最大值 取出d的第一个元素, 如果d(i) 比 MAXnum大就替换成更大值 最后输出.找到d(...

    如何用matlab求解这个方程,试过了牛顿迭代和二分法,没成功,求高手指点...
    包含三角函数的方程是多值解，一般把所有的已知量全代进去，化为变量最少的方程，再求解。上面的方程，经简化为：72*cos(0.1177*cot(x)+0.1963)-sin(x-(0.1177*cot(x)+0.1963))*(2.7766+8.4734\/sin(x))-63.6612=0 则可求解：>> syms x f=72*cos(0.1177*cot(x)+0.1963)-...

    ...一道Matlab的题目,主要是写出可以正确运行的程序,谢谢啦
    “*”是利润，直线是欠款，交点就是还贷完成时，20年总利润-20年总贷款就是未来值，代码如下：clear clc all=20 one=5 p=5.58\/100 for i=1:20 debt if i<=4 own(i)=0;for j=1:i own(i)=one*(1+p)*j+own(i);end end if i>4 own(i)=one*p*4+own(i-1);end profit if...

    悬赏50分,求一段MATLAB程序的解释
    input sequence y = conv(h,x);%函数conv()用于求两个多项式的乘积多项式 n = 0:14;%n是0到14的区间 subplot(2,1,1);%将当前窗口分割成为2行，1列，并在其中的第1个区域绘图 stem(n,y);%Stem函数用来实现离散序列图的绘制 横坐标n个点，，依次值为y xlabel('Time index n'); ...

    求大神调试一下matlab程序,这个程序一会儿可以运行,一会儿又运行不了了...
    子函数combinat_help中，计算结果未返回，这会导致两个问题：函数返回值未定义，会导致报错；即使不报错，由于递归计算结果没有返回，所以也根本无法得到最终的结果。像这样的题目要求，使用递归方法设计程序不算是一个好的选择。最后提个建议：提问时请把代码贴出来，像你这样只贴图片，很少会有人有耐心...

    Matlab gui程序 数据被覆盖的问题,高手帮忙
    先clear后赋值。比如那个变量叫X；xlear X; % 先清除原数据 X = ……; % 赋值为新数据

    相似回答
    大家正在搜