PAGE
PAGE 21
西北工业大学
《基于MATLAB的数字信号处理》实验报告
学 院:
计算机学院
学 号:
姓 名:
专 业:
计算机科学与技术
西北工业大学
2017年 07 月
实验一 MATLAB基本编程实验
一、实验目的及要求
1.熟悉MATLAB运行环境;
2. 掌握MATLAB的基本语法和函数;
3. 掌握MATLAB的基本绘图功能
二、实验设备(环境)及要求
1. 计算机
2. Matlab软件编程实验平台
三、实验内容与步骤
1. 设
求 使
2 . 分别用for或while循环结构编写程序,求出 。并考虑一种避免循环语句的程序设计算法实现同样的运算。
3. 将图形窗口4等分,按顺序绘制出下列函数的曲线图及曲面图并标记。
四、设计思想
1.,则X=B*inv(A)
2. eq \o\ac(○,1)编写循环程序,循环105次,每次循环中计算的值
eq \o\ac(○,2)对做数学运算,可以得到
3.使用subplot将图形窗口分为四个子窗口,在子窗口中进行绘图
五、程序代码及注释
% question 1*************
% 计算矩阵B*inv(A)
A=[ 0 2 1
2 -1 3
-3 3 -4];
B=[ 1 2 3
2 -3 1];
X=B*inv(A)
%*************************
% question2***************
% 分别使用for及while循环计算矩阵K的值
K=0;
for i=1:105
K=K+sqrt(3)/2^(i);
end
K
%
K=0;
i=1;
while i<106
K=K+sqrt(3)/2^(i);
i=i+1;
end
K
%*************************
% question 3*************
% 首先计算x,y1,y2,y3,X,Y,Z的值,然后画图
x=0:0.1:2*pi;
y1=cos(x);
y2=2*sin(2*x);
y3=x.*cos(x);
[X,Y]=meshgrid(-2:0.1:2,-4:0.4:4);
Z=X.^2/2^2+Y.^2/4^2;
subplot(2,2,1)
plot(x,y1)
axis([0 6.5 min(y1)-0.1 max(y1)+0.1]);
xlabel('x')
ylabel('y1')
title('x-y1 graph')
subplot(2,2,2)
plot(x,y2)
axis([0 6.5 min(y2)-0.1 max(y2)+0.1]);
xlabel('x')
ylabel('y2')
title('x-y2 graph')
subplot(2,2,3)
plot(x,y3)
axis([0 6.5 min(y3)-0.1 max(y3)+0.1]);
xlabel('x')
ylabel('y3')
title('x-y3 graph')
subplot(2,2,4)
surf(X,Y,Z);
axis([-2 2 -4 4 0 2]);
xlabel('X')
ylabel('Y')
zlabel('Z')
title('X-Y-Z graph')
%*************************
六、实验结果
1.
2. K=1.7321
3.
实验二 MATLAB基本编程实验
一、实验目的及要求
1. 回顾数字信号处理的主要内容;
2. 掌握利用MATLAB进行信号处理的方法;
3. 了解信号处理工具箱中一些函数的功能;
二、实验设备(环境)及要求
1. 计算机
2. Matlab软件编程实验平台
三、实验内容
1 .任何实数序列x(n)都能分解成为它的偶部分量和奇部分量之和,编写一个MATLAB函数文件,该函数的功能是将一给定序列分解成为它的偶部分量和奇部分量。并通过M文件或命令行调用该函数文件将以下序列分解成为偶部分量和奇部分量。
2. 一数字滤波器由 频率响应函数描述
1)写出其差分方程表示;
2)画出上面滤波器的幅频和相频图;
3)产生信号 的200个样本,通过该滤波器得到输出y(n