东南大学仪器科学和工程学院大三短学期matlab第二次作业.doc

word格式精心整理版 范文范例 学习指导 1、设 其中ω(n)为0均值方差为1的白噪声 分别用xcorr()函数和conv函数，计算x(n)和y(n)的互相关函数 答: k=0:0.0001:1; w=randn(1,length(k)); x=sin(k*pi/10+pi/3); y= sin((k-2)*pi/10+pi/3)+w; xcorr(x,y); 2、2、求系统 频率特性，频率范围为 ，求系统的单位脉冲响应。同时求H(s)的部分分式展开。 答：频率特性 s1=[1 2]; s2=[1 2]; s3=[1 3]; den=conv(conv(s1,s2),s3); num=[4 17 16]; w=logspace(0.1,10); freqs(num,den,w) 系统的单位脉冲响应 s1=[1 2]; s2=[1 2]; s3=[1 3]; den=conv(conv(s1,s2),s3); num=[4 17 16]; sym=tf(num,den); impulse(sym); 部分分式展开 s1=[1 2]; s2=[1 2]; s3=[1 3]; den=conv(conv(s1,s2),s3); num=[4 17 16]; [R P K] = residue(num,den) R = 1.0000 3.0000 -2.0000 P = -3.0000 -2.0000 -2.0000 K = [] 3、系统函数 ，当系统输入为 时，画出系统的输出。同时求系统的频率响应，画出幅频响应和相位响应图。 答： 输出 num=[1 0 0]; den=[1 0.8 0.5]; n=0:1:200; x=5+cos(0.2*pi*n)+2*sin(0.7*pi*n); y=lsim(num,den,x,n); plot(y) 幅频响应和相位响应图 num=[1 0 0]; den=[1 0.8 0.5]; n=0:1:200; x=5+cos(0.2*pi*n)+2*sin(0.7*pi*n); w=logspace(0.1,10); freqs(num,den,w) 4 针对差分方程 编写MATALB程序计算其脉冲响应（即当输入 输出y(n)就是系统的脉冲响应h(n)）.并画图显示波形图 答： b=1;a=[1,-0.7,-0.6,1]; x=[1,zeros(1,200)]; y=filter(b1,a,x); plot(y) 5 用matlab实现由系统的差分方程求出的系统的就冲击响应，并画出其幅度、相位响应，并支出系统的零点和极点。差分方程为： 答： b=[1,-2];a=[1,-1,0.6]; x=[1,zeros(1,200)]; y=filter(b,a,x); plot(y) freqz(y,x) 6、增加数据的长度对谱分析的影响： 模拟信号 ， 以t=0.01*n( )进行采样，求N点的信号频谱，N分别为64、256。比较信号谱结构的不同，并且说明原因。 答： N1=64; n1=0:N1-1; t1=0.01*n1; N2=256; n2=0:N2-1; t2=0.01*n2; x1=2*sin(4*pi*t1)+5*cos(8*pi*t1); x2=2*sin(4*pi*t2)+5*cos(8*pi*t2); subplot(2,2,1); stem(x1) subplot(2,2,2); stem(x2) subplot(2,2,3); X1 = fft(x1,N1); stem(X1) subplot(2,2,4); X2 = fft(x2,N2); stem(X2)