matlab第二次課

1、2.1 2.1 信號的生成與運算信號的生成與運算2.2 matlab2.2 matlab 繪圖之二維繪圖繪圖之二維繪圖2.3 matlab2.3 matlab 程序設(shè)計初步程序設(shè)計初步 主要內(nèi)容主要內(nèi)容 （4學(xué)時）學(xué)時） 直接輸入法直接輸入法 變量表達(dá)式（數(shù)）變量表達(dá)式（數(shù)）規(guī)則規(guī)則： 矩陣元素必須放在矩陣元素必須放在 中；中； 矩陣元素必須用矩陣元素必須用逗號或空格逗號或空格分隔分隔 ； 在在 內(nèi)矩陣的行與行之間必須用內(nèi)矩陣的行與行之間必須用分號分號分隔分隔 矩陣元素可以是任何矩陣元素可以是任何matlab表達(dá)式表達(dá)式 ，可以是實數(shù)，可以是實數(shù) ，也，也 可以是復(fù)數(shù)，復(fù)數(shù)可用特殊函數(shù)可以是復(fù)

2、數(shù)，復(fù)數(shù)可用特殊函數(shù)i，j 輸入輸入 矩陣的創(chuàng)建與修改矩陣的創(chuàng)建與修改2.1 信號的表示與運算信號的表示與運算注意：注意： matlab嚴(yán)格區(qū)分大小寫字母嚴(yán)格區(qū)分大小寫字母，因此，因此a與與A是兩個是兩個不同的變量。不同的變量。2.1 信號的表示與運算信號的表示與運算例：例： 輸入矩陣輸入矩陣 A =1 2 3 4 5 6 7 8 9用用matlab函數(shù)創(chuàng)建矩陣函數(shù)創(chuàng)建矩陣空陣空陣 matlab允許輸入空陣，當(dāng)一項操允許輸入空陣，當(dāng)一項操作無結(jié)果時，返回空陣。作無結(jié)果時，返回空陣。rand 隨機(jī)隨機(jī)矩陣矩陣 y=rand(3)eye 單位矩陣單位矩陣zeros 全部元素都為全部元素都為0的矩陣

3、的矩陣ones 全部元素都為全部元素都為1的矩陣的矩陣diag 對角陣對角陣y=zeros (5,3)y=y=eye(3,3)y=ones(3,3)y=diag(-3:3)y=diag(-2 6 4 12)矩陣的創(chuàng)建與修改矩陣的創(chuàng)建與修改2.1 信號的表示與運算信號的表示與運算用用load函數(shù)裝載矩陣函數(shù)裝載矩陣矩陣的創(chuàng)建與修改矩陣的創(chuàng)建與修改2.1 信號的表示與運算信號的表示與運算如 load sn.txt load a.matsn.txt a.mat為已建立的為已建立的2個文件個文件 直接修改直接修改可用可用 鍵找到所要修改的矩陣，用鍵找到所要修改的矩陣，用鍵移鍵移動到要修改的矩陣元素上即

4、可修改動到要修改的矩陣元素上即可修改。 指令修改指令修改 可以用可以用A( , )= 來修改。來修改。矩陣的創(chuàng)建與修改矩陣的創(chuàng)建與修改2.1 信號的表示與運算信號的表示與運算矩陣的運算矩陣的運算矩陣加、減（矩陣加、減（,）運算）運算規(guī)則：規(guī)則： 相加、減的兩矩陣必須有相同的行和列數(shù)相加、減的兩矩陣必須有相同的行和列數(shù),兩矩陣對應(yīng)元素相加減。兩矩陣對應(yīng)元素相加減。 允許參與運算的兩矩陣之一是標(biāo)量。標(biāo)量允許參與運算的兩矩陣之一是標(biāo)量。標(biāo)量與矩陣的所有元素分別進(jìn)行加減操作。與矩陣的所有元素分別進(jìn)行加減操作。例例: : 已知矩陣已知矩陣a=2 3 4;1 3 5 b=2 5 7;5 3 2 c=3求，

5、求，d=a+b e=c+a2.1 信號的表示與運算信號的表示與運算例例: :2. 矩陣乘（矩陣乘（ ）運算運算規(guī)則：規(guī)則： A矩陣的列數(shù)必須等于矩陣的列數(shù)必須等于B矩陣的行數(shù)矩陣的行數(shù) 標(biāo)量可與任何矩陣相乘。標(biāo)量可與任何矩陣相乘。例：例：d=-1;0;2求：求：f=pi*da=1 2 3;4 5 6;7 8 0b=1;2;3求：求：c=a*b2. 矩陣乘（矩陣乘（ ）運算運算規(guī)則規(guī)則（）（） A矩陣的列數(shù)必須等于矩陣的列數(shù)必須等于B矩陣的行數(shù)矩陣的行數(shù) 標(biāo)量可與任何矩陣相乘標(biāo)量可與任何矩陣相乘。矩陣的運算矩陣的運算2.1 信號的表示與運算信號的表示與運算3 矩陣元素乘（點乘）矩陣元素乘（點乘）

6、a b a，b兩數(shù)組必須有相同的行和列兩數(shù)組必須有相同的行和列 兩數(shù)組相應(yīng)元素相乘兩數(shù)組相應(yīng)元素相乘。a=1 2 3;4 5 6;7 8 9;b=2 4 6;1 3 5;7 9 10;求：求：a.*bans = 2 8 18 4 15 30 49 72 90例例: :矩陣的運算矩陣的運算2.1 信號的表示與運算信號的表示與運算a=1 2 3;4 5 6;7 8 9;b=2 4 6;1 3 5;7 9 10;a.*bans = 2 8 18 4 15 30 49 72 90 a=1 2 3;4 5 6;7 8 9;b=2 4 6;1 3 5;7 9 10;a*bans = 25 37 46 55

7、 85 109 85 133 172 矩陣乘與點乘比較矩陣乘與點乘比較矩陣元素相乘矩陣元素相乘 矩陣相乘矩陣相乘矩陣的運算矩陣的運算2.1 信號的表示與運算信號的表示與運算a./b=b.aa.b=b./aa./b=b.a 都是都是a的元素被的元素被b的對應(yīng)元素除的對應(yīng)元素除a.b=b./a 都是都是 b的元素被的元素被a的對應(yīng)元素除的對應(yīng)元素除例例: a=1 2 3;b=4 5 6; 求：求：c1=a.b c2=b./a c3=b.a c4=a./bc1 = c2= 4.0000 2.5000 2.0000c3 = c4= 0.2500 0.4000 0.5000 給出給出a,b對應(yīng)元素間的商

8、對應(yīng)元素間的商.4 矩陣元素除（點除）矩陣元素除（點除）矩陣的運算矩陣的運算2.1 信號的表示與運算信號的表示與運算 inv 矩陣求逆矩陣求逆 det 行列式的值行列式的值 eig 矩陣的特征值矩陣的特征值 diag 對角矩陣對角矩陣 矩陣轉(zhuǎn)置矩陣轉(zhuǎn)置 sqrt 矩陣開方矩陣開方5. 矩陣的函數(shù)運算矩陣的函數(shù)運算 例例: :a=1,2,3;4,5,6;7,8,9;b=inv(a)c=det(a)d=eig(a)e=diag(a) f=ag=sqrt(a)矩陣的運算矩陣的運算2.1 信號的表示與運算信號的表示與運算信號的生成信號的生成2.1 信號的表示與運算信號的表示與運算p信號可以用關(guān)于時間的

9、函數(shù)來表示；信號可以用關(guān)于時間的函數(shù)來表示；p通過計算信號函數(shù)的采樣值，生成信號序通過計算信號函數(shù)的采樣值，生成信號序列（向量）；列（向量）；p當(dāng)采樣間隔足夠小，可認(rèn)為是連續(xù)信號；當(dāng)采樣間隔足夠小，可認(rèn)為是連續(xù)信號；pMATLAB中信號可用中信號可用 列向量或行向量表示單通道信號列向量或行向量表示單通道信號 矩陣表示多通道信號。矩陣表示多通道信號。信號自變量信號自變量t的常用的常用2種生成方法種生成方法2.1 信號的表示與運算信號的表示與運算p t=linspace(t1,t2,N); 在區(qū)間在區(qū)間t1,t2s上均勻產(chǎn)生上均勻產(chǎn)生N個時間樣點值。個時間樣點值。 pt= t1:1/fs:t2 在

10、區(qū)間在區(qū)間t1,t2s上以采樣率上以采樣率fs取時間樣點值。取時間樣點值。例：例： 生成在生成在0,10s上長度為上長度為30的單位階躍信號，的單位階躍信號，并作圖表示。并作圖表示。2.1 信號的表示與運算信號的表示與運算信號的生成信號的生成t=linspace(0,10,30);y=ones(1,30); Stem(t,y)例：例： 生成幅度為生成幅度為1V，頻率為，頻率為5Hz的正弦信號。的正弦信號。信號長信號長1秒。設(shè)采樣頻率為信號頻率的秒。設(shè)采樣頻率為信號頻率的20倍。并倍。并作圖表示。作圖表示。2.1 信號的表示與運算信號的表示與運算信號的生成信號的生成ts=0;te=1; fs=1

11、00; f1=5; t=(ts:1/fs:te); y1=sin(2*pi*f1*t);plot(t,y1)信號起始和終止時刻信號起始和終止時刻采樣頻率采樣頻率 信號頻率信號頻率生成采樣時刻向量生成采樣時刻向量生成信號向量生成信號向量作圖表示信號作圖表示信號信號的生成信號的生成2.1 信號的表示與運算信號的表示與運算常用周期信號：常用周期信號： 正弦信號正弦信號sin(x) 余弦信號余弦信號cos(x) 方波信號方波信號square(x) 鋸齒信號鋸齒信號sawtooth(x) Sinc信號信號sinc(x)u生成周期三角波生成周期三角波 在在(0,1.5)上生成頻率上生成頻率為為50HZ的三

12、角波信號的三角波信號,采樣率為采樣率為10KHZ.sawtooth(T):對時間向量對時間向量T產(chǎn)生周期三角波產(chǎn)生周期三角波sawtooth(T,WIDTH):WIDTH為為(0,1)的偏離系數(shù)的偏離系數(shù)t=0:0.0001:1.5;y=sawtooth(2*pi*50*t);plot(t,y)00.050.10.150.2-1-0.500.51信號的生成信號的生成2.1 信號的表示與運算信號的表示與運算u生成周期三角波生成周期三角波 在在(0,1.5)上生成頻率上生成頻率為為50HZ的三角波信號的三角波信號,采樣率為采樣率為10KHZ.sawtooth(T):對時間向量對時間向量T產(chǎn)生周期三

13、角波產(chǎn)生周期三角波sawtooth(T,WIDTH):WIDTH為為(0,1)的偏離系數(shù)的偏離系數(shù)t=0:0.0001:1.5;y=sawtooth(2*pi*50*t,0.3);plot(t,y)00.050.10.150.2-1-0.500.5100.050.10.150.2-1-0.500.51信號的生成信號的生成2.1 信號的表示與運算信號的表示與運算u生成周期方波生成周期方波 在在(0,2.5)上生成頻率上生成頻率為為50HZ的三角波信號的三角波信號,采樣率為采樣率為1KHZ.square(T):對時間向量對時間向量T產(chǎn)生峰值為產(chǎn)生峰值為1的的周期周期方方波波square(T,DUT

14、Y):DUTY為占空百分比為占空百分比t=0:0.001:2.5;y=square(2*pi*50*t,30);plot(t,y)00.020.040.060.080.1-2-1.5-1-0.500.511.52信號的生成信號的生成2.1 信號的表示與運算信號的表示與運算信號的生成（信號的生成（離散信號生成離散信號生成）2.1 信號的表示與運算信號的表示與運算例：例：生成單位脈沖序列生成單位脈沖序列 =0)(nn000=1nnnnN=100；n0 =10=10；u=zeros(1,N); 生成生成N N個個0 0序列（序列（u u向量）向量）u(1,n0)=1; 修改修改u u 向量第向量第1

15、010個元素為個元素為1 1 stem(0:N-1,u) 作離散信號圖作離散信號圖信號的生成（信號的生成（離散信號生成離散信號生成）2.1 信號的表示與運算信號的表示與運算例：例：生成單位階躍序列生成單位階躍序列 =0)(nnu0001nnnn N=100;n0 =10;u=ones(1,N); %生成生成N N個個0 0序列（序列（u u向量）向量）u(1,1:n0)=0; %修改修改u u 向量第向量第1010個元素為個元素為1 1 stem(0:N-1,u) %作離散信號圖作離散信號圖信號的生成（信號的生成（信號疊加信號疊加）2.1 信號的表示與運算信號的表示與運算例：例：二個幅度為二個

16、幅度為1V1V，頻率為，頻率為5Hz5Hz和和8Hz8Hz的正弦信號的的正弦信號的疊加。信號長疊加。信號長1 1秒。設(shè)采樣頻率為秒。設(shè)采樣頻率為100Hz100Hz。并將包括。并將包括疊加信號在內(nèi)的三個信號表示在一張圖內(nèi)。疊加信號在內(nèi)的三個信號表示在一張圖內(nèi)。ts=0;te=1; fs=100; f1=5;f2=8; t=(ts:1/fs:te); y1=sin(2*pi*f1*t);y2=sin(2*pi*f2*t);plot(t,y1,t,y2)hold on;plot(t,y1+y2)信號的生成（信號的生成（信號疊加信號疊加）2.1 信號的表示與運算信號的表示與運算2.1 信號的表示與運

17、算信號的表示與運算例：例： 1.均值為均值為0.6、方差為、方差為0.1的的5階正態(tài)分布隨機(jī)矩陣。階正態(tài)分布隨機(jī)矩陣。命令如下：命令如下：y=0.6+sqrt(0.1)*randn(5)2.1 信號的表示與運算信號的表示與運算信號的生成（信號的生成（白噪聲信號白噪聲信號） 2.均值為均值為0、方差為、方差為0.1的白噪聲信號的白噪聲信號 。命令如下：命令如下：t=1:1/100:2; y=sqrt(0.1)*randn(size(t)1. plot-最基本的二維圖形指令最基本的二維圖形指令plot命令自動打開一個圖形窗口命令自動打開一個圖形窗口Figure用直線連接相鄰兩數(shù)據(jù)點來繪制圖形用直線

18、連接相鄰兩數(shù)據(jù)點來繪制圖形根據(jù)圖形坐標(biāo)大小自動縮擴(kuò)坐標(biāo)軸，將數(shù)根據(jù)圖形坐標(biāo)大小自動縮擴(kuò)坐標(biāo)軸，將數(shù)據(jù)標(biāo)尺及單位標(biāo)注自動加到兩個坐標(biāo)軸據(jù)標(biāo)尺及單位標(biāo)注自動加到兩個坐標(biāo)軸上上. plot的功能的功能2.2 matlab 繪圖之繪圖之二維繪圖二維繪圖 plot(x) 缺省自變量繪圖格式缺省自變量繪圖格式，x為向為向量量, 以以x元素值為縱坐標(biāo)，以相應(yīng)元素下標(biāo)元素值為縱坐標(biāo)，以相應(yīng)元素下標(biāo)為橫坐標(biāo)繪圖為橫坐標(biāo)繪圖 plot(x,y) 基本格式基本格式，以，以y(x)的函數(shù)關(guān)的函數(shù)關(guān)系作出直角坐標(biāo)圖，如果系作出直角坐標(biāo)圖，如果y為為nm的矩陣，的矩陣，則以則以x 為自變量，作出為自變量，作出m條曲線條曲

19、線 plot(x1,y1,x2,y2) 多條曲線繪圖格式多條曲線繪圖格式 plot(x,y,s) 開關(guān)格式開關(guān)格式 s設(shè)定曲線顏色和繪圖方式設(shè)定曲線顏色和繪圖方式. 或或plot(x1,y1,s1,x2,y2,s2，) plot的調(diào)用格式的調(diào)用格式2.2 matlab 繪圖之繪圖之二維繪圖二維繪圖S的標(biāo)準(zhǔn)設(shè)定值的標(biāo)準(zhǔn)設(shè)定值 字母字母 顏色顏色 標(biāo)點標(biāo)點 線型線型 y 黃色黃色 點線點線 m 粉紅粉紅 圈線圈線 c 亮藍(lán)亮藍(lán) 線線 r 大紅大紅 字線字線 g 綠色綠色 實線實線 b 藍(lán)色藍(lán)色 星形線星形線 w 白色白色 ： 虛線虛線 k 黑色黑色 (-) 點劃線點劃線2.2 matlab 繪圖之

20、繪圖之二維繪圖二維繪圖 如果已經(jīng)存在一個圖形窗口，如果已經(jīng)存在一個圖形窗口，plot命令則命令則清除當(dāng)前圖形，繪制新圖形清除當(dāng)前圖形，繪制新圖形 單窗口單曲線繪圖；單窗口單曲線繪圖；單窗口多曲線繪圖單窗口多曲線繪圖；單窗口多曲線分圖繪圖單窗口多曲線分圖繪圖；多窗口繪圖多窗口繪圖 可任意設(shè)定曲線顏色和線型可任意設(shè)定曲線顏色和線型 可給圖形加坐標(biāo)網(wǎng)線和圖形加注功能可給圖形加坐標(biāo)網(wǎng)線和圖形加注功能 plot的特點的特點2.2 matlab 繪圖之繪圖之二維繪圖二維繪圖1.單窗口單曲線繪圖單窗口單曲線繪圖例例1：x=0, 0.48,0.84,1,0.91,0.6,0.14 x1, x2, x3, x4

21、, x5, x6, x7,plot (x)2.2 matlab 繪圖之繪圖之二維繪圖二維繪圖2.單窗口多曲線繪圖單窗口多曲線繪圖例例2：t=0:pi/100:2*pi;y=sin(t);y1=sin(t+0.25);y2=sin(t+0.5);plot(t,y,t,y1,t,y2)yy1y22.2 matlab 繪圖之繪圖之二維繪圖二維繪圖例例 3：t=0:pi/100:2*pi; y=sin(t);y1=sin(t+0.25);y2=sin(t+0.5); y3=cos(t);y4=cos(t+0.25);y5=cos(t+0.5); plot(t,y,y1,y2,y3,y4,y5)t=0:

22、pi/100:2*pi;y3=cos(t);y4=cos(t+0.25);y5=cos(t+0.5);plot(t,y3);hold on; plot(t,y4); plot(t,y5);例例 4 4：3. 單窗口多曲線分圖繪圖單窗口多曲線分圖繪圖subplot subplot 子圖分割命令子圖分割命令調(diào)用格式：調(diào)用格式： subplot(m,n,p)subplot(m,n,p) 按從左至右，從上至下排列按從左至右，從上至下排列 行行列列繪圖序號繪圖序號subplot(1,3,1); plot(t,y)subplot(1,3,2); plot(t,y3)subplot(1,3,3); plot

23、(t,y2)subplot(3,1,1);plot(t,y)subplot(3,1,2);plot(t,y3)subplot(3,1,3);plot(t,y2)例：例：x=-5:0.1:5;y=x.2; subplot(2,2,1);plot(x,y); subplot(2,2,2);plot(y,x);t=x subplot(2,2,3);plot(t.*cos(t), t.*sin(t); subplot(2,2,4); plot3(t.*cos(t), t.*sin(t),t); grid on4. 多窗口繪圖多窗口繪圖 figure(n) 創(chuàng)建窗口函數(shù)，創(chuàng)建窗口函數(shù)，n為窗為窗 口順序

24、號?？陧樞蛱?。t=0:pi/100:2*pi;y=sin(t);y1=sin(t+0.25);y2=sin(t+0.5);plot(t,y) 自動出現(xiàn)第一個窗口自動出現(xiàn)第一個窗口figure(2)plot(t,y1) 在第二窗口繪圖在第二窗口繪圖figure(3)plot(t,y2) 在第三窗口繪圖在第三窗口繪圖 5.可任意設(shè)置顏色與線型可任意設(shè)置顏色與線型例例 4：plot(t,y,r-,t,y1,g:,t,y2,b*)6.圖形加注功能圖形加注功能將標(biāo)題、坐標(biāo)軸標(biāo)記、網(wǎng)格線及文字注將標(biāo)題、坐標(biāo)軸標(biāo)記、網(wǎng)格線及文字注釋加注到圖形上，這些函數(shù)為：釋加注到圖形上，這些函數(shù)為： title 給圖形加

25、標(biāo)題給圖形加標(biāo)題 xlable 給給x軸加標(biāo)注軸加標(biāo)注 ylable 給給y軸加標(biāo)注軸加標(biāo)注 text 在圖形指定位置加標(biāo)注在圖形指定位置加標(biāo)注 grid on(off) 打開、關(guān)閉坐標(biāo)網(wǎng)格線打開、關(guān)閉坐標(biāo)網(wǎng)格線 axis 控制坐標(biāo)軸的刻度控制坐標(biāo)軸的刻度例：例：t=0:0.1:10y1=sin(t);y2=cos(t);plot(t,y1,r,t,y2,b-);x=1.7*pi;1.6*pi;y=-0.3;0.8;s=sin(t);cos(t);text(x,y,s);title(正弦和余弦曲線正弦和余弦曲線);legend(正弦正弦,余弦余弦)xlabel(時間時間t),ylabel(正弦

26、、余弦正弦、余弦)gridaxis square例：例：一個幅度為一個幅度為1V1V，頻率為，頻率為5Hz5Hz的正弦信號。信號的正弦信號。信號長長1 1秒。設(shè)采樣頻率為秒。設(shè)采樣頻率為100Hz100Hz。被一個為均值為。被一個為均值為0 0、方、方差為差為0.250.25的白噪聲信號污染。請將原信號，白噪聲的白噪聲信號污染。請將原信號，白噪聲信號和受污染信號在一個圖形框內(nèi)畫出波形。信號和受污染信號在一個圖形框內(nèi)畫出波形。ts=0;te=1; fs=100 ;f1=5; t=ts:1/fs:te; y1=sin(2*pi*f1*t);y2=0.5*randn(size(t);plot(t,y

27、1,t,y2)hold on;plot(t,y1+y2)例題例題 （信號與白噪聲的疊加信號與白噪聲的疊加）綜合例題綜合例題 用戶如想靈活應(yīng)用用戶如想靈活應(yīng)用matlab去解決實際問題，去解決實際問題，充分調(diào)用充分調(diào)用matlab的科學(xué)技術(shù)資源，就需要編的科學(xué)技術(shù)資源，就需要編輯輯m文件文件 包含包含matlab語言代碼的文件稱為語言代碼的文件稱為m文件，其文件，其擴(kuò)展名為擴(kuò)展名為m。 編輯編輯m文件可使用各種文本編輯器。文件可使用各種文本編輯器。 m m文件文件2.3 Matlab程序設(shè)計初步程序設(shè)計初步（教材（教材P128） m文件的類型是普通的文本文件，我們可以文件的類型是普通的文本文件，

28、我們可以使用系統(tǒng)認(rèn)可的文本文件編輯器來建立使用系統(tǒng)認(rèn)可的文本文件編輯器來建立m文件。文件。如如 windows的記事本和的記事本和word等。等。 具體的創(chuàng)建方法：具體的創(chuàng)建方法： 1）在）在matlab命令命令 窗口點擊窗口點擊file 菜單菜單 new m-file m m文件的創(chuàng)建文件的創(chuàng)建2.3 Matlab程序設(shè)計初步程序設(shè)計初步 m m文件的分類文件的分類 m文件的語法類似于文件的語法類似于c語言，但又有其自身特語言，但又有其自身特點。它只是一個簡單的點。它只是一個簡單的ASCII碼文本文件，執(zhí)碼文本文件，執(zhí)行程序時逐行解釋運行程序，行程序時逐行解釋運行程序，matlab是解釋是解

29、釋性的編程語言。性的編程語言。 m文件有兩類文件有兩類: 獨立的獨立的m文件文件 稱命令稱命令(腳本腳本)文件文件 可調(diào)用可調(diào)用m文件文件 稱函數(shù)文件稱函數(shù)文件2.3 Matlab程序設(shè)計初步程序設(shè)計初步 (1). 命令文件命令文件 簡單的簡單的m文件文件 命令文件實際上是一串指令的集合，與在命命令文件實際上是一串指令的集合，與在命 令窗口逐行執(zhí)行文件中的所有指令，其結(jié)果令窗口逐行執(zhí)行文件中的所有指令，其結(jié)果 是一樣的。是一樣的。沒有輸入輸出參數(shù)沒有輸入輸出參數(shù)。2.3 Matlab程序設(shè)計初步程序設(shè)計初步 m m文件的分類文件的分類命令文件舉例命令文件舉例首先建立首先建立命令文件命令文件并以

30、文件名并以文件名mingl.m存盤。存盤。% 切比雪夫切比雪夫II型低通數(shù)字濾波器設(shè)計型低通數(shù)字濾波器設(shè)計clear;close ；clcwp=0.2;ws=0.4;Rp=1;Rs=80;% 輸入指標(biāo)輸入指標(biāo)N,wc=cheb2ord(wp,ws,Rp,Rs)% 求濾波器階次求濾波器階次B,A=cheby2(N,Rs,wc)% 設(shè)計濾波器，得出系數(shù)設(shè)計濾波器，得出系數(shù)freqz(B,A)% 無左端變量時自動畫頻率特性圖無左端變量時自動畫頻率特性圖在在MATLAB的命令窗口中輸入的命令窗口中輸入mingl，將執(zhí)行文件。，將執(zhí)行文件。說明：說明：1）命令文件可以直接運行，不被任何其他）命令文件可以

31、直接運行，不被任何其他m文件調(diào)文件調(diào) 用，用，一般常用于編制項目主程序。一般常用于編制項目主程序。 2）文件名不可以為純數(shù)字。）文件名不可以為純數(shù)字。 3）命令文件的開頭第一行語句建議為：）命令文件的開頭第一行語句建議為：clear;close；clc3.1 Matlab程序設(shè)計初步程序設(shè)計初步 m m文件的分類文件的分類(2). 函數(shù)文件函數(shù)文件u格式：格式： function 返回變量返回變量(輸出變量輸出變量)=函數(shù)名（輸入變量）函數(shù)名（輸入變量） 注釋說明語句段注釋說明語句段 程序語句段程序語句段2.3 Matlab程序設(shè)計初步程序設(shè)計初步Matlab提供了大量的提供了大量的內(nèi)置函數(shù)內(nèi)

32、置函數(shù)文件以供調(diào)用，對于文件以供調(diào)用，對于matlab沒有提沒有提供的、特殊用途的函數(shù)文件，學(xué)習(xí)者可自行按格式編寫供的、特殊用途的函數(shù)文件，學(xué)習(xí)者可自行按格式編寫自定義函自定義函數(shù)文件。數(shù)文件。內(nèi)置函數(shù)文件內(nèi)置函數(shù)文件： matlab庫函數(shù)文件庫函數(shù)文件內(nèi)置函數(shù)的調(diào)用方法：內(nèi)置函數(shù)的調(diào)用方法：使用函數(shù)名并給出相應(yīng)的入口、出口參數(shù)即可。使用函數(shù)名并給出相應(yīng)的入口、出口參數(shù)即可。2.3 Matlab程序設(shè)計初步程序設(shè)計初步內(nèi)置函數(shù)調(diào)用舉例內(nèi)置函數(shù)調(diào)用舉例 sin.m函數(shù)函數(shù) 調(diào)用格式：調(diào)用格式：y=sin(2*x) 實際應(yīng)用中：實際應(yīng)用中： x=0:2*pi/180:2*pi; y=sin(2*x

33、) plot(x,y)01234567-1-0.8-0.6-0.4-0.200.20.40.60.81又如又如 sqrt(x) 求求x的平方根的平方根 abs(x) 求求x的絕對值的絕對值 det(a) 求行列式的值求行列式的值 等等。等等。matlab擁有大量的內(nèi)置數(shù)學(xué)函數(shù)，這些擁有大量的內(nèi)置數(shù)學(xué)函數(shù)，這些數(shù)學(xué)函數(shù)都放在數(shù)學(xué)函數(shù)都放在matlab6p5toolboxmatlab 目錄下，供我們在數(shù)值運算和符號運算中調(diào)用。目錄下，供我們在數(shù)值運算和符號運算中調(diào)用。2.3 Matlab程序設(shè)計初步程序設(shè)計初步 m m文件的分類文件的分類(2). 函數(shù)文件函數(shù)文件u規(guī)則規(guī)則1.函數(shù)函數(shù)m文件第一行

34、必須以單詞文件第一行必須以單詞function作為引導(dǎo)詞，作為引導(dǎo)詞， 遵循如下形式：遵循如下形式： function =()2.m文件的文件名必須是文件的文件名必須是 .m。如：。如：r2b.m3.程序中的變量均為局部變量，不保存在工作空間中。其程序中的變量均為局部變量，不保存在工作空間中。其變量只在函數(shù)運行期間有效。變量只在函數(shù)運行期間有效。2.3 Matlab程序設(shè)計初步程序設(shè)計初步function y=r2b(x)% 該程序用于把該程序用于把 輸入的二進(jìn)制數(shù)轉(zhuǎn)換為十進(jìn)制數(shù)輸入的二進(jìn)制數(shù)轉(zhuǎn)換為十進(jìn)制數(shù)s=0;t=length(x);for i=0:t-1 s=s+x(t-i)*2i;en

35、dy=s;2.3 Matlab程序設(shè)計初步程序設(shè)計初步說明：說明： 1）函數(shù)文件依然以）函數(shù)文件依然以.m作為擴(kuò)展名作為擴(kuò)展名 2）必須在程序開頭以）必須在程序開頭以function作為關(guān)鍵字作為關(guān)鍵字3）函數(shù)文件的函數(shù)保存名必須與）函數(shù)文件的函數(shù)保存名必須與function后的函數(shù)后的函數(shù) 名稱名稱一一 致，如本例的函數(shù)文件必須保存為致，如本例的函數(shù)文件必須保存為r2b.m 4）命令文件的語句中一般）命令文件的語句中一般不出現(xiàn)不出現(xiàn)：clear; close；clc自定義函數(shù)文件自定義函數(shù)文件： 用戶自己編寫、完成所需功能的文件用戶自己編寫、完成所需功能的文件u程序控制結(jié)構(gòu)分類程序控制結(jié)構(gòu)分類 matlab語言的程序結(jié)構(gòu)與其它高級語言是一致的，語言的程序結(jié)構(gòu)與其它高級語言是一致的，主要有順序結(jié)構(gòu)和循環(huán)結(jié)構(gòu)：主要有順序結(jié)構(gòu)和循環(huán)結(jié)構(gòu)：(1).順序結(jié)構(gòu)順序結(jié)構(gòu) 依次順序執(zhí)行程序的各條語句依次順序執(zhí)行程序的各條語句(2).循環(huán)結(jié)構(gòu)循環(huán)結(jié)構(gòu) 被重復(fù)執(zhí)行的一組語句。被重復(fù)執(zhí)行的一組語句。2.3 Matlab程序設(shè)計初步程序設(shè)計初步u程序控制結(jié)構(gòu)分類程序控制結(jié)構(gòu)分類(2).循環(huán)結(jié)構(gòu)循