MATLAB語言基礎(chǔ)(濟(jì)大)

1、 計(jì)算機(jī)仿真技術(shù) 濟(jì)南大學(xué)專版 MATLAB語言基礎(chǔ)語言基礎(chǔ)計(jì)算機(jī)仿真技術(shù)濟(jì)南大學(xué)自動(dòng)化與電氣工程學(xué)院授課教師：謝錦麗E-mail：cse_聯(lián)系方式ecember 13, 20212濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)第5章 MATLAB語言基礎(chǔ)o 5.1 MATLAB概述概述o 5.2 MATLAB的安裝與啟動(dòng)的安裝與啟動(dòng)o 5.3 MATLAB集成環(huán)境集成環(huán)境o 5.4 MATLAB幫助系統(tǒng)幫助系統(tǒng)o 5.5 基本運(yùn)算基本運(yùn)算o 5.6 矩陣函數(shù)矩陣函數(shù)o 5.7 繪圖功能繪圖功能December 13, 20213濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)

2、仿真技術(shù)5.5 基本運(yùn)算基本運(yùn)算o書寫向量書寫向量(vectors)與矩陣與矩陣(matrices)nx=1 2 3 -4 -5 或或 x=1,2,3,-4,-5 %表示行表示行向量向量nx=1;2;3;-4;-5%用分號(hào)表示列向量用分號(hào)表示列向量n輸入矩陣：輸入矩陣：A=1.2 10 15; 3 5.5 2; 4 6.8 7n如果在矩陣如果在矩陣A后面加上分號(hào)后面加上分號(hào); 取消顯示結(jié)果取消顯示結(jié)果nA(i,:) 第第i行，行，A(:,j) 第第j列列December 13, 20214濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.5 基本運(yùn)算基本運(yùn)算nt=5:-1:2 %產(chǎn)生行向量，從

3、產(chǎn)生行向量，從5到到2間隔為間隔為-1 ny=x %表示表示x的反轉(zhuǎn)矩陣為的反轉(zhuǎn)矩陣為y n 寫成寫成 C=1 exp(-0.02); sqrt(2) 3December 13, 20215濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)0.02123eC5.5 基本運(yùn)算基本運(yùn)算n定數(shù)線性采樣定數(shù)線性采樣linespace產(chǎn)生向量：產(chǎn)生向量：x=linespace(n1,n2,n) 從從n1到到n2有有n個(gè)點(diǎn)個(gè)點(diǎn)n例：例：x=linspace(-10,10,5) %從從-10到到10有有5個(gè)點(diǎn)個(gè)點(diǎn)n定數(shù)對(duì)數(shù)采樣定數(shù)對(duì)數(shù)采樣logspace產(chǎn)生向量產(chǎn)生向量: x=logspace(n1,n2,n)

4、 從從10n1到到10n2有有n各點(diǎn)各點(diǎn)n例：例：x=logspace(-1,1,5) %從從0.1到到10有有5個(gè)點(diǎn)個(gè)點(diǎn)December 13, 20216濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.5 基本運(yùn)算基本運(yùn)算o轉(zhuǎn)置與共軛轉(zhuǎn)置轉(zhuǎn)置與共軛轉(zhuǎn)置(transpose and conjugate transpose)n實(shí)矩陣情況：實(shí)矩陣情況： A=1 2 3;4 5 6;7 8 9; B=An復(fù)數(shù)情況：復(fù)數(shù)情況：x=1+sqrt(3)*i 或或 x=1+sqrt(3)*j 用用i或或j表示復(fù)數(shù)表示復(fù)數(shù)n如果如果i和和j已經(jīng)作為變量使用，定義新的已經(jīng)作為變量使用，定義新的復(fù)數(shù)單位：復(fù)數(shù)

5、單位：ii=sqrt(-1) or jj=sqrt(-1)x=1+sqrt(3)*ii or x=1+sqrt(3)*jj n書寫復(fù)數(shù)矩陣書寫復(fù)數(shù)矩陣(complex matrices)： 寫成寫成 X=1 j;-5*j 2December 13, 20217濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)152jXj5.5 基本運(yùn)算基本運(yùn)算n非共軛轉(zhuǎn)置矩陣為：非共軛轉(zhuǎn)置矩陣為：Y=X. 或或 Y=conj(X) n共軛轉(zhuǎn)置矩陣：共軛轉(zhuǎn)置矩陣：Y=XDecember 13, 20218濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.5 基本運(yùn)算基本運(yùn)算o矩陣的旋轉(zhuǎn)矩陣的旋轉(zhuǎn)nrot90(A)

6、%將矩陣將矩陣A旋轉(zhuǎn)旋轉(zhuǎn)90度度o隨機(jī)數(shù)矩陣的產(chǎn)生隨機(jī)數(shù)矩陣的產(chǎn)生nrand(2,2) %產(chǎn)生區(qū)間（產(chǎn)生區(qū)間（0,1）的均勻分）的均勻分布的隨機(jī)數(shù)。布的隨機(jī)數(shù)。nrandn(2,2) %產(chǎn)生正態(tài)分布的隨機(jī)數(shù)矩產(chǎn)生正態(tài)分布的隨機(jī)數(shù)矩陣的函數(shù)陣的函數(shù)December 13, 20219濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.5 基本運(yùn)算基本運(yùn)算o特殊矩陣產(chǎn)生特殊矩陣產(chǎn)生nones(n) 產(chǎn)生產(chǎn)生n乘乘n矩陣，元素值均為矩陣，元素值均為1.nones(m,n) 產(chǎn)生產(chǎn)生m乘乘n矩陣，元素值均為矩陣，元素值均為1.nzeros(n) 產(chǎn)生產(chǎn)生n乘乘n零矩陣零矩陣nzeros(m,n) 產(chǎn)生產(chǎn)

7、生m乘乘n零矩陣零矩陣 o單位矩陣單位矩陣(Identity matrix)neye(n) 產(chǎn)生產(chǎn)生n乘乘n單位矩陣，即對(duì)角線為單位矩陣，即對(duì)角線為1，其它為，其它為0December 13, 202110濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.5 基本運(yùn)算基本運(yùn)算o對(duì)角矩陣對(duì)角矩陣(Diagonal matrix)ndiag(x) 產(chǎn)生對(duì)角矩陣，產(chǎn)生對(duì)角矩陣，x為對(duì)角線元素。比如：為對(duì)角線元素。比如：x=ones(1,4) %產(chǎn)生產(chǎn)生14的行向量，元素值均為的行向量，元素值均為1diag(x) %產(chǎn)生產(chǎn)生44對(duì)角矩陣，對(duì)角線元素為對(duì)角矩陣，對(duì)角線元素為1n如果如果x是是nn矩陣，矩

8、陣，diag(x)是包含是包含x矩陣對(duì)角線元矩陣對(duì)角線元素的列向量：素的列向量：A=1 2 3;4 5 6;7 8 9 %產(chǎn)生產(chǎn)生A矩陣矩陣diag(A) %求求A的對(duì)角線向量的對(duì)角線向量diag( diag(A) )%產(chǎn)生對(duì)角矩陣產(chǎn)生對(duì)角矩陣December 13, 202111濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.5 基本運(yùn)算基本運(yùn)算o對(duì)角矩陣對(duì)角矩陣(Diagonal matrix)ndiag(1:5)o加法與減法加法與減法(addition and subtraction)nA=2 3;4 5;6 7; B=1 0;2 3;0 4; C=A+Bnx=5;4;6; y=x-1

9、o乘法乘法(multiplication)nx=1;2;3; y=4;5;6; z=x*y u=x*yDecember 13, 202112濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.5 基本運(yùn)算基本運(yùn)算o復(fù)數(shù)的幅值與相角復(fù)數(shù)的幅值與相角(magnitude and phase angle of a complex number)n復(fù)數(shù)復(fù)數(shù)z=x+iy=rei 幅值為幅值為r，相角為，相角為，可以表示為可以表示為 r=abs(z) theta=angle(z)o數(shù)組相乘數(shù)組相乘(array multiplication)n矩陣或向量里的元素相乘矩陣或向量里的元素相乘 .* A=1 2 3

10、;0 9 8; B=4 5 6;7 6 5; C=A.*BDecember 13, 202113濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.5 基本運(yùn)算基本運(yùn)算o數(shù)組的冪數(shù)組的冪對(duì)矩陣或向量的元素求冪，表示為對(duì)矩陣或向量的元素求冪，表示為 .x=1 2 3; y=x.2A=1+j 2-2*j;3+4*j 5-j; B=A.2o絕對(duì)值絕對(duì)值(absolute values)abs(A)=sqrt( real(A).2 + imag(A).2 )A=2+2*i 1+3*j; 4+5*j 6-j; abs(A)December 13, 202114濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)

11、5.5 基本運(yùn)算基本運(yùn)算o數(shù)組除法數(shù)組除法(Array division)n數(shù)組相除為數(shù)組相除為 左除左除 ./ 與與 右除右除. nA除以除以B A./B B除以除以A A.Bx=1 2 3; y=4 5 6; u=x./y v=x.yA=1 2 3;1 9 8; B=4 5 6;7 6 5; C=A./B D=A.Bn分母為分母為0時(shí)，時(shí)，MATLAB會(huì)給出警告：會(huì)給出警告：5/0 0/0December 13, 202115濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)第5章 MATLAB語言基礎(chǔ)o 5.1 MATLAB概述概述o 5.2 MATLAB的安裝與啟動(dòng)的安裝與啟動(dòng)o 5.3

12、MATLAB集成環(huán)境集成環(huán)境o 5.4 MATLAB幫助系統(tǒng)幫助系統(tǒng)o 5.5 基本運(yùn)算基本運(yùn)算o 5.6 矩陣函數(shù)矩陣函數(shù)o 5.7 繪圖功能繪圖功能December 13, 202116濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.6 矩陣函數(shù)矩陣函數(shù)o特征值與特征向量特征值與特征向量(Eigenvalues and eigenvectors)nA是一個(gè)是一個(gè)nn矩陣，具有矩陣，具有n個(gè)特征值滿足個(gè)特征值滿足 Ax=x特征值特征值在在MATLAB中的計(jì)算函數(shù)為中的計(jì)算函數(shù)為 eig(A) %產(chǎn)生一個(gè)列向量表示特征值產(chǎn)生一個(gè)列向量表示特征值如果如果A是實(shí)數(shù)并且對(duì)稱的，那么特征值為實(shí)數(shù)，如

13、果是實(shí)數(shù)并且對(duì)稱的，那么特征值為實(shí)數(shù)，如果A不是對(duì)不是對(duì)稱的，特征值經(jīng)常是復(fù)數(shù)：稱的，特征值經(jīng)常是復(fù)數(shù)：A=0 1 0;-1 0 2;3 0 5; eig(A)%A矩陣，求矩陣，求A的特征值的特征值December 13, 202117濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.6 矩陣函數(shù)矩陣函數(shù)n根據(jù)公式根據(jù)公式 AX=XD這里這里D為矩陣形式的特征值，為矩陣形式的特征值，X為矩陣形式的特征向量為矩陣形式的特征向量A=0 1 0;0 0 1;-6 -11 -6; X D=eig(A)December 13, 202118濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.6 矩陣函數(shù)矩

14、陣函數(shù)o特征方程特征方程(caracteristic equation)n矩陣的特征方程的根即為矩陣的特征值。矩陣的特征方程的根即為矩陣的特征值。A=0 1 0;0 0 1;-6 -11 -6; p=poly(A) %矩陣矩陣A的特征方程，系數(shù)列為行向量的特征方程，系數(shù)列為行向量得到得到 p=1 6 11 6表示特征方程的系數(shù)，特征方程為：表示特征方程的系數(shù)，特征方程為：s3+6s2+11s+6 計(jì)算特征方程計(jì)算特征方程p的根：的根：r=roots(p) 計(jì)算矩陣計(jì)算矩陣A的特征值的特征值lameda=eig(A)兩者是相等的。兩者是相等的。December 13, 202119濟(jì)南大學(xué) 自動(dòng)

15、化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.6 矩陣函數(shù)矩陣函數(shù)n已知多項(xiàng)式的根，用函數(shù)已知多項(xiàng)式的根，用函數(shù)poly(r)求多項(xiàng)式。求多項(xiàng)式。r=-3 -2 -1; p=poly(r) %poly函數(shù)內(nèi)為向量，則視為函數(shù)內(nèi)為向量，則視為根，求其多項(xiàng)式根，求其多項(xiàng)式o多項(xiàng)式乘積多項(xiàng)式乘積(卷積卷積, convolution)將兩個(gè)多項(xiàng)式系數(shù)向量相乘將兩個(gè)多項(xiàng)式系數(shù)向量相乘 conv(a,b)a=3 10 25 36 50; b=1 2 10;c=conv(a,b)即即a(s)=3s4+10s3+25s2+36s+50b(s)=s2+2s+10c(s)=(3s4+10s3+25s2+36s+50)*.

16、(s2+2s+10)December 13, 202120濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.6 矩陣函數(shù)矩陣函數(shù)o多項(xiàng)式相除多項(xiàng)式相除(解卷積解卷積, deconvolution)將兩個(gè)多項(xiàng)式系數(shù)向量相除將兩個(gè)多項(xiàng)式系數(shù)向量相除 deconv(a,b)a=3 10 25 36 50; b=1 2 10;q, r=deconv(a,b)q 和和 r 是多項(xiàng)式相除得到的商和余數(shù)是多項(xiàng)式相除得到的商和余數(shù)得到得到 q=3 4 -13, r=0 0 0 22 180即即 a(s)=b(s)*q(s)+r(s)3s4+10s3+25s2+36s+50=(s2+2s+10)*.(3s2+4

17、s-13)+22s+180December 13, 202121濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.6 矩陣函數(shù)矩陣函數(shù)December 13, 202122濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.6 矩陣函數(shù)矩陣函數(shù)o矩陣指數(shù)運(yùn)算矩陣指數(shù)運(yùn)算(matrix exponential)函數(shù)函數(shù)expm(A)指將矩陣指將矩陣A做指數(shù)運(yùn)算：做指數(shù)運(yùn)算：expm(A)=I+A+A2/2!+A3/3!+A=0 1 0;0 0 1;-6 -11 -6;expm(A)expm(eye(3) %33單位矩陣的指數(shù)矩陣單位矩陣的指數(shù)矩陣o矩陣求逆矩陣求逆(inverse of a squ

18、are matrix)A=1 1 2;3 4 0;1 2 5;inv(A)December 13, 202123濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.6 矩陣函數(shù)矩陣函數(shù)o創(chuàng)建三維數(shù)組創(chuàng)建三維數(shù)組n二維數(shù)組第一維稱為二維數(shù)組第一維稱為“行行”，第二維稱為，第二維稱為“列列”，對(duì)，對(duì)于三維數(shù)組，第三維稱為于三維數(shù)組，第三維稱為“頁頁”。MATLAB將三維或?qū)⑷S或者三維以上的數(shù)組統(tǒng)稱為高維數(shù)組。者三維以上的數(shù)組統(tǒng)稱為高維數(shù)組。n使用下標(biāo)引用方法創(chuàng)建三維數(shù)組使用下標(biāo)引用方法創(chuàng)建三維數(shù)組nA(2,2,2)=1; for i = 1:2 for j = 1:2 for k=1:2 A(i,

19、j,k)=i+j+k; end end endDecember 13, 202124濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.6 矩陣函數(shù)矩陣函數(shù)n使用低維數(shù)組創(chuàng)建高維數(shù)組使用低維數(shù)組創(chuàng)建高維數(shù)組nD2=1,2,3;4,5,6;7,8,9; D3(:,:,1)=D2; D3(:,:,2)=2*D2; D3(:,:,3)=3*D2;December 13, 202125濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)第5章 MATLAB語言基礎(chǔ)o 5.1 MATLAB概述概述o 5.2 MATLAB的安裝與啟動(dòng)的安裝與啟動(dòng)o 5.3 MATLAB集成環(huán)境集成環(huán)境o 5.4 MATLAB幫助

20、系統(tǒng)幫助系統(tǒng)o 5.5 基本運(yùn)算基本運(yùn)算o 5.6 矩陣函數(shù)矩陣函數(shù)o 5.7 繪圖功能繪圖功能December 13, 202126濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.7 繪圖功能繪圖功能oMATLAB具有強(qiáng)大的繪圖功能。主要分為二維繪圖具有強(qiáng)大的繪圖功能。主要分為二維繪圖與三維繪圖。要定義不同維度相應(yīng)的坐標(biāo)。與三維繪圖。要定義不同維度相應(yīng)的坐標(biāo)。o二維線型圖形：二維線型圖形：plot . 比如，定義時(shí)間從比如，定義時(shí)間從0到到2，每隔，每隔0.01取一個(gè)時(shí)間點(diǎn)取一個(gè)時(shí)間點(diǎn)，在該時(shí)間序列內(nèi)計(jì)算正弦函數(shù)在該時(shí)間序列內(nèi)計(jì)算正弦函數(shù)y=sin(3*t)的數(shù)值，并的數(shù)值，并以時(shí)間為以時(shí)

21、間為x軸，以函數(shù)值為軸，以函數(shù)值為y軸，畫出時(shí)間響應(yīng)曲線：軸，畫出時(shí)間響應(yīng)曲線：t=0 : 0.01*pi : 2*pi; %生成時(shí)間向量生成時(shí)間向量y=sin(3*t);plot(t,y)o添加網(wǎng)格線，圖形標(biāo)題，添加網(wǎng)格線，圖形標(biāo)題，x軸標(biāo)題，軸標(biāo)題，y軸標(biāo)題：軸標(biāo)題：grid;title(plot of sin(3t);xlabel(t (sec);ylabel(sin(3t);December 13, 202127濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.7 繪圖功能繪圖功能o在圖形上標(biāo)記文字：在圖形上標(biāo)記文字：n text(pi,sin(3*pi),sin3t); %在點(diǎn)在點(diǎn)(

22、,sin(3)出標(biāo)記文字出標(biāo)記文字sin3tn gtext(sin(4) %該命令輸入后，產(chǎn)生該命令輸入后，產(chǎn)生光標(biāo)輸入文字位置光標(biāo)輸入文字位置o限制限制x、y軸區(qū)域軸區(qū)域n 在在MATLAB繪圖中，繪圖中，x軸和軸和y軸的區(qū)域是軸的區(qū)域是自動(dòng)選擇的，有時(shí)需要手動(dòng)選擇區(qū)域自動(dòng)選擇的，有時(shí)需要手動(dòng)選擇區(qū)域. n xlim(0 2*pi); ylim(-2 2); %限定限定x軸和軸和y軸的范圍軸的范圍December 13, 202128濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.7 繪圖功能繪圖功能o選擇圖形形狀選擇圖形形狀n 普通默認(rèn)為實(shí)線圖，即將各點(diǎn)用曲普通默認(rèn)為實(shí)線圖，即將各點(diǎn)用曲

23、線起來，可以選擇其他形狀繪圖線起來，可以選擇其他形狀繪圖.n 線繪圖包括線繪圖包括 “ 實(shí)線實(shí)線”，“- - 短劃短劃線線”，“： 點(diǎn)線點(diǎn)線”，“- . 短劃線與短劃線與點(diǎn)線點(diǎn)線”： plot(t,y,-, t,2*y,-, t,3*y,:, t,4*y,-.)n 點(diǎn)繪圖包括點(diǎn)繪圖包括 “ + ” ，“ o ”，“ * ”，“ x ” 等，詳見等，詳見help plot. plot(t,y,-, t,2*y,-+, t,3*y,:o, t,4*y,-.x)December 13, 202129濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.7 繪圖功能繪圖功能o選擇圖形顏色選擇圖形顏色系統(tǒng)默

24、認(rèn)顏色為藍(lán)色，可以自己選擇圖系統(tǒng)默認(rèn)顏色為藍(lán)色，可以自己選擇圖形顏色：形顏色：redr;greeng; blueb; whitew; yellow y; magenta m(紫色紫色)cyanc;(青色青色)blackkplot(t,y,-r, t,2*y,-+c, t,3*y,:om, t,4*y,-.xg)o選擇線的粗細(xì)大小選擇線的粗細(xì)大小系統(tǒng)默認(rèn)為系統(tǒng)默認(rèn)為0.5，可以自己加粗，可以自己加粗plot(t,y,-r, t,2*y,-+c, t,3*y,:om, t,4*y,-.xg,LineWidth,3)December 13, 202130濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)

25、5.7 繪圖功能繪圖功能o在圖形上顯示希臘字符與符號(hào)在圖形上顯示希臘字符與符號(hào)n 顯示希臘字符：顯示希臘字符： alpha; beta; gamma; Gamma; delta; Delta; omega; Omega; sigma; Sigma; phi; psi; theta; zeta; pi; Pioplot(t,y); text(pi,sin(3*pi),sin(3pi);December 13, 202131濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.7 繪圖功能繪圖功能n 顯示其他數(shù)學(xué)符號(hào)：顯示其他數(shù)學(xué)符號(hào)： infinity leftarrow uparrow right

26、arrow downarrow circ pmplot(t,y); text(pi,sin(3*pi),leftarrowsin(3pi);December 13, 202132濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.7 繪圖功能繪圖功能o變化文字大小與字體變化文字大小與字體plot(t,y); grid; title(plot of sin(3t),Fontsize,20); xlabel(t (sec),Fontangle, italic); ylabel(sin(3t),Fontname, Times New Roman);下標(biāo)：下標(biāo)：y_1 顯示為顯示為y1.上標(biāo)：上標(biāo)：y1

27、顯示為顯示為y1.December 13, 202133濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.7 繪圖功能繪圖功能o復(fù)數(shù)數(shù)據(jù)畫圖復(fù)數(shù)數(shù)據(jù)畫圖n 如果如果z為復(fù)數(shù)，為復(fù)數(shù)，plot(z)自動(dòng)按照實(shí)部和虛部畫圖，自動(dòng)按照實(shí)部和虛部畫圖，即等于即等于plot(real(z),imag(z)n z=2+2*i 3+3*j 4+5*j 6+7j 7+9j 8+6j 9+5j 10+j;plot(z)xlabel(real axis); ylabel(imaginary axis);o極坐標(biāo)圖：極坐標(biāo)圖：polar(theta,rho)n theta 為角度，為角度，rho為半徑。為半徑。x=

28、0: 0.1*pi : pi; %角度為從角度為從0到到，間隔，間隔0.1y=0:1:10;%半徑為從半徑為從0到到10，間隔為，間隔為1.polar(x,y)December 13, 202134濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.7 繪圖功能繪圖功能o對(duì)數(shù)圖(logarithmic plots)：n loglog: x軸和y軸均使用log10數(shù)字變換，即x軸為log10 x,y軸為log10y.x=1:1:100; y=1:10:1000; loglog(x,y)n semilogx: x軸用log10變換，但y軸不變。semilogx(x,y)n semilogy: x軸不

29、變，y軸用log10變換.semilogy(x,y)December 13, 202135濟(jì)南大學(xué) 自動(dòng)化與電氣工程學(xué)院 計(jì)算機(jī)仿真技術(shù)5.7 繪圖功能繪圖功能o使用多圖模式使用多圖模式n 一個(gè)圖形窗口可以分為多個(gè)小圖，每個(gè)小圖單獨(dú)使用一個(gè)圖形窗口可以分為多個(gè)小圖，每個(gè)小圖單獨(dú)使用plot繪圖繪圖subplot(m,n,p) 將窗口分為將窗口分為m乘乘n個(gè)子窗口，個(gè)子窗口，p代表當(dāng)前的子窗口。代表當(dāng)前的子窗口。t=0 : 0.01*pi : 2*pi;y1=sin(t); y2=sin(2*t);y3=sin(t)+sin(2*t); y4=sin(t).*sin(2*t);subplot(2,2,1)plot(t,y1); grid; xlabel(t (sec); ylabel(y_1=sin(t);subplot(2,2,2) plot(t,y2); grid; xlabel(t (sec); ylabel(y_2=sin(2t);subplot(2,2,3) plot(t,y3); grid; xlabel(t (sec); ylabel(y_3=sin(t)+sin(2t);subplot(2,2,4) plot