科技大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)(編輯修改稿)
2024-07-14 21:42 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 mag(p),39。*39。)。hold on。plot(real(z),imag(z),39。o39。)。xlabel(39。t39。)。title(39。零極點(diǎn)圖39。)。grid %運(yùn)行程序發(fā)現(xiàn)不穩(wěn)定，分析沖激響應(yīng) %沖激響應(yīng)t=0::3。h=impulse(sys,t)。subplot(1,2,2)。plot(h)。xlabel(39。t39。)。title(39。沖激響應(yīng)h(t)39。)。Grid(2) 1個(gè)極點(diǎn)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。零極點(diǎn)圖39。)。grid %運(yùn)行程序發(fā)現(xiàn)不穩(wěn)定，分析沖激響應(yīng) %沖激響應(yīng)t=0::3。h=impulse(sys,t)。subplot(1,2,2)。plot(h)。xlabel(39。t39。)。title(39。沖激響應(yīng)h(t)39。)。grid(3) 2個(gè)共軛極點(diǎn)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。零極點(diǎn)圖39。)。grid %運(yùn)行程序發(fā)現(xiàn)不穩(wěn)定，分析沖激響應(yīng) %沖激響應(yīng)t=0::3。h=impulse(sys,t)。subplot(1,2,2)。plot(h)。xlabel(39。t39。)。title(39。沖激響應(yīng)h(t)39。)。grid(4) 2個(gè)共軛極點(diǎn)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。零極點(diǎn)圖39。)。grid %運(yùn)行程序發(fā)現(xiàn)不穩(wěn)定，分析沖激響應(yīng) %沖激響應(yīng)t=0::3。h=impulse(sys,t)。subplot(1,2,2)。plot(h)。xlabel(39。t39。)。title(39。沖激響應(yīng)h(t)39。)。grid(5) 零點(diǎn)在s =0. 5，極點(diǎn)在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。零極點(diǎn)圖39。)。grid %運(yùn)行程序發(fā)現(xiàn)不穩(wěn)定，分析沖激響應(yīng)%沖激響應(yīng)t=0::3。h=impulse(sys,t)。subplot(2,2,2)。plot(h)。xlabel(39。t39。)。title(39。沖激響應(yīng)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) 零點(diǎn)在s =0. 5，極點(diǎn)在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。零極點(diǎn)圖39。)。grid %運(yùn)行程序發(fā)現(xiàn)不穩(wěn)定，分析沖激響應(yīng)%沖激響應(yīng)t=0::3。h=impulse(sys,t)。subplot(2,2,2)。plot(h)。xlabel(39。t39。)。title(39。沖激響應(yīng)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 離散系統(tǒng)分析1. 已知一個(gè)LTI系統(tǒng)的差分方程為： (1)初始條件y(1)=1,y(2)=2，輸入x(n)=u(n)，計(jì)算系統(tǒng)的零輸入響應(yīng)； (2)當(dāng)下面三個(gè)信號(hào)分別通過(guò)系統(tǒng)，分別計(jì)算系統(tǒng)的響應(yīng)： (3)指出這是一個(gè)什么特性的系統(tǒng)。 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。零輸入響應(yīng)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
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1