实验一熟悉环境
一、实验目的
(1)熟悉的主要操作命令。
(2)学会简单的矩阵输入和数据读写。
(3)掌握简单的绘图命令。
(4)用编程并学会创建函数。
(5)观察离散系统的频率响应。
二、实验内容
认真阅读本章附录,在环境下重新做一遍附录中的例子,体会各条命令的含义。在熟悉了基本命令的基础上,完成以下实验。
上机实验内容:
(1)数组的加、减、乘、除和乘方运算。输入[1 2 3 4],[3 4 5 6],求,,.*B,,.^B 并用语句画出A、B、C、D、E、F、G。
实验程序:
[1 2 3 4];
[3 4 5 6];
1:4;
.*.^B;
(4,2,1)(,'') ('时间序列n')('A');
(4,2,2)(,'') ('时间序列n ')('B');
(4,2,3)(,'') ('时间序列n ')('');
(4,2,4)(,'') ('时间序列n ')('');
(4,2,5)(,'') ('时间序列n ')('A.*B');
(4,2,6)(,'') ('时间序列n ')('');
(4,2,7)(,'') ('时间序列n ')('A.^B');
运行结果:
(2)用实现以下序列。
a)x(n)=0.8n0≤n≤15
实验程序:
0:150.8.^n;
(,''); ('时间序列n ')('x(n)=0.8^n');
b)x(n)(0.2+3j)n 0≤n≤15
实验程序:
0:15((0.2+3*j)*n);
(,''); ('时间序列n ')('x(n)((0.2+3*j)*n)'); 运行结果:
a)的时间序列 b)的时间序列c)x(n)=3(0.125π0.2π)+2(0.25π0.1π) 0≤n≤15
实验程序:
0:1:15;
3*(0.125**0.2*)+2*(0.25**0.1*);
(,''); ('时间序列n ');
('x(n)=3*(0.125**0.2*)+2*(0.25**0.1*)');
运行结果:
d)将c)中的x(n)扩展为以16为周期的函数x16(n)(16),绘出四个周期
实验程序:
0:1:63;
3*(0.125**(n,16)+0.2*)+2*(0.25**(n,16)+0.1*);
(,''); ('时间序列n ')('x16(n)');
e)将c)中的x(n)扩展为以10为周期的函数x10(n)(10),绘出四个周期
实验程序:
0:1:39;
3*(0.125**(n,10)+0.2*)+2*(0.25**(n,10)+0.1*);
(,''); ('时间序列n ')('x10(n)');
运行结果:
d )的时间序列
e )的时间序列
(3)x(n)=[11,3,5],产生并绘出下列序列的样本。
a )x 1(n)=2x(2)(1)-2x(n)
实验程序:
0:3;
[1 -1 3 5];
x1(x,[0 -2])2(x,[0 1])3=2*x12-2*x;
(x3,''); ('时间序列n ')('x1(n)=2x(2)(1)-2x(n)');
b )∑=-=5
1
k 2)k n (nx (n) x
实验程序:
0:3;
[1 -1 3 5];
x1(x,[0 1])2(x,[0 2])3(x,[0 3]);
x4(x,[0 4])5(x,[0 5]);
1*x1+2*x2+3*x3+4*x4+5*x5;
(,''); ('时间序列n ');
('x2(n)(1)+2x(2)+3x(3)+4x(4)+5x(5)');
运行结果:
a)的时间序列b)的时间序列
(4)绘出时间函数的图形,对x轴、y轴图形上方均须加上适当的标注。
a) x(t)(2πt) 0≤t≤10s b) x(t)(100πt)(πt) 0≤t≤4s
实验程序:
;
t1=0:0.001:102=0:0.01:4;
(2**t1)(100**t2).*(*t2);
(2,1,1);
(t1) ('t')('x(t)')('x(t)(2**t)');
(2,1,2);
(t2) ('t')('x(t)')('x(t)(100**t2).*(*t2)');
运行结果: