阶跃信号和脉冲信号模拟

阶跃信号模拟

阶跃信号是一种理想的信号模型，通常用于描述系统在某一时刻发生突变或跃变的情况。在数学上，阶跃信号通常表示为一个在某一时刻突然从一个值跃变到另一个值的函数。阶跃信号在信号处理、控制系统等领域具有重要的应用，因为它能够很好地描述系统在某一时刻的快速变化。在实际系统中，很多物理过程都可以近似地用阶跃信号来描述，因此阶跃信号是一个非常重要的信号模型。

以下是Python生成阶跃信号的代码：

import numpy as np
import matplotlib.pyplot as plt

# 定义阶跃函数
def step_function(x):
    return np.where(x >= 0, 1, 0)

# 生成x值
x = np.linspace(-10, 10, 400)

# 计算对应的y值
y = step_function(x)

# 绘制图像
plt.figure(figsize=(8, 4))
plt.plot(x, y, label="Step Function", color='blue')
plt.ylim(-0.1, 1.1)  # 设置y轴的范围
plt.xlabel('x')
plt.ylabel('y')
plt.title('Step Function Graph')
plt.axhline(0, color='black',linewidth=0.5)
plt.axvline(0, color='black',linewidth=0.5)
plt.grid(color = 'gray', linestyle = '--', linewidth = 0.5)
plt.legend()
plt.show()

这里，`condition`是`x >= 0`，这意味着对于输入数组`x`中的每个元素，都会检查这个元素是否大于或等于0。

- 如果某个元素满足条件（即该元素的值大于或等于0），则从`x`（在这个例子中是数字1）中选择元素，因此对应位置的输出将是1。
- 如果某个元素不满足条件（即该元素的值小于0），则从`y`（在这个例子中是数字0）中选择元素，因此对应位置的输出将是0。

这样，`np.where`函数将输入数组`x`转换成一个新的数组，其中所有非负元素都被替换为1，所有负元素都被替换为0。这正是阶跃函数的定义：当输入大于或等于0时，输出为1；否则，输出为0。

 

以下是MATLAB生成阶跃信号的代码：

t = -10:0.1:10;
u = heaviside(t);
plot(t, u);
xlabel('Time (s)');
ylabel('Amplitude');
title('Heaviside Step Function');
grid on;

 

 

脉冲信号模拟

脉冲信号是一种短暂的、非连续变化的信号，通常具有很高的幅度并且只在非常短的时间内存在。在工程和物理领域，脉冲信号经常用于各种测量和通信系统中，因为它们可以携带关于系统性能的重要信息。 

以下是Python生成脉冲信号的代码：

import numpy as np
import matplotlib.pyplot as plt

def pulse_signal(x, pulse_center, pulse_width):
    """
    生成脉冲信号
    :param x: 输入的x值数组
    :param pulse_center: 脉冲中心的位置
    :param pulse_width: 脉冲的宽度
    :return: 对应的脉冲信号数组
    """
    return np.where((x >= pulse_center - pulse_width / 2) & (x <= pulse_center + pulse_width / 2), 1, 0)

# 生成x值
x = np.linspace(-10, 10, 1000)

# 脉冲信号的参数
pulse_center = 0  # 脉冲中心
pulse_width = 2   # 脉冲宽度

# 计算对应的脉冲信号
y = pulse_signal(x, pulse_center, pulse_width)

# 绘制脉冲信号图像
plt.figure(figsize=(10, 4))
plt.plot(x, y, label="Pulse Signal", color='red')
plt.ylim(-0.1, 1.1)  # 设置y轴的范围
plt.xlabel('x')
plt.ylabel('Amplitude')
plt.title('Pulse Signal Graph')
plt.axhline(0, color='black',linewidth=0.5)
plt.axvline(0, color='black',linewidth=0.5)
plt.grid(color = 'gray', linestyle = '--', linewidth = 0.5)
plt.legend()
plt.show()

 

以下是MATLAB生成阶跃信号的代码：

t = -5:0.001:5; % 时间向量，从-5到5，步长为0.001
width = 0.1;    % 脉冲宽度
pulse = rectpuls(t, width); % 生成脉冲信号
plot(t, pulse); % 绘制脉冲信号
xlabel('Time (s)');
ylabel('Amplitude');
title('Rectangular Pulse Signal');
grid on;

 

 

广义脉冲信号模拟

还可以通过高斯钟形函数、双边指数函数、取样函数模拟脉冲函数，以下是模拟代码：

1.高斯钟形函数

高斯脉冲的特点是在中心附近迅速上升并下降，没有硬边界，即没有严格的起始和结束时刻，这使得它在信号处理中具有良好的抗干扰能力。此外，高斯脉冲是完全连续的，没有任何不连续点。

import numpy as np
import matplotlib.pyplot as plt

def delta_function(t, b):
    """
    近似狄拉克δ函数
    :param t: 时间变量
    :param b: 控制尖峰宽度和高度的参数，b越大，尖峰越窄且越高
    :return: δ函数的近似值
    """
    return b * np.exp(-np.pi * (b * t) ** 2)

# 生成t值
t = np.linspace(-0.1, 0.1, 1000)

# 选择一个较大的b值
b_values = [50, 100, 200, 400]

# 绘制图像
plt.figure(figsize=(10, 6))
for b in b_values:
    plt.plot(t, delta_function(t, b), label=f'b={b}')

plt.xlabel('t')
plt.ylabel('Amplitude')
plt.title('Approximation of Dirac Delta Function')
plt.legend()
plt.grid(True)
plt.show()

 

2.双边指数函数

双边指数函数的特点是有明确的起点和终点，形成一个尖锐的峰值。这种形状类似于实际电子设备中的冲激响应，因此常用于模拟理想冲激信号。

 

import numpy as np
import matplotlib.pyplot as plt

def delta_function_approx(t, b):
    """
    近似狄拉克δ函数的另一种表述
    :param t: 时间变量
    :param b: 控制尖峰宽度和高度的参数
    :return: δ函数的近似值
    """
    return (1 / (2 * b)) * np.exp(-np.abs(t) / b)

# 生成t值
t = np.linspace(-0.05, 0.05, 1000)

# 选择几个不同的b值
b_values = [0.01, 0.005, 0.0025]

# 绘制图像
plt.figure(figsize=(10, 6))
for b in b_values:
    plt.plot(t, delta_function_approx(t, b), label=f'b={b}')

plt.xlabel('t')
plt.ylabel('Amplitude')
plt.title('Approximation of Dirac Delta Function')
plt.legend()
plt.grid(True)
plt.show()

 

3.取样函数

取样函数是一个理论上的脉冲序列，通常用来表示离散时间信号。它的表达式可以是一个简单的单位阶跃函数的离散序列，或者是更一般的形式，如狄拉克δ函数的离散序列。取样函数的特点是在离散的时间点上有无限大的振幅，而在其他时间点上振幅为零。这种函数通常用于数字信号处理中，因为实际的数字系统只能处理离散的时间点。

import numpy as np
import matplotlib.pyplot as plt

def sinc_approx(t, b):
    """
    近似狄拉克δ函数的sinc函数形式
    :param t: 时间变量
    :param b: 控制函数形状的参数
    :return: δ函数的近似值
    """
    # 避免除以0的情况
    with np.errstate(divide='ignore', invalid='ignore'):
        result = np.sin(b * t) / (np.pi * t)
        # 当t=0时，手动设置结果为b/pi，这是利用洛必达法则或直接代入sin(bt)/(pi*t)的极限得到的
        result[t == 0] = b / np.pi
    return result

# 生成t值
t = np.linspace(-0.05, 0.05, 1000)

# 选择一个较大的b值
b = 10000  # 使用一个较大的b值来近似δ函数

# 绘制图像
plt.figure(figsize=(10, 6))
plt.plot(t, sinc_approx(t, b), label=f'b={b}')

plt.xlabel('t')
plt.ylabel('Amplitude')
plt.title('Approximation of Dirac Delta Function using Sinc Function')
plt.legend()
plt.grid(True)
plt.show()