信号与系统Matlab作业一.doc

信號與系統Matlab作業一 Matlab程式碼如下: % 電機三 張毓麟 % (a) 做出x1, x2, x3, x4然後畫圖, 用subplot和hold各畫一遍 clear all clc close all % -20=n=100 n=-20:100; % define the functions x1=exp(j*(pi/4)*n); x2=sin(pi*n/8+pi*n/16); x3=(9/10).^n; x4=n+1; % 開始畫圖 % 以下是用Subplot的 % X1實部 subplot(2,3,1) stem(n, real(x1),b); axis([-20,100,-1,1]); title(x1 Real Part); % X1 Image Part subplot(2,3,2); stem(n, imag(x1), b); axis([-20,100,-1,1]); title(x1 Image Part); % X2 subplot(2,3,3); stem(n, x2, b); axis([-20,100,-1,1]); title (x2); % X3 subplot(2,3,4); stem(n, x3, b); axis([-20,100,0,10]); title (x3); % X4 subplot(2,3,5); stem(n, x4, b); axis([-20,100,-50,150]); title (x4); ----------------------------------------------------------------------------------------------------------------------------------------------------------- 所得圖形如下頁 然後接下來是使用hold放在同一張圖上，之後的想用hold都可以用這個改，所以不贅印 % 使用single figure, hold繪製圖形 clc hold on stem(n, real(x1), y); %Yellow stem(n, imag(x1), k); %Black stem(n, x2, b); %Blue stem(n, x3, g); %Green stem(n, x4, r); %Red title (Using hold to combine these vectors in a single figure); legend(x1 Real part, ,x1 Image Part, ,x2, ,x3, ,x4, ); (b) % (b)使用filter B=[1 0.9]; A=[1 -0.25]; y1=filter(B,A,x1); y2=filter(B,A,x2); y3=filter(B,A,x3); y4=filter(B,A,x4); % 開始畫圖 subplot(2,3,1); stem(n, real(y1)); title(y1 Real Part); axis([0 100 -3 3]); subplot(2,3,2); stem(n, imag(y1)); title(y1 Image Part); axis([0 100 -3 3]); subplot(2,3,3); stem(n, y2); title(y2); axis([0 100 -3 3]); subplot(2,3,4); stem(n, y3); title(y3); axis([0 100 0 5]); subplot(2,3,5); stem(n, y4); title(y4); axis([0 100 0 300]); 本題的使用filter即可解決 比較input和output，可得到為eigenfunction (c) % (c)eigenvalue H1=y1./x1; H2=y2./x2; H3=y3./x3; H4=y4./x4; % 將圖繪出 subplot(2,3,1); stem(n, real(H1)); title(H1 Real Part); axis([0 100 0 2]); subplot(2,3,2); stem(n, imag(H1)); title(H1 Image Part); axis([0 100 -1.4 0]); subplot(2,3,3); stem(n, H2); title(H2); axis([0 100 -3 3]); subplot(2,3,4); stem(n, H3); title(H3); axis([0 100 0 3])