一、散點圖

1．1．命令??plot

功能??線性二維圖。在線條多于一條時，若用戶沒有指定使用顏色，則plot循環(huán)使用由當前坐標軸顏色順序屬性（current axes ColorOrder property）定義的顏色，以區(qū)別不同的線條。在用完上述屬性值后，plot又循環(huán)使用由坐標軸線型順序屬性（axes LineStyleOrder property）定義的線型，以區(qū)別不同的線條。

用法??plot(X,Y)?當X,Y均為實數向量，且為同維向量（可以不是同型向量），X=[x(i)]，Y=[y(i)]，則plot(X,Y)先描出點(x(i)，y(i))，然后用直線依次相連；若X，Y為復數向量，則不考慮虛數部分。若X，Y均為同維同型實數矩陣，X = [X(i)]，Y = [Y(i)]，其中X(i),Y(i)為列向量，則plot(X,Y)依次畫出plot(X(i),Y(i))，矩陣有幾列就有幾條線；若X，Y中一個為向量，另一個為矩陣，且向量的維數等于矩陣的行數或者列數，則矩陣按向量的方向分解成幾個向量，再與向量配對分別畫出，矩陣可分解成幾個向量就有幾條線；在上述的幾種使用形式中，若有復數出現，則復數的虛數部分將不被考慮。

plot(Y)?若Y為實數向量，Y的維數為m，則plot(Y)等價于plot(X,Y)，其中x=1：m；若y為實數矩陣，則把y按列的方向分解成幾個列向量，而y?的行數為n，則plot(Y)等價于plot(X,Y)其中x=[1;2;…;n]；在上述的幾種使用形式中，若有復數出現，則復數的虛數部分將不被考慮。

plot(X1,Y1,X2,Y2,…)，其中Xi與Yi成對出現，plot(X1,Y1,X2,Y2,…)將分別按順序取兩數據Xi與Yi進行畫圖。若其中僅僅有Xi或Yi是矩陣，其余的為向量，向量維數與矩陣的維數匹配，則按匹配的方向來分解矩陣，再分別將配對的向量畫出。

plot(X1,Y1,LineSpec1,X2,Y2,LineSpec2…)?將按順序分別畫出由三參數定義Xi,Yi,LineSpeci的線條。其中參數LineSpeci指明了線條的類型，標記符號，和畫線用的顏色。在plot?命令中我們可以混合使用三參數和二參數的形式：

plot(X1,Y1,LineSpec1,X2,Y2,X3,Y3,LineSpec3)

plot(…,'PropertyName',PropertyValue,…)?對所有的用plot生成的line圖形對象中指定的屬性進行恰當的設置。

h = plot(…)?返回line圖形對象句柄的一列向量，一線條對應一句柄值。

說明??參數LineSpec

功能??定義線的屬性。Maltab允許用戶對線條定義如下的特性：

1．線型

表7-1

定義符	-	--	：	-.
線型	實線（缺省值）	劃線	點線	點劃線

2．線條寬度

指定線條的寬度，取值為整數（單位為像素點）

3．顏色

表7-2

定義符	R（red）	G(green)	b(blue)	c(cyan)
顏色	紅色	綠色	蘭色	青色
定義符	M(magenta)	y(yellow)	k(black)	w(white)
顏色	品紅	黃色	黑色	白色

4．標記類型

表7-3

定義符	+	o(字母)	*	.	x
標記類型	加號	小圓圈	星號	實點	交叉號
定義符	d	^	v	>	<
標記類型	棱形	向上三角形	向下三角形	向右三角形	向左三角形
定義符	s	h	P	?	?
標記類型	正方形	正六角星	正五角星	?	?

5．標記大小

指定標記符號的大小尺寸，取值為整數（單位為像素）

6．標記面填充顏色

指定用于填充標記符面的顏色。取值在上表。

7．標記周邊顏色

指定標記符顏色或者是標記符（小圓圈、正方形、棱形、正五角星、正六角星和四個方向的三角形）周邊線條的顏色。取值在上表。

在所有的能產生線條的命令中，參數LineSepc可以定義線條的下面三個屬性：線型、標記符號、顏色進行設置。對線條的上述屬性的定義可用字符串來定義，如：plot(x,y,'-.or')

結合x和y，畫出點劃線（-.），在數據點（x，y）處畫出小圓圈（o），線和標記都用紅色畫出。其中定義符（即字符串）中的字母、符號可任意組合。若沒有定義符，則畫圖命令plot自動用缺省值進行畫圖。若僅僅指定了標記符，而非線型，則plot只在數據點畫出標記符。

1.基本畫圖

程序如下：

view plaincopy to clipboardprint?

x=0:pi/1000:2*pi;??

y1=sin(2*x);??

y2=2*cos(2*x);??

%輸出圖像??

plot(x,y1,'k-',x,y2,'b--');??

title('?Plot?of?f(x)=sin(2x)?and?its?derivative');??

%設置X坐標和Y坐標的標簽??

xlabel('x');??

ylabel('y');??

%制作圖例??

legend('f(x)=sin(2x)','d/dx?f(x)')??

%顯示網格??

grid?on;??

x=0:pi/1000:2*pi; y1=sin(2*x); y2=2*cos(2*x); %輸出圖像 plot(x,y1,'k-',x,y2,'b--'); title(' Plot of f(x)=sin(2x) and its derivative'); %設置X坐標和Y坐標的標簽 xlabel('x'); ylabel('y'); %制作圖例 legend('f(x)=sin(2x)','d/dx f(x)') %顯示網格 grid on;

顯示結果：

2.利用有限個點畫出平滑曲線

view plaincopy to clipboardprint?

x=1:1:10;??

y=[1,4,8,10,11,11.5,12,7,5,1];??

xx=1:0.01:10;??

%使用了函數interp1??

yy=interp1(x,y,xx,'cublic');??

plot(xx,yy,x,y,'.');??

grid?on;??

x=1:1:10; y=[1,4,8,10,11,11.5,12,7,5,1]; xx=1:0.01:10; %使用了函數interp1 yy=interp1(x,y,xx,'cublic'); plot(xx,yy,x,y,'.'); grid on;?

3.繪制三維曲線

view plaincopy to clipboardprint?

t=0:pi/100:20*pi;??

x=sin(t);??

y=cos(t);??

z=t.*sin(t).*cos(t);??

plot3(x,y,z);??

title('Line?in?3-D?Space');??

xlabel('X');ylabel('Y');zlabel('Z');??

grid?on;??

t=0:pi/100:20*pi; x=sin(t); y=cos(t); z=t.*sin(t).*cos(t); plot3(x,y,z); title('Line in 3-D Space'); xlabel('X');ylabel('Y');zlabel('Z'); grid on;

4.繪制三維曲面

view plaincopy to clipboardprint?

[x,y]=meshgrid(-8:0.5:8);??

z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);??

subplot(2,2,1);??

mesh(x,y,z);??

title('mesh(x,y,z)')??

subplot(2,2,2);??

meshc(x,y,z);??

title('meshc(x,y,z)')??

subplot(2,2,3);??

meshz(x,y,z)??

title('meshz(x,y,z)')??

subplot(2,2,4);??

surf(x,y,z);??

title('surf(x,y,z)')??

[x,y]=meshgrid(-8:0.5:8); z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps); subplot(2,2,1); mesh(x,y,z); title('mesh(x,y,z)') subplot(2,2,2); meshc(x,y,z); title('meshc(x,y,z)') subplot(2,2,3); meshz(x,y,z) title('meshz(x,y,z)') subplot(2,2,4); surf(x,y,z); title('surf(x,y,z)')

標準三維曲面：

view plaincopy to clipboardprint?

t=0:pi/20:2*pi;??

[x,y,z]=?cylinder(2+sin(t),30);??

subplot(2,2,1);??

surf(x,y,z);??

subplot(2,2,2);??

[x,y,z]=sphere;??

surf(x,y,z);??

subplot(2,1,2);??

[x,y,z]=peaks(30);???

surf(x,y,z);??

t=0:pi/20:2*pi; [x,y,z]= cylinder(2+sin(t),30); subplot(2,2,1); surf(x,y,z); subplot(2,2,2); [x,y,z]=sphere; surf(x,y,z); subplot(2,1,2); [x,y,z]=peaks(30); surf(x,y,z);

其他函數：

view plaincopy to clipboardprint?

subplot(2,2,1);??

bar3(magic(4))??

subplot(2,2,2);??

y=2*sin(0:pi/10:2*pi);??

stem3(y);??

subplot(2,2,3);??

pie3([2347,1827,2043,3025]);??

subplot(2,2,4);??

fill3(rand(3,5),rand(3,5),rand(3,5),?'y'?)??

subplot(2,2,1); bar3(magic(4)) subplot(2,2,2); y=2*sin(0:pi/10:2*pi); stem3(y); subplot(2,2,3); pie3([2347,1827,2043,3025]); subplot(2,2,4); fill3(rand(3,5),rand(3,5),rand(3,5), 'y' )

繪制多峰函數的瀑布圖和等高線圖：

view plaincopy to clipboardprint?

subplot(1,2,1);??

[X,Y,Z]=peaks(30);??

waterfall(X,Y,Z)??

xlabel('X-axis'),ylabel('Y-axis'),zlabel('Z-axis');??

subplot(1,2,2);??

contour3(X,Y,Z,12,'k');?????%其中12代表高度的等級數??

xlabel('X-axis'),ylabel('Y-axis'),zlabel('Z-axis');??

subplot(1,2,1); [X,Y,Z]=peaks(30); waterfall(X,Y,Z) xlabel('X-axis'),ylabel('Y-axis'),zlabel('Z-axis'); subplot(1,2,2); contour3(X,Y,Z,12,'k'); %其中12代表高度的等級數 xlabel('X-axis'),ylabel('Y-axis'),zlabel('Z-axis');

5.圖形修飾處理

三種圖像著色方式效果顯示：

view plaincopy to clipboardprint?

[x,y,z]=sphere(20);??

colormap(copper);??

subplot(1,3,1);??

surf(x,y,z);??

axis?equal??

subplot(1,3,2);??

surf(x,y,z);shading?flat;??

axis?equal??

subplot(1,3,3);??

surf(x,y,z);shading?interp;??

axis?equal??

[x,y,z]=sphere(20); colormap(copper); subplot(1,3,1); surf(x,y,z); axis equal subplot(1,3,2); surf(x,y,z);shading flat; axis equal subplot(1,3,3); surf(x,y,z);shading interp; axis equal

光照處理后的球面：

view plaincopy to clipboardprint?

[x,y,z]=sphere(20);??

subplot(1,2,1);??

surf(x,y,z);axis?equal;??

light('Posi',[0,1,1]);??

shading?interp;??

hold?on;??

plot3(0,1,1,'p');text(0,1,1,'?light');??

subplot(1,2,2);??

surf(x,y,z);axis?equal;??

light('Posi',[1,0,1]);??

shading?interp;??

hold?on;??

plot3(1,0,1,'p');text(1,0,1,'?light');??

[x,y,z]=sphere(20); subplot(1,2,1); surf(x,y,z);axis equal; light('Posi',[0,1,1]); shading interp; hold on; plot3(0,1,1,'p');text(0,1,1,' light'); subplot(1,2,2); surf(x,y,z);axis equal; light('Posi',[1,0,1]); shading interp; hold on; plot3(1,0,1,'p');text(1,0,1,' light');

圖形的裁剪處理：

view plaincopy to clipboardprint?

[x,y]=meshgrid(-5:0.1:5);??

z=cos(x).*cos(y).*exp(-sqrt(x.^2+y.^2)/4);??

subplot(1,2,1);??

surf(x,y,z);shading?interp;??

title('裁剪之前');??

i=find(x<=0&y<=0);??

z1=z;z1(i)=NaN;??

subplot(1,2,2);??

surf(x,y,z1);shading?interp;??

title('裁剪之后');??

[x,y]=meshgrid(-5:0.1:5); z=cos(x).*cos(y).*exp(-sqrt(x.^2+y.^2)/4); subplot(1,2,1); surf(x,y,z);shading interp; title('裁剪之前'); i=find(x<=0&y<=0); z1=z;z1(i)=NaN; subplot(1,2,2); surf(x,y,z1);shading interp; title('裁剪之后');

從文件載入圖像：

view plaincopy to clipboardprint?

[x,cmap]=imread('flower.jpg');??%讀取圖像的數據陣和色圖陣??

image(x);colormap(cmap);??

axis?image?off????%保持寬高比并取消坐標軸???

[x,cmap]=imread('flower.jpg'); %讀取圖像的數據陣和色圖陣 image(x);colormap(cmap); axis image off %保持寬高比并取消坐標軸

制作動畫：

view plaincopy to clipboardprint?

[X,Y,Z]=peaks(30);???

surf(X,Y,Z)??

axis([-3,3,-3,3,-10,10])??

axis?off;??

shading?interp;??

colormap(hot);??

m=moviein(20);????????????%建立一個20列大矩陣??

for?i=1:20??

view(-37.5+24*(i-1),30)??????%改變視點??

m(:,i)=getframe;????????????%將圖形保存到m矩陣??

end???

movie(m,2);?????????????????%播放畫面2次??

?

1．2．命令?scatter(x1,y,50,c,'o','filled')

?

二、一元線性回歸

2．1．命令?polyfit最小二乘多項式擬合

?[p，S]=polyfit（x，y，m）

多項式y=a1xm+a2xm-1+…+amx+am+1

其中x=（x1，x2，…，xm）x1…xm為（n*1）的矩陣;

y為（n*1）的矩陣；

p=（a1，a2，…，am+1）是多項式y=a1xm+a2xm-1+…+amx+am+1的系數；

S是一個矩陣，用來估計預測誤差.

2．2．命令?polyval多項式函數的預測值

Y=polyval（p，x）求polyfit所得的回歸多項式在x處的預測值Y；

p是polyfit函數的返回值；

x和polyfit函數的x值相同。

2．3．命令?polyconf?殘差個案次序圖

[Y，DELTA]=polyconf（p，x，S，alpha）求polyfit所得的回歸多項式在x處的預測值Y及預測值的顯著性為1-alpha的置信區(qū)間DELTA；alpha缺省時為0.05。

p是polyfit函數的返回值；

x和polyfit函數的x值相同；

S和polyfit函數的S值相同。

2．4?命令?polytool（x，y，m）一元多項式回歸命令

?

2．5．命令regress多元線性回歸（可用于一元線性回歸）

b=regress( Y,??X )

[b, bint,r,rint,stats]=regress(Y,X,alpha)

b?回歸系數

bint?回歸系數的區(qū)間估計

r?殘差

rint?殘差置信區(qū)間

stats?用于檢驗回歸模型的統計量，有三個數值：相關系數R²、F值、與F對應的概率p,相關系數R²越接近1，說明回歸方程越顯著；F > F1-α（k，n-k-1）時拒絕H₀，F越大，說明回歸方程越顯著；與F對應的概率p?時拒絕H₀，回歸模型成立。

Y為n*1的矩陣；

X為（ones(n,1),x1,…,xm）的矩陣；

alpha顯著性水平（缺省時為0.05）。

?

三、多元線性回歸

3．1．命令?regress（見2。5）

3．2．命令?rstool?多元二項式回歸

命令：rstool（x，y，’model’, alpha）

x?為n*m矩陣

y為?n維列向量

model?由下列4個模型中選擇1個（用字符串輸入，缺省時為線性模型）：

linear（線性）：

purequadratic（純二次）：?

interaction（交叉）：

quadratic（完全二次）：

alpha?顯著性水平（缺省時為0.05）

返回值beta?系數

返回值rmse剩余標準差

返回值residuals殘差

?

四、非線性回歸

4．1．命令?nlinfit

[beta,R,J]=nlinfit(X,Y,’’model’,beta0)

X?為n*m矩陣

Y為?n維列向量

model為自定義函數

beta0為估計的模型系數

beta為回歸系數

R為殘差

J

?

4．2．命令?nlintool

nlintool(X,Y,’model’,beta0,alpha)

X?為n*m矩陣

Y為?n維列向量

model為自定義函數

beta0為估計的模型系數

alpha顯著性水平（缺省時為0.05）

?

4．3．命令?nlparci

betaci=nlparci(beta,R,J)

beta為回歸系數

R為殘差

J

返回值為回歸系數beta的置信區(qū)間

?

4．4．命令?nlpredci

[Y,DELTA]=nlpredci(‘model’,X,beta,R,J)

Y為預測值

DELTA為預測值的顯著性為1-alpha的置信區(qū)間；alpha缺省時為0.05。

X?為n*m矩陣

model為自定義函數

beta為回歸系數

R為殘差

J

?

五、其它

命令?grid on

命令?axis（坐標軸）([0??60??0??0.025])

命令?figure?彈出新的畫圖窗口

命令?獲取矩陣的某行某列

x(n,:);?獲取矩陣的第n行

x(:,n);?獲取矩陣的第n列

命令?rcoplot

畫出殘差及其置信區(qū)間：rcoplot（r，rint）

glmfit?一般線性模型擬合??

regstats?回歸統計量診斷?

regstats(responses,DATA,'model')'

stats = regstats(responses,DATA,'model','whichstats')

'Q'Q from the QR Decomposition of X

'R'R from the QR Decomposition of X

'beta'Regression Coefficients 'covb'Covariance of Regression Coefficients

'yhat'Fitted Values of the Response Data

'r'Residuals

'mse'Mean Squared Error

'leverage'Leverage 'hatmat'Hat (Projection) Matrix 's2_i'Delete-1 Variance 'beta_i'Delete-1 Coefficients 'standres'Standardized Residuals 'studres'Studentized Residuals 'dfbetas'Scaled Change in Regression Coefficients 'dffit'Change in Fitted Values 'dffits'Scaled Change in Fitted Values 'covratio'Change in Covariance 'cookd'Cook's Distance 'all'Create all of the above statistics

命令?bar(條圖)

命令?pie(餅圖)

命令?hist(直方圖)

命令?help

命令?mean(平均值)

命令?inv(逆概率分布)

命令?pdf(密度)

命令?cdf(分布函數)

命令?stat(均差與方差)

命令rnd(隨機函數)

命令?std(標準差)

命令?var(方差)

命令?median(中位數)

命令?skewness(偏度)

命令?kurtosis(峰度)

命令?norm（正態(tài)分布）

命令?t（t分布）

命令?f（f分布）

命令?chr2（x2分布）

命令?poiis（泊松分布）

matlab作圖里面如何分別設置雙縱坐標的刻度?
需要設置y軸的刻度，用到以下函數，?
set(gca,'XTick',[0:5:100]) ? ? ? ? ? ? ? ? ? ? ? ? ?%設置x軸間隔?

set(gca,'yTick',[0:10:350]) ? ? ? ? ? ? ? ? ? ? ? ? %設置y軸間隔 ?
雙縱坐標的標注已實現?
[AX]=plotyy(x1,y1,x1,y2);?
set(get(AX(1),'Ylabel'),'string','left Y-axis‘); ? ? ? ? ? %標注左邊y軸?

set(get(AX(2),'Ylabel'),'string','right y-axis'); ? ? ? ? ?%標注右邊y軸 ? ?
雙坐標畫圖實例?
x = 0:0.01:20;?
y1 = 200*exp(-0.05*x).*sin(x); y2 = 0.8*exp(-0.5*x).*sin(10*x); [AX]=plotyy(x1,y1,x1,y2);?
得到兩個axes句柄，AX(1)和AX(2)?
set(AX(1),'yTick',[0:10:350]) ? ? ? ? ? ? ? ? ?%設置左邊Y軸的刻度??

set(AX(2),'yTick',[0:10:350]) ? ? ? ? ? ? ? ? %設置右邊Y軸的刻度?

[AX,H1,H2] = plotyy(x,y1,x,y2,'plot');?

set(AX(1),'XColor','k','YColor','b'); ? ? ? ? ? ? ?%設置左坐標軸顏色?

set(AX(2),'XColor','k','YColor','r');?
HH1=get(AX(1),'Ylabel'); ? ? ? ? ? ? ? ? ? ? %得到左邊坐標軸句柄?

set(HH1,'String','Left Y-axis'); ? ? ? ? ? ? ? ? ?%標注左邊坐標軸?

set(HH1,'color','b'); ? ? ? ? ? ? ? ? ? ? ? ? ??
HH2=get(AX(2),'Ylabel'); ? ? ? ? ? ??
set(HH2,'String','Right Y-axis'); ? ? ? ? ? ? ? ? %標注左邊坐標軸 set(HH2,'color','r');?
set(H1,'LineStyle','-'); ? ? ? ? ? ? ? ? ? ? ? ? %設置曲線特性 set(H1,'color','b');?












set(H2,'LineStyle',':'); set(H2,'color','r');?

總結

以上是生活随笔為你收集整理的MATLAB画图命令zz的全部內容，希望文章能夠幫你解決所遇到的問題。

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