Giao trinh matlab v5.2 P18

Các tham chiếu gián tiếp đến các chương trình khả dụng hay các bộ phận dữ liệu đã được xác định từ trước cho phép nhiều ngôn ngữ định hướng ứng dụng tích hợp được các thao tác khác nhau.
Nội dung trích xuất từ tài liệu:
Giao trinh matlab v5.2 P18Toolbox - Digital signal Processingnhãm tiÕp theo cña lÖnh th× rÊt lín cho bé nhí cña m¸y tÝnh cña b¹n, b¹n cã muèn gi¶mnã thμnh sè nhá, nh− 32 » T = 5; » N = 512; » t = linspace (0, T, N + 1); t = (1 : N); » x1 = 2 * t / T - 1/2 ; x2 = 2 * (T - t) / T - 1/2; » x = min (x1, x2); % tÝn hiÖu tam gi¸c » plot (t, x)Chóng ta tÝnh hÖ sè cña sines vμ cosine. » X = fft (x); » A = real (X) / N; % hÖ sè cosine » B = - imag (X) / N); % hÖ sè sine » sumcos = zeros (N, N); » sumsin = zeros (N, N); » for h = 1 : N sumcos (h, : ) = A(h) * cos (2 * pi * (h - 1) * t/T); sumsin (h, : ) = B (h) * sin (2 * pi * (h - 1) * t/T); end » y = sum (sumcos + sumsin);Chóng ta cã thÓ kiÓm tra c¸c kÕt qu¶ b»ng c¸ch so s¸nh x vμ y, ®å häa cña chóng » plot (t, x, t, y);vμ sè »max (abs (x - y))9. Nh÷ng tÝn hiÖu tÇn sè cao vμ ký hiÖuë h×nh 1.12 ®· chØ ra sù t−¬ng øng gi÷a c«ng suÊt cña tÝn hiÖu vμ biÕn ®æi Fouriercña nã ®èi víi c¸c tÇn sè ®Õn tÇn sè Nyquist. §iÒu nμy trë nªn thó vÞ ®Ó xem ®iÒug× x¶y ra khi chóng ta lÊy mÉu t¹i kho¶ng thêi gian Ts h»ng sè tÝn hiÖu tuÇn hoμnliªn tôc cña tÇn sè cao ®Õn tÇn sè Nyquist Nf = 1/ (2Ts). Nh− chóng ta nh×n thÊyë ®©y, version lÊy mÉu cña tÝn hiÖu ®ång nhÊt víi tÝn hiÖu kh¸c tÇn sè thÊp. HiÖnt−îng nμy gäi lμ dÊu hiÖu tõ C1 , tõ ý nghÜa Latin “other”, nh÷ng c¸i kh¸c. §Ó 176PhÇn 2 - øng dôngToolbox - Digital signal ProcessingnhÊn m¹nh ý nμy chóng ta chän T lμ 5 gi©y, N = 16 lÊy mÉu trong mét chu kú, vμhiÖn ra theo kho¶ng lÊy mÉu víi Ts = T/N vμ tÇn sè mÉu víi fs > 1/Ts. TÝn hiÖu tuÇn hoμn víi chu kú T cã chu kú c¬ b¶n cña T nèi T/k víi k phï hîp.Chóng ta chØ ra tÇn sè cña nã k/T, víi f nhá. Còng nh− tÝn hiÖu, cho kho¶ng c¸ch sin(2πft) vμ cos (2πft). TÇn sè f cã thÓ lu«n viÕt nh− sau f = fapp + nfs Trong ®ã n vμ sè nguyªn vμ 0 ≤ |fapp| < Nf. Nã dÔ dμng kiÓm tra r»ng t¹i c¸c t béisè cña Ts nh− sau t = hTs, sin (2πft) = sin (2πfappt). Thùc tÕ sin(2πft) = sin (2π (fapp + nfs) t) = sin (2π (fapp + nfs) hTs) = sin (2π (fapp hTs + 2πnfshTs) = sin (2π (fapp hTs + 2πnh) = sin (2π fapp)Song tÝn hiÖu x = sin(2π ft), tÇn sè f, khi lÊy mÉu ë tÇn sè fs, lμ kh«ng thÓ ph©n biÖt ®−îctõ tÝn hiÖu x1 = sin(2πfappt) cña tÇn sè thÊp fapp MATLAB cho phÐp chóng ta gi¶i quyÕtvÊn ®Ò vμ biÓu diÔn c¸c dÊu hiÖu. H·y dïng m tÖp sau; alias.m:T=5; % tÇn sè c¬ b¶nNp = 512; %Sè ®iÓm ®Ó chÊmt = linspace(0,T,Np+1;t = t(1:Np); % t×m ®é ph©n gi¶i cña thêi gian %®Ó chÊm ®iÓmN=16; % sè ®iÓm lÊy mÉuTs =T/N; % kho¶ng lÊy mÉufs =1/Ts; % tÇn sè lÊy mÉuts = Ts*(0:(N-1)); % kho¶ng thêi gian lÊy mÉuNf = 1/(2*Ts); % TÇn sè Nyquistf = k/T; % tÇn sè liªn tôc % tÝn hiÖux = sin(2*pi*f*t); % tÝn hiÖu, ®é ph©n gi¶i caoxs = sin(2*pi*f*ts); % tÝn hiÖu, lÊy mÉu ph©n gi¶i% t×m fapp, nh− sau: f =n*fs+fappn = round(f/fs); 177PhÇn 2 - øng dôngToolbox - Digital signal Processingfapp = f-n*fn;xa = sin(2*pi*fapp*t);plot(t,[x;xa],ts,xs,0);str1 = [fs = , num2str(fs), Nf = ,num2str(Nf)];str2 = [k = , num2str(k), f = ,num2str(f)];str3 = [fapp=, num2str(fapp)];str = [str1, ,str2, , str3];title(str);Ch¹y chóng víi lÖnh sau » k= 17; alias H×nh 1.17 tÝn hiÖu tÇn sè cao lÊy mÉu nh− mét tÇn sè thÊp.VÝ dô 1.8: Giao ®éng cña mét tÊm ViÖc tÝnh to¸n ë vÝ dô 1.5 vμ 1.7 cã thÓ cã mét øng dông kü thuËt m« t¶ trong:M¸y kiÓm tra giao ®éng. VÝ dô ®¬n gi¶n cã d¹ng nh− h×nh 1.17. C¸c bé phËn ho¹t ®éngcña m¸y lμ 4 trôc quay, kh«ng cã khèi l−îng giao ®éng m1 ®Õn m4. Nh− m« t¶ ë trªnh×nh1.17 (a), khèi l−îng kh«ng giao ®éng cã thÓ ®o¹n cña vßng trong lμm b»ng s¾t(thÐp) vμ tùa trªn ®Üa quay. Khèi l−îng m1 vμ m2 b»ng nhau, nh−ng quay theo hai h−íng 178PhÇn 2 - øng dôngToolbox - Digital signal Processing®èi nhau, vμ còng nh− vËy ®èi víi khèi l−îng m3 vμ m4. Mét trong nh÷ng bé phËn ®−îcchØ chi tiÕt trªn h×nh 1.17 (b). Cho r»ng kho¶ng c¸ch gi÷a trôc quay qua ®iÓm 0 vμ t©mcña khèi l−îng kh«ng giao ®éng, mi lμ ri ...