matlab_toan_tap_13

Chia sẻ: Kata_8 Kata_8 | Ngày: | Loại File: PDF | Số trang:10

Thêm vào BST

Báo xấu

30
lượt xem 5
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Tham khảo tài liệu 'matlab_toan_tap_13', kỹ thuật - công nghệ, điện - điện tử phục vụ nhu cầu học tập, nghiên cứu và làm việc hiệu quả

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: matlab_toan_tap_13

121 Simpo PDF Merge and Split Unregistered Version - http://www.simpopdf.com >> axis on, grid off % turn the axis on, the grid off H×nh 17.10 >>axis ij % turn the plot upside-down H×nh 17.11
122 Simpo PDF Merge and Split Unregistered Version - http://www.simpopdf.com >> axis square equal % give axis two command at once H×nh 17.12 >> axis xy normal % return to the defaults H×nh 17.13
123 Simpo PDF Merge and Split Unregistered Version - http://www.simpopdf.com 17.6 In h×nh §Ó in c¸c h×nh mµ b¹n võa vÏ hoÆc c¸c h×nh trong ch−¬ng tr×nh cña MATLAB mµ b¹n cÇn, b¹n cã thÓ dïng lÖnh in tõ b¶ng chän hoÆc ®¸nh lÖnh in vµo tõ cöa sæ lÖnh: +) In b»ng lÖnh tõ b¶ng chän: Tr−íc tiªn ta ph¶i chän cöa sæ h×nh lµ cö sæ ho¹t ®éng b»ng c¸ch nhÊn chuét lªn nã, sau ®ã b¹n chän môc b¶ng chän Print tõ b¶ng chän file. Dïng c¸c th«ng sè t¹o lªn trong môc b¶ng chän Print Setup hoÆc Page Setup, ®å thÞ hiÖn t¹i cña b¹n sÏ ®−îc göi ra m¸y in. +) In b»ng lÖnh tõ cöa sæ lÖnh: Tr−íc tiªn b¹n còng ph¶i chän cöa sæ h×nh lµm cöa sæ ho¹t ®éng b»ng c¸ch nhÊn chuét lªn nã hoÆc dïng lÖnh figure(n), sau ®ã b¹n dïng lÖnh in. >> print % prints the current plot to your printer LÖnh orient sÏ thay ®æi kiÓu in: KiÓu mÆc ®Þnh lµ kiÓu portrait, in theo chiÒu ®øng, ë gi÷a trang. KiÓu in landscape lµ kiÓu in ngang vµ kÝn toµn bé trang. KiÓu in tall lµ kiÓu in ®øng nh−ng kÝn toµn bé trang. §Ó thay ®æi kiÓu in kh¸c víi kiÓu mÆc ®Þnh, b¹n dïng lÖnh orient víi c¸c th«ng sè cña nã nh sau: >> orient % What is the current orientation ans= portrait >> orient landscape % print sideways on the page >> orient tall % stretch to fill the vertical page NÕu b¹n muèn t×m hiÓu kü h¬n vÒ chóng th× h·y xem trî gióp trùc tuyÕn vÒ chóng. 17.7 Thao t¸c víi ®å thÞ B¹n cã thÓ thªm nÐt vÏ vµo ®å thÞ ®· cã s½n b»ng c¸ch dïng lÖnh hold. Khi b¹n thiÕt lËp hold on, MATLAB kh«ng bá ®i hÖ trôc ®· tån t¹i trong khi lÖnh plot míi ®ang thùc hiÖn, thay vµo ®ã, nã thªm dêng cong míi vµo hÖ trôc hiÖn t¹i. Tuy nhiªn nÕu nh− d÷ liÖu kh«ng phï hîp víi hÖ trôc to¹ ®é cò, th× trôc ®−îc chia l¹i . ThiÕt lËp hold off sÏ bá ®i cöa sæ figure hiÖn t¹i vµ thay vµo b»ng mét ®å thÞ míi. LÖnh hold mµ kh«ng cã ®èi sè sÏ bËt t¾t chøc n¨ng cña chÕ ®é thiÕt lËp hold tr−íc ®ã. Trë l¹i víi vÝ dô tr−íc: >> x = linspace(0,2*pi,30); >> y = sin(x); >> z = cos(x); >> plot(x,y) H×nh 17.14 B©y giê gi÷ nguyªn ®å thÞ vµ thªm vµo ®−êng cosine
124 Simpo PDF Merge and Split Unregistered Version - http://www.simpopdf.com %Gi÷ nguyªn ®å thÞ ®· vÏ lóc tr−íc >> hold on % hµm logic nµy tr¶ vÒ gi¸ trÞ 1 (true) nÕu hold ë tr¹ng th¸i ON >> ishold ans = 1 >> plot(x,z,'m') >> hold off % hold b©y giê kh«ng cßn ë tr¹ng th¸i ON n÷a. >> ishold ans = 0 Chó ý r»ng ®Ó kiÓm tra tr¹ng th¸i cña hold ta cã thÓ dïng hµm ishold . H×nh 17.15 NÕu b¹n muèn hai hay nhiÒu ®å thÞ ë c¸c cöa sæ figure kh¸c nhau, h·y dïng lÖnh figure trong cöa sæ lÖnh hoÆc chän new figure tõ b¶ng chän file, figure kh«ng cã tham sè sÏ t¹o mét figure míi. B¹n cã thÓ chän kiÓu figure b»ng c¸ch dïng chuét hoÆc dïng lÖnh figure(n) trong ®ã n lµ sè cöa sæ ho¹t ho¹t ®éng. MÆt kh¸c mét cöa sæ figure cã thÓ chøa nhiÒu h¬n mét hÖ trôc. LÖnh subplot(m,n,p) chia cöa sæ hiÖn t¹i thµnh mét ma trËn mxn kho¶ng ®Ó vÏ ®å thÞ, vµ chän p lµ cöa sæ ho¹t ®éng. C¸c ®å thÞ thµnh phÇn ®−îc ®¸nh sè tõ tr¸i qua ph¶i, tõ trªn xuèng d−íi, sau ®ã ®Õn hµng thø hai .v.v. . . VÝ dô: >> x = linspace(0,2*pi,30); >> y = sin(x); >> z = cos(x); >> a = 2*sin(x).*cos(x); >> b = sin(x)./(cos(x)+eps); >> subplot(2,2,1) % pick the upper left of % 2 by 2 grid of subplots >> plot(x,y),axis([0 2*pi -1 1]),title('sin(x)') >> subplot(2,2,2) % pick the upper right of the 4 subplots >> plot(x,z),axis([0 2*pi -1 1]),title('cos(x)') >> plot(x,z),axis([0 2*pi -1 1]),title('cos(x)') >> subplot(2,2,3)% pick the lowwer left of the 4 subplots >> plot(x,a),axis([0 2*pi -1 1]),title('2sin(x)cos(x)') >> subplot(2,2,4)%pick the lowwer right of the 4 subplots >> plot(x,b),axis([0 2*pi -20 20]),title('sin(x)/cos(x)')
125 Simpo PDF Merge and Split Unregistered Version - http://www.simpopdf.com H×nh 17.6 17.8 Mét sè ®Æc ®iÓm kh¸c cña ®å thÞ trong hÖ to¹ ®é ph¼ng • loglog t−¬ng tù nh− plot ngo¹i trõ thang chia lµ logarithm cho c¶ hai trôc. • semilogx t−¬ng tù nh− plot ngo¹i trõ thang chia cña trôc x lµ logarithm cßn thang chia trôc y lµ tuyÕn tÝnh. • semology t−¬ng tù nh− plot ngo¹i trõ thang chia cña trôc y lµ logarithm, cßn thang chia trôc x lµ tuyÕn tÝnh. • area( x, y ) t−¬ng tù nh− plot (x,y) ngo¹i trõ kho¶ng c¸ch gi÷a 0 vµ y ®−îc ®iÒn ®Çy, gi¸ trÞ c¬ b¶n y cã thÓ ®−îc khai b¸o, nh−ng mÆc ®Þnh th× kh«ng. • S¬ ®å h×nh mói tiªu chuÈn ®−îc t¹o thµnh tõ lÖnh pie(a, b), trong ®ã a lµ mét vector gi¸ trÞ vµ b lµ mét vector logic tuú chän. VÝ dô: >> a = [.5 1 1.6 1.2 .8 2.1]; >> pie(a,a==max(a)); >> title('Example Pie Chart') H×nh 17.7 • Mét c¸ch kh¸c ®Ó quan s¸t d÷ liÖu ®ã lµ biªu ®å Pareto, trong ®ã c¸c gi¸ trÞ trong c¸c vector ®−îc vÏ thµnh mét khèi ch÷ nhËt. VÝ dô dïng vector a ®· nãi ë trªn:
126 Simpo PDF Merge and Split Unregistered Version - http://www.simpopdf.com >> pareto(a); >> title('Example Pareto Chart') H×nh 17.18 • §«i khi b¹n muèn vÏ hai hµm kh¸c nhau trªn cïng mét hÖ trôc mµ l¹i sö dông thang chia kh¸c nhau, plotyy cã thÓ lµm ®iÒu ®ã cho b¹n: >> x = -2*pi:pi/10:2*pi; >> y = sin(x);z = 2*cos(x); >> subplot(2,1,1),plot(x,y,x,z), >> title('Two Plots on the same scale'); >> subplot(2,1,2),plotyy(x,y,x,z) >> title('Two plots on difference scale.'); H×nh 17.19 • §å thÞ bar vµ stair cã thÓ sinh ra bëi viÖc dïng lÖnh bar, bar3, barh vµ stairs. D−íi ®©y lµ vÝ dô:
127 Simpo PDF Merge and Split Unregistered Version - http://www.simpopdf.com >> x = -2.9:0.2:2.9; >> y = exp(-x.*x); >> subplot(2,2,1) >> bar(x,y) >> title('Bar chart of bell Curve') >> subplot(2,2,2) >> bar3(x,y) >> title('3-D Bar Chart of a Bell Cuve') >> subplot(2,2,3) >> stairs(x,y) >> title('Stair Chart of a Bell Curve') >> subplot(2,2,4) >> barh(x,y) >> title('Horizontal Bar Chart') H×nh 17.20 • rose(V) vÏ mét biÓu ®å trong to¹ ®é cùc cho c¸c gãc trong vector v, t−¬ng tù ta còng cã c¸c lÖnh rose(v,n) vµ rose(v,x) trong ®ã x lµ mét vector. D−íi ®©y lµ mét vÝ dô: >> v = randn(100,1)*pi; >> rose(v) >> title('Angle Histogram of Random Angle')
128 Simpo PDF Merge and Split Unregistered Version - http://www.simpopdf.com H×nh 17.21 ----------------------oOo---------------------- ch−¬ng 18 §å ho¹ trong kh«ng gian 3 chiÒu MATLAB cung cÊp mét sè hµm ®Ó hiÓn thÞ d÷ liÖu 3 chiÒu nh− c¸c hµm vÏ ®−êng th¼ng trong kh«ng gian 3 chiÒu, c¸c hµm vÏ bÒ mÆt vµ vµ khung d©y vµ mµu cã thÓ ®−îc sö dông thay thÕ cho chiÒu thø t−. 18.1 §å thÞ ®−êng th¼ng. LÖnh plot tõ trong kh«ng gian hai chiÒu cã thÓ më réng cho kh«ng gian 3 chiÒu b»ng lÖnh plot3. Khu«n d¹ng cña plot3 nh− sau: plot3 ( x1, y1, z1, S1, x2, y2, z2, S2, .... ), trong ®ã xn, yn vµ zn lµ c¸c vector hoÆc ma trËn, vµ Sn lµ x©u kÝ tù tuú chän dïng cho viÖc khai b¸o mµu, t¹o biÓu t−îng hoÆc kiÓu ®−êng. Sau ®©y lµ mét sè vÝ dô: >> t = linspace (0, 10*pi); >> plot3(sin(t),cos(t),t) >> title ('Helix'),xlabel('sin(t)') >> ylabel('cos(t)'),zlabel('t')
129 Simpo PDF Merge and Split Unregistered Version - http://www.simpopdf.com H×nh 18.1 Chó ý r»ng: hµm zlabel t−¬ng øng víi hµm hai chiÒu xlabel vµ ylabel. T−¬ng tù nh− vËy, lÖnh axis còng cã khu«n d¹ng: axis ( [xmin xmax ymin ymax zmin zmax ] ) thiÕt lËp giíi h¹n cho c¶ 3 trôc. VÝ dô : >> axis('ij') % thay ®æi h−íng trôc tõ sau ra tr−íc H×nh 18.2 Hµm text còng cã khu«n mÉu nh− sau: text ( x, y, z, string ) sÏ ®Æt vÞ trÝ x©u ‘string ‘ vµo to¹ ®é x, y, z. 18.2 §å thÞ bÒ mÆt vµ l−íi MATLAB ®Þnh nghÜa bÒ mÆt l−íi b»ng c¸c ®iÓm theo h−íng trôc z ë trªn ®−êng kÎ « h×nh vu«ng trªn mÆt ph¼ng x-y. Nã t¹o lªn mÉu mét ®å thÞ b»ng c¸ch ghÐp c¸c ®iÓm gÇn kÒ víi c¸c ®−êng th¼ng. KÕt qu¶ lµ nã tr«ng nh− mét m¹ng l−íi ®¸nh c¸ víi c¸c m¾t l−íi lµ c¸c ®iÓm d÷ liÖu. §å thÞ l−íi nµy th−êng ®−îc sö dông ®Ó quan s¸t nh÷ng ma trËn lín hoÆc vÏ nh÷ng hµm cã hai biÕn.
130 Simpo PDF Merge and Split Unregistered Version - http://www.simpopdf.com B−íc ®Çu tiªn lµ ®a ra ®å thÞ l−íi cña hµm hai biÕn z = f (x, y ), t−¬ng øng víi ma trËn X vµ Y chøa c¸c hµng vµ c¸c cét lÆp ®i lÆp l¹i. MATLAB cung cÊp hµm meshgrid cho môc ®Ých nµy. [ X, Y ] = meshgrid(x, y ), t¹o mét ma trËn X, mµ c¸c hµng cña nã lµ b¶n sao cña vector x, vµ ma trËn Y cã c¸c cét cña nã lµ b¶n sao cña vector y. CÆp ma trËn nµy sau ®ã ®−îc sö dông ®Ó −íc l−îng hµm hai biÕn sö dïng ®Æc tÝnh to¸n häc vÒ m¶ng cña MATLAB. Sau ®©y lµ mét vÝ dô vÒ c¸ch dïng hµm meshgrid. >> x = -7.5:.5:7.5; >> y = x; >> [X,Y] = meshgrid(x,y); X, Y lµ mét cÆp cña ma trËn t−¬ng øng mét l−íi ch÷ nhËt trong mÆt ph¼ng x-y. Mäi hµm z=f(x,y) cã thÓ sö dông tÝnh chÊt nµy. >> R = sqrt(X.^2+Y.^2)+eps; >> % find the distance from the origin (0,0) >> Z = sin(R)./R; % calculate sin(r)/ r Ma trËn R chøa b¸n kÝnh cña mçi ®iÓm trong [X,Y], nã lµ kho¶ng c¸ch tõ mçi ®iÓm ®Õn t©m ma trËn. Céng thªm eps ®Ó kh«ng ®Ó x¶y ra phÐp chia cho 0. Ma trËn Z chøa sine cña b¸n kÝnh chia cho b¸n kÝnh mçi ®iÓm trong s¬ ®å. C©u lÖnh sau vÏ ®å thÞ l−íi: >> mesh(X,Y,Z) H×nh 18.3 §å thÞ trªn lµ ®¬n s¾c. Tuy nhiªn b¹n cã thÓ thay ®æi mµu s¾c víi sù trî gióp cña MATLAB rÊt rÔ dµng nÕu b¹n ®äc ®Õn phÇn colormaps.. Trong vÝ dô nµy, hµm mesh x¾p xÕp gi¸ trÞ cña c¸c phÇn tö cña ma trËn vµo c¸c ®iÓm (XÞ,YÞ,ZÞ) trong kh«ng gian ba chiÒu. mesh còng cã thÓ vÏ mét ma trËn ®¬n t−¬ng tù nh− víi mét ®èi sè; mesh(Z), sö dông c¸c ®iÓm (i,j,ZÞ). Nh− vËy Z ®−îc vÏ ng−îc l¹i víi c¸c chØ sè cña nã, trong tr−êng hîp nµy mesh(Z) chØ ®¬n gi¶n lµ chia l¹i ®é kh¾c c¸c trôc x, y theo c¸c chØ sè cña ma trËn Z. B¹n h·y thö t¹o vÝ dô cho tr−êng hîp nµy?.