2012东南大学MATLAB考试试题及答案.docVIP

MATLAB上机测验题 （考试时间：2:20----4:20） 姓名 ZZZ 学号 04010xyy 考试要求： 1、要求独立完成不得与他人共享，答卷雷同将做不及格处理。 2、答卷用Word文件递交，文件名为学号+姓名.doc，试卷写上姓名及学号。 3、答卷内容包括： 程序； 运行结果及其分析； 图也要粘贴在文档中。 上机考题： 一、系统传递函数为，按照以下要求求解： 求其极零点图，判断系统的稳定性，画出系统的频谱特性； 当系统输入信号为：，时，画出系统的输出。 答：1） num=[1,1,0]; den=[1,-1,0.81]; [a,b,c,d]=tf2ss(num,den); sys=ss(a,b,c,d，-1); pzmap(sys); 分析：系统极零点仅有一个零点在单位圆上，其余均在单位圆内，所以系统稳定。 num=[1,1,0]; den=[1,-1,0.81]; [a,b,c,d]=tf2ss(num,den); sys=ss(a,b,c,d,-1); freqz(num,den,512) 2） num=[1,1,0]; den=[1,-1,0.81]; [a,b,c,d]=tf2ss(num,den); sys=ss(a,b,c,d,-1); n=0:50; e=5+cos(0.2*pi*n)+2*sin(0.7*pi*n); lsim(sys,e,n) 分析：图中灰色波形为输入信号，蓝色为输出信号。 二、系统传递函数为， 画出系统的零极点图，判断稳定性； 给定频率范围为[0,10]，步长为0.1，画出其频率响应； 画出系统的单位脉冲响应。 答：1） num=[1,11,30,0]; den=[1,9,45,97,50]; [a,b,c,d]=tf2ss(num,den); sys=ss(a,b,c,d); pzmap(sys) 分析：系统除一个零点在虚轴上，其余极零点均在虚轴左侧，所以系统稳定； 2） num=[1,11,30,0]; den=[1,9,45,97,50]; [a,b,c,d]=tf2ss(num,den); sys=ss(a,b,c,d); w=0:0.2*pi:20*pi; %2=2*pi*f freqs(num,den,w); 3） num=[1,11,30,0]; den=[1,9,45,97,50]; [a,b,c,d]=tf2ss(num,den); sys=ss(a,b,c,d); impulse(sys) 三、系统传递函数为，初始状态为，激励信号为，求系统的零输入、零状态和全响应。 答： 零输入响应： num=[1,3]; den=[1,3,2]; x0=[4,-1]; [a,b,c,d]=tf2ss(num,den); sys=ss(a,b,c,d); initial(sys,x0) title(零输入相应) 零状态响应： %function epxl function y=epxl(t) if t=0 y=1; else y=0; end %零状态响应 num=[1,3]; den=[1,3,2]; x0=[4,-1]; [a,b,c,d]=tf2ss(num,den); sys=ss(a,b,c,d); t=0.0001:0.1:10; e=exp(-2*t).*(epxl(t)-epxl(t-6)); lsim(sys,e,t) title(零状态相应) 分析：图中蓝色信号为零状态响应波形 全响应： num=[1,3]; den=[1,3,2]; x0=[4,-1]; [a,b,c,d]=tf2ss(num,den); sys=ss(a,b,c,d); t=0.0001:0.1:10; e=exp(-2*t).*(epxl(t)-epxl(t-6)); lsim(sys,e,t,x0) title(全响应) 四、信号，以进行采样，N分别为(1)N＝70、(2)N＝200，给出两种条件下信号x(t)的频谱图(用FFT命令)，并分析差异以及产生差异的原因。 答：(1)%N=70 clc clear all N=70; t=0:0.01:0.01*(N-1); x=2*sin(4*pi*t)+5*cos(8*pi*t); G=fft(x,128)/100; subplot(2,1,1); plot(abs(G));title(N=70幅频特性); subplot(2,1,2); plot(angle(G));title(N=70相频特性); % N=200; clc clear all N=200; t=0:0.01:0.01*(N-1); x=2*sin(4*pi*t)+5*cos(8*pi*t); G=fft(x,128)/100; subplot(