试验一 熟悉MATLAB环境(实验报告)

试验一熟悉MATLAB环境(实验报告)
一、实验目的
(1) 熟悉MA TLAB的主要操作命令；
(2) 学会简单的短阵输入和数据读写；
(3) 掌握简单的绘图命令；
(4) 用MATLAB编程井学会创建函数；
(5) 观察离散系统的频率响应。

二、实验内容
认真阅读本章附录，在MATLAB环境下重新做一遍附录中的例子，体会各条命令的含义。

在熟悉了MATLAB的基本命令的基础上，完成以下实验。

上机实验内容：
（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。

>> 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
>> subplot(3,3,1);stem(A);
>> subplot(3,3,2);stem(B);
>> subplot(3,3,3);stem(C);
>> subplot(3,3,4);stem(D);
>> subplot(3,3,5);stem(E);
>> subplot(3,3,6);stem(F);
>> subplot(3,3,7);stem(G);
（2）用MATLAB 实现下列序列：
a ）()0.8n x n = 015n ≤≤
n=[0:15];
>> a=0.8;
>> x=a.^n
b ）()3cos(0.1250.2)2sin(0.250.1)x n n n ππππ=+++ 015n ≤≤
>> n=[0:15];
>> x=3*cos(0.125*pi*n+o.2*pi)+2*sin(0.25*pi*n+0.1*pi)
（3）用plot 函数绘制二维图形：
a) y=sin(πt)+1/cos(πt)+2
t=0.1:0.001:10;
y=(sin(pi*t)+1)./(cos(pi*t)+2);
plot(t,y);
xlabel('t');
ylabel('y');
(4) 在n=[-10;10]范围内画出以下信号：
a) x1(n)= δ(n+2)
n1=-10;n2=10;n0=-2;
n=n1:n2;
x=[n==n0]; stem(n,x,'filled');
b) x2(n)=δ(n-4)
n1=-10;n2=10;n0=4;
n=n1:n2;
x=[n==n0];
stem(n,x,'filled');
c) x3(n)=2δ(n+2) -δ(n-4)
n1=-10;n2=10;n01=4;n02=-2; n=n1:n2;
x1=[n==n01];x2=[n==n02]; x=2*x2-x1;
stem(n,x,'filled');
（5）产生复信号，并画出其实部与虚部及幅度和相位：a）x1(n)=ej (π/8 )n 0≤n≤32
b) x1(n)=e(-0.1+j0.3)n
n1=[0:32];n2=[-10:10];
x1=exp(j*(pi/8)*n1);
x2=exp((-0.1+j*0.3)*n2);
re1=real(x1);
im1=imag(x1);
a1=abs(x1);
th1=angle(x1);
subplot(3,3,1);stem(re1);
subplot(3,3,2);stem(im1);
subplot(3,3,3);stem(a1);
subplot(3,3,4);stem(th1);
re2=real(x2);
im2=imag(x2);
a2=abs(x2);
th2=angle(x2);
subplot(3,3,5);stem(re2);
subplot(3,3,6);stem(im2);
subplot(3,3,7);stem(a2);
subplot(3,3,8);stem(th2);
(6) 利用conv函数计算序列{8 -2 -1 2 3}和序列{2 3 -1 -3}的卷积，并作图表示卷积结果
f1=[8,-2,-1,2,3];
f2=[2,3,-1,-3];
y=conv(f1,f2);
stem(y,'filled');
(三)、思考题
（1）对于两个子系统级联或并联的系统，如何用matlab计算它们的幅频响应与相频相应？答：级联转换为直接型：cas2dir，并联转换为直接型：par2dir然后应用freqz的函数。