i,a(i)=input(振幅a()F方 wM)=-input(角频率w()=方 end %(2)将各个频谱向量相加合成并画图 %此处应该把各时刻的图形转为动画，此处省略了动画的语句) t=0:0.1:20:It=length(t): %给出时间数组 p=a'*ones(1,It).*exp(j*w*t); %各频谱分量随时间变化的复 数值 q=cumsum(p): %各频谱分量的累加（包括所有节 点) figure(1),plot(real(q(4,:)》,imag(q(4,)),grid on%画合成复信号 的端点轨迹 %(3)将此轨迹的在实轴与虚轴两个方向的投影画成时间信号 figure(2),subplot(2,2,1),plot(real(q(4,:)),imag(q(4,:))),grid on subplot(2,2,3),plot(real(q(4,)》,t),grid on%画出实信号的时间 波形 subplot(2,2,2),plot(t,imag(q(4.:)).grid on %画出虚信号的时 间波形 ◆程序运行的结果 运行程序并按提示输入，如果只取两个幅度相等，频率符号相反的集总频谱， 那么将得到与图911相仿的结果。现取四个集总频谱，输入 a(1F1,w1)=1: a(2)Fl,w2)=1: a(3)=0.5,w(33: a(4)=0.5,w(1)=-4: 在此处，为了显示复信号，我们有意把输入频谱设成不对称的。于是读者将 看到四节杆的运动动画，并得到杆系及其端点在复平面上的轨迹，如图912，改 变了比例尺的轨迹见图9-l3子图()。将它在x,y两方向的投影与时间轴的关系 画在子图b)和（©）中，我们就得到信号与系统课程中常见的实信号曲线。i，a(i)=input('振幅 a(i)= '); w(i)=input('角频率 w(i)='); end % （2）将各个频谱向量相加合成并画图 % 此处应该把各时刻的图形转为动画，此处省略了动画的语句） t=0:0。1:20;lt=length(t); % 给出时间数组 p=a'*ones(1，lt)。*exp(j*w'*t); % 各频谱分量随时间变化的复 数值 q=cumsum(p); % 各频谱分量的累加（包括所有节 点） figure(1)，plot(real(q(4，:))，imag(q(4，:)))，grid on % 画合成复信号 的端点轨迹 % （3）将此轨迹的在实轴与虚轴两个方向的投影画成时间信号 figure(2)，subplot(2，2，1)，plot(real(q(4，:))，imag(q(4，:)))，grid on subplot(2，2，3)，plot(real(q(4，:))，t)，grid on % 画出实信号的时间 波形 subplot(2，2，2)，plot(t，imag(q(4，:)))，grid on % 画出虚信号的时 间波形 ◆ 程序运行的结果 运行程序并按提示输入，如果只取两个幅度相等，频率符号相反的集总频谱， 那么将得到与图 9-11 相仿的结果。现取四个集总频谱，输入 a(1)=1，w(1)=-1； a(2)=1，w(2)=-1； a(3)=0。5，w(3)=3； a(4)=0。5，w(1)=-4； 在此处，为了显示复信号，我们有意把输入频谱设成不对称的。于是读者将 看到四节杆的运动动画，并得到杆系及其端点在复平面上的轨迹，如图 9-12，改 变了比例尺的轨迹见图 9-13 子图(a)。将它在 x，y 两方向的投影与时间轴的关系 画在子图(b)和(c)中，我们就得到信号与系统课程中常见的实信号曲线