Tài liệu matlap toàn tập_3

Tham khảo tài liệu tài liệu matlap toàn tập_3, công nghệ thông tin, kỹ thuật lập trình phục vụ nhu cầu học tập, nghiên cứu và làm việc hiệu quả
Nội dung trích xuất từ tài liệu:
Tài liệu matlap toàn tập_3 51zeros Ma trËn kh«ngVÝ dô VÊn ®Ò: Ta cã m¹ch ®iÖn nh− trong h×nh 7.1 ®−îc m« t¶ b»ng ph−¬ng tr×nh ®iÖn ¸p nót khi nguån®a vµo lµ sãng h×nh sin. H×nh 7.1E = 10o; R1 = 2; L = 10j; C = ; R2 = 10.ë ®©y vi lµ ®iÖn ¸p gi÷a nót thø i vµ ®Êt. Hái ®iÖn ¸p t¹i mçi nót lµ bao nhiªu? Gi¶i ph¸p: §©y lµ vÊn ®Ò vÒ ph©n tÝch pha. Ph−¬ng ph¸p gi¶i bµi nµy lµ gi¶i ph¬−ng tr×nh trªn, vµchuyÓn c¸c kÕt qu¶ vÒ d¹ng thêi gian. Trong MATLAB gi¶i ph¸p sÏ lµ:function circuit % circuit.m script file to solve circuit proplemA(1,1)=1/2; % poke in nonzero values as neededA(1,2)=-1/2;A(2,1)=-1/2;A(2,2)=1/2 + 0.2j + 1/10j;A(2,3)= -1/10j;A(3,2)=-1/10j;A(3,3)=1/10 + 1/10j;y=[-1 0 0]; % right hand side vectorv=Ay % complex solutionvmag=abs(v) % solution magnitudesvphase=angle(v)*180/pi % solution phase in degreestheta=linspace(0,2*pi); % plot results in timev1=vmag(1)*cos(theta-vphase(1));v2=vmag(2)*cos(theta-vphase(2));v3=vmag(3)*cos(theta-vphase(3));thd=theta*180/pi;plot(thd,v1,thd,v2,thd,v3)Sau khi ch¹y ch−¬ng tr×nh trªn, kÕt qu¶ sÏ lµ:v= -4.0000 + 6.0000i -2.0000 + 6.0000i 2.0000 + 4.0000ivmag = 7.2111 52 6.3246 4.4721vphase = 123.6901 108.4349 63.4349 H×nh 7.2 ---------------------oOo---------------------ch−¬ng 8 c¸c phÐp tÝnh LOGIC Vµ QUAN HÖ Thªm vµo nh÷ng to¸n tö ‘truyÒn thèng’, MATLAB cung cÊp to¸n tö logic vµ quan hÖ. B¹n cã thÓquen thuéc víi nh÷ng phÐp to¸n nµy, nÕu b¹n ®· lµm quen víi c¸c ng«n ng÷ lËp tr×nh kh¸c. Môc ®Ýchcña nh÷ng to¸n tö vµ hµm nµy lµ ®Ó tr¶ lêi c©u hái True_False (®óng_sai). §èi víi c¸c sè th× trong to¸n tö logic vµ quan hÖ quy ®Þnh c¸c sè kh¸c kh«ng lµ True cßn sè kh«nglµ False. KÕt qu¶ cña phÐp to¸n logic vµ quan hÖ ®a ra lµ 1 cho True, 0 cho False. 538.1 To¸n tö quan hÖ To¸n tö quan hÖ MATLAB bao gåm tÊt c¶ c¸c phÐp so s¸nh:To¸n tö quan hÖ ý nghÜa < nhá h¬n lín h¬n >= lín h¬n hoÆc b»ng == b»ng ~= kh«ng b»ng To¸n tö quan hÖ MATLAB cã thÓ dïng ®Ó so s¸nh hai m¶ng cã cïng kÝch cì hoÆc so s¸nh métm¶ng víi mét sè ®¬n. Trong tr−êng hîp thø hai, sè ®¬n so s¸nh víi tÊt c¶ c¸c phÇn tö cña m¶ng, kÕtqu¶ tr¶ vÒ gièng nh− kÝch cì cña m¶ng. VÝ dô:>> A = 1:9, B = 9 - AA= 1 2 3 4 5 6 7 8 9B= 8 7 6 5 4 3 2 1 0>> tf = A>4tf= 0 0 0 0 1 1 1 1 1t×m kiÕm c¸c phÇn tö cña A mµ lín h¬n 4. KÕt qu¶ b»ng 0 khi A 4, b»ng 1 khi A>4.>> tf = (A==B)tf= 0 0 0 0 0 0 0 0 0 T×m kiÕm c¸c phÇn tö cña A mµ b»ng víi B. Chó ý sù kh¸c nhau gi÷a = vµ == dïng ®Ó so s¸nh haibiÕn vµ tr¶ vÒ 1 khi chóng b»ng nhau, 0 khi chóng kh¸c nhau; = dïng ®Ó g¸n kÕt qu¶ ®a ra cña to¸n töcho mét biÕn.>> tf = B - (A>2)tf= 8 7 5 4 3 2 1 0 -1 T×m c¸c phÇn tö A>2 vµ bÞ trõ bëi vector B. VÝ dô nµy chØ ra r»ng kÕt qu¶ ®a ra cña to¸n tö logic lµmét m¶ng sè bao gåm c¸c sè kh«ng vµ mét, chóng còng cã thÓ dïng trong c¸c phÐp to¸n sè häc.>> B = B + (B==0)*epsB= Columns 1 through 7 8.0000 7.0000 6.0000 5.0000 4.0000 3.0000 2.0000 Columns 8 through 9 1.0000 0.0000 VÝ dô trªn ®a ra c¸ch thay thÕ c¸c phÇn tö cña B mµ trïng víi kh«ng b»ng sè ®Æc biÖt cña MATLABlµ eps, cã gi¸ trÞ xÊp xØ 2.2e-16. C¸ch thay thÕ nµy ®«i khi cã Ých lµ tr¸nh tr−êng hîp chia cho sèkh«ng nh− vÝ dô sau:>> x = (-3:3)/3 54 -1.0000 -0.6667 -0.3333 0 0.3333 0.6667 1.0000>> sin(x)./xWarning: Divide by zeroans= 0.8415 0.9276 0.9816 NaN 0.9816 0.9276 0.8415 TÝnh to¸n hµm sin(x)/ x ®a ra mét c¶nh b¸o v× phÇn tö thø t− b»ng kh«ng, sin(0)/ 0 kh«ng ®−îc®Þnh nghÜa, MATLAB tr¶ l¹i NaN ( nghÜa lµ kh ...