数学建模学习交流的博客

MATLAB文本数据处理入门篇的课后习题，MATLAB小白到高手的必刷题

木棍随机的折成三段后的概率问题

判断整数是否为回文数

投骰子经过100格的概率

蒙特卡罗模拟求自然常数e

每隔0.2秒在屏幕上随机输出0或者1，当连续出现三次1时，停止输出。

排队论问题的模拟

电子管的更换策略问题

并行计算的parfor循环

违背直觉的三门问题

计算动物能存活的天数的期望值

计算中国式排名

利用蒙特卡罗模拟来计算圆周率π

估计离开家之前能拿到报纸的概率

层次分析法中出现的判断矩阵

randperm函数的练习题

随机请同学吃饭、唱歌和看电影

计算银行存钱的本息

拉马努金的恒等式

怎样将一个方阵A的主对角线元素重新赋值为0

判断A和B是否互为对方的逆矩阵

找到向量中最大值对应的位置的索引

命令加载数据时，可以在工作区看到s3这个字符向量，s3中包含8320个字符，里面包含了西游记一百回的内容梗概。东胜神洲有一花果山，山项上一仙石孕育出一石猴。石猴在所居为“水帘洞”，被群猴拥戴为王。又过三五百年，根据一老猴指点，被菩提祖师收为徒，名曰孙悟空。《西游记》是中国古典四大名著之一，共分为一百回。每一回不仅有独特的标题，还包含了丰富的故事内容。​讲解视频：可以在bilibili搜索《

小王子》（The Little Prince）是法国作家安托万·德·圣-埃克苏佩里（Antoine de Saint-Exupery）于1942年创作的著名儿童文学短篇小说。本题需要大家对《小王子》这本书进行词频分析，即计算书中每个单词出现的频次。​讲解视频：可以在bilibili搜索《

首先观察每一行数据的结构，可以看出每一行要提取的指标由空格隔开，因此我们可以考虑使用split函数在空格处对每一行进行拆分，拆分后每一行变成了多个子字符串；最后，我们会把所有提取出来的数值指标整合到一个8行4列的矩阵M中，每行对应一个催化剂组合的四个指标，从而达成我们的分析目标。现在我们需要从s2的每一行中提取上述介绍的四个指标，并将提取的结果保存到一个8行4列的数值矩阵M中。我们的任务就是从提供的催化剂组合中提取关键的数值指标。​讲解视频：可以在bilibili搜索《乙醇浓度2.1ml/min。

下表左侧给出了四名同学在三次测试中的成绩数据，请对所有同学的测试成绩进行排名，并生成一份详细的排名通知。排名通知是一个带有换行符的字符串标量，每一行包括学生的姓名、第几次测试、具体成绩以及在所有成绩中的排名，具体可参考下表右侧。讲解视频：可以在bilibili搜索“

思路：在第三章的课后习题中，我们见到过类似的题目，当时我们介绍过蒙特卡罗模拟这种方法，蒙特卡罗模拟将所求解的问题同一定的概率模型相联系，用计算机实现统计模拟或抽样来获得问题的近似解。我们可以模拟这个过程N次（N一般要设置的大一点，例如让N等于10万），每次模拟中我们都让一只猫走出山洞，并记录下这只猫所需的时间。一只失明的小猫掉进山洞里，山洞有三个门，其中进入第一个门后走2h后可以回到地面，进入第二个门后走4h会回到原始的出发点，进入第三个门后走6h还是回到原始的出发点。

讲解视频：可以在bilibili搜索“斐波那契数列的递推公式为。

一个三位正整数各位数字的立方和等于该数本身则称该数为水仙花数，例如：1^3 + 5^3+ 3^3 = 153，则153是水仙花数。请你找出所有的水仙花数并将其保存到向量S中。讲解视频：可以在bilibili搜索“

讲解视频：可以在bilibili搜索“

（1）前面我们学习了计算对数的一组函数，分别是log, log2和log10，那么我们应该怎样计算以任意正数为底的对数呢？如果有多的糖，就带回家给媳妇吃。第三步：将该脚本保存到“第二章作业”的文件夹中，命名为“homework.m”。（2）请计算sin(2π )的值，如果你得到的结果不是0，应该怎么给别人解释？第二步：打开MATLAB，新建一个脚本文件，在里面任意输入一些你学过的代码。第四步：更改MATLAB的当前文件夹为“第二章作业”的文件夹。第五步：执行你写的脚本，将计算结果作为注释放到每一行的后面。

请将结果保存到元胞数组c中，c有三列：第一列为各地区的名称，第二列为各地区获得本科组一等奖的队伍数量，第三列为相应的频率（即各地区获得一等奖的队伍数量除以获得本科组一等奖的所有队伍数量299）。此外，请根据获奖数量将元胞数组c按照降序排列，即获奖数量多的地区应排在前面，并将排序后的结果保存在元胞数组cc中。下方给出的代码中用到了一个新的内置函数：tabulate函数，它能实现这个要求。我们的目标是统计获得本科组一等奖的299支队伍中，各地区的数量和出现的频率。（变量s3的下载方法可以看上方的讲解视频）

请你求出任意两名同学共有的兴趣爱好数量，并将结果保存到大小为40行40列的矩阵num中，其中num(ii, jj)表示第ii名同学和第jj名同学的共同爱好数量。（注意：统计不同同学的兴趣爱好时，可以忽略“其他”这个选项，因为它可能涵盖的兴趣范围过于广泛）s2中每个数据都是字符向量，包含由中文顿号分隔的多个兴趣爱好，部分同学的兴趣列表中包含了“其他”作为一个额外的选项。思考：在上面的代码中，每次比较两名同学的兴趣爱好时，都会重复执行strsplit函数来获取兴趣列表。讲解视频：可以在bilibili搜索“

我们使用mat2cell(sd,[3,3,3],[3,3,3])将数独方阵分割成九个宫格，每个宫格都是3×3的方阵，并存储在元胞数组cc中。sort(sd,1)表示沿着行方向对sd的每一列进行排序，如果每一列排序后都是1到9的序列，sort(sd,1) == (1:9)'就会返回一个9行9列全为逻辑值1的方阵，此时使用两次all函数可以得到Condition_1为true。这类似于列的验证，但方向不同。数独是一个9×9的方阵，它由九个宫格构成，每个宫格又由九个小格子构成(图中用颜色区分开的3×3的方阵)。

在上面这段代码中，我们开始时创建一个空的元胞数组SS，它被用来存储那些等级为'A'的学生的信息。这种语法大家可能是首次见到，它是一个非常实用的方法来动态地增加数组的大小，使用end+1可以确保新数据总是被添加到数组的末尾。其中，S的第一列是他们的姓名，假设他们的姓名由4个随机英文字母生成，首字母大写；S的第二列是他们语数外三门科目的成绩，你可以使用长度为3的数值向量保存成绩，假设成绩是区间[0,100]上的随机整数；进一步地，请从元胞数组S中筛选出等级为'A'的学生，并将结果保存到元胞数组SS中。

最后，我们可以输出模拟的k值，表示玩到第k把时手上有炸弹。另外，本题还能继续扩展下去，例如重复上面的模拟过程N次（N可以设置得大一点，例如N等于10万），得到这N次模拟结果的k，并计算这N次k的平均值，这个平均值就能表示你作为农民首次出现炸弹所需的期望轮数。在本题中，既用到了while循环又用到了for循环，且出现了两个不同用途的break，大家课后一定要认真消化，并尝试自己求解这个例题（当然，判断是否存在普通的炸可以不用循环语句，我们在第三章的课后习题中有讲解，详情请看第三章课后习题挑战篇的Q5）。

质数（Prime number），又称素数，指在大于1的自然数中，除了1和该数自身外，无法被其他自然数整除的数（也可定义为只有1与该数本身两个正因数的数）。给定任意一个大于100的自然数n（例如n=135389），请判断n是否为质数。还有下面这个内置函数：p = primes(n) 返回包含所有小于或等于 n 的质数的行向量。思路1：我们可以遍历从2到n-1的所有整数，检查它们是否能够整除n。思路2：缩小搜索的范围，只需要考虑奇数，且搜索的上界调整为ceil(sqrt(n))，即根号n的向上取整。

在每次循环中，如果找到回文数，我们会将构造好的字符向量 tmp 添加到 res 中，左侧代码使用的是已经过时的strvcat函数，右侧代码使用的是char函数。在循环体中：使用cc(ii,:)获取每一行表示的字符向量，将其保存到变量c_ii中，接下来我们使用了deblank函数，这个函数可以删除c_ii末尾的空白字符，最后我们使用disp函数输出处理后的字符向量，这样输出的每一行都不会有末尾的空格。（提示：转换后的二进制数为字符向量类型，将这个字符向量反转顺序，如果和原来的一样就是回文数。

这个字符矩阵有3行、每一行有5个字符（第一列的数字后面跟有一定数量的空格，保证每一行的字符数量相同），而我们这一小节介绍的其他例子要么只有1行（行字符向量），要么只有1列（列字符向量）。在第三章我们介绍过：向量和矩阵都属于数组，只不过向量的维度是一维的，矩阵的维度是二维的。这里的字符向量和字符矩阵也是类似的，字符向量也可以看成字符矩阵的特例，我们下一小节会详细介绍字符矩阵。构建二进制字符向量：在每次循环中，将计算得到的余数使用num2str函数转换为字符，并添加到二进制字符向量c的最前面。

如果是，在对应的位置为逻辑值1，否则为逻辑值0，因此这里的ind 等于[1 0 0 1 0 0 0 0 1 1 0 1 0 0 0 0 1 1 0 0 0 1 0 1 0 0 0 0 1 1 0 0 0 1 1 0 0 0 0 0 1 0 0 0 1 0 0 0 0 0 1 0 0 1 0 0 0] （注意：空格、标点符号等都要被判断哦！使用sum(ind)计算逻辑向量ind中逻辑值1的总数，也就是元音字母在字符向量c中出现的次数，将这个值存储在变量num中，得到num等于17；