用matlab的m文件 求小于10^8的雷劈数，求代码。

以下是一份MATLAB程序，使用筛法求得小于10^8的雷打数：

```matlab

n = 1e8% 设置筛法范围

isPrime = true(n, 1)% 初始化数组为全部为true

isPrime(1) = false% 定义1不是素数

% 首先从2开始，将所有的素数的倍数标记为非素数

for i = 2:sqrt(n)

if isPrime(i)

isPrime(i * i : i : n) = false

end

end

% 输出小于10^8的雷打数

for i = 1:length(isPrime)

if mod(i, 4) == 1 &&isPrime(i)

fprintf('%d\n', i)

end

end

```

在程序中，首先定义了一个筛法范围 `n = 1e8`，然后利用逻辑数组 `isPrime` 来记录每个数是否为素数。初始化时，所有的数默认为素数，将数组全部初始化为 `true`；同时，定义数字1不是素数。之后，从2开始循环，将所有的素数的倍数标记为非素数。具体实现的方法是，如果检测到一个素数i，则将i * i及其之后的所有i的倍数都标记为非素数。

最后，再次循环判断并输出所有小于10^8的雷打数。因为所有的雷打数的形式为4k+1，所以对于满足此条件的素数，都可以视为雷打数。程序使用 `mod(i, 4) == 1` 检查当前数是否满足4k+1的形式，并使用 `isPrime(i)` 检查该数是否为素数。如果两个条件都满足，则输出该数即可。

需要注意的是，在程序中使用了 `sqrt(n)` 函数，以降低时间复杂度。因为n以内的素数最大可能为 `sqrt(n)`，所以只需要进行到 `sqrt(n)` 的筛法，即可得到全部小于n的素数。

1、首先，打开matlab，如图所示将如下所示的程序敲入matlab命令行窗口。

2、然后回车就可以得出matlab软件画出的阶跃函数波形，如图所示。

3、对于这个通过matlab画出来的波形图片的编辑，及其保存，都是可以通过d出窗口的菜单工具来实现的，如图所示。

4、plot(t,y)是以t向量为X轴，u向量为Y轴绘制波形的；axis([-0.5 1 -0.2 1.2])是设置-0.5,1为X轴的最大和最小值，同样的-0.2,1,2为Y轴的最大值和最小值。

---