??在之前的文章

Matlab精读复刻论文

子系统分为虚拟子系统和非虚拟子系统,非虚拟子系统又分为原子子系统和受外部信号触发、外部信号使能的触发子系统和使能子系统，另外还有跟流控制模块lf , Switch,For ，While模块配合使用的动作子系统(Action Subsystem),以及包含多个子系统在同一个时间仅能激活其中一个的可变子系统。

1.子系统

Matlab中的fillmissing函数可以用于填充缺失数据，可以根据相邻的数据进行插值或使用指定的填充方法来填充缺失值。下面是用C语言实现类似功能的代码示例：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

void fill_missing(double *data, int n, char method)
{
    int i, j;
    double x1, x2, y1, y2;
    double a, b;

    for (i = 0; i < n; i++) {
        if (isnan(data[i])) {
            if (method == 'l') {  // 线性插值
                j = i + 1;
                while (j < n && isnan(data[j])) {
                    j++;
                }
                if (j == n) {
                    data[i] = data[i - 1];
                } else {
                    x1 = i - 1;
                    x2 = j;
                    y1 = data[i - 1];
                    y2 = data[j];
                    a = (y2 - y1) / (x2 - x1);
                    b = y1 - a * x1;
                    data[i] = a * i + b;
                }
            } else if (method == 'm') {  // 平均值填充
                j = i - 1;
                while (j >= 0 && isnan(data[j])) {
                    j--;
                }
                x1 = j + 1;
                j = i + 1;
                while (j < n && isnan(data[j])) {
                    j++;
                }
                x2 = j - 1;
                data[i] = (data[x1] + data[x2]) / 2;
            } else if (method == 'p') {  // 前向填充
                j = i - 1;
                while (j >= 0 && isnan(data[j])) {
                    j--;
                }
                if (j < 0) {
                    data[i] = 0;
                } else {
                    data[i] = data[j];
                }
            } else if (method == 'b') {  // 后向填充
                j = i + 1;
                while (j < n && isnan(data[j])) {
                    j++;
                }
                if (j == n) {
                    data[i] = 0;
                } else {
                    data[i] = data[j];
                }
            }
        }
    }
}

int main()
{
    double data[] = {1, 2, 3, NAN, 5, 6, NAN, 8, NAN};
    int n = sizeof(data) / sizeof(double);
    int i;

    printf("Original data: ");
    for (i = 0; i < n; i++) {
        printf("%.1f ", data[i]);
    }
    printf("\n");

    fill_missing(data, n, 'l');  // 线性插值
    printf("Filled data (linear): ");
    for (i = 0; i < n; i++) {
        printf("%.1f ", data[i]);
    }
    printf("\n");

    fill_missing(data, n, 'm');  // 平均值填充
    printf("Filled data (mean): ");
    for (i = 0; i < n; i++) {
        printf("%.1f ", data[i]);
    }
    printf("\n");

    fill_missing(data, n, 'p');  // 前向填充
    printf("Filled data (forward): ");
    for (i = 0; i < n; i++) {
        printf("%.1f ", data[i]);
    }
    printf("\n");

    fill_missing(data, n, 'b');  // 后向填充
    printf("Filled data (backward): ");
    for (i = 0; i < n; i++) {
        printf("%.1f ", data[i]);
    }
    printf("\n");

    return 0;

因为之前用生成对抗网络及众多变体生成诸如心电信号，肌电信号，脑电信号，微震信号，机械振动信号，雷达信号等，但生成的信号在频谱或者时频谱上表现很差，所以暂时先不涉及到这些复杂信号，仅仅以手写数字图像为例进行说明，因为Python相关的资源太多了，我就不凑热闹了，使用的编程环境为MALAB R2021B。

我们生活中经常会问遇到手机拍照内存不够该怎么办？朋友圈经常看到各种运动轨迹 APP记录的跑步路线图，我们感到好玩，但是我们不禁还想问如何实现的呢？其实这都和手机中的传感器有关。手机是一个高度集成的电子产品，除了能通讯外，他还集中了很多的传感器，包括加速度、角速度、磁场、方向、位置。那么上面有关拍照上传或者读取传感器数据的问题都可以用matlab 软件实现，以下讲进行讲解。

该程序运行环境为MATLAB R2018A。

%%%arctanx=x-x^3/3+x^5/5-...

%%%pi/4=1-1/3+1/5-1/7+1/9-....

s=0; n=1;

ps=pi*0.25;