Octave教程五：控制语句if/while/for以及函数的定义和使用
Octave教程五：控制语句if/while/for以及函数的定义和使用
在这篇文章中，我将告诉你该怎样为你的Octave程序写控制语句，诸如if/while/for这些语句，以及如何定义和使用方程。
>> 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
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
>> a = 1:10;
>> a
a =
1 2 3 4 5 6 7 8 9 10
>> v
v =
2 4 8 16 32 64 128 256 512 1024
>> V
V =
0
0
0
0
0
0
0
0
0
0
>> v = v'
v =
2
4
8
16
32
64
128
256
512
1024
>> 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 ++;
> 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
>>现在我们来说如何在Octave里定义函数：
这个文件只有三行，第一行写着function y = squareThisNumber(x)，这就告诉octave，我想返回一个值，并且这个返回的值，我想放在变量y中；另外，它告诉了octave这个函数有一个参数，就是参数x，还有定义的函数体，也就是y等于x的平方。
现在让我们尝试调用这个函数。
>> squareThisNumber(5)
error: 'squareThisNumber' undefined near line 1 column 1
>>出错：这个函数未被定义
这是因为octave不知道在哪里找这个函数
>> pwd
ans = C:\fakepath\octave4
>> cd 'D:\soft\Octave\user'
>> squareThisNumber(5)
parse error near line 1 of file D:\soft\Octave\user\squareThisNumber.m
syntax error
>>>
^出错，修改SquareThisNumber.m文件，如下所示：
此时，程序便可正常运行：
>> SquareThisNumber(5)
ans = 25
>> SquareThisNumber(5)
ans = 25
>> SquareThisNumber(8)
ans = 64
>>>> % Octave search path (advanced/optional)
>> addpath ('C:\fakepath\octave4')
>> SquareThisNumber(5)
ans = 25
>> pwd
ans = D:\soft\Octave\user
>> addpath ('D:\soft\Octave\user')
>> cd 'C:\'
>> SquareThisNumber(5)
ans = 25
>> cd
>> pwd
ans = C:\Users\huyang
>>我们再来看另外一种情况：
我们在D:/soft/Octave/user目录下创建squareAddCubeThisNumber.m文件，内容为：
>>
>> [a,b] = squareAddCubeThisNumber(5)
a = 25
b = 125
>>>> [a,b] = squareAddCubeThisNumber(5);
>> a
a = 25
>> b
b = 125
>>接下来，我想设计一个复杂一点的函数：
比方说，我有一个数据集，像这样，数据点为[1,1],[2,2],[3,3]
我想做的事是定义一个octave函数来计算代价函数J（θ），就是计算不同θ值所对应的代价函数值J。
首先让我们把数据放在octave里，如下所示：
>> 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];
>>接下来，我定义一个costFunctionJ.m函数，如下图所示：
源代码如下：
function J = costFunctionJ(X,Y,theta)
%X is the "design matrix" containing our training examples
%Y is the class labels
m = size(X,1); %number of training examples
predictions = X * theta; %predictions of hypothesis on all m examples
sqrErrors = (predictions - Y) .^ 2; % squared errors
J = 1/(2*m) * sum(sqrErrors);接下来，我们写如下代码：
>> P = Y .^ 2
P =
1
4
9
>> sum(P)
ans = 14
>> J = costFunctionJ(X,Y,theta)
J = 0
>>屏幕输出J = 0，这是因为，如果我的数据集x为[1;2;3]，y也为[1;2;3]，然后设置θ0为0，θ1 = 1，这就恰好是45°的斜线，这条线是可以完美拟合我的数据集的；而相反地，如果我设置θ = [0;0]，那么这个假设所计算出的J就如下图所示：
>>
>> 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
>>