信号与系统实验3 连续信号的频域分析

问题1 求实验图所示的周期信号(T=2 =1)的傅里叶级数,做出其前3、9、21、45项谐波的合成波形并与原信号做比较,作出其单边幅度谱和相位谱

t = -6:0.01:6;
d = -6:2:6;
fxx = pulstran(t,d,'tripuls'); %三角脉冲
f1 = fourierseries(3,t);    %fourierseries 自定义函数
f2 = fourierseries(9,t);
f3 = fourierseries(21,t);
f4 = fourierseries(45,t);
N = length(fxx);
f5 = fft(fxx,N);
f = (0:N-1)'*100/N;
f6 = angle(f5);
subplot(3,2,1);plot(t,fxx,'r',t,f1,'b');xlabel('t');ylabel('f(t)');title('前3项谐波');
grid on
axis([-6 6 -0.1 1.1]);
subplot(3,2,2);plot(t,fxx,'r',t,f2,'b');xlabel('t');ylabel('f(t)');title('前9项谐波');
grid on
axis([-6 6 -0.1 1.1]);
subplot(3,2,3);plot(t,fxx,'r',t,f3,'b');xlabel('t');ylabel('f(t)');title('前21项谐波');
grid on
axis([-6 6 -0.1 1.1]);
subplot(3,2,4);plot(t,fxx,'r',t,f4,'b');xlabel('t');ylabel('f(t)');title('前45项谐波');
grid on
axis([-6 6 -0.1 1.1]);
subplot(3,2,5);stem(f(1:600),abs(f5(1:600)),'.');xlabel('w');ylabel('F(jw)');title('单边频谱图');
grid on
subplot(3,2,6);stem(f6,'.');xlabel('w');ylabel('f6');title('相位谱');
grid on
axis([0 1200 -200 200]);
print(gcf,'-dpng','3-1.png');

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值