信号与系统实习

1信號與系統實習實習二：Matlab簡介（指令介紹）

指導教授：沈岱範助教：葉思佑

簡瑞成劉建良賴欣茹2常用的數學函數3計算向量元素統計量的常用函數4矩陣的基本演練A=[123;456;789]B=[357;246;789]C=25矩陣乘法A*Bans=2837466488112 100139178A=123456789B=3572467896取出矩陣的值

A(2,2)ans=5A(1:2,1:2)ans=1245A=123456789B=3572467897純量乘法

A.*C

ans=24681012141618A=123456789C=28矩陣冪次方A.^Cans=149162536496481C.*Bans=610144812141618A=123456789B=357246789C=29刪除列或行A(2,:)=[]A=123789A(:,3)=[]A=124578A=12345678910轉置A'ans=147258369A=12345678911練習題1.利用指令取出[671011]2.利用指令取得A的轉置3.利用刪除的指令使得矩陣變為4.取出A的第四列與第四行做矩陣的相乘運算

Ans:600A=12345678

910111213141516A=12456813141612基本的繪圖指令13基本繪圖2-D繪圖給定X-data以及(或)Y-data

定義圖形顏色以及線條樣式…etc

繪製單一線條圖形

plot(xdata,ydata,‘color_linestyle_marker')14其他指令應用圖形範例TITLEXLABELYLABELLEGENDCURVES15x=1:0.1:2*pi;y=sin(x);z=cos(x);

subplot(2,1,1),plot(x,y,’r*’);

subplot(2,1,2),plot(x,z,’g--’);Subplot指令16subplot(2,1,1)subplot(2,1,2)12Subplot說明plot(x,y,’r*’);plot(x,z,’g--’);17範例1/5畫一個cosine的連續訊號假設有一函數y=4cos(6πt+π/6)時間點介於0與1之間，每隔0.01取一點求值繪出連續圖形則其MATLAB的程式為t=0:0.01:1;y=4*cos(6*pi*t+pi/6);plot(t,y)18假設有一函數y=4cos(6πt+π/6)時間點介於0與1之間，每隔0.01取一點求值，試繪出離散圖形則其MATLAB的程式為n=0:0.01:1;y=4*cos(6*pi*n+pi/6);stem(n,y)範例2/5註：stem用法同plot19

假設有一函數

x(t)=Asin(ω0t+ψ)exp(-at)

A=15,ω0=10π,ψ=0,a=6t介於0到1之間，每隔0.01取一點求值繪出連續與離散圖形其MATLAB的方程式為clear,clc;t=0:0.01:1.0;y=15*sin(10*pi*t).*exp(-6*t);plot(t,y,'g')holdon,stem(t,y,'r')範例3/5註：1.clear:清除workspace中的變數2.clc:清空commandwindow3.holdon:將目前描繪的圖形繪在上一個圖形中20holdonplot(t,y,'g')stem(t,y,'r')stem(t,y,'r')plot(t,y,'g')stem(t,y,'r')plot(t,y,'g')holdonplot(t,y,'g')stem(t,y,'r')之前的圖形被覆蓋保留之前的圖形21clear,clct=0:0.05:2*pi;f1=1,phase1=0;f2=2,phase2=0;x1=cos(2*pi*f1*t+phase1);x2=cos(2*pi*f2*t+phase2);subplot(311),plot(t,x1);title('frequency=1Hz');subplot(312),plot(t,x2);title('frequency=2Hz');subplot(313),plot(t,x1,'g',t,x2,'r');legend('frequency=1','frequency=2')頻率為1與2的弦波圖形範例4/522弦波相乘做調變1t=0:0.01:2*pi;f1=1,f2=2;x1=sin(2*pi*f1*t)x2=cos(2*pi*f2*t);x3=sin(2*pi*f1*t).*cos(2*pi*f2*t);subplot(311),plo