用MATLAB實(shí)現(xiàn)常用的離散時(shí)間信號(hào)及其時(shí)域運(yùn)算

1、用MATLAB實(shí)現(xiàn)常用的離散時(shí)間信號(hào)及其時(shí)域運(yùn)算摘要:1. 運(yùn)用MATLAB勺繪圖指令繪制離散時(shí)間信號(hào)；2. 用MATLA實(shí)現(xiàn)離散時(shí)間信號(hào)的時(shí)域運(yùn)算。關(guān)鍵詞：MATLAB 離散時(shí)間信號(hào)時(shí)域運(yùn)算實(shí)驗(yàn)?zāi)康模?. 復(fù)習(xí)MATLA語(yǔ)言以及常用指令2. 學(xué)習(xí)和掌握用MATLA語(yǔ)言產(chǎn)生離散時(shí)間信號(hào)的編程方法3. 通過(guò)編程繪制出離散時(shí)間信號(hào)波形，加深理解信號(hào)的時(shí)域運(yùn)算實(shí)驗(yàn)原理：1.1單位階躍序列和單位樣值序列離散時(shí)間信號(hào)只在某些離散瞬間時(shí)給出信號(hào)的值，因此，它是時(shí)間上不連續(xù)的序列。單位階躍序列和單位樣值序列在離散信號(hào)與系統(tǒng)的分析中是兩個(gè)非常典型的序列，分別記為u(n)和3(n),他們的定義分別如下：0 (n

2、 0)u(n)=U (n 工 0)(n)二1 (n = 0)0 (n = 0)若單位節(jié)約序列的起始點(diǎn)為 no,單位樣值序列出現(xiàn)在no時(shí)刻，則表達(dá)式為:u(n - no)0 (n no)1 (n 一 no)1 (n = no)(n - no)二io (n 式 no)2離散時(shí)間信號(hào)的時(shí)域運(yùn)算與連續(xù)時(shí)間信號(hào)的研究類(lèi)似，在離散系統(tǒng)分析中，經(jīng)常遇到離散時(shí)間的運(yùn)算，包括 兩信號(hào)的相加、相乘以及序列本身的移位、反褶、尺度等等，也需要了解在運(yùn)算過(guò)程 中序列的以及相應(yīng)的波形的變化。序列x(n)的反褶:x( _ n)序列x(n)的移位：x( n - no)序列x(n)的尺度：x ( a n )兩序列x(n)和xn

3、)相乘:x(n) x£n)兩序列x(n)和xn)相乘:x(n) x£n)兩序列x(n)和x(n)相加減:x(n) _ X2(n)兩序列x(n)和xn)相乘:x(n) x£n)試驗(yàn)任務(wù):1繪制用于產(chǎn)生以下信號(hào)的通用程序1 x( n)二 A ( nn o)2 x( n) = Au( n - mi) - u( n - m2)1.1程序和圖形見(jiàn)附錄11.2程序和圖形見(jiàn)附錄22已知波形，求表達(dá)式并繪出各信號(hào)波形f ( k - 2)u( k)f(k2)f ( k - 2)u( k - 2)2.程序和圖形見(jiàn)附錄3附錄1:% n1:起始時(shí)刻；n0:跳變時(shí)刻；n2 :終止時(shí)刻;t1

4、=i nput('請(qǐng)輸入起始時(shí)刻n1=');t0=i nput('請(qǐng)輸入跳變時(shí)刻n0=');t2=i nput('請(qǐng)輸入終止時(shí)刻n2=');A=i nput('請(qǐng)輸入幅度A=');t=t1:t2; n=len gth(t); tt=t1:t0;n 1=le ngth(tt); x=zeros(1, n); x(n 1)=A*1; stem(t,x),gridoncleartitle('x(n)')axis(t1 t2 -0.2 A+1)n1=-2, n0=1, n2=3, A=2附錄2:Clear%n1:起始時(shí)刻；

5、m1第一次跳變時(shí)刻；m2:第二次跳變時(shí)刻；n2 :終止時(shí)刻;t1=input('請(qǐng)輸入起始時(shí)刻n1=');t0=in put('請(qǐng)輸入第一次跳變時(shí)刻m仁');t3=input('請(qǐng)輸入第二次跳變時(shí)刻m2=');t2=input('請(qǐng)輸入終止時(shí)刻n2=');A=input('請(qǐng)輸入幅度A=');t=t1:t2;n=len gth(t);tt=t0:t2;n 1=le ngth(tt);x=zeros(1, n-n 1),o nes(1, n1);ttt=t3:t2;n 2=le ngth(ttt);y=zeros(1

6、, n-n 2),o nes(1, n2);z=A*(x-y);stem(t,z),gridontitle( 'x ( n)')axis(t1, t2 ,min( z)-1, max(z)+1)n1=-3， n2=3， m仁-1，m2=1，A=2附錄3:1. 離散序列反褶函數(shù)(Isfz.m)fun ctionf,k=lsfz(f1,k1)f=fliplr(f1);k=-fliplr(k1);2. 離散序列移位函數(shù)(lspy.m)fun cti onf,k=lspy(ff,kk,kO)k=kk+kO;f=ff;3. 離散序列相乘函數(shù)(lsxc.m)fun ctionf,k=lsx

7、c(f1,f2,k1,k2)k=mi n(min (k1),mi n(k2):max(max(k1),max(k2); s仁 zeros(1,le ngth(k);s2=s1;s1(fi nd(k>=mi n(k1)&(k<=max(k1)=1)=f1;s2(fi nd(k>=mi n( k2)&(k<=max(k2)=1)=f2; f=s1.*s2;Retur n4.主程序:'f(k)')'f(k-2)u( n)')'f(-k+2)')cleari=1;for m=-3:4if m<0f(i)=m+

8、3;elseif (m=4)f(i)=0；elsef(i)=3；endi=i+1;end%f(k)m=-3:4subplot(221),stem(m,f,'filled' ),grid on ,title(axis(mi n( m)-1,max(m)+1,mi n(f)-0.5,max(f)+1)i=1;for n=-4:6if n<0u(i)=0;elseu(i)=1;endi=i+1;endn=-4:6m=-3:4%f(k-2)y1,k1=lspy(f,m,2);%f(-k)y2,k2=lsfz(f,m);%u(k-2)y3,k3=lspy(u, n,2);%f(k-2)u(k)z1,k=lsxc(y1,u,k1, n)subplot(222),ste m( k,z1,'filled' ),grid on ,title(axis(mi n(k)-1,max(k)+1,mi n(z1)-0.5,max(z1)+1)%f(-k+2)z2,w2=lspy(y2,k2,2)subplot(223),stem(w2,z2,'filled' ),grid on ,title(axis(mi n(w2)-1,max(w2)+1,mi n(z2)-1,max(z2)+1)%f(k-2)u(k-2)z3,w3=lsxc(y1,y