Matlab上机模拟试题.docVIP

Matlab上机模拟试题 计算当时，的值。 计算 求解矩阵方程： 求解矩阵方程： 计算序列{-2 0 1 –1 3}和序列{1 2 0 -1}的离散卷积。 用MATLAB计算差分方程 P138 157 当输入序列为时的输出结果。 b=[1 0.7 -0.45 -0.6]; x=[1 zeros(1,N-1)]; k=0:1:N-1; y=filter(a,b,x); stem(k,y) xlabel(n);ylabel(幅度) 7、用MATLAB计算差分方程 P156 所对应的系统函数的 A=[1 0.7 -0.45 -0.6]; freqz(B,A); 用FFT计算序列的频谱。 xn=[1 3 -1 1 2 3 3 1]; Xn=fft(xn); figure(1); plot(abs(Xn)); 9、求下列直接型系统函数的零、极点。 b=[1,-0.1,-0.3,-0.3,-0.2]; a=[1,0.1,0.2,0.2,0.5]; disp(H(z)); printsys(b,a,z); [z,p,k]=tf2zp(b,a) 10、画出序列 的图形并补零，序列长度为20. xn=[1 2 3 4 5]; b=zeros(1,15); xn=[xn,b]; stem(xn) 产生一个1024点的正弦加白噪声序列，并绘其频谱。 x=[0:1023]; y=sin(x); y1=awgn(y,10); freqz(y1) 已知信号，求其FFT变换的幅值和相位值，并绘出结果。 y=2.*(cos(2.*x)).*(cos(100.*x)); y1=fft(y); freqz(y1) 13、给定， 。分别求两序列相加、相卷积。用4个子图显示。 n1=0:4; x1=[1 2 3 4 5]; subplot(4,1,1),stem(n1,x1);ylabel(x1); n2=0:3; x2=[ 1 2 1 2]; subplot(4,1,2),stem(n2,x2); ylabel(x2); n=0:4; x1=[x1 ,zeros(1,5-length(n1))]; x2=[ x2,zeros(1,5-length(n2))]; x=x1+x2; y=conv(x1,x2); subplot(4,1,3),stem(n,x); ylabel(x); subplot(4,1,4),stem(y); ylabel(y) 生成一个长度为50的汉宁窗，并观察期频率特性。 n=51; window=hanning(n); [h,w]=freqz(window,1); subplot(1,2,1) stem(window); subplot(1,2,2); plot(w/pi,20*log(abs(h)/abs(h(1)))); 15、已知有限长序列x(n)=[7,6,5,4,3,2],求DFT和IDFT，要求：画出序列傅立叶变换对应的幅度谱和相位谱；画出原信号与傅立叶逆变换IDFT[X(k)]的图形进行比较。 xn=[7,6,5,4,3,2]; N=length(xn); n=0:(N-1);k=0:(N-1); Xk=xn*exp(-j*2*pi/N).^(n*k); x=(Xk*exp(j*2*pi/N).^(n*k))/N; subplot(2,2,1),stem(n,xn); title(x(n)); subplot(2,2,2),stem(n,abs(x)); title(IDFT|X(k)|); subplot(2,2,3),stem(k,abs(Xk)); 已知周期序列的主值x(n)=[7,6,5,4,3,2],求x(n)周期重复次数为3次时的DFS和IDFS。要求：画出原信号序列的主值和周期序列的图形；画出离散傅立叶变换对应的幅度谱和相位谱.