《实验一 熟悉matlab环境》由会员分享，可在线阅读，更多相关《实验一 熟悉matlab环境（9页珍藏版）》请在金锄头文库上搜索。
1、实验一 熟悉 MATLAB 环境一、实验目的（1） 熟悉 MATLAB 的主要操作命令。（2） 学会简单的矩阵输入和数据读写。（3） 掌握简单的绘图命令。（4） 用 MATLAB 编程并学会创建函数。（5） 观察离散系统的频率响应。二、实验内容（1） 数组的加、减、除和乘方运算。输入 A=1 2 3 4，B=3 4 5 6，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.B 并用 stem 语句画出A、B、C、D、E、F、G。（2） 用 MATLAB 实现下列序列：a) nx8.0)(15b) je)32.(c) )1.025.sin().0.cos) nnnx 5d) 将 c)
2、中的 扩展以 16 为周期的函数 绘出四个)(x ),6(6x周期。e) 将 c)中的 扩展以 10 为周期的函数 绘出四个)(n ),10()10n周期。（3） ，产生并绘出下列序列的样本5,31)(nx)(21()2(1 nx（4） 绘出下列时间函数的图形，对 x 轴、y 轴以及图形上方均须加上适当的标注：a) )2sin()ttxst10b) )si(10cot4（5） 编写函数 stepshift(n0,n1,n2)实现 u(n-n0)，n1n0n2 ，绘出该函数的图形，起点为 n1，终点为 n2。 （6） 给定一因果系统 求出并绘)9.067.1/()21() 2121 zzzH制
3、的幅频响应和相频响应。)(zH（7） 计算序列8 -2 -1 2 3和序列2 3 -1 -3的离散卷积，并作图表示卷积结果。（8） 求以下差分方程所描述系统的单位脉冲响应 h(n), 50n)1(2)(06.)1(.0)( xnyny三、思考题（1） 比较实验内容第（2）题中 d)和 e)两小题的结果，试说明对于周期性信号，应当如何采样，才能保证周期扩展后与原信号保持一致？（2） 对于有限长序列，如何用 MATLAB 计算其 DTFT?四、实验报告要求（1） 简述实验目的及原理。（2） 按实验步骤附上实验程序。（3） 按实验步骤附上有关离散系统的频率特性曲线。（4） 简要回答思考题。（1）A=
4、1 2 3 4A =1 2 3 4B=3 4 5 6B =3 4 5 6C=A+BC =4 6 8 10D=A-BD =-2 -2 -2 -2E=A.*BE =3 8 15 24F=A./BF =0.3333 0.5000 0.6000 0.6667G=A.BG =1 16 243 40960 2 4024 A0 2 40246 B0 2 40510 C0 2 4-2-10 D0 2 402040 E0 2 400.51 F0 2 405000 GA=1 2 3 4; B=3 4 5 6;C=A+B;D=A-B;E=A.*B; F=A./B;G=A.B;subplot(3,3,1);stem(
5、A); title(A);subplot(3,3,2);stem(B); title(B);subplot(3,3,3);stem(C); title(C);subplot(3,3,4);stem(D); title(D);subplot(3,3,5);stem(E); title(E);subplot(3,3,6);stem(F); title(F);subplot(3,3,7);stem(G); title(G);（2）a）n=0:15;x=0.8.nx =Columns 1 through 14 1.0000 0.8000 0.6400 0.5120 0.4096 0.3277 0.26
6、21 0.2097 0.1678 0.1342 0.1074 0.0859 0.0687 0.0550Columns 15 through 16 0.0440 0.0352b) n=0:15;x=exp(0.2+3j)*n)x =Columns 1 through 8 1.0000 -1.2092 + 0.1724i 1.4324 - 0.4168i -1.6602 + 0.7509i 1.8780 - 1.1942i -2.0650 + 1.7677i 2.1923 - 2.4934i -2.2212 + 3.3928iColumns 9 through 16 2.1010 - 4.4854
7、i -1.7673 + 5.7857i 1.1398 - 7.3006i -0.1198 + 9.0242i -1.4106 -10.9326i 3.5900 +12.9763i -6.5776 -15.0719i 10.5514 +17.0909ic) n=0:15;x=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi)x =Columns 1 through 14 3.0451 3.3495 2.3714 0.2076 -2.3814 -4.3399 -4.8652 -3.8251 -1.8090 0.2145 1.4328 1.6083 1.
8、1453 0.7759Columns 15 through 16 1.0610 2.0091d) n=0:15;x=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi);y=x x x xy =Columns 1 through 14 3.0451 3.3495 2.3714 0.2076 -2.3814 -4.3399 -4.8652 -3.8251 -1.8090 0.2145 1.4328 1.6083 1.1453 0.7759Columns 15 through 28 1.0610 2.0091 3.0451 3.3495 2.3714 0
9、.2076 -2.3814 -4.3399 -4.8652 -3.8251 -1.8090 0.2145 1.4328 1.6083Columns 29 through 42 1.1453 0.7759 1.0610 2.0091 3.0451 3.3495 2.3714 0.2076 -2.3814 -4.3399 -4.8652 -3.8251 -1.8090 0.2145Columns 43 through 56 1.4328 1.6083 1.1453 0.7759 1.0610 2.0091 3.0451 3.3495 2.3714 0.2076 -2.3814 -4.3399 -4
10、.8652 -3.8251Columns 57 through 64 -1.8090 0.2145 1.4328 1.6083 1.1453 0.7759 1.0610 2.0091e) n=0:15;x=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi);z=x(1:10);x=z z z zx =Columns 1 through 14 3.0451 3.3495 2.3714 0.2076 -2.3814 -4.3399 -4.8652 -3.8251 -1.8090 0.2145 3.0451 3.3495 2.3714 0.2076Col
11、umns 15 through 28 -2.3814 -4.3399 -4.8652 -3.8251 -1.8090 0.2145 3.0451 3.3495 2.3714 0.2076 -2.3814 -4.3399 -4.8652 -3.8251Columns 29 through 40 -1.8090 0.2145 3.0451 3.3495 2.3714 0.2076 -2.3814 -4.3399 -4.8652 -3.8251 -1.8090 0.21450 20 40 60 80-5-4-3-2-101234 y0 10 20 30 40-5-4-3-2-101234 xn=0:
12、15;x=0.8.n;n=0:15;x=exp(0.2+3j)*n);n=0:15;x=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi);n=0:15;x=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi);y=x x x x;n=0:15;x=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi);z=x(1:10);x=z z z z;subplot(1,2,1);stem(y);title(y);subplot(1,2,2);stem(x);title(
13、x);（3）x=1 -1 3 5;z1=x 0 0 0;z2=0 0 x 0;z3=0 0 0 x;x1=2*z1-z3-2*z2x1 = 2 -2 4 11 -5 -13 -5x=1 -1 3 5;z=x 0 0 0 0 0;z1=0 x 0 0 0 0;z2=0 0 x 0 0 0;z3=0 0 0 x 0 0;z4=0 0 0 0 x 0;z5=0 0 0 0 0 x;Z=z1;z2;z3;z4;z5;n=1:length(Z);x2=n.*sum(Z)x2 = 0 2 0 12 40 48 49 64 45（4）t=0:0.1:10;x=sin(2*pi*t);plot(t,x);ax
14、is(0 10 -1.2 1.2); xlabel(时间);ylabel(幅值)0 1 2 3 4 5 6 7 8 9 10-1-0.8-0.6-0.4-0.200.20.40.60.81值值值值t=0:0.05:4;x=sin(2*pi*t).*cos(100*pi*t);plot(t,x);axis(0 4 -1.2 1.2); xlabel(时间);ylabel(幅值)0 0.5 1 1.5 2 2.5 3 3.5 4-1-0.8-0.6-0.4-0.200.20.40.60.81值值值值（5） n0=5;n1=1;n2=10;n=n1:n2;u=zeros(1,n0-n1) ones(1,n2-n0+1);stem(n,u);1 2 3 4 5 6 7 8 9 1000.10.20.30.40.50.60.70.80.91（6）t=0:1:600*pi/600;b=1 20.5 1;a=1 -0.67 0.9;freqz(b,a,t)a =1.0000
