f(t)?g(t)cos(100t)，在接收端又被解调为g0(t)?f(t)cos(100t)，并通过低
通滤波器
?1H(j?)???0??30其他
得到信号g1(t)。试用MATLAB绘制上述各个信号的时域波形和频谱。
程序、实验结果及解释说明：
clear;
[t,omg,FT,IFT] = prefourier([-5,5],length(-5:0.01:5),[-250,250],1000); g=3*cos(10*t)+2*cos(20*t); f=g.*cos(100*t); g0=f.*cos(100*t); G0=FT*g0;
H=(omg<30&omg>-30); G1=G0.*H; g1=IFT*G1; G=FT*g; F=FT*f;
subplot(4,2,1)
plot(t,g);ylabel('g'); subplot(4,2,2)
plot(t,f);ylabel('f'); subplot(4,2,3)
plot(t,g0);ylabel('g0'); subplot(4,2,4)
plot(t,g1);ylabel('g1'); subplot(4,2,5)
plot(omg,abs(G)),ylabel('G'); subplot(4,2,6)
plot(omg,abs(F)),ylabel('F'); subplot(4,2,7)
plot(omg,abs(G0)),ylabel('G0'); subplot(4,2,8)
plot(omg,abs(G1)),ylabel('G1');
function [t,omg,FT,IFT]=prefourier(Trg,N,OMGrg,K) %Trg 为时域起止范围 %N 为时域抽样点数
%OMGrg 为频域起止范围 %K 为频域抽样点数 %t 为抽样时间点 %omg 为抽样频率点
%FT 为傅里叶变换矩阵（将该矩阵左乘信号即得信号的傅里叶变换）
10
%IFT 为傅里叶反变换矩阵（将该矩阵左乘信号即得信号的傅里叶反变换） T=Trg(2)-Trg(1);
t=linspace(Trg(1),Trg(2)-T/N,N)'; OMG=OMGrg(2)-OMGrg(1);
omg=linspace(OMGrg(1),OMGrg(2)-OMG/K,K)'; FT=T/N*exp(-j*kron(omg,t'));
IFT=OMG/2/pi/K*exp(j*kron(t,omg'));
11
12
13
实验四 连续时间LTI系统的复频域分析
8. 已知某二阶系统的零极点分别为p1??100，p2??200，z1?z2?0（二重零点），试用MATLAB绘出该系统在0 ~ 1kHz频率范围内的幅频特性曲线（要求用频率响应的几何求解方法实现），并说明该系统的作用。
程序、实验结果及解释说明：
clear;
f=0:0.01:1000; w=2*pi*f; p1=-100; p2=-200;
A1=abs(j*w-p1); A2=abs(j*w-p2); B1=abs(j*w);
H=(B1.^2)./(A1.*A2); plot(w,H)
xlabel('频率（rad/s）') ylabel('幅频响应')
14