方波的合成与分解.docxVIP

一、方波的合成与分解

1、傅立叶级数分析的原理:

任何周期信号都可以用一组三角函数｛sin(n3ot),cos(nw°t)｝的组合表示：

x(t)=￥aejn00=￥a[cos(n①t)+jsin(n①t)]n n 0 0

n=一3 n=一3

这表明傅立叶级数可以表示为连续时间的周期信号，也即是连续时间周期信号可以分解

2兀

为无数多个复指数谐波分量。在这里。孔为傅立叶级数的系数，o0=—称为基波频率。

2、建立方波信号的模型：

思考：如何建立连续周期方波信号？①预置一个周期内的方波信号：

--A(-T/2t0)

一个完整周期内的信号表达式：x0(t)日

La (0tT/2)

②对方波信号以周期T进行平移：x(t)=党x(t-nT)

0

n=s

通过以上的两个步骤我们可以建立一个连续周期方波信号，为降低方波信号分解与合成的复杂程度，可以预置方波信号为奇谐信号，此连续时间周期方波信号如下：

x=sign(冗一t),te(0,2兀)

x(t)=芝x(t一2兀n),te(-8,8)

0

n=-8

3、方波信号分解：

1.510.50

1.5

1

0.5

0

-0.5

T=

dx

II1111111111111111

「

1111

11

11

I11111111

T T1 1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

「 「

1 1

1 1

1 1

1 1

1 1

1 1

1 1

1 1

「

1

1

1

1

1

1

1

1

r t ■

■

■

■

■

■

■

■

■

T T -

i i

i i

i i

i i

i i

i i

i i

i i

1 1

1 1

2冗，①

:, 0

=?1

1

1

11111111

1 1

iiiiiiiiiiiiiiiiii

dt

jk^①易—

p]Figure1

-口，回]S3.

FileEditViewInsertToolsDesktop

Window

Help

nsg捻

□

B

■Q

1

[1-(-1)k]

jk0

jk

0,

为偶数

ejkt=24sin(kt),k为奇数

k=一3乙k丸

k=一3

k=1

由以上可知道，此方波信号可以分解为各奇次谐波。

方波的各奇次谐波可由如下Matlab语句得到:

t=0:0.001:2*pi;

w=1;

N=[1;3;5;7;9;11];

k=1;

whilek=6n=N(k,:); %基次谐波的次数

b=4./(pi*n);x=b*sin(w*n*t);subplot(3,2,k),plot(t,x);xlabelt),ylabel(x(t));axis([0,2*pi,-2,2]),legend(N次谐波])k=k+1;

end

输出图像如下:

4、方波信号合成：

对连续周期方波信号各谐波分量(基波分量、三次波分量、五次波分量??…)分别进行求和运算，步骤如下：

考查一个完整周期(0?2n)这段时间内的信号，画出结果，并显示。

画出基波分量匕，并显示，观察匕与原周期方波信号的误差大小。

将三次谐波f(t)加到第二步之上，画出结果，并显示，观察J与原周期方波信号的

3

误差大小。

J(t)=f(t)+f(t)

TOC\o1-5\h\z1 3

将五次谐波f5(t)加到第三步之上，画出结果，并显示，观察七与原周期方波信号的误差大小。

J(t)=f(t)+f(t)+f(t)

13 5

将七次谐波f(t)与九次谐波f(t)加到第四步之上，画出结果，并显示，观察J与

7 9 3

原周期方波信号的误差大小。

可由如Tmatlab程序进行合成：

function[P]=hecheng(N)

ifnargin==0,N=1;end

t=0:0.001:2火pi;

w=1;

k=1:2:N;

b=4./(pi*k);

x=b*sin(w*k*t);

plot(t,x);

xlabel(t),ylabel(x(t));

axis([0,2火pi,-2,2]),legend(第叫次合成);gridon

输出图像如下：

当n继续变大时，分别取n27,81,201,得到对应合成图像如下: