Matlab作圖一

1.曲線圖


PLOT(X,Y,S):

(命令行輸入help plot可以看到用法,快捷方式f1選中也行)

X,Y是向量,分別表示點(diǎn)集的橫坐標(biāo)和縱坐標(biāo)

PLOT(X,Y)//畫實(shí)線

同一個(gè)圖畫多個(gè)曲線:

PLOT(X,Y1,S1(第一條曲線),X,Y2,S2,……,X,Yn,Sn)//將多條線畫在一起

例:在[0,2*pi]用紅線畫sin(x),用綠圈畫cos(x).

x=linspace(0,2*pi,30);//起點(diǎn)是0,終點(diǎn)是2pi,中間均勻取30個(gè)點(diǎn)

y=sin(x);

z=cos(x);

plot(x,y,'r',x,z,’go'(綠色的圈畫出cosx曲線))



2.符號(hào)函數(shù)(顯函數(shù)、隱函數(shù)和參數(shù)方程)畫圖


(1) ezplot

ezplot(‘f(x)’,[a,b])

? ? 表示在a<x<b繪制顯函數(shù)f=f(x)的函數(shù)圖

ezplot(‘f(x,y)’,[xmin,xmax,ymin,ymax])

? ? 表示在區(qū)間xmin<x<xmax和 ymin<y<ymax繪制隱函數(shù)f(x,y)=0的函數(shù)圖

ezplot(‘x(t)’,’y(t)’,[tmin,tmax])

? ? ? 表示在區(qū)間tmin<t<tmax繪制參數(shù)方程 x=x(t),y=y(t)的函數(shù)圖

例:在[0,pi]上畫y=cos(x)的圖形

輸入命令:ezplot(‘sin(x)’,[0,pi])


輸入命令:ezplot(‘cos(t)^3’,’sin(t)^3’,[0,2*pi])


輸入命令:ezplot('exp(x)+sin(x*y)'(右邊等于0可以省略不寫,'exp(x)=-sin(x*y)'也行),[-2,0.5,0,2])

(2) fplot


教材上的fplot用在程序會(huì)有問(wèn)題

fplot(‘fun’,lims)

? 表示繪制字符串fun指定的函數(shù)在lims=[xmin,xmax]的圖形.

注意:

[1] fun必須是M文件的函數(shù)名或是獨(dú)立變量為x的字符串.?

[2] fplot函數(shù)不能畫參數(shù)方程和隱函數(shù)圖形,但在一個(gè)圖上可以畫多個(gè)圖形。


先建M文件myfun1.m:

? function y=myfun1(x)

? ? ? y=exp(2*x)+sin(3*x^2)

再輸入命令:

fplot(@(myfun1)myfun1,[-1,2])

例:在[-2,2]范圍內(nèi)繪制函數(shù)tanh的圖形

解 fplot(@(x) tanh(x),[-2,2])


輸入命令:

fplot(@(x)[tanh(x),sin(x),cos(x)],2*pi*[-1,1])

3. 對(duì)數(shù)坐標(biāo)圖

loglog(Y) 表示 x、y坐標(biāo)都是對(duì)數(shù)坐標(biāo)系

semilogx(Y)? ? ? 表示 x坐標(biāo)軸是對(duì)數(shù)坐標(biāo)系

semilogy(…)? ? ? 表示y坐標(biāo)軸是對(duì)數(shù)坐標(biāo)系

plotyy? ? 有兩個(gè)y坐標(biāo)軸,一個(gè)在左邊,一個(gè)在右邊

例:用方形標(biāo)記創(chuàng)建一個(gè)簡(jiǎn)單的loglog

輸入命令:

? x=logspace(-1,2);//默認(rèn)取50個(gè)點(diǎn)

? loglog(x,exp(x),'-s')

? grid on

例:創(chuàng)建一個(gè)簡(jiǎn)單的半對(duì)數(shù)坐標(biāo)圖

?輸入命令:

? ? x=0:.1:10;

? ? semilogy(x,10.^x)

//subplot(2,3,x)產(chǎn)生兩行三列6個(gè)圖,x代表第幾個(gè)圖

4.空 間 曲 線

1、 一條曲線

PLOT3(x,y,z(xyz是n維向量,分別表示曲線上點(diǎn)集的橫坐標(biāo)、縱坐標(biāo)、函數(shù)值),s(指定顏色、線形等))

例 在區(qū)間[0,10*pi]畫出參數(shù)曲線x=sin(t),y=cos(t),

? ? ? z=t.

t=0:pi/50:10*pi;

? ? ? plot3(sin(t),cos(t),t)

? ? ? rotate3d? %旋轉(zhuǎn),有501個(gè)點(diǎn)

2、多條曲線

PLOT3(x,y,z)
其中x,y,z是都是m*n矩陣,其對(duì)應(yīng)的每一列表示一條曲線.

例 畫多條曲線觀察函數(shù)Z=(X+Y).^2.

x=-3:0.1:3;y=1:0.1:5;

? [X,Y]=meshgrid(x,y);//返回2個(gè)矩陣

Z=(X+Y).^2;

plot3(X,Y,Z)

(這里meshgrid(x,y)的作用是產(chǎn)生一個(gè)以向量x為行、向量y為列的矩陣)

5.空 間 曲 面

(1) surf(x,y,z)?畫出數(shù)據(jù)點(diǎn)(x,y,z)表示的曲面

數(shù)據(jù)矩陣。分別表示數(shù)據(jù)點(diǎn)的橫坐標(biāo)、縱坐標(biāo)、函數(shù)值

例:畫函數(shù)Z=(X+Y).^2的圖形.

解? x=-3:0.1:3;

? ? y=1:0.1:5;

? ? [X,Y]=meshgrid(x,y);

? ? Z=(X+Y).^2;

? ? surf(X,Y,Z)

? ? shading? flat? ? %將當(dāng)前圖形變得平滑

(2) Mesh(x,y,z) 畫網(wǎng)格曲面

例 畫出曲面Z=(X+Y).^2在不同視角的網(wǎng)格圖.

解 x=-3:0.1:3; y=1:0.1:5;

? ? [X,Y]=meshgrid(x,y);

? ? Z=(X+Y).^2;

? ? mesh(X,Y,Z)?

(3)meshz(X,Y,Z) 在網(wǎng)格周圍畫一個(gè)curtain圖(如,參考平面)?

例 繪peaks的網(wǎng)格圖

解 輸入命令:

? [X,Y]=meshgrid(-3:.125:3);//步長(zhǎng)為0.125

? Z=peaks(X,Y);

? Meshz(X,Y,Z)

6.處理圖形

1、在圖形上加格柵、圖例和標(biāo)注

(1)GRID ON: 加格柵在當(dāng)前圖上

? ? ? ? GRID OFF: 刪除格柵

(2)hh = xlabel(string): 在當(dāng)前圖形的x軸上加圖例string

hh = ylabel(string): 在當(dāng)前圖形的y軸上加圖例string

hh = zlabel(string): 在當(dāng)前圖形的z軸上加圖例string

hh = title(string): 在當(dāng)前圖形的頂端上加圖例string

例 在區(qū)間[0,2*pi]畫sin(x)的圖形,并加注圖例“自變量

? ? X”、“函數(shù)Y”、“示意圖”, 并加格柵.

解 x=linspace(0,2*pi,30);

y=sin(x);

plot(x,y)

xlabel('自變量X')

ylabel('函數(shù)Y')

title('示意圖')

grid on


(3) hh = gtext(‘string’)

命令gtext(‘string’)用鼠標(biāo)放置標(biāo)注在現(xiàn)有的圖上.運(yùn)行命令gtext(‘string’)時(shí),屏幕上出現(xiàn)當(dāng)前圖形,在

圖形上出現(xiàn)一個(gè)交叉的十字,該十字隨鼠標(biāo)的移動(dòng)移動(dòng),當(dāng)按下鼠標(biāo)左鍵時(shí),該標(biāo)注string放在當(dāng)前十交叉的位置.

例 在區(qū)間[0,2*pi]畫sin(x),并分別標(biāo)注“sin(x)”

? ? ”cos(x)”.

解 x=linspace(0,2*pi,30);

y=sin(x);

z=cos(x);

plot(x,y,x,z)

gtext(‘sin(x)’);gtext(’cos(x)’)

2、定制坐標(biāo)

axis([xmin xmax ymin ymax zmin zmax])

axis auto? ??將坐標(biāo)軸返回到自動(dòng)缺省值

例 在區(qū)間[0.005,0.01]顯示sin(1/x)的圖形。

解 x=linspace(0.0001,0.01,1000);

y=sin(1./x);

plot(x,y)

axis([0.005 0.01 –1 1])

3、圖形保持

(1) hold on? ?保持當(dāng)前圖形, 以便繼續(xù)畫圖到當(dāng)前圖上

? ? ? hold? off? ?釋放當(dāng)前圖形窗口

例 將y=sin(x)、y=cos(x)分別用點(diǎn)和線畫出在同一屏幕上。

解 x=linspace(0,2*pi,30);

y=sin(x);

z=cos(x)

plot(x,z,’:’)

hold on

plot(x,y)

(2) figure(h)新建h窗口,激活圖形使其可見,并把它置于其它圖形之上


例 區(qū)間[0,2*pi]新建兩個(gè)窗口分別畫出y=sin(x);

? ? z=cos(x)。

解 x=linspace(0,2*pi,100);

? ? y=sin(x);z=cos(x);

? ? plot(x,y);

? ? title('sin(x)');

? ? pause//暫停,新建一個(gè)窗口畫圖

figure(2);

? ? plot(x,z);

? ? title('cos(x)');

4、分割窗口

h=subplot(mrows,ncols,thisplot)

劃分整個(gè)作圖區(qū)域?yàn)閙rows*ncols塊(逐行對(duì)塊訪問(wèn))并激活第thisplot塊,其后的作圖語(yǔ)句將圖形畫在該塊上。

subplot(mrows,ncols,thisplot)

激活已劃分為mrows*ncols塊的屏幕中的第thisplot塊,其后的作圖語(yǔ)句將圖形畫在該塊上。

subplot(1,1,1)

命令subplot(1,1,1)返回非分割狀態(tài)。

例 將屏幕分割為四塊,并分別畫出y=sin(x),z=cos(x),a=sin(x)*cos(x),b=sin(x)/cos(x)。

解x=linspace(0,2*pi,100);

y=sin(x); z=cos(x);?

a=sin(x).*cos(x);b=sin(x)./(cos(x)+eps)

subplot(2,2,1);plot(x,y),title(‘sin(x)’)

subplot(2,2,2);plot(x,z),title(‘cos(x)’)

subplot(2,2,3);plot(x,a),title(‘sin(x)cos(x)’)

subplot(2,2,4);plot(x,b),title(‘sin(x)/cos(x)’)

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡(jiǎn)書系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

友情鏈接更多精彩內(nèi)容