Giao trinh matlab v5.2 P15

Tùy theo ngôn ngữ, các vùng trống gián tiếp có thể bao gồm các tham chiếu, mà thật ra, chúng là các con trỏ (pointer) chỉ đến những vùng chứa khác của bộ nhớ, được cài đặt trong các biến hay nhóm các biến.
Nội dung trích xuất từ tài liệu:
Giao trinh matlab v5.2 P15ToolBox - Digital Signal Processing » V= 0.01; %V » R= 15; % Ohm » omega= 2*pi*50*10^6; %rad/s » C= 100*10^(-12); %F » Zc= -j/(omega*C); %Ohm Ta tÝnh ®−îc tæng trë cña m¹ch » Z= R+ Zc Z= 16.0000 -31.8310i Vμ dßng ®iÖn » I = V/Z I= 1.2606e-04+ 2.5097e-04i §iÖn ¸p r¬i trªn ®iÖn trë » Vr = R*I Vr= 0.0018 + 0.0039i Vμ qua tô ®iÖn chóng ta cã thÓ ®å thÞ biªn pha . Sö dông MATLAB subplot » subplot (1, 2, 1) C¸c th«ng sè 1, 2, 1 cã ý nghÜa lμ 1 vμ 2 ma trËn cña h×nh vÏ sinh ra vμ chÊm ®iÓmtheo mét ®iÓm ®Çu. Trong MATLAb 3.5 c¸c th«ng sè nμy sÏ ®−îc viÕt lμ (121) Gåm c¸c gi¸ trÞ V, Vr, Vc nh− pha cuèi cïng. §Ó ®Þnh nghÜa pha chóng ta cÇn ®−a®iÓm gèc cña chóng » VV= [0 V]; VVc = [0 Vc]; VVr = [0 Vr] TiÕp ®Õn ta chÊm ®iÓm ®iÖn thÕ pha » plot (real(VV) , imag(VV) §Ó quan s¸t ®−îc thang gãc pha chÝnh x¸c cña trôc thËt cÇn ph¶i nh− ®èi víi trôc¶o . Cã nghÜa lμ khung chÊm ®iÓm lμ h×nh vu«ng » axis(square) Khung cña c¸c trôc sÏ ®−îc sinh ra nh− sau » axis ([0 0.012 -0.006 0.006]) 144PhÇn 2 - øng dôngToolBox - Digital Signal Processing PhÇn më réng cña c¶ hai trôc lμ 0.012. C¸c gi¸ trÞ nμy cã thÓ ®−îc nhËn bëi thönghiÖm vμ sai sè, nh− ®−îc hiÖn thÞ ë ®å thÞ , khi ®¸nh gi¸ trôc tèt nhÊt vμ dïng c©ulÖnh axis .HoÆc ta cã thÓ sö dông mét hμm kh¸c lμ abs ®Ó t×m gi¸ trÞ biªn pha lín nhÊtvμ ®Ó ®iÒu chØnh l¹i trôc Chóng ta còng cã thÓ dïng hμm hold ®Ó gi÷ cho ®å thÞ chång lªn hai h×nhkh¸c,Vr vμ Vc » hold on » plot(real(VVr), imag(VVr)) » plot(real(VVc), imag(VVc)) Toμn bé ®å thÞ sÏ ®−îc viÕt nh− sau: » title((a)) » xlabel(Real), ylabel(Imaginary) B¹n còng cã thÓ x¸c ®Þnh c¸c h×nh b»ng c¸ch tù ®¸nh dÊu c¸c ®iÓm tøc lμ dïnghμm gtext » text(real(V), imag(V), V) » text(real(Vr), imag(Vr), Vr) » text(real(Vc), imag(Vc), Vc) Sau ®ã chÊm ®iÓm toμn bé » hold off Ta còng cã thÓ tÝnh tæng cña ®iÖn ¸p nh− sau(xem thªm h×nhvÏ): » Vr + Vc ans = 0.0100 - 0.0000i » 180*(angle(Vr) - -angle(Vc))/pi ans = 90 145PhÇn 2 - øng dôngToolBox - Digital Signal Processing H×nh 1.3. DiÖn ¸p trªn m¹ch nãi tiÕp RC §Ó cã ®−îc ®å thÞ vμ tÝnh to¸n ®iÖn ¸p ta lμm theo nh÷ng b−íc sau » f= 50*10^6; % tÇn sè, Hz » T = 1/f; % chu kú, s » omega= 2*pi*f; %tÇn sè gãc, rad/s » t=0: T/50 : 2*T; % chuçi c¸c gi¸ trÞ » v = V*sin(omega*t); » vr = abs(Vr)*sin(omega*t + angle(Vr)); » vc = abs(Vc)*sin(omega*t + angle(Vc)); §Ó cã ®−îc ®iÓm nh− phÇn tö thø hai cña ma trËn ®å thÞ . » subplot(1, 2, 2) vμ hcÊm ®iÓm theo nh÷ng lÖnh sau: » plot (t, 1000*v, t, 1000*vr, t, 1000*vc) 146PhÇn 2 - øng dôngToolBox - Digital Signal Processing » grid, title((b)) » xtabel(t,s), ylabel(mV) »text(t(5), 1000*v(5), V) »text(t(20), 1000*v(20), Vr) »text(t(50), 1000*v(5), Vc) Nh− ta ®· biÕt m¹ch RC lμ mét m¹ch läc th«ng thÊp. Gi¶ sö ta cho vμo ®Çu vμo10mV, b©y giê ta xem ®Çu ra cña tô ®iÖn. Ký hiÖu ®iÖn ¸p vμo lμ Vin vμ ®iÖn ¸p ra Vout . 1 Z = R + Rc = R + jωC Vin 1 I= = Vin R − j / ωC Z 1 1 Vout = ZcI = Vin jωC R − j / ωC Hμm ®¸p øng tÇn (dïng hμm biÕn ®æi) cña bé läc x¸c ®Þnh bëi tØ sè Vout 1 H ( jω ) = = Vin 1 + jωRC Gi¸ trÞ ωc = 1/RC gäi lμ tÇn sè c¾t ë d¹ng kh«ng thø nguyªn 1 H ( jω ) = 1 + jω / ω c hoÆc ω/ωc = 2πf/2πfc =f/fc Chóng ta nhËn ®−îc d¹ng kh«ng thø nguyªn cña hμm biÕn ®æi 1 H( f ) = 1 + jf / f cCh ...