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

總結

以上是生活随笔為你收集整理的西北工业大学matlab安装,西北工业大学Matlab实验报告.doc的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。