Octave教程

Windows安装Octave

http://wiki.octave.org/Octave_for_Microsoft_Windows

基本操作(Basic Operations)

octave:1> PS1('>> ');
>>
>> path
>> pause
>> clear ; close all; clc
>> fprintf(['test print %%s : %s \n' ...
         'test print %%d : %d \n' ...
         'test print %%g : %g \n' ...
         'test print %%f : %f \n' ...  
        ], 'String', 10, 2.33553e-11, 0.12345678);
test print %s : String 
test print %d : 10 
test print %g : 2.33553e-11 
test print %f : 0.123457 
>>
>> 1+2
ans =  3
>>  5*8
ans =  40
>>  1/2
ans =  0.50000
>>  2^6
ans =  64
>>  2^3
ans =  8
>>  1 == 2 % false
ans = 0
>>  1 == 1 % true
ans = 1
>>  1 ~= 2 % true
ans = 1
>>  1 ~= 2 % true
ans = 1
>>  1 ~= 1 % false
ans = 0
>>  1 ~= 1 % false  not equal
ans = 0
>>  1 && 0
ans = 0
>>  1 && 0 % AND
ans = 0
>>  1 || 0 % or
ans = 1
>>  xor(1,0)
ans = 1
>>  xor(1,0) % yi or
ans = 1
>> PS1('>> ');
>> a=3
a =  3
>> a=4;
>> a
a =  4
>> c=(3>=4)
c = 0
>> c
c = 0
>> a=pi
a =  3.1416
>> a=pi;
>> a
a =  3.1416
>> b=e;
>> b
b =  2.7183
>> disp(a);
 3.1416
>> disp(sprintf('2 decimals: %0.2f', a))
2 decimals: 3.14
>> disp(sprintf('10 decimals: %0.10f', a))
10 decimals: 3.1415926536
>> disp(sprintf('10 decimals: %0.10f', b))
10 decimals: 2.7182818285
>> format long
>> a
a =  3.141592653589793
>> b
b =  2.718281828459045
>>  format short
>> a
a =  3.1416
>> b
b =  2.7183
>> A = [1 2; 3 4; 5 6];
>>  A
A =

   1   2
   3   4
   5   6

>>  A
A =

   1   2
   3   4
   5   6

>>  v=[1 2 3];
>>  v =[1; 2; 3]
v =

   1
   2
   3

>>  v=1:0.1:2;
>>  v
v =

 Columns 1 through 7:

    1.0000    1.1000    1.2000    1.3000    1.4000    1.5000    1.6000

 Columns 8 through 11:

    1.7000    1.8000    1.9000    2.0000

>>  v=1:6; 
>>  v
v =

   1   2   3   4   5   6

>>  ones(2,3)
ans =

   1   1   1
   1   1   1

>>  2 * ones(2,3)
ans =

   2   2   2
   2   2   2

>>  w=ones(1,3)
w =

   1   1   1

>>  w=zeros(1,3)
w =

   0   0   0

>>  w=rand(1,3)
w =

   0.066378   0.828597   0.569681

>>  rand(3,3)
ans =

   0.2443954   0.3866398   0.2332407
   0.6611382   0.6353541   0.0066938
   0.2465706   0.2734220   0.6748704

>>  rand(3,3)
ans =

   0.71983   0.47680   0.30126
   0.19673   0.29252   0.57751
   0.99891   0.57380   0.31022

>>  rand(3,3)
ans =

   0.26842   0.82579   0.22678
   0.63709   0.67602   0.11723
   0.83274   0.48216   0.26778

>>  rand(3,3)
ans =

   0.790568   0.666798   0.084361
   0.343994   0.482338   0.860989
   0.908285   0.967817   0.935013

>>  randn(1,3) % gaosi fenbu
ans =

  -0.098394  -0.147581   0.085802

>>  randn(1,3) % gaosi fenbu
ans =

   1.10694  -1.99933  -0.17373

>>  w= -6 + sqrt(10) * (randn(1,10000));
>>  hist(w)
>>  hist(w,50)
>> 
>>  I = eye(4)
I =

Diagonal Matrix

   1   0   0   0
   0   1   0   0
   0   0   1   0
   0   0   0   1

>>  help eye
>>  help help

移动数据(Moving Data Around)

>>  A = [1 2; 3 4; 5 6]
A =

   1   2
   3   4
   5   6

>>  size(A)
ans =

   3   2

>>  size(A, 1)
ans =  3
>>  size(A, 2)
ans =  2
>>  size(A, 3)
ans =  1
>>  v = [1 2 3 4]
v =

   1   2   3   4

>>  length(v)
ans =  4
>>  length(A)
ans =  3

>> load featuresX.dat
>> load ('priceY.dat')
>> who
Variables in the current scope:

A          ans        featuresX  priceY     v

>> whos
Variables in the current scope:

   Attr Name           Size                     Bytes  Class
   ==== ====           ====                     =====  =====
        A              3x2                         48  double
        ans            1x23                        23  char
        featuresX      3x1                         24  double
        priceY         3x1                         24  double
        v              1x4                         32  double

Total is 39 elements using 151 bytes

>> priceY
priceY =

   1000
   2000
   3000

>> featuresX
featuresX =

   1234
   5678
   1235

>> clear featuresX
>>
>> whos
Variables in the current scope:

   Attr Name        Size                     Bytes  Class
   ==== ====        ====                     =====  =====
        A           3x2                         48  double
        ans         1x23                        23  char
        priceY      3x1                         24  double
        v           1x4                         32  double

Total is 36 elements using 127 bytes

>> clear
>> whos
>> load ('priceY.dat')
>> v = priceY(2:3)
v =

   2000
   3000

>> save hello.mat v;
>> clear
>> whos
>> load hello.mat
>> whos
Variables in the current scope:

   Attr Name        Size                     Bytes  Class
   ==== ====        ====                     =====  =====
        v           2x1                         16  double

Total is 2 elements using 16 bytes

>> v
v =

   2000
   3000

>> save hello.txt -ascii % save as text (ASCII)
>> A = [1 2;3 4;5 6]
A =

   1   2
   3   4
   5   6

>> A(3,2)
ans =  6
>>
>> A(2,:) % ":" means every element along that row/column
ans =

   3   4

>> A(:,2)
ans =

   2
   4
   6

>> A([1 3], :)
ans =

   1   2
   5   6

>> A
A =

   1   2
   3   4
   5   6

>> A([1 2 3], :)
ans =

   1   2
   3   4
   5   6

>> A(:,2) = [10; 20; 30]
A =

    1   10
    3   20
    5   30

>> A = [A, [100; 101; 102]] % append another column vector to right
A =

     1    10   100
     3    20   101
     5    30   102

>> A = [A; [100 101 102]] % append another column vector to right
A =

     1    10   100
     3    20   101
     5    30   102
   100   101   102

>> A(:)
ans =

     1
     3
     5
   100
    10
    20
    30
   101
   100
   101
   102
   102

>> A = [1 2;3 4;5 6]
A =

   1   2
   3   4
   5   6

>> B = [11 12; 13 14; 15 16]
B =

   11   12
   13   14
   15   16

>> C = [A B]
C =

    1    2   11   12
    3    4   13   14
    5    6   15   16

>> [B A]
ans =

   11   12    1    2
   13   14    3    4
   15   16    5    6

>> [A B C]
ans =

    1    2   11   12    1    2   11   12
    3    4   13   14    3    4   13   14
    5    6   15   16    5    6   15   16

>> [A; B]
ans =

    1    2
    3    4
    5    6
   11   12
   13   14
   15   16

>> [A, B]
ans =

    1    2   11   12
    3    4   13   14
    5    6   15   16

计算数据(Computing on Data)

>> A = [1 2; 3 4; 5 6]
A =

   1   2
   3   4
   5   6

>> B = [11 12; 13 14; 15 16]
B =

   11   12
   13   14
   15   16

>> C = [1 1; 2 2]
C =

   1   1
   2   2

>> A * C
ans =

    5    5
   11   11
   17   17

>> A .* B
ans =

   11   24
   39   56
   75   96

>> 3 * A
ans =

    3    6
    9   12
   15   18

>> A*3
ans =

    3    6
    9   12
   15   18

>> A
A =

   1   2
   3   4
   5   6

>> B
B =

   11   12
   13   14
   15   16

>> A.^2
ans =

    1    4
    9   16
   25   36

>> v = [1; 2; 3]
v =

   1
   2
   3

>> 1./v
ans =

   1.00000
   0.50000
   0.33333

>> 1./A
ans =

   1.00000   0.50000
   0.33333   0.25000
   0.20000   0.16667

>> log(v)
ans =

   0.00000
   0.69315
   1.09861

>> exp(v)
ans =

    2.7183
    7.3891
   20.0855

>> abs(v)
ans =

   1
   2
   3

>> abs([-1; 2; -3])
ans =

   1
   2
   3

>> -v
ans =

  -1
  -2
  -3

>> -A
ans =

  -1  -2
  -3  -4
  -5  -6

>> -1*v
ans =

  -1
  -2
  -3

>> v+ones(length(v),1)
ans =

   2
   3
   4

>> ones(length(v),1)
ans =

   1
   1
   1

>> ones(length(v),2)
ans =

   1   1
   1   1
   1   1

>> zeros(length(v),2)
ans =

   0   0
   0   0
   0   0

>> v.+1
ans =

   2
   3
   4

>> v.-1
ans =

   0
   1
   2

>> 1.+v
ans =

   2
   3
   4

>> v+1
ans =

   2
   3
   4

>> v-1
ans =

   0
   1
   2

>> A
A =

   1   2
   3   4
   5   6

>> A'
ans =

   1   3   5
   2   4   6

>> A''
ans =

   1   2
   3   4
   5   6

>> A''''
ans =

   1   2
   3   4
   5   6

>> a = [1 15 2 0.5]
a =

    1.00000   15.00000    2.00000    0.50000

>> val = max(a)
val =  15
>> [val, ind] = max(a)
val =  15
ind =  2
>> max(A)
ans =

   5   6

>> a < 3
ans =

  1  0  1  1

>> A < 3
ans =

  1  1
  0  0
  0  0

>> find(a<3)
ans =

   1   3   4

>> a
a =

    1.00000   15.00000    2.00000    0.50000

>> magic(3)
ans =

   8   1   6
   3   5   7
   4   9   2

>> magic(4)
ans =

   16    2    3   13
    5   11   10    8
    9    7    6   12
    4   14   15    1

>> magic(3)
ans =

   8   1   6
   3   5   7
   4   9   2

>> magic(3)
ans =

   8   1   6
   3   5   7
   4   9   2

>> magic(3)
ans =

   8   1   6
   3   5   7
   4   9   2

>> magic(3)
ans =

   8   1   6
   3   5   7
   4   9   2

>> [r,c] = find(A >= 7)
r = [](0x1)
c = [](0x1)
>> A
A =

   1   2
   3   4
   5   6

>> A=magic(3)
A =

   8   1   6
   3   5   7
   4   9   2

>> [r,c] = find(A >= 7)
r =

   1
   3
   2

c =

   1
   2
   3

>> [r,c] = find(A >= 2)
r =

   1
   2
   3
   2
   3
   1
   2
   3

c =

   1
   1
   1
   2
   2
   3
   3
   3

>> a
a =

    1.00000   15.00000    2.00000    0.50000

>> sum(a)
ans =  18.500
>> prod(a)
ans =  15
>> floor(a)
ans =

    1   15    2    0

>> ceil(a)
ans =

    1   15    2    1

>> rand(3)
ans =

   0.43298   0.41511   0.89855
   0.21145   0.15275   0.96643
   0.67367   0.31471   0.15306

>> max(rand(3), rand(3))
ans =

   0.99326   0.98554   0.83527
   0.78799   0.71117   0.92512
   0.99326   0.98996   0.34543

>> A
A =

   8   1   6
   3   5   7
   4   9   2

>> max(A, [], 1)
ans =

   8   9   7

>> max(A, [], 2)
ans =

   8
   7
   9

>> max(A)
ans =

   8   9   7

>> max(max(A))
ans =  9
>> A(:)
ans =

   8
   3
   4
   1
   5
   9
   6
   7
   2

>> max(A(:))
ans =  9
>>
>>
>>
>> A = magic(9)
A =

   47   58   69   80    1   12   23   34   45
   57   68   79    9   11   22   33   44   46
   67   78    8   10   21   32   43   54   56
   77    7   18   20   31   42   53   55   66
    6   17   19   30   41   52   63   65   76
   16   27   29   40   51   62   64   75    5
   26   28   39   50   61   72   74    4   15
   36   38   49   60   71   73    3   14   25
   37   48   59   70   81    2   13   24   35

>> sum(A,1)
ans =

   369   369   369   369   369   369   369   369   369

>> sum(A,2)
ans =

   369
   369
   369
   369
   369
   369
   369
   369
   369

>> eye(9)
ans =

Diagonal Matrix

   1   0   0   0   0   0   0   0   0
   0   1   0   0   0   0   0   0   0
   0   0   1   0   0   0   0   0   0
   0   0   0   1   0   0   0   0   0
   0   0   0   0   1   0   0   0   0
   0   0   0   0   0   1   0   0   0
   0   0   0   0   0   0   1   0   0
   0   0   0   0   0   0   0   1   0
   0   0   0   0   0   0   0   0   1

>> A .* eye(9)
ans =

   47    0    0    0    0    0    0    0    0
    0   68    0    0    0    0    0    0    0
    0    0    8    0    0    0    0    0    0
    0    0    0   20    0    0    0    0    0
    0    0    0    0   41    0    0    0    0
    0    0    0    0    0   62    0    0    0
    0    0    0    0    0    0   74    0    0
    0    0    0    0    0    0    0   14    0
    0    0    0    0    0    0    0    0   35

>> sum(A .* eye(9))
ans =

   47   68    8   20   41   62   74   14   35

>> sum(sum(A .* eye(9)))
ans =  369
>> flipud(eye(9))
ans =

Permutation Matrix

   0   0   0   0   0   0   0   0   1
   0   0   0   0   0   0   0   1   0
   0   0   0   0   0   0   1   0   0
   0   0   0   0   0   1   0   0   0
   0   0   0   0   1   0   0   0   0
   0   0   0   1   0   0   0   0   0
   0   0   1   0   0   0   0   0   0
   0   1   0   0   0   0   0   0   0
   1   0   0   0   0   0   0   0   0

>> sum(sum(A .* flipud(eye(9))))
ans =  369
>>
>>
>>
>>
>> A = magic(3)
A =

   8   1   6
   3   5   7
   4   9   2

>> pinv(A) % pseudoinverse
ans =

   0.147222  -0.144444   0.063889
  -0.061111   0.022222   0.105556
  -0.019444   0.188889  -0.102778

>> inv(A)
ans =

   0.147222  -0.144444   0.063889
  -0.061111   0.022222   0.105556
  -0.019444   0.188889  -0.102778

>> temp = pinv(A)
temp =

   0.147222  -0.144444   0.063889
  -0.061111   0.022222   0.105556
  -0.019444   0.188889  -0.102778

>> temp * A
ans =

   1.00000   0.00000  -0.00000
  -0.00000   1.00000   0.00000
   0.00000   0.00000   1.00000

>> abs(temp * A)
ans =

   1.00000   0.00000   0.00000
   0.00000   1.00000   0.00000
   0.00000   0.00000   1.00000

数据绘制(Plotting Data)

>>  t=[0:0.01:0.98];
>>  t
t =

 Columns 1 through 11:

   0.00000   0.01000   0.02000   0.03000   0.04000   0.05000   0.06000   0.07000   0.08000   0.09000   0.10000

 Columns 12 through 22:

   0.11000   0.12000   0.13000   0.14000   0.15000   0.16000   0.17000   0.18000   0.19000   0.20000   0.21000

 Columns 23 through 33:

   0.22000   0.23000   0.24000   0.25000   0.26000   0.27000   0.28000   0.29000   0.30000   0.31000   0.32000

 Columns 34 through 44:

   0.33000   0.34000   0.35000   0.36000   0.37000   0.38000   0.39000   0.40000   0.41000   0.42000   0.43000

 Columns 45 through 55:

   0.44000   0.45000   0.46000   0.47000   0.48000   0.49000   0.50000   0.51000   0.52000   0.53000   0.54000

 Columns 56 through 66:

   0.55000   0.56000   0.57000   0.58000   0.59000   0.60000   0.61000   0.62000   0.63000   0.64000   0.65000

 Columns 67 through 77:

   0.66000   0.67000   0.68000   0.69000   0.70000   0.71000   0.72000   0.73000   0.74000   0.75000   0.76000

 Columns 78 through 88:

   0.77000   0.78000   0.79000   0.80000   0.81000   0.82000   0.83000   0.84000   0.85000   0.86000   0.87000

 Columns 89 through 99:

   0.88000   0.89000   0.90000   0.91000   0.92000   0.93000   0.94000   0.95000   0.96000   0.97000   0.98000

>>  y1=sin(2*pi*4*t);
>>  plot(t,y1);
>>  y2=cos(2*pi*4*t);
>>  plot(t,y2);

>>  plot(t,y1);
>>  hold on;
>>  plot(t,y2,'r');
>>  xlabel('time');
>>  ylabel('value');
>>  legend('sin','cos');
>>  legend('cos','sin');
>>  legend('sin','cos');
>>  title('my plot');

>>  print -dpng 'myPolt.png'
>>  close

>>  figure(1); plot(t,y1);
>>  figure(2); plot(t,y2,'r');

>>  subplot(1,2,1);
>>  plot(t,y1)
>>  subplot(1,2,2)
>>  plot(t,y2,'r')
>>  title myPlot
>>  axis([0.5 1 -1 1])
>>  clf;

>>  A=magic(5);
>>  imagesc(A);
>>  imagesc(A), colorbar, colormap gray;
>>  imagesc(magic(15)), colorbar, colormap gray;

控制语句:for,while,if语句(Control Statements)

>>  v=zeros(10,1)
v =

   0
   0
   0
   0
   0
   0
   0
   0
   0
   0

>>  for i=1:10
>   v(i) = 2^i;
> end
>>  v
v =

      2
      4
      8
     16
     32
     64
    128
    256
    512
   1024

>>  indices=1:10
indices =

    1    2    3    4    5    6    7    8    9   10

>>  for i=indices
>     disp(i)
>   end
 1
 2
 3
 4
 5
 6
 7
 8
 9
 10
>> 
>>  v
v =

      2
      4
      8
     16
     32
     64
    128
    256
    512
   1024

>>  i=1
i =  1
>>  while i <=5
>     v(i) = 100;
>     i = i +1;
>   end;
>>  v
v =

    100
    100
    100
    100
    100
     64
    128
    256
    512
   1024

>>  i=1;
>>  while true
>     v(i)=999;
>     i = i + 1;
>     if i==6,
>       break;
>     end;
>   end;
>>  v
v =

    999
    999
    999
    999
    999
     64
    128
    256
    512
   1024

>>  v(1)
ans =  999
>>  v(1)=2;
>>  if v(1)==1,
>     disp('The value is one');
>   elseif v(1) == 2,
>     disp('The value is two');
>   else
>     disp('The value is not one or two');
>   end;
The value is two

函数的定义和使用(Defining and using functions)

>>  function y = squareThisNumber(x);
>   y = x^2;
>   end;
>> 
>>  squareThisNumber(5)
ans =  25

>>  % octave serch path (advanced/optional)
>>  addpath('C:\users\');

>>  % you can also put the function into the 'suareThisNumber.m' file under addpath
>>  % then you can invoke this function use 'squareThisNumber(5)'

>>  function [y1,y2] = squareAndCubeThisNumber(x);
>     y1 = x^2;
>     y2 = x^3;
>   end;
>>  squareAndCubeThisNumber(3)
y1 =  9
y2 =  27
ans =  9
>>  [a,b]=squareAndCubeThisNumber(3)
y1 =  9
y2 =  27
a =  9
b =  27
>>  a
a =  9
>>  b
b =  27

>>  X = [1 1; 1 2; 1 3]
X =

   1   1
   1   2
   1   3

>>  y = [1; 2; 3]
y =

   1
   2
   3

>>  theta = [0;1]
theta =

   0
   1

>>  function J = costFunctionJ(X, y, theta)
>   % X is the "design matrix" containing out training examples.
>   % y is the class labels
>
>   m = size(X, 1);                     % number of traning examples
>   predictions = X * theta;            % predictions of hypothesis on all m examples
>   sqrErrors = (predictions - y).^2;   % squared errors
>
>   J = 1/(2*m) * sum(sqrErrors);
>   end;
>>  j=costFunctionJ(X,y,theta)
j = 0

>>  size(X,1)
ans =  3
>>  X*theta
ans =

   1
   2
   3

>>  sum(X)
ans =

   3   6

>>  A = magic(9)
A =

   47   58   69   80    1   12   23   34   45
   57   68   79    9   11   22   33   44   46
   67   78    8   10   21   32   43   54   56
   77    7   18   20   31   42   53   55   66
    6   17   19   30   41   52   63   65   76
   16   27   29   40   51   62   64   75    5
   26   28   39   50   61   72   74    4   15
   36   38   49   60   71   73    3   14   25
   37   48   59   70   81    2   13   24   35

>>  sum(A,1)
ans =

   369   369   369   369   369   369   369   369   369

>>  sum(A,2)
ans =

   369
   369
   369
   369
   369
   369
   369
   369
   369

>>  sum(A)
ans =

   369   369   369   369   369   369   369   369   369

>>  sum(sum(A))
ans =  3321

>>  theta=[0;0]
theta =

   0
   0

>>  j=costFunctionJ(X,y,theta)
j =  2.3333
>>  (1^2 + 2^2 + 3^2)/(2*3)
ans =  2.3333

转载于:https://www.cnblogs.com/keyshaw/p/10654226.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值