《数字信号处理》实验报告 03051235 陆义良
实验一 傅立叶级数展开法生成 x(t)
一、Matlab 源代码
/*work1.m*/
%傅立叶级数展开合成方波
t=0:0.01:2*pi;
f1=2/pi*sin(t); % 基波
f3=2/pi*(sin(3*t)/3); %三次谐波
f5=2/pi*(sin(5*t)/5);
f7=2/pi*(sin(7*t)/7);
y1=f1+f3;
y2=f1+f3+f5;
y3=f1+f3+f5+f7;
subplot(2,2,1) %在第一个子窗口画 基波 分量
plot(t,f1),hold on
y=1*sign(pi-t); %画方波信号
plot(t,y,'c:')
title('方波的合成-基波')
subplot(2,2,2) %在第二个子窗口画 (基波+3 次谐波) 分量
plot(t,y1),hold on
y=1*sign(pi-t);
plot(t,y,'c:')
title('方波的合成-基波+3 次谐波')
subplot(2,2,3) %在第三个子窗口画 (基波+3 次谐波+ 5 次谐波) 分量
plot(t,y2),hold on
y=1*sign(pi-t);
plot(t,y,'c:')
title('方波的合成-基波+3 次谐波+5 次谐波')
subplot(2,2,4) %第四个子窗口画 (基波+3 次谐波+5 次谐波+7 次谐波) 分量
plot(t,y3),hold on
y=1*sign(pi-t);
plot(t,y,'c:')
title('方波的合成-基波+3 次谐波+5 次谐波+7 次谐波')
二、实验结果
1
评论4