发布时间 : 星期二 文章长江大学信号与系统matlab实验答案更新完毕开始阅读ca1cf7e9f242336c1eb95eca
幅度响应108642000.10.20.30.50.6频率?(x? rad/sample)0.40.70.80.9相位响应0-20-40-6000.10.20.30.50.6频率?(x? rad/sample)0.40.70.80.9实验6周期信号的分解与合成
上机内容
? 用MATLAB画出图1的频谱和谐波次数N=20时合成的波形。观察N变化时合成波形和吉布斯现象。其中图1是一个周期的波形。
32.521.510.5001234图1
Eg6
代码如下：
function [An]=ZQXHFS(x,T0,T,N) %UNTITLED 计算连续周期信号的频谱
% x为周期信号的第一周期的字符串表达式，T0为第一周期的起始时间，T为周期，N为最高谐波次数
t=linspace(T0,T0+T,10000);w=2*pi/T; fm=eval([x]); for k=0:N
a(k+1)=2/T*trapz(t,fm.*cos(k*w*t));
b(k+1)=2/T*trapz(t,fm.*sin(k*w*t)); end
An=a-1i*b;An(1)=a(1)/2; f=abs(An);
A=angle(An)*180/pi; A(1)=A(1)/2;
fm_max=max(f);fm_min=min(f);
f_max=max(f);f_min=min(f);df=(f_max-f_min)/10; A_max=max(A);A_min=min(A);dA=(A_max-A_min)/10; t1=[t,t+T,t+2*T];ft=[fm,fm,fm];
subplot(3,1,1);plot(t1,ft,'LineWidth',2);ylabel('f(t)'); title('周期信号的波形','FontSize',8)
axis([t1(1) t1(end) 0 3]);grid,set(gca,'FontSize',8) subplot(3,1,2);h=stem(0:N,f,'.');
title('单边幅度频谱','FontSize',8);ylabel('An的模','FontSize',8);
set(h(2),'color','r','LineWidth',2),set(h(1),'color','r','LineWidth',2)
axis([0 N -0.2 1.5]);grid,set(gca,'FontSize',8)
subplot(3,1,3);h=stem(0:N,A,'.');title('单边相位频谱','FontSize',8) set(h(2),'color','r','LineWidth',2),set(h(1),'color','r','LineWidth',2)
grid,set(gca,'FontSize',8);ylabel('An的相位（度）');
w=2/T;xw=strcat('n\\Omega, \\Omega=',num2str(w),'\\pi');xlabel(xw); end
T=4;N=20;
y='t.*(heaviside(t)-heaviside(t-3))'; figure(1)
A_n=ZQXHFS(y,0,T,N); t=linspace(0,3*T,400); n=0:N;omega_0=2*pi/T; y1=A_n*exp(1i*omega_0*n'*t); yt=real(y1);
figure(2),plot(t,yt,'linewidth',2); title('周期信号合成（前20项的部分和）',8)
结果：
周期信号的波形3210024681012f(t)单边幅度频谱2An的模100246810单边相位频谱1214161820An的相位（度）2000-200024681012n?, ?=0.5?14161820
实验7滤波器频域分析与卷积分析
上机内容
滤波器为RC低通滤波器，输入信号为周期锯齿波如图1所示。当RC=0.5、2时，用MATLAB分析系统的输入频谱、输出频谱以及系统的时域响应。再用卷积分析的方法画出滤波器输出波形。
Eg7
部分代码如下（不会做，没做完） t=linspace(0,4,400); w=linspace(-15,15,30);
f=sym('t*(Heaviside(t)-Heaviside(t-3))') F=fourier(f), F=simple(F) f1=subs(f); Fv=subs(F);
F1=abs(Fv); P1=angle(Fv)*180/pi; subplot(3,1,1),plot(t,f1,'linewidth',2); grid;ylabel('f(t)');
subplot(3,1,2),stem(w,F1,'.'); grid;ylabel('|F(j\\omega)|');
结果：