这个函数的频谱图怎么画呀

这个函数的频谱图怎么画呀

问题描述:

img

                                     2A
                           4 十
                                    T。
                                                                       2
x(t)=
                                   2A
                                                                                            0
                                                                0≤ts
                                   To
                                                                                          2

你好,这是要进行傅里叶展开啊

clc;clear
syms T0 A t n
x1 = A+2*A/T0*t;
x2 = A-2*A/T0*t;
F = int(x1*cos(2*n*pi*t/T0), t, -T0/2, 0)+...
    int(x2*cos(2*n*pi*t/T0), t, 0, T0/2);
a0 = 2/T0*(int(x1,t,-T0/2, 0) + int(x2,t,0,T0/2))
an = 2/T0*F

得到

a0 =
A

an =

(4*A*sin((pi*n)/2)^2)/(n^2*pi^2)

我们令A=1;n从1到10取值,有

A = 1;
n = [1:10];
Amp = (4*A*sin((pi*n)/2).^2)./(n.^2*pi^2);
stem([0,n],[A,Amp])

频谱图

img

用matlab自动运行