方波信号产生公式

方波信号产生公式

    int step=100;
    int k=1;

    for(int i = 0; i < DATASIZE; i ++)
    {
            t_int_data=10;
            if((i>(k*step -step))&&(i<k*step))
            {
               t_int_data=0;
            }

            if(i>k*step)
            {
                k +=2;
            }

}

波形图:

在这里插入图片描述

  • 0
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
方波信号的合成与分解是一个非常经典的问题,下面给出Python代码实现。 合成方波信号:假设我们要合成一个频率为$f$的方波信号,其基波为$A_0\sin(2\pi ft)$,而其谐波为$n=3,5,7,\cdots$的正弦波,其幅值为$\frac{A_0}{n}$(即基波幅值的$\frac{1}{n}$倍),相位为$0$。则合成的方波信号为: $$ x(t)=\frac{4A_0}{\pi}\sum_{n=0}^{\infty}\frac{\sin((2n+1)2\pi ft)}{2n+1} $$ 下面给出Python代码实现: ```python import numpy as np import matplotlib.pyplot as plt # 生成方波信号 def square_wave(f, T, A=1): t = np.linspace(0, T, 1000) x = np.zeros_like(t) for n in range(1, 100, 2): x += (4*A/(n*np.pi))*np.sin(2*np.pi*f*n*t) return x # 绘制方波信号 f = 10 T = 1/f x = square_wave(f, T) plt.plot(x) plt.show() ``` 分解方波信号:假设我们已经获得了一个方波信号$x(t)$,我们要将其分解成基波和谐波的叠加形式。假设方波信号的周期为$T_0$,则基波的频率为$f_0=\frac{1}{T_0}$,基波幅值为$x_0=\frac{1}{T_0}\int_0^{T_0}x(t)dt$,而谐波的频率为$n\cdot f_0$($n=2,3,\cdots$),谐波幅值为$x_n=\frac{2}{T_0}\int_0^{T_0}x(t)\sin(2\pi n\frac{t}{T_0})dt$。则分解后的方波信号为: $$ x(t)=x_0+\sum_{n=2}^{\infty}x_n\sin(2\pi n\frac{t}{T_0}) $$ 下面给出Python代码实现: ```python import numpy as np import matplotlib.pyplot as plt # 读取方波信号 x = np.loadtxt('square_wave.txt') # 分解方波信号 T0 = 100 f0 = 1/T0 x0 = (1/T0)*np.trapz(x, dx=T0/len(x)) xn = [] for n in range(2, 100, 2): xn.append((2/T0)*np.trapz(x*np.sin(2*np.pi*n*np.arange(len(x))*T0/len(x)), dx=T0/len(x))) xn = np.array(xn) # 合成方波信号 t = np.linspace(0, T0, len(x)) y = x0 + np.sum(np.outer(xn, np.sin(2*np.pi*np.arange(2, 100, 2)*t/T0)), axis=0) # 绘制原始方波信号和分解/合成后的信号 plt.subplot(211) plt.plot(x) plt.title('Original Square Wave') plt.subplot(212) plt.plot(y) plt.title('Decomposed and Composed Square Wave') plt.show() ``` 注意:在分解方波信号时,我们使用了NumPy中的`trapz`函数来进行积分,其实现方式与传统的复合梯形公式相同。此外,我们还使用了NumPy中的`outer`函数来实现外积的计算。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

W001123456789

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值