信号合成与分解-三角波、锯齿波程序
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
clc;
clear all;
close all;
n=7;
T0=2;A=2;
NofT0=2;
%周期信号的时域描述
tn_i=1;
for tn=0:0.01:NofT0*T0
if(rem (tn,T0)<=T0/4)
y_t(tn_i)=4*A*rem(tn,T0)/T0;
elseif (T0/4<=rem(tn,T0))&(rem(tn,T0)<=3*T0/4) y_t(tn_i)=-4*A*(rem(tn,T0)-T0/2)/T0;
elseif(3*T0/4<= rem (tn,T0)<=T0)
y_t(tn_i)=4*A*( rem (tn,T0)-T0)/T0;
end;
t_t(tn_i)=tn;
tn_i=tn_i+1;
end;
%周期信号的频域描述
t=0:0.01:NofT0*T0;
x=0;
pi=3.1415926;
w0=2*pi/T0;
for i=1:n
fw(i)=(2*i-1)*w0;
a(i)=(-1)^(i-1)*(8*A/(pi^2*(2*i-1)^2));
fai(i)=0;
y(i,:)=a(i)*sin(fw(i)*t);
x=x+y(i,:);
end;
subplot(2,2,1);
plot(t_t,[y_t;x]);
title('三角波、三角波合成图')
subplot(2,2,2);
plot(t,[x; y]);
title('0-n次谐波及合成图')
subplot(2,2,3);
stem(fw,a);
title('三角波幅频谱图')
subplot(2,2,4);
stem(fw,fai);
title('三角波相频谱图')
clc;
clear all;
close all;
n=7;
T0=2;A=2;
NofT0=2;
%周期信号的时域描述
tn_i=1;
for tn=0:0.01:NofT0*T0
y_t(tn_i)=A* rem (tn,T0)/T0;
t_t(tn_i)=tn;
tn_i=tn_i+1;
end;
%周期信号的频域描述
t=0:0.01:NofT0*T0;
x=A/2;
pi=3.1415926;
w0=2*pi/T0;
for i=1:n
fw(i)=i*w0;
a(i)=-A/(pi*i);
fai(i)=0;
y(i,:)=a(i)*sin(fw(i)*t);
x=x+y(i,:);
end;
subplot(2,2,1);
plot(t_t,[y_t;x]);
title('锯齿波、锯齿波合成图') subplot(2,2,2);
plot(t,[x; y]);
title('0-n次谐波及合成图') subplot(2,2,3);
stem(fw,a);
title('锯齿波幅频谱图')
subplot(2,2,4);
stem(fw,fai);
title('锯齿波相频谱图')。