东华大学信号课程设计2015第三次.docVIP

实验内容 1. 绘制幅频响应和相频响应曲线： 代码: b=[-1 1];a=[1 1]; freqs(b,a); 图形： (2) 绘制激励e(t)波形： t=0:0.001:50; et=sin(t)+sin(3*t); plot(t,et); grid on; title(e(t)); xlabel(t); 求r(t): syms t; et=sin(t)+sin(3*t); E=fourier(et); syms w; H=(1-i*w)/(1+i*w); R=E*H; rt=ifourier(R,t) 结果：rt=-6/5*cos(t)*cos(2*t)-2/5*cos(t)-4/5*sin(3*t) 绘制响应r(t)的波形： t=0:0.001:50; rt=-6/5*cos(t).*cos(2*t)-2/5*cos(t)-4/5*sin(3*t); plot(t,rt); grid on; title(r(t)); xlabel(t); 2. (1)调制与解调过程： t1=0:0.1:20; syms t v; x=sin(t-10)/(t-10); subplot(321); ezplot(x,[0,20]); title(原函数); Fx=fourier(x,v); subplot(322) ezplot(Fx,[0,20]); axis([-20,20,-8,8]); title(频谱); y=sin(t1-10)./(t1-10).*cos(500*t1); subplot(323); plot(t1,y); title(调制后); y1=x*cos(500*t); Fy=fourier(y1,v); subplot(324); ezplot(Fy,[0 20]); axis([-20,20,-8,8]); title(频谱); z=y1*cos(500*t); Fz=fourier(z,v); G=-heaviside(v-1)+heaviside(v+1); Fx1=Fz*G; x1=ifourier(Fx1,v,t); subplot(325); ezplot(2*x1,[0 20]); title(解调后); subplot(326); ezplot(2*Fx1,[0 20]); axis([-20,20,-8,8]); title(频谱); 图形： (2) t1=0:0.1:20; syms t v; x=sin(t-10)/(t-10); subplot(321); ezplot(x,[0,20]); title(原函数); Fx=fourier(x,v); subplot(322) ezplot(Fx,[0,20]); axis([-20,20,-8,8]); title(频谱); y=sin(t1-10)./(t1-10).*cos(500*t1); subplot(323); plot(t1,y); title(调制后); y1=x*cos(500*t); Fy=fourier(y1,v); subplot(324); ezplot(Fy,[0 20]); axis([-20,20,-8,8]); title(频谱); z=y1*cos(499*t+12); Fz=fourier(z,v); G=-heaviside(v-1)+heaviside(v+1); Fx1=Fz*G; x1=ifourier(Fx1,v,t); subplot(325); ezplot(2*x1,[0 20]); title(解调后); subplot(326); ezplot(2*Fx1,[0 20]); axis([-20,20,-8,8]); title(频谱); 3. 采样： dt=0.005; t=-10:dt:10; ff=0.1; f=sin(t)./t; subplot(221); plot(t,f); xlabel(t); title(f(t)); Ts=1; ts=-10:Ts:10; fs=sin(ts)./ts; subplot(222); plot(t,f,:); hold on; stem(ts,fs); xlabel(t); title(f_s(t)); w=linspace(-10*pi,10*pi,4001); F=f.*exp(-j*t.*w)*dt; subplot(223); plot(w/pi,abs(F)); xlabel(\omega/\pi); title(|F(\omega)|); wc=linspace(-10*pi,10*p