Octave 控制语句

深藏阁楼爱情的钟 2022-04-18 05:12 288阅读 0赞

首先将v值设为一个10行1列的零向量。

![20181108204202937.png][]

现在写一个for循环，让i等于1到10，写出来就是 for i = 1:10 然后设v(i)的值等于2的i次方 循环的最后记得写上end; 这里的空格没有什么关系，加了一些空格只是为了让缩进后的代码看起来结构更清晰，但其实这里的空格没有意义，然后打印结果：

![20181108204954738.png][]

另一种方式是设置indices（索引）等于1到10，这时indices就是一个从1到10的序列。

![20181108205239891.png][]

那么也可以写for i = indices 这实际上和for i = 1:10是一样的 也可以写disp(i) 也能得到一样的结果。所以这就是一个for循环。

![20181108205548727.png][]

如果你对break语句和continue语句比较熟悉，也可以在Octave环境的循环里使用。但是首先需要先告诉你while循环是如何工作的？

那么先来写个while循环：

![20181108210512952.png][]

所以这是什么意思呢？让i的取值从1开始，然后让v(i)等于100，再让i递增1直到i大于5停止。

![20181108210843921.png][]

现在来看一下结果，原来的向量v是2的次方，而现在已经把向量v的前5个元素用100覆盖掉了，这就是一个while循环的句法结构。

现在来分析另外一个例子，这里将展示如何使用break语句：

![20181108211727403.png][]

这个逻辑是让i等于1，然后开始下面的增量循环，while语句重复设置v(i)等于999，不断让i加1，然后当i达到6时，就用break中止while循环，所以最好的效果是取向量v的前5个元素，把它们设置为999。

![20181108212034680.png][]

这就是if语句的句法结构和while语句，并且注意要有end，这里是有两个end的，第一个end结束的是if语句，第二个end结束的是while语句。

现在来看使用if-else语句时更一般的句法结构。

v(1)现在是等于999的，在下面的例子中，令v(1)等于2。

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMzExNjA0_size_16_color_FFFFFF_t_70][]

如果需要退出Octave，可以键入exit命令，然后回车，就会退出Octave，或者命令quit也可以。

最后来说一说函数，如何定义和调用函数？

这是我Octave的当前路径，在里面存了一个预先定义的文件，名为"squarethisnumber.m"。在Octave环境下定义函数，你需要创建一个文件，用你的函数名来命名，然后以.m的后缀结尾。当Octave发现这文件，它知道应该在什么位置查找squarethisnumber.m这个函数的定义。

这里只是想建议你，如果你使用的是微软的Windows系统，那么可以使用写字板程序而不是记事本来打开这些文件。如果你有别的什么文本编辑器也是可以的。记事本有时会把代码的间距弄得很乱。如果你只有记事本程序，那也能用。但如果你有写字板的话，建议你使用写字板或者其他的文本编辑器来编辑函数。

现在来说如何在Octave里定义函数

![20181110110511523.png][]

文件有三行，第一行写了`function y = squareThisNumber(x) 这就告诉Octave我想返回一个值把它存放在变量y里，另外它告诉了Octave，这个函数有一个自变量x。定义的函数主体是y = x^2`

`现在我们尝试来调用这个函数`

![20181110103016286.png][]

这样实际上是行不通的，Octave说这个函数未被定义，这是因为Octave不知道在哪里找这个文件。所以像之前一样，使用pwd，现在不在正确目录下(这个函数在"C:\\Octave\\program"里面)。

![20181110103202653.png][]

因此我们把路径设为"C:\\Octave\\program"。

![20181110105029522.png][]

出错，修改squareThisNumber.m文件。**注意，如果直接在word文档里改，可能还是会出现上张图的那种错误，可以进入Octave的squareThisNumber.m的这个文件里改。**如下图所示：

![20181110110545310.png][]

现在如果键入squareThisNumber(5)，返回值是25：

![2018111010530695.png][]

还有一种更高级的功能：

如果你知道search path(搜索路径)这个术语，如果你想要修改Octave的搜索路径。你可以把下面这部分作为一个进阶知识。仅适用于那些熟悉编程语言中搜索路径概念的人。

你可以使用addpath命令，添加路径"C:\\Octave\\program"，将该目录添加到Octave的搜索路径，这样即使你跑到了其他路径下面，Octave依然知道会在"C:\\Octave\\program"目录下寻找函数。这样，即使现在在不同的目录下，它仍然知道在哪里可以找到"squareThisNumber"这个函数。

![20181110110343703.png][]

但是，如果你不熟悉搜索路径的概念，也没关系，只要确保在执行函数之前先用cd命令设置到函数所在的目录下，实际上也是一样的效果。

Octave还有一个其他许多编程语言都没有的概念。那就是它可以允许你定义一个函数能够返回多个值。这里有一个例子，定义一个函数叫"squareAndCubeThisNumber(x)"，这说的就是该函数返回值是两个值y1和y2，y1是x^2，y2是x^3。这就是说函数会真的返回两个值。所以你们可能的熟悉C、C++或Java通常情况下认为作为函数返回值只能是一个值，但Octave的语法结构就不一样，可以返回多个值。

![2018111011145440.png][]

现在如果键入\[a,b\] = squareAndCubeThisNumber(5);

![20181110112106273.png][]

所以这样往往很方便，如果你需要定义一个函数并且返回多个值。

最后来演示一个更复杂一点的函数的例子。

例如，现在有一个数据集，如下图所示，数据点为(1,1)，(2,2)，(3,3)。

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMzExNjA0_size_16_color_FFFFFF_t_70 1][]

要做的就是定义一个Octave函数来计算代价函数J(θ)，不同θ值所对应的代价函数值J，首先把数据放到Octave里。

把设计矩阵设置为 X = \[1 1;1 2;1 3\]

![20181110113321219.png][]

第一项表示x\_0项，矩阵的第一列。第二列表示训练样本的三个x值。

再来设置y = \[1;2;3\]

![20181110113600845.png][]

这三个数字是y轴的对应值。

现在我们设定θ为\[0;1\]

![20181110113843737.png][]

在我的Octave的当前路径里已经有我预先定义的代价函数J，如果打开函数的定义是下面这样的，里面有一些注释用于解释输入变量。接下来几步，设定m为训练样本的数量，因此行数为X。然后计算预测值X\*theta。再下面就是计算平方误差，预测值减去y值，然后每一项进行平方。最后，就可以计算代价函数J，并且Octave知道J是要返回的值，因为J出现在了函数的定义里。

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMzExNjA0_size_16_color_FFFFFF_t_70 2][]

然后键入 J = costFunctionJ(X,y,theta) 就会计算出J = 0。

![20181110155122413.png][]

这是因为，如果数据集X为\[1;2;3\]，y也为\[1;2;3\]，然后设置theta\_0 = 0，theta\_1 = 1，这样形成的假设函数就恰好是倾斜度为45度的直线，这样就可以完美拟合数据集。而相反地，如果设置theta = \[0;0\]，也就是theta\_0 = 0，theta\_1 = 0，那么这个假设预测所有的值都为0，然后再来计算代价函数。

![20181110155548635.png][]

实际上就等于1的平方（也就是第一个样本的平方误差），加上2的平方，再加上3的平方。然后除以2\*m（也就是训练样本数的2倍）。

![20181110155856569.png][]

答案也是2.3333，因此这也就验证了，这里的函数计算出了正确的代价函数。至少用这几个简单的训练样本做的验证是没问题的。完整性检查证明定义的代价函数J确实是可以计算出正确的代价函数的，至少基于这里的训练集成立。也就是这里的x和y，这几个简单的训练样本。

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMzExNjA0_size_16_color_FFFFFF_t_70 3][]

[20181108204202937.png]: /images/20220418/0436d78469b4480b838d329edcda9989.png
[20181108204954738.png]: /images/20220418/b67bda86d1df44d788654048b4ce4dff.png
[20181108205239891.png]: /images/20220418/400e8506a01e47718c480426f7329451.png
[20181108205548727.png]: /images/20220418/d11e75322ff143a4936df53a6ff1ed7e.png
[20181108210512952.png]: /images/20220418/005a7b53ec5c4d558c1a7b2f728b7927.png
[20181108210843921.png]: /images/20220418/e57b807453a64ace8ff15947087a0b59.png
[20181108211727403.png]: /images/20220418/0784cbfc95ee4dd4b82e240f9afab98a.png
[20181108212034680.png]: /images/20220418/d4b919f781284d5f9f069a54e062127c.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMzExNjA0_size_16_color_FFFFFF_t_70]: /images/20220418/5c7a026e24184f55bb35196eb62ee8f8.png
[20181110110511523.png]: /images/20220418/89992f3c3e704751ab6351840c6f7c4a.png
[20181110103016286.png]: /images/20220418/09549c62aa96435b95c25b68ac99ab5c.png
[20181110103202653.png]: /images/20220418/df062e7c67784b4099cf88e718b24da3.png
[20181110105029522.png]: /images/20220418/8bd2de41ceb94431ba20346fb7c15dcf.png
[20181110110545310.png]: /images/20220418/71d0d8cce56e4ea09fd3b38a77d80ce3.png
[2018111010530695.png]: /images/20220418/815e997eea28481c942686f6842a71bd.png
[20181110110343703.png]: /images/20220418/afcb07a37a3b44d3adc953b689728841.png
[2018111011145440.png]: /images/20220418/c49a832a1f744faf85d2e05cca82371f.png
[20181110112106273.png]: /images/20220418/0e88d24a325b463299472fdf5b062316.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMzExNjA0_size_16_color_FFFFFF_t_70 1]: /images/20220418/bd5b859e10da4359808a39882785bf0e.png
[20181110113321219.png]: /images/20220418/ef7f1c69f0ac4644963346ffcc822ac6.png
[20181110113600845.png]: /images/20220418/5fc1dee59ee249139111f7f245bbd781.png
[20181110113843737.png]: /images/20220418/87cc2e204a8144589c1613a4ae6e6692.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMzExNjA0_size_16_color_FFFFFF_t_70 2]: /images/20220418/a856e220a2074ece8a73ff5a4ea7e2be.png
[20181110155122413.png]: /images/20220418/59ea6f8eb69f482fa79df7c430a757d3.png
[20181110155548635.png]: /images/20220418/918032da63c244b19d6d60405c5dc9fc.png
[20181110155856569.png]: /images/20220418/bfe855e67d3b474b9652eb2a3fda4b6b.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzQxMzExNjA0_size_16_color_FFFFFF_t_70 3]: /images/20220418/c566e07770464c419e52389bc10dbe8d.png