周期信号的分解与合成.pptx

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
4
xlabel('n');
ylabel('Fn'); title('周期矩形脉冲的频谱图'); 3. %使用不同频率的谐波合成方波信号;注意观察随着谐波数的增加合成的波形 %发生的变化。
t=-1:0.001:1;
omega=2*pi;
y=square(2*pi*t,50);
plot(t,y);grid on
学海无 涯
实验一 周期信号的分解与合成
一、实验目的
1. 用同时分析法观测50Hz 非正弦周期信号的频谱。 2. 观测基波和其谐波的合成。
二、实验原理
1. 一个非正弦周期函数可以用一系列频率成整数倍的正弦函数来表示,其中与非正弦 具 有相同频率的成分称为基波或一次谐波,其它成分则根据其频率为基波频率的2、3、4、...、 n 等倍数分别称二次、三次、四次、...、n 次谐波,其幅度将随谐波次数的增加而减小, 直至无穷小。
t=-2:0.001:2;
omega=pi;
y=rectpuls(t)+rectpuls(t-2)+rectpuls(t+2);
axis([-2.2 2.2 0 1.2]);
plot(t,y);grid on
xlabel('t'); ylabel('周期矩形脉冲信号');
axis([-2.2 2.2 -1.5 1.5]);
五、思考题
1. 什么样的周期性函数没有直流分量和余弦项?
附: 1. 50HZ 方波信号的频谱。
>> w1=
; %基波角频率
>> n=0:1:30; >>bn=
; %三角级数中系数bn,参考书 p122
>> stem(n*w1,bn),grid on
>> xlabel('\omega(rad/s)'),ylabel('bn')
2.三角波
1
学海无 涯
3.半波 4.全波 5.矩形波
三、预习要求
在做实验前必须认真复习教材中关于周期性信号傅利叶级数分解的有关内容。
四、实验内容
1. 50HZ 方波信号的频谱。 2. 周期矩形脉冲的频谱;脉冲宽度为 1;周期为 4;则基波角频率为 0.5pi 3. 使用不同频率的谐波合成方波信号;注意观察随着谐波数的增加合成的波形发生的变化。 4. 使用不同频率的谐波合成矩形脉冲信号;注意观察随着谐波数的增加合成的波形。
>> title('方波信号频谱分析图') 2. 周期矩形脉冲的频谱;脉冲宽度为 1;周期为 4;则基波角频率为 0.5pi
tao=
;
w1=
;
n=-15:1:15;
2
学海无涯
fn=
; %矩形脉冲级数系数 fn,参考书 p130,用 matlab 自带函数 sinc
stem(n,fn),grid on
plot(t,y)
hold on;
plot(t,x);
hold off;
xlabel('t'); ylabel('部分和的波形');
axis([-1 1 -1.5 1.5]),grid on
title(['最大谐波数=',num2str(n_max(k))])
end 4. %使用不同频率的谐波合成矩形脉冲信号;注意观察随着谐波数的增加合成的波形 %发生的变化。
n_max=[1 3 5 11 47];
N=length(n_max);
for k=1:N
n=1:1:n_max(k);
来自百度文库
3
学海无 涯
a=2*sin(n*omega/2)./(n*pi); x=a*cos(omega*n'*t); x=x+1/2; figure; plot(t,y) hold on; plot(t,x); hold off; xlabel('t'); ylabel('部分和的波形'); axis([0 2 -1.5 1.5]),grid on title(['最大谐波数=',num2str(n_max(k))]) end
2. 不同频率的谐波可以合成一个非正弦周期波,反过来,一个非正弦周期波也可以分 解为无限个不同频率的谐波成分。 3. 一个非正弦周期函数可用傅里叶级数来表示,级数各项系数之间的关系可用一各个 频 谱来表示,不同的非正弦周期函数具有不同的频谱图,各种不同波形及其傅氏级数表达式 见 表1-1
表1-1 各种不同波形的傅里叶级数表达式(下) 1.方波
xlabel('t'); ylabel('周期方波信号');
axis([-1 1 -1.5 1.5]);
n_max=[1 3 5 11 47];
N=length(n_max);
for k=1:N
n=1:2:n_max(k); b=4./(pi*n);
x=b*sin(omega*n'*t);
figure;
相关文档
最新文档