Matlab编程 入门(五)

                  Matlab(五)

1.希尔伯特(Hilbert)矩阵

Matlab由 A = hilb(3)生成

Trial>> A = hilb(3)

A =

    1.0000    0.5000    0.3333

    0.5000    0.3333    0.2500

    0.3333    0.2500    0.2000

Trial>> [n,d] = rat(A)

n =

     1     1     1

     1     1     1

     1     1     1

d =

     1     2     3

     2     3     4

     3     4     5

[n,d] = rat(A)//将A中元素转化成分数将分子放入矩阵n,分母放入d

2.基础数据变换和数论函数

(1) floor()

  n = floor(X)

将X中元素向－∞取整

(2) ceil()

  n = ceil(X)

将X中元素向＋∞取整

(3)round()

   n = round(X)

将X中元素四舍五入

(4)fix()

   n = fix(X)

将X中的元素向0取整

(5)rem()

  N = rem(A,C)

将A中元素对C取模求余

(6)gcd()

  A = gcd(m,n)

求整数m与n的最大公约数

Trial>> m = sym(1856120);n = sym(1483720);

Trial>> A = [gcd(m,n),lcm(m,n)]

A =

[ 1960, 1405082840]

(7)lcm()

  A = lcm(m,n)

求整数m与n的最小公倍数

(8)factor()

 M = factor(n)

对n进行质因数分解

Trial>> M = factor(456)

M =

     2     2     2     3    19

(9)isprime()

V1 = isprime(v)判断v的各个值是否质数,若是则置1否则置0

Trial>> v = 1:1:9;

Trial>> v1 = isprime(v)

v1 =

  1×9 logical 数组

   0   1   1   0   1   0   1   0   0

问:求0~100的全部质数

Trial>> v = 1:100;

Trial>> B = v(isprime(v))

B =

  1 至 13 列

     2     3     5     7    11    13    17    19    23    29    31    37    41

  14 至 25 列

    43    47    53    59    61    67    71    73    79    83    89    97