rref代码（REF和RREF）

今天给各位分享rref代码的知识，其中也会对REF和RREF进行解释，如果能碰巧解决你现在面临的问题，别忘了关注本站，现在开始吧！

本文目录一览：

• 1、如何利用MatlAB软件用矩阵求逆、矩阵除法、矩阵分解求下面线性方程组的解？
• 2、通达信函数 REF 带的参数太少是什么意思？
• 3、帮忙做一个MATLAB的用来算阶梯型矩阵RREF()的代码
• 4、C语言错误：left operand must be l-value

如何利用MatlAB软件用矩阵求逆、矩阵除法、矩阵分解求下面线性方程组的解？

A=[5,1,-1,1;1,3,-1,-1;-1,-3,-1,5;1,1,-1,4]

C=inv(A) /*这是求矩阵 A的逆矩阵

B=[1 2 3 4]'

C*B得解

D=[A B]

rref(D) /*这是求矩阵D的行简化阶梯矩阵，得出未知数的解

矩阵的除法：A\B命令

通达信函数 REF 带的参数太少是什么意思？

这段通达信函数代码错误句可以修改,修改后正确效果如下图

帮忙做一个MATLAB的用来算阶梯型矩阵RREF()的代码

在matlab命令行输入

 type rref

就会看到系统自带的rref的源代码

function [A,jb] = rref(A,tol)

%RREF   Reduced row echelon form.

%   R = RREF(A) produces the reduced row echelon form of A.

%

%   [R,jb] = RREF(A) also returns a vector, jb, so that:

%       r = length(jb) is this algorithm's idea of the rank of A,

%       x(jb) are the bound variables in a linear system, Ax = b,

%       A(:,jb) is a basis for the range of A,

%       R(1:r,jb) is the r-by-r identity matrix.

%

%   [R,jb] = RREF(A,TOL) uses the given tolerance in the rank tests.

%

%   Roundoff errors may cause this algorithm to compute a different

%   value for the rank than RANK, ORTH and NULL.

%

%   Class support for input A:

%      float: double, single

%

%   See also RANK, ORTH, NULL, QR, SVD.

%   Copyright 1984-2005 The MathWorks, Inc. 

%   $Revision: 5.9.4.3 $  $Date: 2006/01/18 21:58:54 $

[m,n] = size(A);

% Does it appear that elements of A are ratios of small integers?

[num, den] = rat(A);

rats = isequal(A,num./den);

% Compute the default tolerance if none was provided.

if (nargin  2), tol = max(m,n)*eps(class(A))*norm(A,'inf'); end

% Loop over the entire matrix.

i = 1;

j = 1;

jb = [];

while (i = m)  (j = n)

   % Find value and index of largest element in the remainder of column j.

   [p,k] = max(abs(A(i:m,j))); k = k+i-1;

   if (p = tol)

      % The column is negligible, zero it out.

      A(i:m,j) = zeros(m-i+1,1);

      j = j + 1;

   else

      % Remember column index

      jb = [jb j];

      % Swap i-th and k-th rows.

      A([i k],j:n) = A([k i],j:n);

      % Divide the pivot row by the pivot element.

      A(i,j:n) = A(i,j:n)/A(i,j);

      % Subtract multiples of the pivot row from all the other rows.

      for k = [1:i-1 i+1:m]

         A(k,j:n) = A(k,j:n) - A(k,j)*A(i,j:n);

      end

      i = i + 1;

      j = j + 1;

   end

end

% Return "rational" numbers if appropriate.

if rats

    [num,den] = rat(A);

    A=num./den;

end

C语言错误：left operand must be l-value

error C2106: '=' :left operand must be l-value 表示等号左边必须为一个左值，该处错误为等号左边为*p+i\j是表达式，不是一个左值。

问题分析

C语言的符号优先级表格(部分)如下表：

由该表格可知取值符 * 运算优先级高于加法运算符 +，因此 *p+i 和 *p+j 的运算均为先取值再加上变量，结果为一个数值，无法作为左值被赋值，只能用于赋值或比较故产生了该错误。

解决方法

该设计思路为先计算p+i作为地址再进行取值操作，故此处可以使用括号修改代码，修改结果如下：

此时赋值符 = 左边均为左值即可被赋值的变量或存储空间。故使用赋值符 = 时应对其左边的内容稍加关注，避免此类由于优先级导致的语法错误。

扩展资料

值得注意的是，赋值符正确使用时会附带隐含的类型转换，如果赋值运算符两边的数据类型不相同，系统将自动进行类型转换，即把赋值号右边的类型换成左边的类型。具体规定如下：

1、实型赋予整型，舍去小数部分。此处应注意C语言的舍去一般采用截断舍去，对整数部分不影响。

2、整型赋予实型，数值不变，但将以浮点形式存放，即增加小数部分，一般小数部分的值为0，但由于浮点数精度问题可能产生误差。

3、字符型赋予整型，由于字符型为一个字节，而整型为二个字节，故将字符的ASCII码值放到整型量的低八位中，高八位为0。整型赋予字符型，只把低八位赋予字符量。此处实际赋值时与选择的编译器有关，具体编译器需要具体分析。

关于rref代码和REF和RREF的介绍到此就结束了，不知道你从中找到你需要的信息了吗 ？如果你还想了解更多这方面的信息，记得收藏关注本站。