以下是补充完整的程序:
clear all;
fs = 10000;
t = 0:1/fs:0.1;
f0 = 50;
sum = 0;
subplot(2,1,1)
for n = 1:2:9
fn = (1/n)*sin(2*pi*n*f0*t);
sum = sum + fn;
plot(t,fn)
hold on
end
plot(t,sum,'k','LineWidth',2)
xlabel('时间')
ylabel('振幅')
title('方波近似合成')
subplot(2,1,2)
plot(t,sum,'k','LineWidth',2)
xlabel('时间')
ylabel('振幅')
title('合成方波')