科技大學數字信號處理課程設計(編輯修改稿)

　

【文章內容簡介】 mag(p),39。*39。)。hold on。plot(real(z),imag(z),39。o39。)。xlabel(39。t39。)。title(39。零極點圖39。)。grid %運行程序發現不穩定，分析沖激響應 %沖激響應t=0::3。h=impulse(sys,t)。subplot(1,2,2)。plot(h)。xlabel(39。t39。)。title(39。沖激響應h(t)39。)。Grid(2) 1個極點s =0，增益k =1； a=[1,0]。k=1。b=[1]*k。sys=tf(39。s39。)。sys=tf(b,a)。p=roots(a)。z=roots(b)。subplot(1,2,1)。plot(real(p),imag(p),39。*39。)。hold on。plot(real(z),imag(z),39。o39。)。xlabel(39。t39。)。title(39。零極點圖39。)。grid %運行程序發現不穩定，分析沖激響應 %沖激響應t=0::3。h=impulse(sys,t)。subplot(1,2,2)。plot(h)。xlabel(39。t39。)。title(39。沖激響應h(t)39。)。grid(3) 2個共軛極點s=177。j5，增益 k =1 ；a=conv([1 5j],[1 5j])。k=1。b=[1]*k。sys=tf(39。s39。)。sys=tf(b,a)。p=roots(a)。z=roots(b)。subplot(1,2,1)。plot(real(p),imag(p),39。*39。)。hold on。plot(real(z),imag(z),39。o39。)。xlabel(39。t39。)。title(39。零極點圖39。)。grid %運行程序發現不穩定，分析沖激響應 %沖激響應t=0::3。h=impulse(sys,t)。subplot(1,2,2)。plot(h)。xlabel(39。t39。)。title(39。沖激響應h(t)39。)。grid(4) 2個共軛極點s=0 .1177。 j5，增益k =1； a=conv([1 +5j],[1 ])。k=1。b=[1]*k。sys=tf(39。s39。)。sys=tf(b,a)。p=roots(a)。z=roots(b)。subplot(1,2,1)。plot(real(p),imag(p),39。*39。)。hold on。plot(real(z),imag(z),39。o39。)。xlabel(39。t39。)。title(39。零極點圖39。)。grid %運行程序發現不穩定，分析沖激響應 %沖激響應t=0::3。h=impulse(sys,t)。subplot(1,2,2)。plot(h)。xlabel(39。t39。)。title(39。沖激響應h(t)39。)。grid(5) 零點在s =0. 5，極點在s=0 .1177。j5，增益k =1； a=conv([1 +5j],[1 ])。k=1。b=[1 ]*k。sys=tf(39。s39。)。sys=tf(b,a)。p=roots(a)。z=roots(b)。subplot(2,2,1)。plot(real(p),imag(p),39。*39。)。hold on。plot(real(z),imag(z),39。o39。)。xlabel(39。t39。)。title(39。零極點圖39。)。grid %運行程序發現不穩定，分析沖激響應%沖激響應t=0::3。h=impulse(sys,t)。subplot(2,2,2)。plot(h)。xlabel(39。t39。)。title(39。沖激響應h(t)39。)。grid%幅頻特性w=0::5。h=freqs(b,a,w)。subplot(2,2,3)。plot(w,abs(h))。xlabel(39。w39。)。title(39。幅頻特性39。)。grid(6) 零點在s =0. 5，極點在s=0 .1177。j5，增益k =1； a=conv([1 +5j],[1 ])。k=1。b=[1 ]*k。sys=tf(39。s39。)。sys=tf(b,a)。p=roots(a)。z=roots(b)。subplot(2,2,1)。plot(real(p),imag(p),39。*39。)。hold on。plot(real(z),imag(z),39。o39。)。xlabel(39。t39。)。title(39。零極點圖39。)。grid %運行程序發現不穩定，分析沖激響應%沖激響應t=0::3。h=impulse(sys,t)。subplot(2,2,2)。plot(h)。xlabel(39。t39。)。title(39。沖激響應h(t)39。)。grid%幅頻特性w=0::5。h=freqs(b,a,w)。subplot(2,2,3)。plot(w,abs(h))。xlabel(39。w39。)。title(39。幅頻特性39。)。grid 離散系統分析1. 已知一個LTI系統的差分方程為： (1)初始條件y(1)=1,y(2)=2，輸入x(n)=u(n)，計算系統的零輸入響應； (2)當下面三個信號分別通過系統，分別計算系統的響應： (3)指出這是一個什么特性的系統。 N=50。n=0:N1。a=[1,]。b=[,]。x=zeros(N)。 zi=filtic(b,a,[1,2])。y=filter(b,a,x,zi)。x1=cos(pi/10.*n)。x2=cos(pi/5.*n)。x3=cos(pi/10*7.*n)。y1=filter(b,a,x1,zi)。y2=filter(b,a,x2,zi)。y3=filter(b,a,x3,zi)。figure(1)subplot(4,1,1)。stem(n,y)。xlabel(39。n39。)。ylabel(39。y(n)39。)。title(39。零輸入響應39。)。subplot(4,1,2)。stem(n,y1)。xlabel(39。n39。)。ylabel(39。y(n)39。)。title(39。x(n)=cos(pi/10.*n)39。)。subplot(4,1,3)。stem(n,y2)。xlabel(39。n39。)。ylabel(39。y(n)39。)。title(39。x(n)=cos(pi/5.*n)39。)。subplot(4,1,4)。stem(n,y3)。xlabel(39。n39。)。ylabe