信号matlab题目.doc

1.离散时间系统： stem：画数字序列针状图。参数：画线处坐标t；画线处数值x。 例：t=[0:0.1:2]; x=cos(pi*t+0.6); stem(t,x); conv：实现两个序列的卷积。参数：序列向量a,b。输出参数：a,b卷积的结果向量y。 例：a=[1 1 1 1 1]; b=[1 2 3 4 5 6 7 8 9]; c=conv(a,b); stem(c); 设单位脉冲响应h(n)等于sin(0.5n),n≥0，输入x(n)等于sin(0.2n),n≥0。计算当n=0,1,…,40时的响应y(n).（2-5） n=0:40; x = sin(.2*n); h = sin(.5*n); y = conv(x,h); stem(n,y(1:length(n))) 有连续时间信号：画出信号的波形。 t = 0:20/400:20; w1 = 1; w2 = 4; w3 = 8; A1 = input(Input the amplitude A1 for w1 = 1: ); A2 = input(Input the amplitude A2 for w2 = 4: ); A3 = input(Input the amplitude A3 for w3 = 8: ); x = A1*cos(w1*t)+A2*cos(w2*t+pi/3)+A3*cos(w3*t+pi/2); 5.画出矩形脉冲序列。 t = -3:6/1000:3; N = input(Number of harmonics ); c0 = 0.5; w0 = pi; xN = c0*ones(1,length(t)); % dc component for k=1:2:N, % even harmonics are zero theta = ((-1)^((k-1)/2) - 1)*pi/2; xN = xN + 2/k/pi*cos(k*w0*t + theta); end plot(t,xN) title([Example 3.4, N = ,num2str(N)]) xlabel(Time (sec)) ylabel([x,num2str(N),(t)]) 6.画出信号的幅频谱和相频谱3.8 w = 0:0.2:50; b = 10; X = (1)./(b+j*w); clf subplot(211),plot(w,abs(X)); % plot magnitude of X title(Example 3.8) xlabel(Frequency (rad/sec)) ylabel(|X|) subplot(212),plot(w,angle(X)*180/pi); % plot angle of X in deg. xlabel(Frequency (rad/sec)) ylabel(Angle(X), degrees) subplot(111) 7.利用符号运算计算拉氏变换求的拉氏变换：6.3 syms x b t x = exp(-b*t); X = laplace(x) % x = sym(1); % creates a symbolic expression for 1 X = laplace(x) 8.用residue求拉氏反变换。；6.17 num = [ 1 2]; den = [1 4 3 0]; [r,p] = residue(num,den) 如果是复数极点（或重极点等），类似。6.18 num = [ 1 -2 1]; den = [1 3 4 2]; [r,p] = residue(num,den) 9.设传函为，输入为，求y(t)6.34 syms s H X Y H = 2/(s^2+4*s+4); X = 1/s; Y = H*X; y = ilaplace(Y) % use simplify(y) to simplify result ezplot(y,[0,10]); axis([0 10 0 .75]) 10.求z反变换，求解x[n]7.16 syms X x z X = (8*z^3+2*z^2-5*z)/(z^3-1.75*z+.75); x = iztrans(X) % 上面部分求得了x. % 下面部分求得了x[n]，第一个是x[0]. q = 10; num = [8 2 -5 0]; den = [1 0 -1.75 .75]; x = filter(num,den,[1 zeros(1,q-1)]) 11.求系统的响应7.20 根据理论计算，可以求得系统传函为 N0时，x=0，定义一初始向量，用filter命令可以得到阶跃输入响