西工大信号与系统_实验1

1、、实验目的西北工业大学信号与系统实验报告西北工业大学专业word可编辑N=L2;M=7；n=I0:2*H-L;Xm'£ln(2+pi+Kin)/1T)；%tith='pg网 of sin( (2+pi+M+nJ/V) subplot <2, 2, 3);stem (n, Xro);ElabelC n');ylabeLC k; (n/ );N=L2;N=10;n=lO:2*H-L;Km=sln(.2*pl*K*rii/T);%ti t le = ' phase cf sin.( (2+pi+M+nJ /N) subplot 12, 2, 4J ;s

2、tem in, Xm);slabeLC ;ylabelC zlOCn)1 );学习使用MATLAB,巩固课本知识二、实验要求认真完成三、实验设备(环境)电脑，带MATLAB;四、实验内容与步骤认真看实验册，完成后面习题五、实验结果1.2离散时间正弦信号N=12;M=4;n=f0:2ttf-l;%title="phase of sinff的哪n)制'subplot (2, 2,1);steifi(n, Kui 1xlabele');ylabelt' (n)');N=L2;M=5;n=0:2*tf-l;小号IM (2*国*Wni看”%title='

3、 phase of sinftLtWinj/N)1 subplot 之外；5tem(n,加;xlabfil /n );ylal>elC x5 in)');a.I 也 Figure 1- X却阳 0物E 西侬 ffiAffi IHCD 蒯I 口出 精即出）上图分别是 0<n<2N-1,M=4,5,7,10 时，Xmn的图像。由上图可看出，当M=4时，基波周期T=3 ; M=5时，基波周期T=12; M=7时，基波 周期T=12 ; M=10时，基波周期 T=6 ;所以当M=4时，得到的最小整数周期为3。题中彳t号 Xm(n尸sin(2 jMn/N)的频率 w=2 ttM

4、/N ,由公式得周期 T=2k兀/w ,有T=kN/M(k=1,2,)。当N/M 为正整数时，最小周期T=N/M ;当N/M 为有理数时，不管M>N 或者M<N ,都有最小周期T=N ;当N/M为无理数时，该序列不是周期序列b.k-2wk=75:. o g ：?：;：=三inH.m ： subplot (2,2, 2)7| st t" n, "?1 : ?：_alel1 a' jylabel (- x2 ()');»k=(2*pi*k)75Kk=sin(wk*n) ubplot 2, 4) stem (n, xk) ?label (?

5、i label (J x6(n)')以上是代码，下图是运行结果才 Figure 1 X文性由丽F）爸看（Vl 1ftAm IWD总面电胃口州）ffiKCH）j j j k j,:卷*“ 0r盟 口由上图可看出，图上一共有3个唯一的信号。当k=1和k=6的时候的图像是一样的。因为档 k= 1 时，wk=(2*PI)/5,k=6 时，wk=2PI+(2*PI)/5,即 w6 =2PI+w1 ,因为sin函数的周期是2PI,所以他俩的图像是一样的c.代码如下：N=6;subplot (.3, 1T 1) nl-0:24;xl-cos (2*pl+nl. /N) +sin(3+pi*nt. /

6、N);stem(.nl, xl) grid;subplot (3, 1, 2) n2=0:24;工2=2联c口与24口2. /N)(3>hi2. /N)二stem(n2, x2)耳 rid;ubplet(3, 1T 3)n3=0;48;茂3=cq5(2+pi+n3. /N)+3+5in(3卡口i卡n3. / (2+N): stem(n3, x3)arid:图像如下：* Figure 1文F)强辑间查看% 搐入工具(T)里面。)窗口的帮助(H)信3餐、。©05包口国国可得出结论：如果2*pi/w0不是有理数，则该信号不是周期的1.3离散时间信号时间变量的变换a. nx=zeros

7、(1,3) 2 0 1 -1 3 zeros(1,3); 图像如下:+ Figure 1 U X空伸门潴锵旧有看(V)国入工具(D定面国口 (W)精助(H) a j同'a0qh 口困fsq1J'cI'711CC£111-3*2-101234567b.代码如下：x=zeros(1,11);x(4)=2;x(6)=1;x(7)=-1；x(8)=3;n=-3:7;n1=n-2;n2=n+1;n3=-n;n4=-n+1;y1=x ;y2=x ;y3=x ;y4=x ;c:代码和结果如下结果n3-D：48,工3=曰5 ：加 i和1 /li)*3*s1ti i.5*p L

8、tnS, /:sten(n3.s3)/id：» s=:erP5(L, Ll);i(4)=2;i (6) =1；10)-3;hf-3：7;nl=jr2;施irLMi;n4=-ti+l;jL=i; y即黯 y>i; Rr; aibplDt(2, 1) steolnL it); title (' yr);grid; sihpletli?)stenM祖;titleC y.：');二力；比 士口 他母：titlef );grid;5ifc(pct(i2,3)title (' 1);srH»底7 Figure 1iW 身触酹M 歌® Ilaj陋苗

9、口副粗®3 OO 血卜图是结果图A Figure 1 X文件旧事旧百看应)注人(1)工具工)卓面身口幽用助(H)日U i Q二气。3小 园卫1-10-505专业word可编辑X超前2得到y1,; x延日1得到y2; x倒置再延时1得到y3; x倒置再延时2得到y4.1.4 离散时间系统a:证明：当 n=0 时，若输入 x1n =6n,则 x1(0) =1,y10)=1;若 x2n=26n,则 x2(0)=2 , y2(0)=sin(pi/4)= (根号 2) /2,结果不等于2倍y1(0)，因此是非线形的。b ：证明：yn=xn+xn+1=un+un+1, 当 n=-1 时，y-1=

10、u-1+u0=1, 即当n<0时，存在y(n)不等于0,因此原系统不是因果的MATLAB：(使用了 heaviside 函数)for n=-5:9x(n)=heaviside(n);endfor n=-6:9y(n)=heaviside(n)+heaviside(n+1);endxy在实验过程中，以上代码总是显示不知如何解决。卡标索引必须为正整数类型或逻辑类型.1.5 实现一阶差分方程1.6 连续时间复指数信号» x=sym(r sinJ)警告： Support ct' strings that are not valid variable names or defin

11、e a number will be removed in a future release. To create syrobolic expressions, first create symbolic variables and then use operations on them. > In sracons510n (line 1536)In sym>convert；Chiir (Line 14411Ind line 1198)In sii (line 177)sin(2*pi+t) /T)>> x5=subs (xh *115)=sln(2+pi+t)/5)&

12、#187; t=0：10；>> ezplot (x5, t);a:代码如图:过 程 中 发 现 了 课 本 中 的就创建MA LLAB中的符号表达式表示工)oX的变量是单一的字符串T和'T:函数日小日用于1 变最，所以必须将yn的基波周期设置为某一监殊值 用 SLtJS必二 at>s(X,5,'T)f于是,X5 tt是髭口(2n才市)的二个符号表达式口试对心创建符号表达式，并用电p山画出网个周期内I 果正无误的话，应该与图L2所示一致.fhi对下面侑已创律格甘表法式就是函数x5的两个参数 5和T”的位置放反了 。代码的结果如下图： Figure 1 X文件旧

13、胸指出 查看M 插入0)工具卬 桌面。窗口(W)帮助国) d d Q W 堤乂 口国 和书上的图1.2是一致的。b:正余弦函数分别定义如下：T=4>> i=syiiC sin(2*pi*t/T)T);y=synC cDS(2*pi*t/I)');f, 4);中subsxxl 招*1;t=0:32;ezplot(zr t);T=8x=sym C sin(2*pi*t/T)J ；y=symC cos(2*pi*t/T)1 );xl=subs (x, ' T' 1 8);yl=subs (y, ' T' , 8);z=xl+j1;t=0:32;ez

14、plot (zr t);文件（D 筑福田 有看M 插入出 工具 京面）窗口皿 帮助但）口日d）玲a-以/豆:口园niT=16>> x=symC sln(2*pi*t/T)J);y=sym(7 cos(2*pi*t/T)");xl=subs (x,7 T'，16);yl=subs (y, Il , 16):z=xltyl;t=0:32;ezplot (z, t);专业word可编辑次 Fiqure 1 X文件(F)箱走日查看(V) A(|J工具m桌面(D)窗口帮助E)口1口0|4口用口cos(t力助顺打8)X(t)的基波周期为：T1= T/21.8连续时间信号的能量与功率xl=sym(? cos (pi*t/5) *);x2=symC sin(pit/5)1);x3=syn)(, (exp(i*(2*pi*t/3) + (exp(i*(pi*t)');a: b: x1的图像 T=10专业word可编辑"X Figure 1 X文件麻辑查看(Y)播入3 工具CD桌面Q)窗口物帮助(WX2的图像：T=10盯盲 H J7"-1